Численные методы вычисления корней

[XeN0N]

вот оно

[BarT(48rus)]

нет, спасибо канечно, но нето....ладно, что нибудь придумаем)

[XeN0N]

табулирование функции для вычисление отрезков изоляции:

Общая идея в том, что если на концах отрезка функция имеет разные значения, значит в нем она проходит через 0 и там нужно искать корень

	for (i = a; i < b; i+=h) {
		n=d=0;
		if(f(i)==0&&i==a) out <<i<<endl; // если начало промежутка - решение, то вывести его
		if((f(i)>0&&f(i+h)<0) || (f(i)<0&&f(i+h)>0)) // если на концах промежутка функция имеет разные знаки, то найти в нем корень и вывести в файл
		{
			out <<solve(i,i+h,e,&n,&d);
			out <<" n="<<n<<" d="<<d<<endl;
		}
		if (i+h)
			if(f(i+h)==0) out <<i<<endl; // если узловая точка или конец промежутка - решение, то вывести его
	}

Это кусок программы, которая решает гибридным методом ньютона-половинных делений.
функция solve - найти корень на данном промежутке, n и d - переменные показывающие сколько было ньютоновских операций и сколько из них не попали в отрезок изоляции
a - начало промежутка
b - конец
h - шаг
f(x) - заданная функция
out <<n; - вывести в файл out переменную n

[BarT(48rus)]

Друзья, не подскажете команду, которая обрывает цикл в паскале, а то чо то не найду никак))))заранее спасибо)))

[XeN0N]

break .

[BarT(48rus)]

что то походу он и внутренний repeat и внешний вайл рвет)))что то не то))))может еще варианты есть?

[Sazary]

Сообщение от BarT(48rus)

что то походу он и внутренний repeat и внешний вайл рвет)))

break выходит только из внутреннего цикла (из которого он был вызван)

[BarT(48rus)]

хм))))ну может тогда глянете код???может что то не правильно???

writeln('xmax=',xmax);
  while (xmin<=xmax) do
  begin
   n:=0;
    repeat
     xc:=xmin;
     xn:=8*sin(xc)/p[i];
     dx:=abs(xn-xc);
     n:=n+1;
     if (n>10000) then break
    until dx<e;
  xmin:=xmin+1;
  end;
 writeln ('xn=', xn:12:5);
 writeln('dx=', dx:12:7);
 writeln('N=', n);
 writeln('p[',i,']=', p[i]);
 end;
end.

а выходит он когда первое только более 10000, а дальше не хочет сравнивать(

[Sazary]

Сложно искать ошибки, когда код нельзя скомпилировать. Где остальная часть? И что она должна делать? (приведенное задание в первом посте врядли подойдет, т.к. не в курсе что такое интервал изоляции). А также: что не работает?

То, что break выходит только из внутреннего цикла может проиллюстрировать следующая программа:

uses crt;
var
a,b : integer;
begin
clrscr;
b:=0;
while b<3 do
 begin
 writeln('b= ',b);
 a:=0;
 repeat
 inc(a);
 if a=2 then break;
 until a=3;
 inc(b);
 end;
readln;
end.

[BarT(48rus)]

хм, странно, ну тогда глянь всю прогу, она вот)

program laba2;
type
Massiv=array[1..4] of real;
var
p:massiv;
n, i:integer;
x0, e, xc, xn, dx, xmin, xmax:real;
begin
write ('e=');
readln (e);
for i:=1 to 4 do
 begin
  write('p[',i,']=');
  readln(p[i]);
  xmin:=-8/p[i];
  writeln('xmin=',xmin);
  xmax:=8/p[i];
  writeln('xmax=',xmax);
  while (xmin<=xmax) do
  begin
   n:=0;
    repeat
     xc:=xmin;
     xn:=8*sin(xc)/p[i];
     dx:=abs(xn-xc);
     n:=n+1;
     if (n>10000) then break
    until dx<e;
  xmin:=xmin+1;
  end;
 writeln ('xn=', xn:12:5);
 writeln('dx=', dx:12:7);
 writeln('N=', n);
 writeln('p[',i,']=', p[i]);
 end;
end.

она компилится, но рабит не так, как задумано, она должна выходить из внутреннего цикла при условии n>10000, а она похоже выходит из всей, ды не похоже, а точно)