implementation
//PROGRAM D11R1
//Driver for routine MOMENT
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%11.4f';
  NPTS = 10000; NBIN = 100;  NDAT = NPTS + NBIN;  PIO2=1.5707963;
var
  F:TextFile;   I,J,K,NLIM:Integer;
  X,AVE,ADEV,SDEV,VAR1,SKEW,CURT:real;
  DATA:array[0..NDAT] of real;
begin
  I:=1;
  For J:= 1 To NBIN do
  begin
    X:= PI * J / NBIN;
    NLIM:= Round(Sin(X) * PI / 2 * NPTS / NBIN);
    For K:= 1 To NLIM do
    begin
      DATA[I]:= X;
      I:= I + 1;
    end;
  end;
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p11\d11r1.dat');
  Rewrite(F);
  Writeln(F, 'Moments of a sinusoidal distribution');
  MOMENT(DATA, NPTS, AVE, ADEV, SDEV, VAR1, SKEW, CURT);
  Writeln(F);
  Writeln(F, '                       Calculated  Expected');
  Writeln(F, 'Mean :              ',Format(s1,[AVE]),Format(s1,[PI / 2]));
  Writeln(F, 'Average Deviation : ',Format(s1,[ADEV]),'     0.5708');
  Writeln(F, 'tandard Deviation : ',Format(s1,[SDEV]),'     0.6837');
  Writeln(F, 'Variance :          ',Format(s1,[VAR1]),'     0.4674');
  Writeln(F, 'Skewness :          ',Format(s1,[SKEW]),'     0.0000');
  Writeln(F, 'Kurtosis :          ',Format(s1,[CURT]),'    -0.8062');
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p11\d11r1.dat');
end;