implementation
//PROGRAM D11R8
//Driver for routine CHSONE
uses
  unit2;
  {$R *.DFM}
Function EXPDEV:real;
begin
    EXPDEV:= -Ln(Random);
end;
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%12.2f'; s2 = '.###0e+00';  NPTS = 2000;  NBINS = 10;
var
  F:TextFile;
  BINS,EBINS:array[0..NBINS] of real;
  I,J,IBIN:integer;  X,DF,CHSQ,PROB:real;
begin
  Randomize;
  For J:=1 To NBINS do
    BINS[J]:=0;
  For I:=1 To NPTS do
  begin
    X:=EXPDEV;
    IBIN:=Trunc(X * NBINS / 3) + 1;
    If IBIN <= NBINS Then BINS[IBIN]:=BINS[IBIN] + 1;
  end;
  For I:=1 To NBINS do
    EBINS[I]:=3 * NPTS / NBINS * Exp(-3 * (I - 0.5) / NBINS);
  CHSONE(BINS, EBINS, NBINS, -1, DF, CHSQ, PROB);
  //Generate Gaussian distributed data
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p11\d11r8.dat');
  Rewrite(F);
  Writeln(F);
  Writeln(F,'     Expected    Observed');
  For I:=1 To NBINS do
    Writeln(F,Format(s1,[EBINS[I]]),Format(s1,[BINS[I]]));
  Writeln(F);
  Writeln(F,'    Chi-squared: ', FormatFloat(s2,CHSQ));
  Writeln(F,'    Probability: ', FormatFloat(s2,PROB));
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p11\d11r8.dat');
end;