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

Program n;
uses crt;
const e=0.0001;
var
   x0,y0,x,y,gy,I,Ix,Iy,dx,dy:real;
   d:boolean;
Function F1(x,y:real):real;
begin
     F1:=((sin(x*y+0.1))/cos(x*y+0.1))-x*x;
end;
function F2(x,y:real):real;
begin
     F2:=(0.5*x*x)+(2*y*y)-1;
end;
function F1x(x:real):real;
begin
     F1x:=(y/cos(x*y+0.1))-2*x;
end;
function F2x(x:real):real;
begin
     F2x:=4*y;
end;
function F1y(y:real):real;
begin
     F1y:=(x/cos(x*y+0.1));
end;
function F2y(y:real):real;
begin
     F2y:=x;
end;
Begin  ClrScr;
       x0:=0; y0:=0.1; gy:=1;
  Repeat
        d:=false;
        I:=F1x(x0)*gy-F2x(x0)*F1y(y0);
        Ix:=-F1(x,y)*gy+F2(x,y)*F1y(y0);
        Iy:=-F1x(x0)*F2(x,y)+F2x(x0)*F1(x,y);
     dx:=Ix/I;
     dy:=Iy/I;
     x:=x0+dx;
     y:=y0+dy;
   if (abs(x-x0)<e) and (abs(y-y0)<e) then d:=true;
      x0:=x;
      y0:=y;
  Until d=true;
        writeln('Reshenie sistemy');
        writeln('x1=',x:4:3,' x2=',y:4:3);
 Readln;
   End.