99 FILES FILO,FILO
100 PRINT "&&&L100&&&STARTING"
105 SCRATCH #2
110 READ D$,O$
111 WRITE #2,D$
112 WRITE #2,O$
120 READ S1,S2,S3,N,F1,F2
121 WRITE #2,S1
122 WRITE #2,S2
123 WRITE #2,S3
124 WRITE #2,N
125 WRITE #2,F1
126 WRITE #2,F2
130 LET S4=S2+S3
140 LET S5=S4+S1
150 READ M,V1
160 FOR I=1 TO S4
170 READ S(I)
171 WRITE #2,S(I)
180 NEXT I
190 FOR I=1 TO S4
200 READ A$(I)
201 WRITE #2,A$(I)
210 NEXT I
220 FOR I=1 TON 
230 READ L(I),R(I)
231 WRITE #2,L(I)
232 WRITE #2,R(I)
240NEXT I
250 READ T1,T2,T3,T4
251 WRITE #2,T1
252 WRITE #2,T2
253 WRITE #2,T3
254 WRITE #2,T4
260 LET E2=0
262 LET E3=0
270 LET E2=FNT(T2)-FNT(T1)
280 IF S1=1 THEN 300
290 LET E3=FNT(T4)-FNT(T3)
300 LET S6=S4+2
310 FOR I=1 TO N
320 READ P(I)
321 WRITE #2,P(I)
330 NEXT I
400 GOSUB 1000
500 GOSUB 2000
600 GOTO 99999
1000 PRINT "THIS STATEMENT SHOULD BE FILES FILE1,FILE1"
1010 DIM A$(18),S(18),L(100),P(100),R(100),Y(100)
1020 DIM C(20,20),Z(20,20)   
1030 DIM Q(20),X(20)
1033 DIM D(100),N(100)
1038 RESTORE #1
1040 READ #1,D$,O$
1050 READ #1,S1,S2,S3,N,F1,F2
1060 LET S4 = S2 + S3
1070 LET S5=S4+1
1080 READ #1,M,V1
1090 FOR I = 1 TO S4 
1100 READ #1,S(I)
1110 NEXT I
1120 FOR I = 1 TO S4 
1130 READ #1,A$(I)
1140 NEXT I
1150 FOR I = 1 TO N
1160 READ #1,L(I),R(I)
1170 NEXT I
1180 READ #1,T1,T2,T3,T4
1190 LET E2=E3=0   
1200 LET E2 = FNT(T2)-FNT(T1)
1210 IF S1=1 THEN 1230
1220 LET E3 = FNT(T4)-FNT(T3)
1230 LET S6=S4+2   
1240 FOR I = 1 TO N
1250 READ #1,P(I)
1260 NEXT I
1270 MAT Z = ZER(S6,S6)      
1280 MAT C = ZER(S6,S6)      
1290 FOR I=1 TO S6 
1300 LET X(I)=Q(I)=O 
1310 NEXT I
1320 FOR I = 1 TO N
1330 LET Y(I) = F1*(P(I)-F2) 
1340 IF S1=1 THEN 1430
1350 IF L(I)>S2 THEN 1380
1360 LET Y(I) = Y(I)-E2      
1370 GO TO 1390
1380 LET Y(I) = Y(I)-E3      
1390 IF R(I)>S2 THEN 1420
1400 LET Y(I) = Y(I)+E2      
1410 GO TO 1430
1420 LET Y(I) = Y(I)+E3      
1430 LET L = L(I)  
1440 LET R = R(I)  
1450 LET Q(L) = Q(L)+Y(I)    
1460 LET Q(R) = Q(R)-Y(I)    
1470 LET Z(L,L) = Z(L,L)+1   
1480 LET Z(R,R) = Z(R,R)+1   
1490 LET Z(L,R) = Z(L,R)-1   
1500 LET Z(R,L) = Z(R,L)-1   
1510 LET Z(S5,L) = Z(S5,L)+1 
1520 LETZ(L,S5) =Z(S5,L)     
1530 LET Z(S5,R) = Z(S5,R)-1 
1540 LET Z(R,S5) = Z(S5,R)   
1550 LET Z(S5,S5) = Z(S5,S5)+1 
1560 LET Q(S5) = Q(S5)+Y(I)  
1570 NEXT I
1580 FOR I = 1 TO S4 
1590  LET Z(S6,I) = S(I)     
1600 LET Z(I,S6) = S(I)      
1610 NEXT I
1620 MAT C = INV(Z)
1630 FOR I = 1 TO S6 
1640 FOR J = 1 TO S6 
1650 LET X(I) = X(I) +C(I,J)*Q(J)      
1660 NEXT J
1670 NEXT I
1680 SCRATCH#2     
1690 WRITE #2,D$
1691 WRITE #2,O$
1700 WRITE #2,S1
1701 WRITE #2,S2
1702 WRITE #2,S3
1703 WRITE #2,S4
1704 WRITE #2,S5
1705 WRITE #2,S6
1710 WRITE #2,T1
1711 WRITE #2,T2
1712 WRITE #2,T3
1713 WRITE #2,T4
1720 WRITE #2,N
1721 WRITE #2,M
1722 WRITE #2,V1
1723 WRITE #2,F1
1724 WRITE #2,F2
1725 WRITE #2,E2
1726 WRITE #2,E3
1730 FOR I=1 TO S4 
1740 WRITE #2,A$(I)
1741 WRITE #2,S(I)
1750 NEXT I
1760 FOR I=1 TO N  
1770 WRITE #2,L(I)
1771 WRITE #2,P(I)
1772 WRITE #2,R(I)
1773 WRITE #2,Y(I)
1780 NEXT I
1790 FOR I=1 TO S6 
1800 WRITE #2,X(I)
1810 NEXT I
1820 DEF FNT(T)=(T-20)*(-40.6+(T-20)*(-.95+.01*(T-20)))    
1830 PRINT 
1840 PRINT 
1850 PRINT "*****THISSHOULD BE A CHAIN*****"
1855 RETURN
2000 PRINT
2010 PRINT
2020 PRINT
2030 PRINT TAB(25),"S A N D I A"
2040 PRINT
2050 PRINT"          LEAST SQUARE ANALYSIS OF STANDARD CELL DATA"     
2060 PRINT
2070 PRINT 
2080 PRINT 
2090 PRINT 
2130 READ #1,D$
2131 READ #1 O$
2140 READ #1, S1, S2, S3, S4, S5, S6   
2150 READ #1, T1, T2, T3, T4 
2160 READ  #1, N,M,V1,F1,F2,E2,E3      
2170 FOR I=1TO S4  
2180 READ #1, A$(I), S(I)    
2190 NEXT I
2200 FOR I=1 TO N  
2210 READ #1, L(I),P(I),R(I),Y(I)      
2220 NEXT I
2230 FOR I=1TO S6  
2240 READ #1, X(I) 
2250 NEXT I
2260 LET V2 = 0    
2270 FOR I = 1 TO N
2280 LET L = L(I)  
2290 LET R = R(I)  
2300 LET D(I) = Y(I)-(X(L)-X(R)+X(S5)) 
2310 LET V2 = D(I)*D(I)+V2   
2320 NEXT I
2330 LET V3=N-S4   
2340 LET W3=0      
2350 LET V4=SQR(V2/V3)       
2360 LET W2 = 0    
2370 PRINT "DATE-----" D$ TAB(30) "OBSERVER-----" O$       
2380 PRINT" DESIGN----" S2 "X" S3 ",",N "OBSERVATIONS"     
2390 PRINT 
2400 PRINT " STD DEV =" V4"ACC STD DEV =" V1     
2410 PRINT " DEG OF FREEDOM =" V3      
2420 PRINT 
2430 LET W1 = 0    
2440 IF W2 = 0 THEN 2460
2450 GO SUB 3640
2460 FOR I = 1 TO N 
2470 LET L = L(I) 
2480 LET R = R(I) 
2490 IF D(I)=0 THEN 2520
2500 LET X1=ABS(D(I))       
2510 LET D(I)=(D(I)/X1)*(INT(X1*100+.5))/100    
2520 IF W2 = 1 THEN 2560
2530 IF ABS(D(I))<3*V1 THEN 2580
2540 LET W1 = W1 + 1
2550 IF W1 = 1 THEN 2600
2560 PRINT I,TAB(6),A$(L),TAB(16),A$(R),TAB(27),
2561 PRINT P(I),TAB(35),Y(I)
2570 PRINT TAB(49) D(I)     
2580 NEXT I       
2590 GO TO 2630
2600 PRINT"       THE FOLLOWING DEVIATIONS EXCEED 3 SIGMA"
2610 GO SUB 3640
2620 GO TO 2560
2630 IF W2 = 1 THEN 2650
2640 IF W1 = 0 THEN 2670
2650 GO SUB 3670
2660 GO TO 2680
2670 PRINT" NO DEVIATIONS EXCEED 3 SIGMA."      
2680 IF W2 = 1 THEN 2740
2690 PRINT" FOR A COMPLETE LISTING OF DEV. ENTER 1, OTHERWISE ENTER 0"
2700 INPUT W2     
2710 PRINT
2720 IF W2 = 1 THEN 2430
2730 IF W1 = 0 THEN 2800
2740 PRINT" ENTER NO. OF OBS. TO BE DELETED"    
2750 INPUT N5     
2760 IF N5 <> O THEN 3260
2770 PRINT
2780 PRINT
2790 PRINT
2800 PRINT
2810 PRINT
2820 IF M <>0 THEN 2850
2830 PRINT" MEAN OF GROUP WAS NOT GIVEN."       
2840 GO TO 2860
2850 PRINT "REF GRP MEAN REDUCED BY 1.01, IN MICROVOLTS IS" M       
2860 PRINT
2870 PRINT" LEFT-RIGHT COMPONENT =" X(S5)       
2880 PRINT" TEMP GRP1 =" T2,"TEMP CORR GRP1 =" E2 
2890 IF S1 = 1 THEN 2910
2900 PRINT" TEMP GRP2 =" T4,"TEMP CORR GRP2 =" E3 
2910 PRINT
2920 PRINT "GRP1 CORR TO" T1 "DEG"    
2930 IF S1=1 THEN 2960
2940 PRINT
2950 PRINT "GRP2 CORR TO" T3 "DEG"    
2960 PRINT
2970 PRINT" ESTIMATES OF DIFF FROM MEAN AND TEMP CORRECTED EMFS"    
2980 PRINT
2990 PRINT "CELL NO"," ","EST VALUE","EMF 101+" 
3000 GO SUB 3670
3010 LET P5=0     
3020 FOR I = 1 TO S4
3030 PRINT A$(I), 
3040 IF S(I) = 1 THEN 3090
3050 PRINT" UNK ",
3060 LET P5=P5+1  
3070 LET W3=W3+X(I) 
3080 GO TO 3100
3090 PRINT" STD ",
3100 PRINT X(I),  
3110 IF M = 0 THEN 3140
3120 PRINT INT((X(I)+M)*100+.5)/100   
3130 GO TO 3150
3140 PRINT" "     
3150 NEXT I       
3160 IF M = 0 THEN 3220
3170 LET W3=W3/P5 
3180 LET W3 = INT((W3+M)*10+.5)       
3190 GO SUB 3670
3200 PRINT
3210 PRINT" MEAN OF UNK GRP REDUCED BY 1.01 IS" W3
3220 FOR I = 1 TO 5 
3230 PRINT
3240 NEXT I       
3250 STOP 
3260 PRINT N5 "LINE(S) ARE TO BE DELETED"       
3270 PRINT
3280 LET H=0      
3290 FOR K=1 TO N5
3300 PRINT "ENTER ONE OBSERVATION NUMBER TO BE DELETED"   
3310  INPUT N(K)  
3320 PRINT
3330 NEXT K       
3340 FOR I=1 TO (N5-1)      
3350 FOR J=1 TO (N5-I)      
3360 IF N(J)<N(J+1) THEN 3400
3370 LET T5=N(J)  
3380 LET N(J)=N(J+1)
3390 LET N(J+1)=T5
3400 NEXT J       
3410 NEXT I       
3420 PRINT "THE FOLLOWING OBSERVATIONS WILL BE DELETED"   
3430 PRINT
3440 FOR K=1 TO N5
3450 LET N6=N(K)  
3460 LET L = L(N6)
3470 LET R = R(N6)
3480 PRINT N6 TAB(6) A$(L) TAB(16) A$(R) TAB(27) P(N6)    
3490 NEXT K       
3500 PRINT
3510 FOR K=1 TO N5
3520 LET N6=N(K)  
3530 LET N6=N6-H  
3540 LET N = N-1  
3550 FOR I = N6 TO N
3560 LET P(I) = P(I+1)      
3570 LET L(I) = L(I+1)      
3580 LET R(I) = R(I+1)      
3590 NEXT I       
3600 LET H=H+1    
3610 NEXT K       
3620 PRINT
3630 GO TO 3690
3640 PRINT" OBS" TAB(7) "LEFT" TAB(17) "RIGHT" TAB(27) "OBS."       
3650 PRINTTAB(38) "OBS" TAB(49) "DEV."
3660 PRINT" NO." TAB(7) "CELL" TAB(17) "CELL" TAB(38) "T-CORR"      
3670 PRINT"----------------------------------------------------------"
3680 RETURN       
3690 SCRATCH #2   
3700 WRITE #2,D$
3701 WRITE #2,O$
3710 WRITE #2,S1
3711 WRITE #2,S2
3712 WRITE #2,S3
3713 WRITE #2,N
3714 WRITE #2,F1
3715 WRITE #2,F2
3720 WRITE #2,M
3721 WRITE #2,V1
3730 FOR I=1 TO S4
3740 WRITE #2, S(I) 
3750 NEXT I       
3760 FOR I=1 TO S4
3770 WRITE #2, A$(I)
3780 NEXT I       
3790 FOR I= 1 TO N
3791 WRITE #2,L(I)
3792 WRITE #2,R(I)
3793 NEXT I
3799 FOR I=1 TO N
3800 WRITE #2,T1
3801 WRITE #2,T2
3802 WRITE #2,T3
3803 WRITE #2,T4
3820 WRITE #2, T1  T2  T3  T4 
3840 WRITE #2, P(I) 
3850 NEXT I       
3860 PRINT
3870 PRINT "THIS SHOULB BE A CHAIN PILOT1"
9000 DATA "1/28/70",M8/GU-RUN1
9030 DATA 2,6,6
9045 DATA 24,1,0
9060 DATA 8240.70,1
9075 DATA 1,1,1,1,1,1,0,0,0,0,0,0
9090 DATA G1-1,G1-2,G1-3,G1-4,G1-5,G1-6
9105 DATA M8-1,M8-2,M8-3,M8-4,M8-5,M8-6
9120 DATA 7,2,7,4,11,4,11,6,9,6,9,2
9135 DATA 8,1,8,3,12,3,12,5,10,5,10,1
9150 DATA 1,11,1,9,5,9,5,7,3,7,3,11
9165 DATA 2,12,4,12,4,8,6,8,6,10,2,10
9180 DATA 28.,28.001,28.,28.001
9195 DATA 2.7,5.3,6.,2.1,12.8,4.0
9210 DATA 15.8,7.4,5.7,12.6,12.6,14.1
9225 DATA -14.2,-14.9,-13.4,-14.2,-5.0,-5.7
9240 DATA -3.5,-5.9,-7.6,-13.8,-12.1,-3.3
99999 END
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        