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 D4R19
//Driver for routine BESSY
uses
  unit2;
{$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
  ACTUAL:array[1..5] of real;X:array[1..5] of real;
  N:array[1..5] of Integer;
  F:TextFile;  j:Integer;
const
  s1='%8.2f';s2='.#####0E+00';s3='###';
begin
  ACTUAL[1]:=-1.65068;     ACTUAL[2]:=-0.005868084;
  ACTUAL[3]:=-0.453695;    ACTUAL[4]:=0.135403;
  ACTUAL[5]:=-0.359814;
  N[1]:=2;   N[2]:=2;
  N[3]:=5;   N[4]:=5;
  N[5]:=10;
  X[1]:=1;   X[2]:=10;
  X[3]:=5;   X[4]:=10;
  X[5]:=10;
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p4\d4r19.dat');
  Rewrite(F);
  Writeln(F, ' N      X         Actual         BESSY(N,X)');
  For j:= 1 To 5 do
  begin
    Writeln(F,Formatfloat(s3,N[j]),Format(s1,[X[j]]),'    ',
      FormatFloat(s2,Actual[j]),'    ',FormatFloat(s2,BESSY(N[J],X[J])));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p4\d4r19.dat');
end;

end.