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 D4R26
//Driver for routine BESSI
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='.000000e+00';S3='####';
begin
  ACTUAL[1]:=0.00501669;      ACTUAL[2]:=0.0221684;
  ACTUAL[3]:=10.3312;         ACTUAL[4]:=0.00982568;
  ACTUAL[5]:=0.00458;
  N[1]:=2;     N[2]:=3;
  N[3]:=3;     N[4]:=5;
  N[5]:=10;
  X[1]:=0.2;   X[2]:=1;
  X[3]:=5;     X[4]:=2;
  X[5]:=5;
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p4\d4r26.dat');
  Rewrite(F);
  Writeln(F, '   N       X       Actual       BESSI(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,BESSI(N[j],X[J])));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p4\d4r26.dat');
end;

end.