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

MaxErr1:=0; VidnErr1:=0; MaxErr2:=0; VidnErr2:=0;

dx:=(Xb-Xa)/Cx; dy:=(Yb-Ya)/Cy;

For i:=0 To Cx DoFor j:=0 To Cy Do

Begin

p1:=Func(Xa+i*dx,Ya+j*dy);

p2:=Drib(Xa+i*dx,Ya+j*dy);

p3:=Polinom(Xa+i*dx,Ya+j*dy);

If Abs(p1-p3)>MaxErr1 Then

Begin

MaxErr1:=Abs(p1-p3); VidnErr1:=Abs((p1-p3)/p1);

end;

If Abs(p1-p2)>MaxErr2 Then

Begin

MaxErr2:=Abs(p1-p2); VidnErr2:=Abs((p1-p2)/p1);

End;

End;

Assign(f,'mix.txt'); Append(f);

WriteLn(f,nx:4,ny:4,MaxErr2:19:12,VidnErr2:19:12,MaxErr1:19:12,VidnErr1:19:12);

Close(f);

End;

Begin

For cc:=1 To 10 Do For cc1:=1 To 5 Do

Begin

DataInput; cx:=33; cy:=33;

WriteLn('Nx=',nx,' Ny=',ny);