home *** CD-ROM | disk | FTP | other *** search
/ C!T ROM 2 / ctrom_ii_b.zip / ctrom_ii_b / PROGRAM / PASCAL / NRPAS13 / AVEVAR.DEM < prev    next >
Text File  |  1991-04-29  |  776b  |  37 lines
  1. PROGRAM d13r4(input,output);
  2. (* driver for routine AVEVAR *)
  3. CONST
  4.    npts=1000;
  5.    eps=0.1;
  6. TYPE
  7.    glnparray = ARRAY [1..npts] OF real;
  8. VAR
  9.    gliset : integer;
  10.    glgset : real;
  11.    glinext,glinextp : integer;
  12.    glma : ARRAY [1..55] OF real;
  13.    i,idum,j : integer;
  14.    ave,shift,vrnce : real;
  15.    data : glnparray;
  16.  
  17. (*$I MODFILE.PAS *)
  18. (*$I RAN3.PAS *)
  19.  
  20. (*$I GASDEV.PAS *)
  21.  
  22. (*$I AVEVAR.PAS *)
  23.  
  24. BEGIN
  25. (* generate gaussian distributed data *)
  26.    gliset := 0;
  27.    idum := -5;
  28.    writeln('shift':9,'average':11,'variance':12);
  29.    FOR i := 1 to 11 DO BEGIN
  30.       shift := (i-1)*eps;
  31.       FOR j := 1 to npts DO
  32.          data[j] := shift+i*gasdev(idum);
  33.       avevar(data,npts,ave,vrnce);
  34.       writeln(shift:8:2,ave:11:2,vrnce:12:2)
  35.    END
  36. END.
  37.