File KLMGRV.FC (FOCAL source file)

Directory of image this file is from
This file as a plain text file

C-PS/8 FOCAL, 1971
 1.10 T "KOLMOGOROV-SMIRNOV TEST FOR NORMALITY",!!
 1.15 A "NUMBER DATA POINTS ",NP;T !
 1.20 T "INPUT DATA POINT,SPACE,NO. OF OCCURENCES,SPACE",!!
 1.25 S I=1;S J=0
 1.30 A X(I),NX(I);S J=J+NX(I);S I=I+1
 1.35 I (J-NP)1.3,1.38,1.38
 1.38 T "WAIT"
 1.40 S ND=I-1;S SU=0.;F I=1,ND;D 3.5
 1.45 S MN=SU/NP;S SU=0.;F I=1,ND;D 3.6
 1.50 S SD=FSQT(SU/(NP-1.));S TP=10000;S Q=-TP;S SX=0
 1.53 T !!,"    X(T)      F(I)      FN(I)      ARG"
 1.54 T "      F(ARG)  F(ARG)-FN(I)",!!
 1.55 F IJ=1,ND;D 5
 1.80 T "PICK LARGEST NO. IN RH COLUMN",!
 1.85 T "DF = ",%2,NP,!
 1.90 Q
 3.50 S SU=SU+NX(I)*X(I);R
 3.60 S SU=SU+NX(I)*((X(I)-MN)^2);R
 4.02 S KEY=10.
 4.03 I (B)4.05,4.05,4.04
 4.04 S B=-B;S KEY=-10.
 4.05 S C=B+6.;S D=B-6.
 4.10 S A1=.238619*C+D;S A2=-.238619*C+D
 4.20 S A3=.661209*C+D;S A4=-.661209*C+D
 4.30 S A5=.93247*C+D;S A6=-.93247*C+D
 4.35 S A1=FEXP(-A1*A1/8.);S A2=FEXP(-A2*A2/8.)
 4.40 S A3=FEXP(-A3*A3/8.);S A4=FEXP(-A4*A4/8.)
 4.60 S A5=FEXP(-A5*A5/8.);S A6=FEXP(-A6*A6/8.)
 4.70 S XX=.467914*(A1+A2)+.360762*(A3+A4)+.171324*(A5+A6)
 4.80 S XX=.19947*C*XX-.00027
 4.85 I (KEY)4.86,4.9,4.9
 4.86 S XX=1.-XX
 4.90 R
 5.01 F I=1,ND;D 6
 5.20 S Q=TP;S B=(TP-MN)/SD;D 4
 5.30 S SX=SX+NX(KL)
 5.40 S NF=SX/NP;S DF=FABS(XX-NF)
 5.43 T %7.05,TP,NX(KL),NF,B,XX,DF,!
 5.45 S TP=100000
 5.90 R
 6.10 I (Q-X(I))6.2,6.9,6.9
 6.20 I (TP-X(I))6.9,6.9,6.3
 6.30 S TP=X(I);S KL=I
 6.90 R



Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with.  I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation. 

PDP-8 Home Page   PDP-8 Site Map   PDP-8 Site Search