File STAT13.FC (FOCAL source file)

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

X()0N L$C-PS/8 FOCAL, 1971?MB	LJS X=X(100*BA+I); S Y=X(100*BE+I)?MB4LOI (X-SA)6.3,6.9,6.3?MCAL^I (Y-SA)6.4,6.9,6.4?MCIL(S N1=N1+1?MCXL2S X(N1)=X;S X(N1-100)=Y?MC\MZR?MC"P^D 4.32?MC/PS N1=0;F I=1,N; D 6?MDTP(F I=1,N1;S X=X(I);S Y=X(I-100);S X(I)=FABS(X-Y);
 X(I-100)=FSGN(X-Y)?MD^P+S VT=0;S VP=-1?MD(P-F I=1,N1;D 10?MD0P/D 15;D 14?MD=P2S Y(1)=0; S Y(-1)=0;?METP4F I=N0,N1;S Y(X(I-100))=Y(X(I-100))+X(I)?ME-P7S N1=N1-N0+1;I (Y(1)-Y(-1)) 8.6;S Y(1)=Y(-1)?MMVP<T %6.02;S T=Y(1);I (N1-25) 8.65,8.65,8.7?MMQFT !,  ," 
 ",(T-N1*(N1+1)/4)/FSQT((N1*(N1+1)*(2*N1+1))/24)?MF*TES J=N1?MF>TJS X1=X(100*VT+I);S X2=X(100*VT+J)?MGMTTI (X1-X2)10.4,10.3,10.3?MG!T^S X(100*VT+I)=X2;S X(100*VT+J)=X1?MHATS X(100*VP+I)=X(100*VP+I);S X(100*VP+J)=X(100*VP+J);D 10.6?MHST(S J=J-1; I (J-I)10.
,10.9,10.1?MH4T<S X1=X(100*VP+I);S X(100*VP+I)=X(100*VP+J);S X(100*VP+J)=X1?MH8UZR?MIF\FS N0=N0;S I=N0-1;S T=0?MIX\JS I=I+1;I (I-N1)14.15,14.15;R?MI7\OS X(VT*100+I+1)=X(VT*100+I+1);S X(VT*100+I)=X(VT*100+I)?MNO\PI (X(VT*100+I)-X(VT*100+I+1)) 14.17,14.2?MJ
\TS NT=0;S J=I; S M=0?MJ/\WS J=J+1;S X(VT*100+J)=X(VT*100+J);?MKG\YI (X(VT*100+J)-X(VT*100+I))14.3,14.23,14.3?MM<\^S NT=J-I;F K=I-N0+1,J-N0;;S M=M+K?MK/\(S I=J-1;S T=T+(NT^3-NT)/12;G 14.1?MK8^JS I=0;S J=0?MLK^TS J=J+1;I (J-N1)15.3;S N0=I+1;R?MLZ^^I (X(J))
5.2,15.35,15.2?M^S I=I+1;G 15.2?MMLDJA !!!!,"ANALYSIS ",LI;IF (LI-13)2.2,8.30,2.2?MBTH A "VARIABLES  ",BA,BE,!?ML=DTL R SORT 2.15?MFDQAT !," T ",Y(1);G 8.8?MF$QPT !,"NUMBER OF NONZERO DIFF=",%4.00,N1;G 2.1?MK[\F K=I,J-1;S X(VT*100+K)=M/(J-I)?MJN\QS X(I+VT
100)=I-N0+1;G 14.1?MX(L S STAT13?M?M)=I-N0+1;G 14.1?!?TE TO OUTPUT FILE 
",LI;I (LI-0YES)2.01,2.3,2.01?MDQDBSET U=10;G 1.96?MMEB I (U)1.5,1.5?MB1B<A IN(1);I (FABS(IN(1)-SL)-1E-4)1.82,1.82,1.65?MCBCFS N=N+1;F I=PECTED FREQ.",,%6.02?ML%T2F I=1,RA;T !;F J=1,KO;D 7;I (EX)10.6,10.6;T EX?ML.T<T !!!,%6.02?MNLVJT !;F J=1,KO;T N(10*I+J
?MBTBDE?MM>BLA "LIST OF EXP VAL. 
",LI?MB7BWD 10?MMQBVI (LI-0YES)1.24,1.23,1.24?MMABJA !, "NUMBER OF ROWS ",RA,"COLUMNS  ",KO,!?MB BM

Note: Lines longer than 256 characters were wrapped



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