unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls;
Function FUNC(X:real):real;
type
  TForm1 = class(TForm)
    Button1: TButton;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation
//PROGRAM D9R3
//Driver for routine BRENT
uses
  unit2;
  {$R *.DFM}
Function FUNC(X:real):real;
begin
    FUNC:= BESSJ0(X);
end;

procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%14.6f'; s2 = '#0';  TOL = 0.000001;  EQL = 0.0001;
var
  F:TextFile;
  AMIN:array[0..20] of real;
  I,J,NMIN,IFLAG:integer;  AX,BX,CX,FA,FB,FC,B,XMIN:real;
begin
  //輸出計(jì)算結(jié)果到文件
  AssignFile(F, 'd:\delphi_shu\p9\d9r3.dat');
  Rewrite(F);
  Writeln(F, 'Minima of the function BESSJ0');
  NMIN:= 0;
  Writeln(F, '   Min. #      X          BESSJ0(X)     BESSJ1(X)');
  For I:=1 To 100 do
  begin
    AX:= I;
    BX:=I + 1;
    MNBRAK(AX, BX, CX, FA, FB, FC);
    B:= BRENT(AX, BX, CX, TOL, XMIN);
    If NMIN = 0 Then
    begin
      AMIN[1]:= XMIN;
      NMIN:= 1;
      Writeln(F,'    ',FormatFloat(s2,NMIN),Format(S1,[XMIN]),
                    Format(S1,[BESSJ0(XMIN)]),Format(S1,[BESSJ1(XMIN)]));
    end
    Else
    begin
      IFLAG:= 0;
      For J:= 1 To NMIN do
        If Abs(XMIN - AMIN[J]) <= EQL * XMIN Then IFLAG:= 1;
      If IFLAG = 0 Then
      begin
        NMIN:= NMIN + 1;
        AMIN[NMIN]:= XMIN;
        Writeln(F,'    ',FormatFloat(s2,NMIN),Format(S1,[XMIN]),
                    Format(S1,[BESSJ0(XMIN)]),Format(S1,[BESSJ1(XMIN)]));
      end;
    end;
  end;
  CloseFile(F);
  //屏幕顯示計(jì)算結(jié)果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p9\d9r3.dat');
end;
end.