C-PS/8 FOCAL, 1971
 1.01 C     FROM CHISQR BY M.J. MC KEOWN (DECUS: FOCAL8-40)
 1.02 C     MODIFIED BY JOHN HUISMAN
 1.03 C     OMSI SOFTWARE DEVELOPMENT GROUP
 1.04 C     LAST CHANGE: 8/19/71
 1.05 
 1.06 C     THE PROGRAM WILL PRINT THE DATA MATRIX WITH TWO 
 1.07 C VALUES IN EACH CELL.  ONE IS THE OBSERVED VALUE AND
 1.08 C THE OTHER IS THE EXPECTED VALUE AS CALCULATED.  YOU MAY HAVE UP TO
 1.09 C 10 DIGIT PRECISION, BUT IT IS VARIABLE.  YOU MAY ENTER A
 1.10 C MATRIX NO LARGER THAN 25 ELEMENTS.  THE WIDTH THAT THE MATRIX
 1.11 C MAY BE DEPENDS ON THE PRECISION YOU CHOOSE, AND THE WIDTH
 1.12 C OF YOUR OUTPUT DEVICE.
 1.13 
 1.14 C ** RS IS THE ROW SUM.
 1.15 C ** CS IS THE COLUMN SUM.
 1.16 
 1.17 I (FITR(/5)-(/5)) 1.18,1.3,1.18
 1.18 T "YOU DO NOT HAVE ENOUGH FEATURES.  "
 1.19 T "RESTART FOCAL WITH FEATURE 3.",!;Q
 1.20 
 1.30 S X=FOUT(24)+FOUT(1);F K=1,100;E
 1.31 T "HOW MANY PLACES ACCURATE
";A AC;S X1=(100*AC+3)/100;T %(X1)
 1.32 A "ROWS:"R,"COLUMNS:"C;S DF=(R-1)*(C-1);S L=R*C
 1.33 S K=0;F I=1,R;S RS(I-1)=RP;S RP=0;D 1.4
 1.34 S RS(I-1)=RP
 1.35 T %(X1),!"NO MORE INPUT!"!;G 2.01
 1.40 F J=1,C;D 1.41
 1.41 S K=K+1;T !,%2 "SAMPLE"I,","J,": ";G 1.42
 1.42 A N(K);I (10^AC-N(K)) 1.43,1.43;S RP=RP+N(K);S CS(J)=CS(J)+N(K)
 1.43 T "TOO MANY DIGITS - INPUT AGAIN OR INCREASE ACCURACY",!,"
";G 1.42
 2.01 F J=1,1,R;S T=T+RS(J)
 2.02 S J=1;S I=1;S K=1
 2.03 S E(K)=(RS(I)/T)*CS(J);S J=J+1;S K=K+1
 2.04 I (C-J)2.05,2.03,2.03
 2.05 S J=1;S I=I+1
 2.06 I (R-I)3.01,2.03,2.03
 3.01 F K=1,1,L;S X2=X2+((N(K)-E(K))^2)/E(K)
 4.01 T !!;S K=1;S Y=1;S N=1;S I=1;D 10;D 8.0;T !
 4.02 F J=1,1,C;D 6.0
 4.03 T "*",!,"*";F J=1,C;D 5;T "*"
 4.04 T "*  ","RS:"RS(N);S N=N+1;T !
 4.05 F J=1,1,C;D 7.0
 4.06 T "*";T !;D 8.0;T !;S I=I+1;I (R-I)4.07,4.02,4.02
 4.07 T !;F J=1,1,C;D 9.0
 4.08 T %,"  TOTAL:"T,!!,"CHI-SQUARED:"X2
 4.09 T %(X1),"   DEGREES OF FREEDOM:"DF,!!
 4.10 A "MORE
"Z;I (Z-0YES) 4.11,1.3,4.11
 4.11 Q
 5.01 T :(1+J*(9+AC))
 6.01 T "*  O:"N(K),;D 5.01;S K=K+1
 7.01 T "*  E:"E(Y),;D 5.01;S Y=Y+1
 8.01 F M=1,(C*(AC+9)+1);T "*"
 9.01 T " CS:"CS(J)," "
10.01 T "E = EXPECTED VALUE     O = OBSERVED VALUE",!!