implementation
//PROGRAM D11R5
//Driver for routine AVEVAR
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%9.2f';  NPTS = 1000;  EPS = 0.1;
var
  F:TextFile;   I,J:Integer;
  DATA:array[0..NPTS] of real;    AVE,SHIFT,VAR1:real;
begin
  New(ISET); New(GSET);
  Randomize;
  //Generate Gaussian distributed data
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p11\d11r5.dat');
  Rewrite(F);
  Writeln(F, '    Shift   Average   Variance');
  For I:= 1 To 11 do
  begin
    SHIFT:= (I - 1) * EPS;
    For J:= 1 To NPTS do
      DATA[J]:= SHIFT + I * GASDEV;
    AVEVAR(DATA, NPTS, AVE, VAR1);
    Writeln(F,Format(s1,[SHIFT]), Format(s1,[AVE]),Format(s1,[VAR1]));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p11\d11r5.dat');
end;