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

Delta:=Xij(i,k)*Yij(j,k)/

( D1^[i,j]+

Teta(k,j)*D1^[i,k]+

Teta(k,i)*D1^[k,j]+

Teta(k,i)*Teta(k,j)*D1^[k,k]

);

End;

Var i,j,s,k,Mx:Integer;

Begin

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

D1^[i,j]:=Func(X[i],Y[j]);

End;

k:=0;

D2^:=D1^;

If Nx>Ny Then Mx:=Nx Else Mx:=Ny;

While k
Begin

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

If i>j Then s:=i Else s:=j;

If s=k Then B[i,j]:=D2^[i,j];

End;

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

D2^[i,j]:=Delta(k,i,j);

End;