Составить подпрограмму определения максимума функции, заданной в дискретных точках

var
  nx,ny,i,j:integer;
  ax,bx,ay,by,x0,y0,max,maxx,maxy,hx,hy,x,y,z:real;
begin
  writeln('Введите диапазон изменения X от a до b (два числа через пробел)');
  readln(ax,bx);
  writeln('Введите диапазон изменения Y от a до b (два числа через пробел)');
  readln(ay,by);
  writeln('Введите число разбиений по X (целое)');
  readln(nx);
  writeln('Введите число разбиений по Y (целое)');
  readln(ny);
  x0 := ax;
  y0 := ay;
  hx := (bx-ax)/Nx;
  hy := (by-ay)/Ny;
  x := x0;
  y := y0;
  max := 6*sqr(x)+sqr(y)-2*x*y+2;
  for i:=0 to Nx do begin
    x := x0 + i * hx;
    for j:=0 to Ny do begin
      y := x0 + i * hx;
      z := 6*sqr(x)+sqr(y)-2*x*y+2;
      if z > max then begin
        max := z;
        maxX := x;
        maxY := y;
      end;
    end;
  end;
  writeln('maxZ = ',max);
  writeln('x = ',maxX);
  writeln('y = ',maxY);
end.

uses crt;
type func=function(x,y:real):real;
{$F+}
function z1(x,y:real):real;
begin
z1:=6*x*x+y*y+2*x*y+2;
end;
function z2(x,y:real):real;
begin
z2:=2*x*x+3*y*y+2*x+3*y-3*x*y;
end;
{$F-}
procedure max(z:func;x0,y0,hx,hy:real;n:integer;var max:real);
var x,y,t:real;
    i,j:integer;
begin
max:=z(x0,y0);
x:=x0;
for i:=1 to n do
 begin
  y:=y0;
  for j:=1 to n do
   begin
    t:=z(x,y);
    if t>max then max:=t;
    y:=y+hy;
   end;
  x:=x+hx;
 end;
end;
var x0,y0,hx,hy,mx1,mx2:real;
    n:integer;
begin
clrscr;
write('Введите начало диапазона по Х x0=');
readln(x0);
write('Введите шаг по Х hx=');
readln(hx);
write('Введите начало диапазона по Y y0=');
readln(y0);
write('Введите шаг по Y hy=');
readln(hy);
writeln('Функция z1=6x^2+y^2-2xy+2');
max(z1,x0,y0,hx,hy,n,mx1);
writeln('Максимум=',mx1:0:3);
writeln('Функция z2=2x^2+3y^2+2x+3y-3xy');
max(z2,x0,y0,hx,hy,n,mx2);
writeln('Максимум=',mx2:0:3);
readln
end.

unit modul;
interface
type func=function(x,y:real):real;
function z1(x,y:real):real;
function z2(x,y:real):real;
procedure max(z:func;x0,y0,hx,hy:real;n:integer;var max:real);
implementation
{$F+}
function z1(x,y:real):real;
begin
z1:=6*x*x+y*y+2*x*y+2;
end;
function z2(x,y:real):real;
begin
z2:=2*x*x+3*y*y+2*x+3*y-3*x*y;
end;
{$F-}
procedure max;
var x,y,t:real;
    i,j:integer;
begin
max:=z(x0,y0);
x:=x0;
for i:=1 to n do
 begin
  y:=y0;
  for j:=1 to n do
   begin
    t:=z(x,y);
    if t>max then max:=t;
    y:=y+hy;
   end;
  x:=x+hx;
 end;
end;
end.

uses crt,modul;
var x0,y0,hx,hy,mx1,mx2:real;
    n:integer;
begin
clrscr;
write('Введите начало диапазона по Х x0=');
readln(x0);
write('Введите шаг по Х hx=');
readln(hx);
write('Введите начало диапазона по Y y0=');
readln(y0);
write('Введите шаг по Y hy=');
readln(hy);
writeln('Функция z1=6x^2+y^2-2xy+2');
max(z1,x0,y0,hx,hy,n,mx1);
writeln('Максимум=',mx1:0:3);
writeln('Функция z2=2x^2+3y^2+2x+3y-3xy');
max(z2,x0,y0,hx,hy,n,mx2);
writeln('Максимум=',mx2:0:3);
readln
end.