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

nx:=cc; ny:=cc1*2-1;

{ For i:=0 To Nx Do X[i]:=(Xa+Xb)/2+(Xb-Xa)*Cos(Pi*i/Nx)/2;

For i:=0 To Ny Do Y[i]:=(Ya+Yb)/2+(Yb-Ya)*Cos(Pi*i/Ny)/2;}

For i:=0 To Nx Do X[i]:=Xa+(Xb-Xa)*i/Nx;

For i:=0 To Ny Do Y[i]:=Ya+(Yb-Ya)*i/Ny;

End;

Procedure BuildCoefTable;

Function Xij(i,j:Integer):Real;

Begin

If i>j Then Xij:=X[i]-X[j] Else Xij:=1;

End;

Function Yij(i,j:Integer):Real;

Begin

If i>j Then Yij:=Y[i]-Y[j] Else Yij:=1;

End;

Function Teta(t,s:Integer):Integer;

Begin

If s>t Then Teta:=-1 Else Teta:=0;

End;

Function Delta(k,i,j:Integer):Real;

Begin

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

( D1^[i,j]+

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

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