Решение нелинейных уравнений методом Ньютона

procedure root(g,h,dg,dh:FT;a,b,eps:real;var x:real);
var d,c:real; sw:boolean;
function f(x:real):real;
begin f:=g(x)-h(x) end;
function df(x:real):real;
begin df:=dg(x)-dh(x) end;
begin
 sw:=((f(a)>0) and ((f((a+b)/2))>((f(a)+f(b))/2))) or
     ((f(a)<0) and ((f((a+b)/2))<((f(a)+f(b))/2)));
 if sw then
   begin
    repeat
     c:=b-f(b)/df(b); b:=c
    until f(c)*f(c-eps)<0;
   x:=c-eps/2
   end  else
   begin
    repeat
     c:=a-f(a)/df(a); a:=c
    until f(c)*f(c+eps)<0;
   x:=c+eps/2
   end
end;