Const
 x01=5.0;
 y01=20.0;
 x02=1.0;
 y02=0.0;
 Eps=0.00001;

Type
 Matr=Array[1..2,1..2] of Real;

Var
 x,y,dx,dy,Det_0,Det_X,Det_Y:Real;
 M:Matr;

Function Determ(D:Matr):Real;
begin
 Determ:=D[1,1]*D[2,2]-D[1,2]*D[2,1];
end;

Function f1(z1,z2:real):Real;
begin
 f1:=Sqr(z1)-z2-1;
end;

Function f2(z1,z2:real):Real;
begin
 f2:=z1*6-z2-7;
end;

Function Der_x_f1(z1,z2:real):Real;
begin
  Der_x_f1:=z1*2;
end;

Function Der_y_f1(z1,z2:real):Real;
begin
  Der_y_f1:=-1.0;
end;

Function Der_x_f2(z1,z2:real):Real;
begin
 Der_x_f2:=6.0;
end;

Function Der_y_f2(z1,z2:real):Real;
begin
 Der_y_f2:=-1.0;
end;

Begin
 x:=x01;
 y:=y01;
 Repeat
  M[1,1]:=Der_x_f1(x,y);
  M[1,2]:=Der_y_f1(x,y);
  M[2,1]:=Der_x_f2(x,y);
  M[2,2]:=Der_y_f2(x,y);
  Det_0:=Determ(M);
  M[1,1]:=-f1(x,y);
  M[2,1]:=-f2(x,y);
  Det_X:=Determ(M);
  M[1,1]:=Der_x_f1(x,y);
  M[1,2]:=-f1(x,y);
  M[2,1]:=Der_x_f2(x,y);
  M[2,2]:=-f2(x,y);
  Det_Y:=Determ(M);
  dx:=Det_X/Det_0;
  dy:=Det_Y/Det_0;
  x:=x+dx;
  y:=y+dy;
 Until (Abs(dx)<Eps) and (Abs(dy)<Eps);
 Writeln('x1 = ',x:0:6);
 Writeln('y1 = ',y:0:6);
 Writeln;
 x:=x02;
 y:=y02;
 Repeat
  M[1,1]:=Der_x_f1(x,y);
  M[1,2]:=Der_y_f1(x,y);
  M[2,1]:=Der_x_f2(x,y);
  M[2,2]:=Der_y_f2(x,y);
  Det_0:=Determ(M);
  M[1,1]:=-f1(x,y);
  M[2,1]:=-f2(x,y);
  Det_X:=Determ(M);
  M[1,1]:=Der_x_f1(x,y);
  M[1,2]:=-f1(x,y);
  M[2,1]:=Der_x_f2(x,y);
  M[2,2]:=-f2(x,y);
  Det_Y:=Determ(M);
  dx:=Det_X/Det_0;
  dy:=Det_Y/Det_0;
  x:=x+dx;
  y:=y+dy;
 Until (Abs(dx)<Eps) and (Abs(dy)<Eps);
 Writeln('x2 = ',x:0:6);
 Writeln('y2 = ',y:0:6);
 Readln
End.