@@@@@@@@@@@@@@@@@@@@
C-8K FOCAL @1969

01.01 C***** DEMONSTRATION PROGRAM FOR SUBROUTINE 'NORAND' *****
01.02 C-----COMPUTES MEAN, STANDARD DEVIATION, SKEWNESS, AND KURTOSIS
01.03 C-----   FOR SAMPLES OF SIZE N AND PLOTS A BAR GRAPH OF THE
01.04 C-----   FREQUENCY DISTRIBUTION.
01.10 E
01.15 S IC=8/25
01.20 A !!"POPULATION STANDARD DEVIATION: "S,"  SAMPLE SIZE: "N
01.30 T !"  MEAN      SD      SKEW     KURT.
01.35 T !%5.04,0,%7.04,S,0,3,"  <<<  EXPECTED VALUES
01.40 S X1=0;S X2=0;S X3=0;S X4=0;F I=1,N;D 30;D 1.9;D 2
01.50 S M1=X1/N;S M2=X2/N;S M3=X3/N;S M4=X4/N
01.52 S U2=M2-M1^2;S U3=M3-3*M1*M2+2*M1^3
01.54 S U4=M4-4*M1*M3+6*M1^2*M2-3*M1^4
01.58 T !%5.04,M1,%7.04,FSQT(U2);S SK=1;I (0-U3)1.62,1.62,1.6
01.60 S SK=-1
01.62 T "  "%5.04,SK*FSQT(U3^2/U2^3),%7.04,U4/U2^2
01.64 T "  <<<  SAMPLE ESTIMATES
01.70 I (N-200)1.72;F I=1,25;S H(I)=H(I)*200/N
01.72 D 3
01.80 T !!!;QUIT
01.90 S X1=X1+X;S X2=X2+X^2;S X3=X3+X^3;S X4=X4+X^4

02.01 C----- FREQUENCY TABULATOR
02.10 S T=-4+IC;S X=X/S
02.20 F J=1,25;I (X-T)2.3;S T=T+IC
02.30 S H(J)=H(J)+1;S J=26

03.01 C------PLOT
03.05 T !!"..........FREQUENCY DISTRIBUTION
03.10 F J=1,25;T !".";I (H(J)-1)3.01;F I=1,H(J);T "X

30.01 C***** SUBROUTINE 'NORAND' *****
30.02 C-----GENERATES A NORMALLY DISTRIBUTED RANDOM NUMBER
30.03 C-----   FOR A GIVEN POPULATION STANDARD DEVIATION (S)
30.10 S X=100*FRAN( );S X=X-FITR(X);S Y=100*FRAN( );S Y=Y-FITR(Y)
30.20 S X=(X-.5)*10;I (1/FEXP(X*X/2)-Y)30.1;S X=X*S
*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@Hx~