home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / PROG / PASCAL / NRPAS13.ZIP / MOMENT.DEM < prev    next >
Encoding:
Text File  |  1991-04-29  |  1.1 KB  |  41 lines
  1. PROGRAM d13r0(input,output);
  2. (* driver for routine MOMENT *)
  3. CONST
  4.    pi=3.14159265;
  5.    npts=5000;
  6.    nbin=100;
  7.    nppnb=5100;   (* nppnb=npts+nbin *)
  8. TYPE
  9.    narray = ARRAY [1..nppnb] OF real;
  10. VAR
  11.    adev,ave,curt,sdev,skew : real;
  12.    vrnce,x : real;
  13.    i,j,k,nlim : integer;
  14.    data : narray;
  15.  
  16. (*$I MODFILE.PAS *)
  17. (*$I MOMENT.PAS *)
  18.  
  19. BEGIN
  20.    i := 1;
  21.    FOR j := 1 to nbin DO BEGIN
  22.       x := pi*j/nbin;
  23.       nlim := round(sin(x)*pi/2.0*npts/nbin);
  24.       FOR k := 1 to nlim DO BEGIN
  25.          data[i] := x;
  26.          i := i+1
  27.       END
  28.    END;
  29.    writeln('moments of a sinusoidal distribution');
  30.    writeln;
  31.    moment(data,npts,ave,adev,sdev,vrnce,skew,curt);
  32.    writeln('calculated':39,'expected':11);
  33.    writeln;
  34.    writeln('Mean :',' ':19,ave:12:4,pi/2.0:12:4);
  35.    writeln('Average Deviation :',' ':6,adev:12:4,(pi/2.0)-1.0:12:4);
  36.    writeln('Standard Deviation :',' ':5,sdev:12:4,0.683667:12:4);
  37.    writeln('Variance :',' ':15,vrnce:12:4,0.467401:12:4);
  38.    writeln('Skewness :',' ':15,skew:12:4,0.0:12:4);
  39.    writeln('Kurtosis :',' ':15,curt:12:4,-0.806249:12:4)
  40. END.
  41.