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;
  ISET:^integer; GSET:^real;
implementation
//PROGRAM D11R11
//Driver for routine KSTWO
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%16.6f'; S2='.###0E+00';  N1 = 2000;  N2 = 1000;  EPS = 0.1;
var
  F:TextFile;
  DATA1:array[0..N1] of real;  DATA2:array[0..N2] of real;
  I,J:integer;  Var1,FACTR,D,PROB:real;
begin
  New(ISET); New(GSET);
  Randomize;
  ISET^:=0;
  //Generate Gaussian distributed data
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p11\d11r11.dat');
  Rewrite(F);
  Writeln(F);
  Writeln(F,'     Variance Ratio    K-S Statistic    Probability');
  For J:= 1 To N1 do
    DATA1[J]:= GASDEV;
  For I:= 1 To 11 do
  begin
    Var1:= 1 + (I - 1) * EPS;
    FACTR:= Sqrt(Var1);
    For J:= 1 To N2 do
      DATA2[J]:= FACTR * GASDEV;
    KSTWO(DATA1, N1, DATA2, N2, D, PROB);
        Writeln(F,Format(s1,[Var1]),Format(s1,[D]),'        ',
                                           FormatFloat(s2,PROB));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p11\d11r11.dat');
end;

end.