Метод Ньютона для систем нелин. уравнений(Pascal)

[Мельхиора]

Добрый день. Необходимо написать программу для решения системы нелинейных уравнений методом Гаусса. Паскаль находит ошибку "Duplicate identifier" на выделенной красным строке. Не могу понять, что именно следует исправить, подскажите пожалуйста.

Program z3v;
uses crt;
const e=0.0001;
var
   x0,y0,x,y,gy,I,Ix,Iy,dx,dy:real;
   d:boolean;
Function f(x,y:real):real;
begin
     f:=2*sqr(x)+sqr(y)-1;
end;
function g(x,y:real):real;
begin
     g:=y+(sqr(x)*x);
end;
function fx(x:real):real;
begin
     fx:=2*x+2;
end;
function gx(x:real):real;
begin
     gx:=-1/(2*(x-1));
end;
function fy(y:real):real;
begin
     fy:=2*y+2;
end;
function gy(y:real):real;
     gy:=1;
end;
Begin  ClrScr;
       x0:=0; y0:=0.1; gy:=1;
  Repeat
        d:=false;
        I:=fx(x0)*gy-gx(x0)*fy(y0);
        Ix:=-f(x,y)*gy+g(x,y)*fy(y0);
        Iy:=-fx(x0)*g(x,y)+gx(x0)*f(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.

[Poma][a]

у тебя переменная gy и функция gy что-то надо переименовать

[Мельхиора]

Исправила методом тыка

end;
function qy(y:real):real;
     gy:=1;

g на q, больше ошибки здесь не выдает, Но
Новая ошибка в строчке

begin
     g:=y+(sqr(x)*x);
end;

"205: Floating point overflow"

[Poma][a]

function gy(y:real):real;
     gy:=1;
end;

begin'а нет

205 Floating point overflow (Переполнение при операции с плавающей запятой).
В результате выполнения операции с плавающей запятой получено слишком большое вещественное число.

[Мельхиора]

Begin заметила, дописала. Про плавающую точку нашла в интернете, но как, где и что конкретно здесь следует исправить, не пойму

Program z3v;
uses crt;
const e=0.0001;
var
   x0,y0,x,y,gy,I,Ix,Iy,dx,dy:real;
   d:boolean;
Function f(x,y:real):real;
begin
     f:=2*sqr(x)+sqr(y)-1;
end;
function g(x,y:real):real;
begin
     g:=y+(sqr(x)*x);
end;
function fx(x:real):real;
begin
     fx:=2*x+2;
end;
function gx(x:real):real;
begin
     gx:=-1/(2*(x-1));
end;
function fy(y:real):real;
begin
     fy:=2*y+2;
end;
function qy(y:real):real;
begin
     gy:=1;
end;
...

[Poma][a]

странно у меня ошибку не выдаёт...

[Мельхиора]

Сообщение от Poma][a

странно у меня ошибку не выдаёт...

Не могли бы вы, пожалуйста, написать ответ?
Задание уже решено методами Зейделя и Якоби, ответ
Решение системы: x[1]=0.001 x[2]=0.000