d2r11.txt

《Delphi常用数值算法集》的配书源码

implementation
//PROGRAM D2R11
//Driver for routine SPLIN2
uses
  unit2;
  {$R *.DFM}

procedure TForm1.Button1Click(Sender: TObject);
var
  F:TextFile;
  X1X2,XX1,XX2,F1,FF:real;  I,J:Integer;
  X1,X2:array[0..10] of real;
  Y,Y2:MATRX2;
const
  M = 10;
  N = 10;
  s1='%12.6f';
begin
  SetLength(Y,11,11);
  SetLength(Y2,11,11);
  For I:=1 To M DO
    X1[I]:= 0.2 * I;
  For I:=1 To N DO
    X2[I]:= 0.2 * I;
  For I:=1 To M DO
  begin
    For J:=1 To N DO
    begin
      X1X2:= X1[I] * X2[J];
      Y[I,J]:= X1X2 *Exp(-X1X2);
    end;
  end;
  SPLIE2(X1,X2,Y,M,N,Y2);
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p2\d2r11.dat');
  Rewrite(F);
  Writeln(F, '       X1          X2        SPLIN2      Actual');
  For I:= 1 To 10 do
  begin
    XX1:= 0.1 * I;
    XX2:= Sqr(XX1);
    SPLIN2(X1,X2,Y,Y2,M,N,XX1,XX2,F1);
    X1X2:= XX1 * XX2;
    FF:= X1X2 * Exp(-X1X2);
    Writeln(F,Format(s1,[XX1]),Format(s1,[XX2]),
                      Format(s1,[F1]),Format(s1,[FF]));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p2\d2r11.dat');
end;