implementation
//PROGRAM D11R2
//Driver for routine MDIAN1
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%12.6f';    NPTS = 50;
var
  F:TextFile;   I,J:Integer;
  DATA:array[0..50] of real;    XMED:real;
begin
  New(ISET); New(GSET);
  Randomize;
  For I:= 1 To NPTS do
    DATA[I]:= GASDEV;
  MDIAN1(DATA, NPTS, XMED);
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p11\d11r2.dat');
  Rewrite(F);
  Writeln(F, 'Gaussian distrib., zero mean, unit variance');
  Writeln(F);
  Writeln(F, 'Median of data set is   ', Format('%8.6f',[XMED]));
  Writeln(F);
  Writeln(F, 'Sorted data');
  Writeln(F);
  For I:= 1 To 10 do
  begin
    J:= (I - 1) * 5;
    Writeln(F, Format(s1,[DATA[J + 1]]),Format(s1,[DATA[J + 2]]),
     Format(s1,[DATA[J+3]]),Format(s1,[DATA[J+4]]),Format(s1,[DATA[J+5]]));
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p11\d11r2.dat');
end;