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-13)2.2,8.30,2.2
 2.20 L R SORT 2.15
 4.32 A "VARIABLES  ",BA,BE,!
 6.10 S X=X(100*BA+I); S Y=X(100*BE+I)
 6.15 I (X-SA)6.3,6.9,6.3
 6.30 I (Y-SA)6.4,6.9,6.4
 6.40 S N1=N1+1
 6.50 S X(N1)=X;S X(N1-100)=Y
 6.90 R
 8.30 D 4.32
 8.35 S N1=0;F I=1,N; D 6
 8.40 F I=1,N1;S X=X(I);S Y=X(I-100);S X(I)=FABS(X-Y);S X(I-100)=FSGN(X-Y)
 8.43 S VT=0;S VP=-1
 8.45 F I=1,N1;D 10
 8.47 D 15;D 14
 8.50 S Y(1)=0; S Y(-1)=0;
 8.52 F I=N0,N1;S Y(X(I-100))=Y(X(I-100))+X(I)
 8.55 S N1=N1-N0+1;I (Y(1)-Y(-1)) 8.6;S Y(1)=Y(-1)
 8.60 T %6.02;S T=Y(1);I (N1-25) 8.65,8.65,8.7
 8.65 T !," T ",Y(1);G 8.8
 8.70 T !,  ," Z ",(T-N1*(N1+1)/4)/FSQT((N1*(N1+1)*(2*N1+1))/24)
 8.80 T !,"NUMBER OF NONZERO DIFF=",%4.00,N1;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
15.10 S I=0;S J=0
15.20 S J=J+1;I (J-N1)15.3;S N0=I+1;R
15.30 I (X(J))15.2,15.35,15.2
15.35 S I=I+1;G 15.2