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

GetG:=(xx-X[k-1])*(yy-Y[k-1])/

(B[k,k]+GetH(k+1,k)+GetL(k+1,k)+GetG(k+1));

End;

Begin

If Nx
Drib:=B[0,0]+GetH(1,0)+GetL(1,0)+GetG(1);

End;

Function Polinom(xx,yy:Real):Real;

Var p,q,s,s1,p1,q1:Real; i,j,k:Integer;

Begin

s:=0;

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

p:=1; q:=1;

For k:=0 To Nx Do If k<>i Then p:=p*(xx-X[k])/(X[i]-X[k]);

For k:=0 To Ny Do If k<>j Then q:=q*(yy-Y[k])/(Y[j]-Y[k]);

s1:=p*q*Func(X[i],Y[j]);

s:=s+s1;

End;

Polinom:=s;

End;

Procedure GetMaxError;

Var i,j:Integer; dx,dy,MaxErr1,p1,p2,p3,VidnErr1,MaxErr2,VidnErr2:Real; F:Text;

Begin