Вычисление корней любой функции методом Ньютона(касательных)

@Острослав · Регистрация: 17.03.2013

Author24 — интернет-сервис помощи студентам

Здравствуйте. Прошу помочь мне с программой по поиску корней методом касательных. Программа у меня присутствует, но она по каким-либо причинам не работает. Выводит все время 0. Что дано: есть некая функция, мы можем только подавать ей любое х и взамен получать у. И так вот надо найти корень методом касательных. Я написал функцию для поиска касательной(diff) и процедуру (kasat) для поиска корня. У нас еще есть некоторый отрезок [n;m]. И собственно он задается так, чтобы функция пересекала ось Ох только один раз.

Pascal

function diff(a,b,c:longint;x,h:extended):extended;
var dif,diff2,d:extended;
begin
 d:=0.001;
 dif:=-11111111;
 while abs(dif - diff2) > h do begin
  diff2:=dif;
  dif:=( func(a,b,c,(x+d)) - func(a,b,c,x) ) / d;
  d:= d / 2;
 end;
 diff:=dif;
end;

Pascal

procedure kasat(a,b,c,n,m:longint;h:extended);
var x,x2,xk,st:extended;
begin
 x:=n;
 x2:=m;
 xk:= h+h;
 if func(a,b,c,x)*diff(a,b,c,x,h)>0 then st:=x2
  else st:=x;
 while (func(a,b,c,xk)>-h) and (func(a,b,c,xk)<h) do begin
  xk:=st - func(a,b,c,st) / diff(a,b,c,st,h);
  st:=xk;
 end;
 write(xk:0:3);
end;

	17.03.2013, 12:03