d4r28.txt

《Delphi常用数值算法集》的配书源码

implementation
//PROGRAM D4R28
//Driver for routine BESSIAN
uses
  unit2;
{$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
  ACTUAL:array[0..4] of real; FF:array[0..4] of real;
  F:TextFile;  j:Integer;
  A,X:real;IH:integer;
const
  s1='%6.2f';s2='.000000e+00';S3='###0';
begin
    ACTUAL[0]:=0.359208;       ACTUAL[1]:=0.0238836;
    ACTUAL[2]:=0.000954254;    ACTUAL[3]:=0.0000272471;
    A:=0.5;
    IH:=1;
    X:=0.2;
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p4\d4r28.dat');
  Rewrite(F);
  Writeln(F, 'N   AX   IH    X       Actual        BESSIAN');
  BESIAN(X,A,4,IH,FF);
  For j:= 0 To 3 do
  begin
      Writeln(F,FormatFloat(s3,j),Format(s1,[A]),'   ',
          FormatFloat(s3,IH),Format(s1,[X]),'    ',
          FormatFloat(s2,Actual[j]),'    ',FormatFloat(s2,FF[J]));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p4\d4r28.dat');
end;