построение графика

uses GraphABC,PointRect;
 
 
function y(x: real;z:real): real;       //заданная функция
begin
  y:=(x+z)*(x+z)-(2*z-1)*(2*z-1);
end;
 
function xp(x:real;z:real;Alpha:real):real; //косоуг проекция
begin
     xp:=z*(1/2)*cos(Alpha)+x;
end;
 
function yp(x:real;z:real;Alpha:real):real; //косоуг проекция
begin
     yp:=z*(1/2)*sin(Alpha)+y(x,z);
end;
 
var
nx,nz,i,j:integer;
x1,x2,z1,z2,hx,hz,max_xp,max_yp,min_xp,min_yp,tec_x,tec_y,kx,ky,Alpha:real;
xp1,yp1,xp2,yp2,xp3,yp3,xp4,yp4:real;
a:array [1..4] of Point;
 
begin      //основная программа
writeln('Введите нижнюю границу x');
readln(x1);
writeln('Введите верхнюю границу x');
readln(x2);
writeln('Введите нижнюю границу z');
readln(z1);
writeln('Введите верхнюю границу z');
readln(z2);
writeln('Введите чиcло разбиения Nx');
readln(nx);
writeln('Введите чиcло разбиения Nz');
readln(nz);
 
hx:=(x2-x1)/nx;
hz:=(z2-z1)/nz;
Alpha:=Pi/3;
 
max_xp:=xp(x1,z1,Alpha);  //для коэф как кау
min_xp:=xp(x1,z1,Alpha);
max_yp:=yp(x1,z1,Alpha);
min_yp:=yp(x1,z1,Alpha);
 
 
for i:=0 to nx do
for j:=0 to nz do
begin
tec_x:=xp(x1+i*hx,z1+j*hz,Alpha);
tec_y:=yp(x1+i*hx,z1+j*hz,Alpha) ;
 
if tec_x>max_xp then max_xp:=tec_x;
if tec_y>max_yp then max_yp:=tec_y;
if tec_x<min_xp then min_xp:=tec_x;
if tec_y<min_yp then min_yp:=tec_y;
end;
 
 
kx:=WindowWidth/(max_xp-min_xp);
ky:=WindowHeight/(max_yp-min_yp);
 
 
for i:=0 to (nx-1) do
for j:=0 to (nz-1) do
 
begin
xp1:=xp(x1+i*hx,z1+j*hz,Alpha);
yp1:=yp(x1+i*hx,z1+j*hz,Alpha);
a[1]:=PointF(round(kx*(xp1-min_xp)), round(ky*(max_yp-yp1)) );
xp2:=xp(x1+(i+1)*hx,z1+j*hz,Alpha);
yp2:=yp(x1+(i+1)*hx,z1+j*hz,Alpha);
a[2]:=PointF(round(kx*(xp2-min_xp)), round(ky*(max_yp-yp2)) );
xp3:=xp(x1+(i+1)*hx,z1+(j+1)*hz,Alpha) ;
yp3:=yp(x1+(i+1)*hx,z1+(j+1)*hz,Alpha);
a[3]:=PointF(round(kx*(xp3-min_xp)), round(ky*(max_yp-yp3)) );
xp4:=xp(x1+i*hx,z1+(j+1)*hz,Alpha);
yp4:=yp(x1+i*hx,z1+(j+1)*hz,Alpha);
a[4]:=PointF(round(kx*(xp4-min_xp)), round(ky*(max_yp-yp4)) );
polygon(a,4);
end;
end.