unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls;

type
  TForm1 = class(TForm)
    Button1: TButton;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
type
  matrx2=array of array of real;
var
  Form1: TForm1;

implementation
//PROGRAM D13R2
//Driver for routine ADI
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%7.2f';   JMAX =11;  PI =3.1415926;
var
  F:TextFile;   A,B,C,D,E,FF,G,U:matrx2;   STR1:String;
  I,J,K,MIDL:integer;    ALPHA,BETA,EPS,ALIM,AAA:real;
begin
  SetLength(A,12,12);
  SetLength(B,12,12);
  SetLength(C,12,12);
  SetLength(D,12,12);
  SetLength(E,12,12);
  SetLength(FF,12,12);
  SetLength(G,12,12);
  SetLength(U,12,12);
  For I:=1 To JMAX do
  begin
    For J:=1 To JMAX do
    begin
      A[I, J]:=-1;
      B[I, J]:=2;
      C[I, J]:=-1;
      D[I, J]:=-1;
      E[I, J]:=2;
      FF[I, J]:=-1;
      G[I,J]:=0;
      U[I, J]:=0;
    end;
  end;
  MIDL:=JMAX div 2 + 1;
  G[MIDL, MIDL]:=2;
  ALPHA:= 2 * (1 - Cos(PI / JMAX));
  BETA:= 2 * (1 - Cos((JMAX - 1) * PI / JMAX));
  ALIM:= Ln(4 * JMAX / PI);
  K:= 0;
  repeat
    K:= K + 1;
  until Exp(K*Ln(2)) > ALIM;
  EPS:= 0.0001;
  ADI(A, B, C, D, E, FF, G, U, JMAX, K, ALPHA, BETA, EPS);
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p13\d13r2.dat');
  Rewrite(F);
  Writeln(F);
  Writeln(F,'ADI Solution:');
  Writeln(F);
  For I:=1 To JMAX do
  begin
    STR1:=' ';
    For J:=1 To JMAX do
      STR1:=STR1 + '  '+ Copy(FloatToStr(U[I,J]),1,5);
    Writeln(F,STR1);
  end;
  Writeln(F);
  Writeln(F,'Test that sulotion satisfies Difference Eqns:');
  Writeln(F);
  For I:=2 To JMAX - 1 do
  begin
    For J:=2 To JMAX - 1 do
    begin
      AAA:=U[I + 1, J] + U[I - 1, J] + U[I, J + 1] + U[I, J - 1];
      G[I, J]:= -4 * U[I, J] + AAA;
    end;
    Writeln(F,Format(s1,[G[I,2]]),Format(s1,[G[I,3]]),
      Format(s1,[G[I,4]]),Format(s1,[G[I,5]]),Format(s1,[G[I,6]]),
      Format(s1,[G[I,7]]),Format(s1,[G[I,8]]),Format(s1,[G[I,9]]),
      Format(s1,[G[I,10]]));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p13\d13r2.dat');
end;

end.