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

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);}

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