unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls;
type
    matrx2=array of array of real;
type
  TForm1 = class(TForm)
    Button1: TButton;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

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

procedure TForm1.Button1Click(Sender: TObject);
var
  F:TextFile;
  X1X2: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]:= Sqr(X1X2);
    end;
  end;
  SPLIE2(X1,X2,Y,M,N,Y2);
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p2\d2r10.dat');
  Rewrite(F);
  Writeln(F, 'Second derivatives from SPLIE2');
  Writeln(F);
  Writeln(F, 'Natural spline assumed');
  For I:= 1 To 5 do
    Writeln(F,Format(s1,[Y2[I,1]]),Format(s1,[Y2[I,2]]),
      Format(s1,[Y2[I,3]]),Format(s1,[Y2[I,4]]),Format(s1,[Y2[I,5]]));
  Writeln(F);
  Writeln(F, 'Actual second derivatives');
  For I:= 1 To 5 do
  begin
    For J:= 1 To 5 do
      Y2[I,J]:= 2 * Sqr(X1[I]);
    Writeln(F,Format(s1,[Y2[I,1]]),Format(s1,[Y2[I,2]]),
        Format(s1,[Y2[I,3]]),Format(s1,[Y2[I,4]]),Format(s1,[Y2[I,5]]));
  end;

  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p2\d2r10.dat');
end;

end.