Інтерполяція функції

Uses Crt;

Const MaxX=50;

MaxY=50;

Type MyArr=Array[0..MaxX,0..MaxY] Of Real;

Var Nx,Ny,Cx,Cy:Integer;

X:Array[0..MaxX] Of Real;

Y:Array[0..MaxY] Of Real;

B:MyArr;

Xa,Xb,Ya,Yb:Real;

D1,D2:^MyArr;

cc,cc1:Integer;

Function Func(x,y:Real):Real;

Begin

Func:=1/(x*x+y*y+x*y);

End;

Procedure DataInput;

Var i,j:Integer;

Begin

{ Write('Input Xa : '); ReadLn(Xa);

Write('Input Xb : '); ReadLn(Xb);

Write('Input Ya : '); ReadLn(Ya);

Write('Input Yb : '); ReadLn(Yb);}

Xa:=1; Xb:=2; Ya:=1; Yb:=2;

{ Write('Input Nx : '); ReadLn(Nx);

Write('Input Ny : '); ReadLn(Ny);}