home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / pascal / library / dos / nrpas / mdian2.dem < prev    next >
Encoding:
Text File  |  1994-04-11  |  849 b   |  42 lines
  1. PROGRAM d13r2(input,output);
  2. (* driver for routine MDIAN2 *)
  3. CONST
  4.    npts=50;
  5. TYPE
  6.    glsarray = ARRAY [1..npts] OF real;
  7.    narray = glsarray;
  8. VAR
  9.    glinext,glinextp : integer;
  10.    glma : ARRAY [1..55] OF real;
  11.    gliset : integer;
  12.    glgset : real;
  13.    i,idum : integer;
  14.    xmed : real;
  15.    data : narray;
  16.  
  17. (*$I MODFILE.PAS *)
  18. (*$I RAN3.PAS *)
  19.  
  20. (*$I GASDEV.PAS *)
  21.  
  22. (*$I SORT.PAS *)
  23.  
  24. (*$I MDIAN1.PAS *)
  25.  
  26. (*$I MDIAN2.PAS *)
  27.  
  28. BEGIN
  29.    gliset := 0;
  30.    idum := -5;
  31.    FOR i := 1 to npts DO BEGIN
  32.       data[i] := gasdev(idum)
  33.    END;
  34.    mdian2(data,npts,xmed);
  35.    writeln('Data drawn from a gaussian distribution');
  36.    writeln('with zero mean, unit variance');
  37.    writeln;
  38.    writeln('median according to mdian2 is',xmed:9:6);
  39.    mdian1(data,npts,xmed);
  40.    writeln('median according to mdian1 is',xmed:9:6)
  41. END.
  42.