home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / KSTWO.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.0 KB  |  50 lines
  1. PROGRAM d13r11(input,output);
  2. (* driver for routine KSTWO *)
  3. CONST
  4.    n1=500;
  5.    n2=100;
  6.    eps=0.2;
  7. TYPE
  8.    glarray1 = ARRAY [1..n1] OF real;
  9.    glarray2 = ARRAY [1..n2] OF real;
  10.    glsarray = glarray1;
  11. VAR
  12.    glinext,glinextp : integer;
  13.    glma : ARRAY [1..55] OF real;
  14.    gliset : integer;
  15.    glgset : real;
  16.    i,idum,j : integer;
  17.    d,factr,prob,varnce : real;
  18.    data1 : glarray1;
  19.    data2 : glarray2;
  20.  
  21. (*$I MODFILE.PAS *)
  22. (*$I PROBKS.PAS *)
  23.  
  24. (*$I SORT.PAS *)
  25.  
  26. (*$I RAN3.PAS *)
  27.  
  28. (*$I GASDEV.PAS *)
  29.  
  30. (*$I KSTWO.PAS *)
  31.  
  32. BEGIN
  33.    gliset := 0;
  34.    idum := -1357;
  35.    FOR i := 1 to n1 DO BEGIN
  36.       data1[i] := gasdev(idum)
  37.    END;
  38.    writeln('variance ratio':18,'k-s statistic':15,'probability':14);
  39.    idum := -2468;
  40.    FOR i := 1 to 11 DO BEGIN
  41.       varnce := 1.0+(i-1)*eps;
  42.       factr := sqrt(varnce);
  43.       FOR j := 1 to n2 DO BEGIN
  44.          data2[j] := factr*gasdev(idum)
  45.       END;
  46.       kstwo(data1,n1,data2,n2,d,prob);
  47.       writeln(varnce:15:6,d:15:6,prob:15:6)
  48.    END
  49. END.
  50.