implementation
//PROGRAM D2R9
//Driver for routine POLIN2
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
  X1A,X2A:array[0..5] of real;  YA:matrx2;
  X1,X2,FF,Y,DY:real;
  F:TextFile;  I,J:Integer;
const
  N=5; PI=3.1415926; s1='%12.6f'; 
begin
  //輸出計(jì)算結(jié)果到文件
  SetLength(YA,6,6);
  AssignFile(F, 'd:\delphi_shu\p2\d2r9.dat');
  Rewrite(F);
  For I:=1 To N do
  begin
      X1A[I]:=I * PI / N;
      For J:=1 To N do
      begin
          X2A[J]:=J / N;
          YA[I,J]:=Sin(X1A[I])*EXP(X2A[J]);
      end;
  end;
  //Test 2-dimensional interpolation
  Writeln(F,'      ');
  Writeln(F,'       x1          x2         f(x)    interpolated    error');
  For I:=1 To 4 do
  begin
      X1:=(-0.1 + I/5)*PI;
      For J:=1 To 4 do
      begin
          X2:=-0.1 + J / 5;
          FF:=Sin(X1) * EXP(X2);
          POLIN2(X1A,X2A,YA,N,N,X1,X2,Y,DY);
          Writeln(F,Format(s1,[X1]),Format(s1,[X2]),
               Format(s1,[FF]),Format(s1,[Y]),Format(s1,[DY]));
      end;
      Writeln(F,'*********************************************************');
  end;
  CloseFile(F);
  //屏幕顯示計(jì)算結(jié)果
  memo1.ScrollBars:=ssVertical;
  memo1.Lines.LoadFromFile('d:\delphi_shu\p2\d2r9.dat');
end;