C-PS/8 FOCAL, 1971
 2.10 A !!!!,"ANALYSIS ",LI;IF (LI-11)2.2,4.14,2.2
 2.20 L R SORT 2.15
 4.10 I (X^2-.01)4.12;S X=X/2; D 4; S X=X^2;R
 4.12 S X=1+X+X^2/2+X^3/6+X^4/24+X^5/120+X^6/720;R
 4.14 A "VARIABLE ",BA,!,"FOR THE NORMAL POPULATION STATE MEAN ",
 4.15 A ME,"  STANDARD DEVIATION ",SD,!;S N1=0;S VT=0;F I=1,N;D 11
 4.16 S VP=-1;S DF=0
 4.18 F I=1,N1;D 10
 4.20 S MA=0;F I=1,N1;D 7
 4.22 T %8.03
 4.24 T !,"MEDIAN"
 4.26 S X=(X(FITR((N1+1)/2))+X(FITR((N1+2)/2)))/2
 4.28 T X,!,
MAX.DEVIATION
,!,"NUMBER ",%4.00,N1;G 2.1
 7.10 S Z=(X(I)-ME)/SD; S X1=FABS(Z)
 7.15 S T=1/(1+0.2316419*X1)
 7.20 S X=-X1*X1/2; D 4
 7.30 S P=0.3989423*X*T
 7.35 S X=T*(1.78147794+T*(-1.82125598+T*1.33027443))
 7.37 S P=1-P*(0.31938153+T*(-0.35656378+X))
 7.40 I (Z)7.6,7.6,7.65
 7.60 S P=1-P
 7.65 I (MA-FABS(P-I/N1))7.7,7.7,7.9
 7.70 S MA=FABS(P-I/N1)
 7.90 C
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.40 S J=J-1; I (J-I)10.9,10.9,10.1
10.90 R
11.10 I (X(BA*100+I)-SA)11.2,11.9,11.2
11.20 S N1=N1+1;S X(100*VT+N1)=X(BA*100+I)
11.90 R