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

C-PS/8 FOCAL, 1971
 2.10 A !!!!,"ANALYSIS ",LI;IF (LI-14)2.2,5.10,2.2
 2.20 L R SORT 2.15
 4.32 A "VARIABLES  ",BA,BE,!;I (BA-BE)4.33,4.34,4.34
 4.33 R
 4.34 T !,"FIRST MUST BE SMALLER",!;G 4.32
 5.05 G 5.5
 5.10 D 4;S AN=2;S IN(1)=BA;S IN(2)=BE;G 5.14
 5.12 S AN=KB;T !;F I=1,KB;S IN(I)=I
 5.14 S N1=0;F I=1,AN;F J=1,N;D 5
 5.16 S VT=0;S VP=-1
 5.18 F I=1,N1;D 10
 5.20 S N0=1;D 14
 5.22 F I=1,KB;S N(I)=0;S IN(I)=0
 5.24 F I=N0,N1;S J=X(I-100);S N(J)=N(J)+1;S IN(J)=IN(J)+X(I)
 5.26 I (AN-2)5.58,5.58,5.28
 5.28 S DF=0;F I=1,AN;S DF=DF+(IN(I)^2)/N(I)
 5.30 S DF=((12/(N1*(N1+1)))*DF-3*(N1+1))/(1-T/(N1^3-N1))
 5.32 I (N1-15)5.34,5.34,5.38
 5.34 T !,%8.03,"H",DF,%4.00;F I=1,AN;T N(I)
 5.36 I (AN-3)5.5,2.1;T !,"TOO FEW D.F.";G 2.1
 5.38 T %8.03,"H",DF,"  DF  ",%4.00,AN-1;G 2.1
 5.50 S X=X(IN(I)*100+J)
 5.52 I (X)5.54,5.56,5.54
 5.54 S N1=N1+1;S X(N1)=X;S X(N1-100)=IN(I)
 5.56 R
 5.58 S N(1)=N(BA);S N(2)=N(BE);S DF=N(1)*N(2)+N(1)*(N(1)+1)/2-IN(BA)
 5.60 I (2*DF-N(1)*N(2))5.62,5.62;S DF=N(1)*N(2)-DF
 5.62 I (20-N(1))5.66;I (20-N(2))5.66
 5.64 T !,%4.00,"U",DF,
  N(1)  
,
  N(2)  
;G 2.1
 5.66 S I=(DF-N(1)*N(2)/2)/FSQT(N(1)*N(2)*(N(1)+N(2)+1)/12)
 5.68 T %8.03,!,"Z UNCORRECTED ",I," Z CORRECTED
 5.70 S DF=(DF-N(1)*N(2)/2)/FSQT(N(1)*N(2)/(N1*(N1-1)))
 5.72 T DF/FSQT((N1^3-N1)/12-T);G 2.1
10.05 S J=N1
10.10 S X1=X(100*VT+I);S X2=X(100*VT+J)
10.20 I (X1-X2)10.4,10.3,10.3
10.30 S X(100*VT+I)=X2;S X(100*VT+J)=X1
10.35 S X(100*VP+I)=X(100*VP+I);S X(100*VP+J)=X(100*VP+J);D 10.6
10.40 S J=J-1; I (J-I)10.9,10.9,10.1
10.60 S X1=X(100*VP+I);S X(100*VP+I)=X(100*VP+J);S X(100*VP+J)=X1
10.90 R
14.06 S N0=N0;S I=N0-1;S T=0
14.10 S I=I+1;I (I-N1)14.15,14.15;R
14.15 S X(VT*100+I+1)=X(VT*100+I+1);S X(VT*100+I)=X(VT*100+I)
14.16 I (X(VT*100+I)-X(VT*100+I+1)) 14.17,14.2
14.17 S X(I+VT*100)=I-N0+1;G 14.1
14.20 S NT=0;S J=I; S M=0
14.23 S J=J+1;S X(VT*100+J)=X(VT*100+J);
14.25 I (X(VT*100+J)-X(VT*100+I))14.3,14.23,14.3
14.30 S NT=J-I;F K=I-N0+1,J-N0;;S M=M+K
14.35 F K=I,J-1;S X(VT*100+K)=M/(J-I)
14.40 S I=J-1;S T=T+(NT^3-NT)/12;G 14.1