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;

var
  Form1: TForm1;

implementation
//PROGRAM D4R8
//Driver for routine GSER
uses
  unit2;
{$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
  ACTUAL:array[1..5] of real;X,A:array[1..5] of real;
  F:TextFile;  j:Integer;  GAMSER,GLN:real;
const
  s1='%6.2f';s2='%11.6f';
begin
  ACTUAL[1]:=0.911975;    ACTUAL[2]:=0.095163;
  ACTUAL[3]:=0.009105;    ACTUAL[4]:=0.038732;
  ACTUAL[5]:=0.940427;
  A[1]:=0.1;              A[2]:=1;
  A[3]:=2;                A[4]:=6;
  A[5]:=11;
  X[1]:=0.316228;         X[2]:=0.1;
  X[3]:=0.141421;         X[4]:=2.44949;
  X[5]:=16.58312;
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p4\d4r8.dat');
  Rewrite(F);
  Writeln(F, '   A         X       Actual    GSER(A,X)  GAMMLN(A)    GLN');
  For j:= 1 To 5 do
  begin
    GSER(GAMSER, A[j], X[j], GLN);
    Writeln(F,Format(s1,[A[j]]),Format(s2,[X[j]]),
                    Format(s2,[Actual[j]]),Format(s2,[GAMSER]),
                    Format(s2,[GAMMLN(A[j])]),Format(s2,[GLN]));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p4\d4r8.dat');
end;
end.