unit1.pas

Numeric Programs

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 D4R20
//Driver for routine BESSJ
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.114903;     ACTUAL[2]:=0.25463;
  ACTUAL[3]:=-0.234062;    ACTUAL[4]:=0.263062E-9;
  ACTUAL[5]:=0.207486;
  N[1]:=2;   N[2]:=2;
  N[3]:=5;   N[4]:=10;
  N[5]:=10;
  X[1]:=1;   X[2]:=10;
  X[3]:=10;  X[4]:=1;
  X[5]:=10;
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p4\d4r20.dat');
  Rewrite(F);
  Writeln(F, ' N      X        Actual       BESSJ(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,BESSJ(N[J],X[J])));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p4\d4r20.dat');
end;

end.