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 D4R18
//Driver for routine BESSY1
uses
  unit2;
{$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
  ACTUAL:array[1..5] of real;X:array[1..5] of real;
  F:TextFile;  j:Integer;  
const
  s1='%8.2f';s2='%14.7f';
begin
  ACTUAL[1]:=-0.7812128;     ACTUAL[2]:=0.3246744;
  ACTUAL[3]:=0.1478631;      ACTUAL[4]:=0.2490154;
  ACTUAL[5]:=-0.2100814;
  X[1]:=1;   X[2]:=3;
  X[3]:=5;   X[4]:=10;
  X[5]:=13;
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p4\d4r18.dat');
  Rewrite(F);
  Writeln(F, '     X         Actual       BESSY1(X)');
  For j:= 1 To 5 do
  begin
    Writeln(F,Format(s1,[X[j]]),Format(s2,[Actual[j]]),
        Format(s2,[BESSY1(X[J])]));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p4\d4r18.dat');
end;

end.