Найти экстремум функции 2-х переменных методом градиента

var E,x1,x2,xk1,xk2,dy1,dy2,d2y1,d2y2,d2y12,h,Hh,Norma,Yk,Yx: real;
    k,k1: integer;
function Y(x1,x2:real): real;
begin
    Y:=100*sqr(X1-sqr(X2))+sqr(1-X1);
end;
procedure dy(x1,x2: real; var dy1, dy2: real);
begin
    dy1:=200*X1-200*X2*X2-2*X1+2;
    dy2:=-400*X1*X2+400*X2*X2*X2;
end;
function N(dy1,dy2: real): real;
begin
    N:= sqrt(sqr(dy1)+sqr(dy2));
end;
begin
writeln('e -' + 'na4alnie pribl:');
write('E= '); readln(E);
write('X0 = (x0;y0) = '); readln(x1,x2);
write('shag h = '); readln(Hh);
writeln('to4nost ? = ',E:7:6);
k:=0;
yx:=y(x1,x2);
dy(x1,x2,dy1,dy2);
NORMA:=N(dy1,dy2);
if NORMA>E then begin
repeat
    h:=Hh;
    k1:=0;
    repeat
      if k1>0 then h:=h/2;
      xk1:=x1-h*dy1; xk2:=x2-h*dy2;
      x1:=xk1; x2:=xk2;
      yk:=y(xk1,xk2);
      inc(k1);
    until yk<=yx;
    yx:=yk;
    dy(x1,x2,dy1,dy2);
    NORMA:=N(dy1,dy2);
    inc(k);
until NORMA<=E; end
else
writeln('reshenie zada4i: minF(x) = ',yx:8:5,',');
writeln('X = [',x1:5:3,';',x2:5:3,'].');
end.