mirror of
https://github.com/retro-software/B5500-software.git
synced 2026-03-02 17:44:40 +00:00
2c72f7fd1d
1. Commit library tape images, directories, and extracted text files. 2. Commit additional utilities under Unisys-Emode-Tools.
1585 lines
125 KiB
Fortran
1585 lines
125 KiB
Fortran
FILE 2 = TAPE2, UNIT = TAPE, RECORD = 10, BUFFER = 2 00000100
|
|
FILE 3 = TAPE3, UNIT = TAPE, RECORD = 10, BUFFER = 2 00000200
|
|
FILE 5 = CRA, UNIT = READER, RECORD = 10, BUFFER = 2 00000300
|
|
FILE 6 = LPA, UNIT = PRINTER, RECORD = 15, BUFFER = 2 00000400
|
|
CID 0901HS 15 150 $BMD02V ANALYSIS OF VARIANCE FOR FACTORIAL DESIGN00000500
|
|
CBMD02V ANALYSIS OF VARIANCE FOR FACTORIAL DESIGN JULY 22, 1965 00000600
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00000700
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00000800
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[,FMT%60[ 00000900
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00001000
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00001100
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00001200
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00001300
|
|
C 00001400
|
|
904 FORMAT%78H1BMD02V - ANALYSIS OF VARIANCE FOR FACTORIAL DESIGN - VE00001500
|
|
XRSION OF JULY 22, 1965/ 00001600
|
|
X41H HEALTH SCIENCES COMPUTING FACILITY, UCLA//[ 00001700
|
|
C 00001800
|
|
MTAPE#5 00001900
|
|
CALL CONSTN 00002000
|
|
100 READ (5, 900) OK,PROB,NR,NV,NDF,MEAN,MT,NO,%NORD%I[,I#1,00002100
|
|
1 4[,%L%J[,J#1,8[,NTAPE,KL 00002200
|
|
ONR#%&6HPROBLM[ 00002300
|
|
C IF%OK*%-ONR[&ONR*%-OK[[ 290, 3, 290 00002400
|
|
IF (OK-ONR) 290,3,290 00002500
|
|
290 IF%OK-%&6HFINISH[[ 291, 400, 291 00002600
|
|
291 WRITE ( 6, 929 ) 00002700
|
|
GO TO 400 00002800
|
|
3 IF%NV*%NV-9[[ 6001,291,291 00002900
|
|
6001 IF%KL*%KL-6[[ 6002,291,291 00003000
|
|
6002 IF%%NTAPE-2[*%NTAPE-3[[ 8,291,8 00003100
|
|
8 CALL TPWD%NTAPE,MTAPE[ 00003200
|
|
2 WRITE (6,904 ) 00003300
|
|
WRITE (6, 905) PROB,NV 00003400
|
|
REWIND 2 00003500
|
|
REWIND 3 00003600
|
|
IERROR#0 00003700
|
|
IF%NR[ 41, 41, 42 00003800
|
|
41 NR#1 00003900
|
|
42 WRITE (6, 919) NR 00004000
|
|
WRITE (6, 924 ) 00004100
|
|
DO 43 I#1,NV 00004200
|
|
43 WRITE (6, 925) I,L%I[ 00004300
|
|
DO 4 I#1,NV 00004400
|
|
LF%I[#L%I[ 00004500
|
|
4 LO%I[#I 00004600
|
|
IF%NO[ 7, 7, 5 00004700
|
|
5 IF%NO-4[ 6005,6005,291 00004800
|
|
6005 DO 6 I#1,NO 00004900
|
|
I1#NORD%I[ 00005000
|
|
NSIZE%I[#L%I1[ 00005100
|
|
6 LO%I1[#LO%I1[&100 00005200
|
|
7 ONR#NR 00005300
|
|
NK#1 00005400
|
|
DO 10 I#1,NV 00005500
|
|
10 NK#NK*L%I[ 00005600
|
|
OK#NK 00005700
|
|
DF#MT 00005800
|
|
OK1#OK/DF 00005900
|
|
NK1#NK/MT 00006000
|
|
OK2#NK1 00006100
|
|
IF%OK1-OK2[ 11,11,12 00006200
|
|
11 NCARD#NK1 00006300
|
|
GO TO 13 00006400
|
|
12 NCARD#NK1&1 00006500
|
|
13 TMEAN#0.0 00006600
|
|
ONNN#OK*ONR&1.0 00006700
|
|
TOTAL#0.0 00006800
|
|
IF%NDF[ 303, 303, 300 00006900
|
|
300 READ (5,926)FISH,NT,%N%I[,TEMP%I[,I#1,NT[ 00007000
|
|
SUM#%&6HSPECTG[ 00007100
|
|
C IF%FISH*%-SUM[&SUM*%-FISH[[ 578, 579, 578 00007200
|
|
IF (FISH-SUM) 578,579,578 00007300
|
|
578 WRITE (6, 931 ) 00007400
|
|
GO TO 400 00007500
|
|
579 WRITE (6,930 ) 00007600
|
|
WRITE (6,927 ) 00007700
|
|
DO 301 I#1,NT 00007800
|
|
IF%N%I[*%N%I[-11[[ 301,578,578 00007900
|
|
301 WRITE (6,928)I,N%I[,TEMP%I[ 00008000
|
|
303 CALL VFCHCK%KL[ 00008100
|
|
KL#KL*12 00008200
|
|
READ (5, 932) %FMT%I[, I#1,KL[ 00008300
|
|
IF%NTAPE[479,479,480 00008400
|
|
479 NTAPE#5 00008500
|
|
480 IF%NR-1[14,14,17 00008600
|
|
14 READ (NTAPE,FMT)%X%I[,I#1,NK[ 00008700
|
|
IF%NDF[ 580, 580, 591 00008800
|
|
580 DO 585 I#1,NK 00008900
|
|
585 TMEAN#TMEAN&X%I[ 00009000
|
|
GO TO 15 00009100
|
|
591 DO 595 I#1,NK 00009200
|
|
DO 305 II#1,NT 00009300
|
|
DF#X%I[ 00009400
|
|
NDF#N%II[ 00009500
|
|
CONST#TEMP%II[ 00009600
|
|
CALL TRANS %DF,NDF,CONST,ONNN,IERROR,1,II[ 00009700
|
|
IF%IERROR[ 100, 305, 305 00009800
|
|
305 X%I[#DF 00009900
|
|
595 TMEAN#TMEAN&X%I[ 00010000
|
|
15 IF%MEAN[ 695, 695, 690 00010100
|
|
690 IF%MEAN-2[ 6004,6004,291 00010200
|
|
6004 WRITE (3) %X%I[, I#1,NK[ 00010300
|
|
695 TMEAN#TMEAN/OK 00010400
|
|
DO 16 I#1,NK 00010500
|
|
X%I[#X%I[-TMEAN 00010600
|
|
16 TOTAL#TOTAL&X%I[*X%I[ 00010700
|
|
GO TO 30 00010800
|
|
17 DO 18 I#1,NK 00010900
|
|
18 X%I[#0.0 00011000
|
|
IF%NDF[ 630, 630, 631 00011100
|
|
630 ASSIGN 317 TO NPT 00011200
|
|
GO TO 633 00011300
|
|
631 ASSIGN 635 TO NPT 00011400
|
|
633 DO 21 I#1,NR 00011500
|
|
NX#1 00011600
|
|
DO 20 J#1,NCARD 00011700
|
|
READ (NTAPE,FMT)%TSUM%K[,K#1,MT[ 00011800
|
|
DO 19 K#1,MT 00011900
|
|
GO TO NPT, %317,635[ 00012000
|
|
635 DO 318 II#1,NT 00012100
|
|
DF#TSUM%K[ 00012200
|
|
NDF#N%II[ 00012300
|
|
CONST#TEMP%II[ 00012400
|
|
CALL TRANS%DF,NDF,CONST,ONNN,IERROR,I,II[ 00012500
|
|
IF%IERROR[350,318,318 00012600
|
|
318 TSUM%K[#DF 00012700
|
|
19 CONTINUE 00012800
|
|
317 WRITE (2)%TSUM%K[,K#1,MT[ 00012900
|
|
DO 20 K#1,MT 00013000
|
|
X%NX[#X%NX[&TSUM%K[ 00013100
|
|
IF%NK-NX[ 21, 21, 20 00013200
|
|
20 NX#NX&1 00013300
|
|
21 CONTINUE 00013400
|
|
END FILE 2 00013500
|
|
REWIND 2 00013600
|
|
DO 22 I#1,NK 00013700
|
|
TMEAN#TMEAN&X%I[ 00013800
|
|
22 X%I[#X%I[/ONR 00013900
|
|
IF%MEAN[ 720, 720, 710 00014000
|
|
710 WRITE (3) %X%I[, I#1,NK[ 00014100
|
|
720 TMEAN#TMEAN/%OK*ONR[ 00014200
|
|
WITHIN#0.0 00014300
|
|
DO 25 I#1,NR 00014400
|
|
NX#1 00014500
|
|
DO 24 J#1,NCARD 00014600
|
|
READ (2)%TSUM%K[,K#1,MT[ 00014700
|
|
DO 24 K#1,MT 00014800
|
|
WITHIN#WITHIN&%TSUM%K[-X%NX[[**2 00014900
|
|
TOTAL #TOTAL &%TSUM%K[-TMEAN[**2 00015000
|
|
IF%NK-NX[ 25, 25, 24 00015100
|
|
24 NX#NX&1 00015200
|
|
25 CONTINUE 00015300
|
|
REWIND 2 00015400
|
|
23 DO 26 I#1,NK 00015500
|
|
26 X%I[#X%I[-TMEAN 00015600
|
|
30 WRITE (6, 915) 00015700
|
|
WRITE (6, 910) TMEAN 00015800
|
|
TEMP%1[#TMEAN 00015900
|
|
IF%NV-1[ 33, 33, 39 00016000
|
|
33 SUM#TOTAL-WITHIN 00016100
|
|
NDF#L%1[-1 00016200
|
|
DF#NDF 00016300
|
|
SMEAN#SUM/DF 00016400
|
|
WRITE (6, 921) NDF,SUM,SMEAN 00016500
|
|
NDFZ#L%1[*NR-1 00016600
|
|
NDF#NDFZ-NDF 00016700
|
|
DF#NDF 00016800
|
|
SMEAN#WITHIN/DF 00016900
|
|
WRITE (6, 922) NDF,WITHIN,SMEAN 00017000
|
|
WRITE (6, 917) NDFZ,TOTAL 00017100
|
|
IF%NO[ 730, 730, 34 00017200
|
|
730 IF%MEAN[ 100, 100, 255 00017300
|
|
34 NDF#L%1[ 00017400
|
|
DO 35 I#1,NDF 00017500
|
|
35 ST1%1,I[#X%I[ 00017600
|
|
GO TO 251 00017700
|
|
39 NVM#NV-1 00017800
|
|
FMT%1[#%&1H%[ 00017900
|
|
FMT%2[#%&6H1H [ 00018000
|
|
DO 60 I#3,10 00018100
|
|
60 FMT%I[#%&6H1X, [ 00018200
|
|
FMT%11[#%&6H9X,I6,[ 00018300
|
|
FMT%12[#%&6HF22.5,[ 00018400
|
|
FMT%13[#%&6HF17.5,[ 00018500
|
|
FMT%14[#%&1H[[ 00018600
|
|
DO 31 I#1,255 00018700
|
|
TSUM%I[#0.0 00018800
|
|
31 PQ%I[#0.0 00018900
|
|
I5#%8-NV[&1 00019000
|
|
DO 32 I#1,NVM 00019100
|
|
LL%I[#1 00019200
|
|
I1#I&1 00019300
|
|
DO 32 J#I1,NV 00019400
|
|
32 LL%I[#LL%I[*L%J[ 00019500
|
|
LL%NV[#1 00019600
|
|
DO 40 I#1,NV 00019700
|
|
I1#I-1 00019800
|
|
I2#NV-I1 00019900
|
|
I3#8-I1 00020000
|
|
L%I3[#L%I2[ 00020100
|
|
40 LO%I3[#LO%I2[ 00020200
|
|
K9#0 00020300
|
|
NDFZ#0 00020400
|
|
GSUM#0.0 00020500
|
|
WRITE (6, 911) 00020600
|
|
WRITE (6, 912 ) 00020700
|
|
DO 48 I#1,8 00020800
|
|
48 MI%I[#0 00020900
|
|
MT#1 00021000
|
|
MT1#1 00021100
|
|
DO 231 IJ#1,NVM 00021200
|
|
FMT%IJ&2[#%&6HI1, [ 00021300
|
|
MARK#0 00021400
|
|
IJK#IJ 00021500
|
|
JW#NVM-%IJ-1[ 00021600
|
|
50 DO 51 I#I5,8 00021700
|
|
51 M%I[#0 00021800
|
|
J1#I5&%NVM-JW[ 00021900
|
|
DO 52 I#I5,J1 00022000
|
|
52 M%I[#1 00022100
|
|
Z#1.0 00022200
|
|
70 CALL TABLE 00022300
|
|
IF%MARK[ 205,205,231 00022400
|
|
205 CALL XLOCAT %MEAN[ 00022500
|
|
NDF#1 00022600
|
|
MDF#1 00022700
|
|
I6#I5&%IJ-1[ 00022800
|
|
DO 206 I#I5,I6 00022900
|
|
NDF#NDF*MM%I[ 00023000
|
|
206 MDF#MDF*%MM%I[&1[ 00023100
|
|
FNDF#NDF 00023200
|
|
DF#MDF 00023300
|
|
SUM#ONR*DF*SUM/OK 00023400
|
|
K9#K9&1 00023500
|
|
DO 204 I#1,IJK 00023600
|
|
I4#I5&%I-1[ 00023700
|
|
ROSE#K5%I4[ 00023800
|
|
204 PQ%K9[#%PQ%K9[&ROSE[*10.0 00023900
|
|
PQ%K9[#PQ%K9[/10.0 00024000
|
|
I7#I6-I5 00024100
|
|
IF%I7-1[ 229, 207, 207 00024200
|
|
207 TTSUM#0.0 00024300
|
|
DO 220 JK#1,I7 00024400
|
|
ICE#0 00024500
|
|
KL#JK 00024600
|
|
DO 209 I#I5,8 00024700
|
|
209 MD%I[#0 00024800
|
|
I8#8-I7 00024900
|
|
I9#I8&%JK-1[ 00025000
|
|
DO 210 I#I8,I9 00025100
|
|
210 MD%I[#1 00025200
|
|
Y#1.0 00025300
|
|
211 CALL PART %MD,I8,Y,ICE,KL[ 00025400
|
|
IF%ICE[ 212, 212, 220 00025500
|
|
212 FISH#0.0 00025600
|
|
K#0 00025700
|
|
DO 215 I#I8,8 00025800
|
|
K#K&1 00025900
|
|
IF%MD%I[[ 215, 215, 213 00026000
|
|
213 I4#I5&%K-1[ 00026100
|
|
ROSE#K5%I4[ 00026200
|
|
FISH#%FISH&ROSE[*10.0 00026300
|
|
215 CONTINUE 00026400
|
|
FISH#FISH/10.0 00026500
|
|
DO 217 I#1, 255 00026600
|
|
IF%PQ%I[-FISH[ 217, 216, 211 00026700
|
|
216 TTSUM#TTSUM&TSUM%I[ 00026800
|
|
GO TO 211 00026900
|
|
217 CONTINUE 00027000
|
|
220 CONTINUE 00027100
|
|
SUM#SUM-TTSUM 00027200
|
|
229 SMEAN#SUM/FNDF 00027300
|
|
WRITE (6, FMT) %K5%I[,I#I5,I6[,NDF,SUM,SMEAN 00027400
|
|
NDFZ#NDFZ&NDF 00027500
|
|
GSUM#GSUM&SUM 00027600
|
|
TSUM%K9[#SUM 00027700
|
|
GO TO 70 00027800
|
|
231 CONTINUE 00027900
|
|
IF%NR-1[ 235, 235, 240 00028000
|
|
235 SUM#TOTAL-GSUM 00028100
|
|
NKDF#NK-1 00028200
|
|
NDF#NKDF-NDFZ 00028300
|
|
DF#NDF 00028400
|
|
SMEAN#SUM/DF 00028500
|
|
WRITE (6, 916) NDF,SUM,SMEAN 00028600
|
|
WRITE (6, 917) NKDF,TOTAL 00028700
|
|
GO TO 252 00028800
|
|
240 NDF#1 00028900
|
|
DO 241 I#I5,8 00029000
|
|
241 NDF#NDF*MM%I[ 00029100
|
|
DF#NDF 00029200
|
|
SUM#TOTAL-WITHIN-GSUM 00029300
|
|
SMEAN#SUM/DF 00029400
|
|
DO 250 I#1,NV 00029500
|
|
250 N%I[#I 00029600
|
|
I#NV&2 00029700
|
|
FMT%I[#%&6HI1, [ 00029800
|
|
WRITE (6, FMT) %N%I[,I#1,NV[,NDF,SUM,SMEAN 00029900
|
|
DF#ONR*OK-OK 00030000
|
|
NDF#DF 00030100
|
|
SMEAN#WITHIN/DF 00030200
|
|
WRITE (6, 920) NDF,WITHIN,SMEAN 00030300
|
|
DF#ONR*OK-1.0 00030400
|
|
NDF#DF 00030500
|
|
WRITE (6, 917) NDF,TOTAL 00030600
|
|
252 IF%NO[ 254, 254, 251 00030700
|
|
251 CALL FINISH 00030800
|
|
254 IF%MEAN[ 100, 100, 255 00030900
|
|
255 CALL MEANS%PROB,TMEAN,MEAN[ 00031000
|
|
GO TO 100 00031100
|
|
350 IF%J-NCARD[351,353,353 00031200
|
|
351 J#J&1 00031300
|
|
DO 352 JXX#J,NCARD 00031400
|
|
352 READ (NTAPE,FMT)%TSUM%K[,K#1,MT[ 00031500
|
|
353 IF%I-NR[354,100,100 00031600
|
|
354 I#I&1 00031700
|
|
DO 355 IXX#I,NR 00031800
|
|
DO 355 JXX#1,NCARD 00031900
|
|
355 READ (NTAPE,FMT)%TSUM%K[,K#1,MT[ 00032000
|
|
GO TO 100 00032100
|
|
900 FORMAT%A6,A2,I3,3I1,I2,5I1,8I3,23X,I2,I2[ 00032200
|
|
905 FORMAT%13H0PROBLEM NO. A2///20H NUMBER OF VARIABLESI6[ 00032300
|
|
909 FORMAT%I6,F21.5,F17.5[ 00032400
|
|
910 FORMAT%11H0GRAND MEANF20.5///[ 00032500
|
|
911 FORMAT%59H0SOURCE OF DEGREES OF SUMS OF M00032600
|
|
1EAN[ 00032700
|
|
912 FORMAT%61H VARIATION FREEDOM SQUARES SQ00032800
|
|
1UARES//[ 00032900
|
|
915 FORMAT%1H [ 00033000
|
|
916 FORMAT%18H RESIDUAL I6,F22.5,F17.5[ 00033100
|
|
917 FORMAT%18H TOTAL I6,F22.5[ 00033200
|
|
919 FORMAT%21H NUMBER OF REPLICATESI5[ 00033300
|
|
920 FORMAT%18H WITHIN REPLICATESI6,F22.5,F17.5[ 00033400
|
|
921 FORMAT%18H MEANS I6,F22.5,F17.5[ 00033500
|
|
922 FORMAT%18H WITHIN I6,F22.5,F17.5[ 00033600
|
|
924 FORMAT%27H0VARIABLE NO. OF LEVELS[ 00033700
|
|
925 FORMAT%I6,11X,I4[ 00033800
|
|
926 FORMAT%A6,I1,8%I2,F6.0[[ 00033900
|
|
927 FORMAT%35H0CARD NO. TRANS CODE CONSTANT[ 00034000
|
|
928 FORMAT%1H I5,11X,I2,8X,F10.5[ 00034100
|
|
929 FORMAT%22H0ERROR ON PROBLEM CARD[ 00034200
|
|
930 FORMAT%1H06X,21HTRANS-GENERATION CARD[ 00034300
|
|
931 FORMAT%31H0ERROR ON TRANS-GENERATION CARD[ 00034400
|
|
932 FORMAT%12A6[ 00034500
|
|
400 IF%MTAPE-5[402,402,401 00034600
|
|
401 CONTINUE 00034700
|
|
402 CALL EXIT 00034800
|
|
END 00034900
|
|
CCONSTN SUBROUTINE CONSTN FOR BMD02V 00035000
|
|
SUBROUTINE CONSTN 00035100
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00035200
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00035300
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[ 00035400
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00035500
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00035600
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00035700
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00035800
|
|
CO%2,1,1[# 0.70710680 00035900
|
|
CO%2,1,2[#-0.70710680 00036000
|
|
CO%3,1,1[#-0.70710680 00036100
|
|
CO%3,1,2[# 0.0 00036200
|
|
CO%3,1,3[# 0.70710680 00036300
|
|
CO%3,2,1[# 0.40824825 00036400
|
|
CO%3,2,2[#-0.81649650 00036500
|
|
CO%3,2,3[# 0.40824825 00036600
|
|
CO%4,1,1[#-0.67082040 00036700
|
|
CO%4,1,2[#-0.22360680 00036800
|
|
CO%4,1,3[# 0.22360680 00036900
|
|
CO%4,1,4[# 0.67082040 00037000
|
|
CO%4,2,1[# 0.5 00037100
|
|
CO%4,2,2[#-0.5 00037200
|
|
CO%4,2,3[#-0.5 00037300
|
|
CO%4,2,4[# 0.5 00037400
|
|
CO%4,3,1[#-0.22360680 00037500
|
|
CO%4,3,2[# 0.67082040 00037600
|
|
CO%4,3,3[#-0.67082040 00037700
|
|
CO%4,3,4[# 0.22360680 00037800
|
|
CO%5,1,1[#-0.63245560 00037900
|
|
CO%5,1,2[#-0.31622780 00038000
|
|
CO%5,1,3[# 0.0 00038100
|
|
CO%5,1,4[# 0.31622780 00038200
|
|
CO%5,1,5[# 0.63245560 00038300
|
|
CO%5,2,1[# 0.53452254 00038400
|
|
CO%5,2,2[#-0.26726127 00038500
|
|
CO%5,2,3[#-0.53452254 00038600
|
|
CO%5,2,4[#-0.26726127 00038700
|
|
CO%5,2,5[# 0.53452254 00038800
|
|
CO%5,3,1[#-0.31622780 00038900
|
|
CO%5,3,2[# 0.63245560 00039000
|
|
CO%5,3,3[# 0.0 00039100
|
|
CO%5,3,4[#-0.63245560 00039200
|
|
CO%5,3,5[# 0.31622780 00039300
|
|
CO%5,4,1[# 0.11952286 00039400
|
|
CO%5,4,2[#-0.47809146 00039500
|
|
CO%5,4,3[# 0.71713718 00039600
|
|
CO%5,4,4[#-0.47809146 00039700
|
|
CO%5,4,5[# 0.11952286 00039800
|
|
CO%6,1,1[#-0.59761432 00039900
|
|
CO%6,1,2[#-0.35856859 00040000
|
|
CO%6,1,3[#-0.11952286 00040100
|
|
CO%6,1,4[# 0.11952286 00040200
|
|
CO%6,1,5[# 0.35856859 00040300
|
|
CO%6,1,6[# 0.59761432 00040400
|
|
CO%6,2,1[# 0.54554475 00040500
|
|
CO%6,2,2[#-0.10910895 00040600
|
|
CO%6,2,3[#-0.43643579 00040700
|
|
CO%6,2,4[#-0.43643579 00040800
|
|
CO%6,2,5[#-0.10910895 00040900
|
|
CO%6,2,6[# 0.54554475 00041000
|
|
CO%6,3,1[#-0.37267794 00041100
|
|
CO%6,3,2[# 0.52174911 00041200
|
|
CO%6,3,3[# 0.29814235 00041300
|
|
CO%6,3,4[#-0.29814235 00041400
|
|
CO%6,3,5[#-0.52174911 00041500
|
|
CO%6,3,6[# 0.37267794 00041600
|
|
CO%6,4,1[# 0.18898222 00041700
|
|
CO%6,4,2[#-0.56694667 00041800
|
|
CO%6,4,3[# 0.37796445 00041900
|
|
CO%6,4,4[# 0.37796445 00042000
|
|
CO%6,4,5[#-0.56694667 00042100
|
|
CO%6,4,6[# 0.18898222 00042200
|
|
CO%6,5,1[#-0.06299407 00042300
|
|
CO%6,5,2[# 0.31497035 00042400
|
|
CO%6,5,3[#-0.62994070 00042500
|
|
CO%6,5,4[# 0.62994070 00042600
|
|
CO%6,5,5[#-0.31497035 00042700
|
|
CO%6,5,6[# 0.06299407 00042800
|
|
CO%7,1,1[#-0.56694667 00042900
|
|
CO%7,1,2[#-0.37796445 00043000
|
|
CO%7,1,3[#-0.18898222 00043100
|
|
CO%7,1,4[# 0.0 00043200
|
|
CO%7,1,5[# 0.18898222 00043300
|
|
CO%7,1,6[# 0.37796445 00043400
|
|
CO%7,1,7[# 0.56694667 00043500
|
|
CO%7,2,1[# 0.54554475 00043600
|
|
CO%7,2,2[# 0.0 00043700
|
|
CO%7,2,3[#-0.32732685 00043800
|
|
CO%7,2,4[#-0.43643580 00043900
|
|
CO%7,2,5[#-0.32732685 00044000
|
|
CO%7,2,6[# 0.0 00044100
|
|
CO%7,2,7[# 0.54554475 00044200
|
|
CO%7,3,1[#-0.40824825 00044300
|
|
CO%7,3,2[# 0.40824825 00044400
|
|
CO%7,3,3[# 0.40824825 00044500
|
|
CO%7,3,4[# 0.0 00044600
|
|
CO%7,3,5[#-0.40824825 00044700
|
|
CO%7,3,6[#-0.40824825 00044800
|
|
CO%7,3,7[# 0.40824825 00044900
|
|
CO%7,4,1[# 0.24174696 00045000
|
|
CO%7,4,2[#-0.56407624 00045100
|
|
CO%7,4,3[# 0.08058232 00045200
|
|
CO%7,4,4[# 0.48349392 00045300
|
|
CO%7,4,5[# 0.08058232 00045400
|
|
CO%7,4,6[#-0.56407624 00045500
|
|
CO%7,4,7[# 0.24174696 00045600
|
|
CO%7,5,1[#-0.10910895 00045700
|
|
CO%7,5,2[# 0.43643580 00045800
|
|
CO%7,5,3[#-0.54554475 00045900
|
|
CO%7,5,4[# 0.0 00046000
|
|
CO%7,5,5[# 0.54554475 00046100
|
|
CO%7,5,6[#-0.43643580 00046200
|
|
CO%7,5,7[# 0.10910895 00046300
|
|
CO%7,6,1[# 0.03289758 00046400
|
|
CO%7,6,2[#-0.19738549 00046500
|
|
CO%7,6,3[# 0.49346374 00046600
|
|
CO%7,6,4[#-0.65795166 00046700
|
|
CO%7,6,5[# 0.49346374 00046800
|
|
CO%7,6,6[#-0.19738549 00046900
|
|
CO%7,6,7[# 0.03289758 00047000
|
|
CO%8,1,1[#-0.54006178 00047100
|
|
CO%8,1,2[#-0.38575842 00047200
|
|
CO%8,1,3[#-0.23145505 00047300
|
|
CO%8,1,4[#-0.07715168 00047400
|
|
CO%8,1,5[# 0.07715168 00047500
|
|
CO%8,1,6[# 0.23145505 00047600
|
|
CO%8,1,7[# 0.38575842 00047700
|
|
CO%8,1,8[# 0.54006178 00047800
|
|
CO%8,2,1[# 0.54006178 00047900
|
|
CO%8,2,2[# 0.07715168 00048000
|
|
CO%8,2,3[#-0.23145505 00048100
|
|
CO%8,2,4[#-0.38575843 00048200
|
|
CO%8,2,5[#-0.38575843 00048300
|
|
CO%8,2,6[#-0.23145505 00048400
|
|
CO%8,2,7[# 0.07715168 00048500
|
|
CO%8,2,8[# 0.54006178 00048600
|
|
CO%8,3,1[#-0.43082013 00048700
|
|
CO%8,3,2[# 0.30772862 00048800
|
|
CO%8,3,3[# 0.43082013 00048900
|
|
CO%8,3,4[# 0.18463720 00049000
|
|
CO%8,3,5[#-0.18463720 00049100
|
|
CO%8,3,6[#-0.43082013 00049200
|
|
CO%8,3,7[#-0.30772867 00049300
|
|
CO%8,3,8[# 0.43082013 00049400
|
|
CO%8,4,1[# 0.28203800 00049500
|
|
CO%8,4,2[#-0.52378487 00049600
|
|
CO%8,4,3[#-0.12087343 00049700
|
|
CO%8,4,4[# 0.36260294 00049800
|
|
CO%8,4,5[# 0.36260294 00049900
|
|
CO%8,4,6[#-0.12087343 00050000
|
|
CO%8,4,7[#-0.52378487 00050100
|
|
CO%8,4,8[# 0.28203800 00050200
|
|
CO%8,5,1[#-0.14978626 00050300
|
|
CO%8,5,2[# 0.49215485 00050400
|
|
CO%8,5,3[#-0.36376663 00050500
|
|
CO%8,5,4[#-0.32097056 00050600
|
|
CO%8,5,5[# 0.32097056 00050700
|
|
CO%8,5,6[# 0.36376663 00050800
|
|
CO%8,5,7[#-0.49215485 00050900
|
|
CO%8,5,8[# 0.14978626 00051000
|
|
CO%8,6,1[# 0.06154573 00051100
|
|
CO%8,6,2[#-0.30772866 00051200
|
|
CO%8,6,3[# 0.55391160 00051300
|
|
CO%8,6,4[#-0.30772866 00051400
|
|
CO%8,6,5[#-0.30772866 00051500
|
|
CO%8,6,6[# 0.55391160 00051600
|
|
CO%8,6,7[#-0.30772866 00051700
|
|
CO%8,6,8[# 0.06154573 00051800
|
|
CO%8,7,1[#-0.01706972 00051900
|
|
CO%8,7,2[# 0.11948809 00052000
|
|
CO%8,7,3[#-0.35846429 00052100
|
|
CO%8,7,4[# 0.59744050 00052200
|
|
CO%8,7,5[#-0.59744050 00052300
|
|
CO%8,7,6[# 0.35846429 00052400
|
|
CO%8,7,7[#-0.11948809 00052500
|
|
CO%8,7,8[# 0.01706972 00052600
|
|
CO%9,1,1[#-0.51639775 00052700
|
|
CO%9,1,2[#-0.38729832 00052800
|
|
CO%9,1,3[#-0.25819888 00052900
|
|
CO%9,1,4[#-0.12909943 00053000
|
|
CO%9,1,5[# 0.0 00053100
|
|
CO%9,1,6[# 0.12909943 00053200
|
|
CO%9,1,7[# 0.25819888 00053300
|
|
CO%9,1,8[# 0.38729832 00053400
|
|
CO%9,1,9[# 0.51639775 00053500
|
|
CO%9,2,1[# 0.53181614 00053600
|
|
CO%9,2,2[# 0.13295404 00053700
|
|
CO%9,2,3[#-0.15194747 00053800
|
|
CO%9,2,4[#-0.32288837 00053900
|
|
CO%9,2,5[#-0.37986867 00054000
|
|
CO%9,2,6[#-0.32288837 00054100
|
|
CO%9,2,7[#-0.15194747 00054200
|
|
CO%9,2,8[# 0.13295404 00054300
|
|
CO%9,2,9[# 0.53181614 00054400
|
|
CO%9,3,1[#-0.44494913 00054500
|
|
CO%9,3,2[# 0.22247458 00054600
|
|
CO%9,3,3[# 0.41316705 00054700
|
|
CO%9,3,4[# 0.28603873 00054800
|
|
CO%9,3,5[# 0.0 00054900
|
|
CO%9,3,6[#-0.28603873 00055000
|
|
CO%9,3,7[#-0.41316705 00055100
|
|
CO%9,3,8[#-0.22247458 00055200
|
|
CO%9,3,9[# 0.44494913 00055300
|
|
CO%9,4,1[# 0.31288283 00055400
|
|
CO%9,4,2[#-0.46932425 00055500
|
|
CO%9,4,3[#-0.24583650 00055600
|
|
CO%9,4,4[# 0.20113896 00055700
|
|
CO%9,4,5[# 0.40227793 00055800
|
|
CO%9,4,6[# 0.20113896 00055900
|
|
CO%9,4,7[#-0.24583650 00056000
|
|
CO%6,4,8[#-0.46932425 00056100
|
|
CO%9,4,9[# 0.31288283 00056200
|
|
CO%9,5,1[#-0.18490004 00056300
|
|
CO%9,5,2[# 0.50847512 00056400
|
|
CO%9,5,3[#-0.18490004 00056500
|
|
CO%9,5,4[#-0.41602510 00056600
|
|
CO%9,5,5[# 0.0 00056700
|
|
CO%9,5,6[# 0.41602510 00056800
|
|
CO%9,5,7[# 0.18490004 00056900
|
|
CO%9,5,8[#-0.50847512 00057000
|
|
CO%9,5,9[# 0.18490004 00057100
|
|
CO%9,6,1[# 0.08989331 00057200
|
|
CO%9,6,2[#-0.38204659 00057300
|
|
CO%9,6,3[# 0.49441324 00057400
|
|
CO%9,6,4[# 0.02247332 00057500
|
|
CO%9,6,5[#-0.44946658 00057600
|
|
CO%9,6,6[# 0.02247332 00057700
|
|
CO%9,6,7[# 0.49441324 00057800
|
|
CO%9,6,8[#-0.38204659 00057900
|
|
CO%9,6,9[# 0.08989331 00058000
|
|
CO%9,7,1[#-0.03413943 00058100
|
|
CO%9,7,2[# 0.20483660 00058200
|
|
CO%9,7,3[#-0.47795207 00058300
|
|
CO%9,7,4[# 0.47795207 00058400
|
|
CO%9,7,5[# 0.0 00058500
|
|
CO%9,7,6[#-0.47795207 00058600
|
|
CO%9,7,7[# 0.47795207 00058700
|
|
CO%9,7,8[#-0.20483660 00058800
|
|
CO%9,7,9[# 0.03413943 00058900
|
|
CO%9,8,1[# 0.00881476 00059000
|
|
CO%9,8,2[#-0.07051811 00059100
|
|
CO%9,8,3[# 0.24681340 00059200
|
|
CO%9,8,4[#-0.49362680 00059300
|
|
CO%9,8,5[# 0.61703350 00059400
|
|
CO%9,8,6[#-0.49362680 00059500
|
|
CO%9,8,7[# 0.24681340 00059600
|
|
CO%9,8,8[#-0.07051811 00059700
|
|
CO%9,8,9[# 0.00881476 00059800
|
|
RETURN 00059900
|
|
END 00060000
|
|
CFINISH SUBROUTINE FINISH FOR BMD02V 9-5-6300060100
|
|
SUBROUTINE FINISH 00060200
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00060300
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00060400
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[ 00060500
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00060600
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00060700
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00060800
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00060900
|
|
WRITE (6, 935) 00061000
|
|
WRITE (6, 936) %NORD%J[, J#1,NO[ 00061100
|
|
OR#NR 00061200
|
|
DO 270 I#1,NO 00061300
|
|
WRITE (6, 935) 00061400
|
|
K#NORD%I[ 00061500
|
|
KL#LF%K[ 00061600
|
|
KI#KL-1 00061700
|
|
DO 253 M1#1,KI 00061800
|
|
C%M1,1[#0.0 00061900
|
|
DO 253 M2#1,KL 00062000
|
|
253 C%M1,1[#C%M1,1[&CO%KL,M1,M2[*ST1%I,M2[ 00062100
|
|
NDF#NK/KL 00062200
|
|
DF#NDF 00062300
|
|
DO 257 J#1,KI 00062400
|
|
257 C%J,1[#C%J,1[*C%J,1[*OR/DF 00062500
|
|
DO 267 J#1,KI 00062600
|
|
GO TO %261, 262, 263, 264[, J 00062700
|
|
261 WRITE (6, 930) K, C%1,1[ 00062800
|
|
IF%KI-2[ 268, 267, 267 00062900
|
|
262 WRITE (6, 931) K,C%2,1[ 00063000
|
|
IF%KI-3[ 268, 267, 267 00063100
|
|
263 WRITE (6, 932) K,C%3,1[ 00063200
|
|
IF%KI-4[ 268, 267, 267 00063300
|
|
264 REM#0.0 00063400
|
|
DO 265 M1#4,KI 00063500
|
|
265 REM#REM&C%M1,1[ 00063600
|
|
WRITE (6, 933) K, REM 00063700
|
|
GO TO 268 00063800
|
|
267 CONTINUE 00063900
|
|
268 RES#0.0 00064000
|
|
DO 269 J#1,KI 00064100
|
|
269 RES#RES&C%J,1[ 00064200
|
|
WRITE (6, 934) K,RES 00064300
|
|
270 CONTINUE 00064400
|
|
IF%NO-2[ 100, 360, 369 00064500
|
|
360 I2#NSIZE%1[ 00064600
|
|
I3#NSIZE%2[ 00064700
|
|
NO1#NORD%1[ 00064800
|
|
NO2#NORD%2[ 00064900
|
|
KB1#1 00065000
|
|
KB2#2 00065100
|
|
IN#1 00065200
|
|
IF%NV-2[ 362, 362, 367 00065300
|
|
362 K#0 00065400
|
|
DO 365 I#1,I2 00065500
|
|
DO 365 J#1,I3 00065600
|
|
K#K&1 00065700
|
|
365 ST%1,I,J[#X%K[ 00065800
|
|
367 CALL MATRIX 00065900
|
|
CALL REPORT 00066000
|
|
GO TO 100 00066100
|
|
369 IF%NO-3[ 100, 370, 379 00066200
|
|
370 I2#NSIZE%1[ 00066300
|
|
I3#NSIZE%2[ 00066400
|
|
NO1#NORD%1[ 00066500
|
|
NO2#NORD%2[ 00066600
|
|
KB1#1 00066700
|
|
KB2#2 00066800
|
|
IN#1 00066900
|
|
CALL MATRIX 00067000
|
|
CALL REPORT 00067100
|
|
I3#NSIZE%3[ 00067200
|
|
NO2#NORD%3[ 00067300
|
|
KB2#3 00067400
|
|
IN#2 00067500
|
|
CALL MATRIX 00067600
|
|
CALL REPORT 00067700
|
|
I2#NSIZE%2[ 00067800
|
|
NO1#NORD%2[ 00067900
|
|
KB1#2 00068000
|
|
IN#3 00068100
|
|
CALL MATRIX 00068200
|
|
CALL REPORT 00068300
|
|
GO TO 100 00068400
|
|
379 I2#NSIZE%1[ 00068500
|
|
I3#NSIZE%2[ 00068600
|
|
NO1#NORD%1[ 00068700
|
|
NO2#NORD%2[ 00068800
|
|
KB1#1 00068900
|
|
KB2#2 00069000
|
|
IN#1 00069100
|
|
CALL MATRIX 00069200
|
|
CALL REPORT 00069300
|
|
I3#NSIZE%3[ 00069400
|
|
NO2#NORD%3[ 00069500
|
|
KB2#3 00069600
|
|
IN#2 00069700
|
|
CALL MATRIX 00069800
|
|
CALL REPORT 00069900
|
|
I3#NSIZE%4[ 00070000
|
|
NO2#NORD%4[ 00070100
|
|
KB2#4 00070200
|
|
IN#3 00070300
|
|
CALL MATRIX 00070400
|
|
CALL REPORT 00070500
|
|
I2#NSIZE%2[ 00070600
|
|
I3#NSIZE%3[ 00070700
|
|
NO1#NORD%2[ 00070800
|
|
NO2#NORD%3[ 00070900
|
|
KB1#2 00071000
|
|
KB2#3 00071100
|
|
IN#4 00071200
|
|
CALL MATRIX 00071300
|
|
CALL REPORT 00071400
|
|
I3#NSIZE%4[ 00071500
|
|
NO2#NORD%4[ 00071600
|
|
KB2#4 00071700
|
|
IN#5 00071800
|
|
CALL MATRIX 00071900
|
|
CALL REPORT 00072000
|
|
I2#NSIZE%3[ 00072100
|
|
NO1#NORD%3[ 00072200
|
|
KB1#3 00072300
|
|
IN#6 00072400
|
|
CALL MATRIX 00072500
|
|
CALL REPORT 00072600
|
|
930 FORMAT%I6,8H LINEARF20.5[ 00072700
|
|
931 FORMAT%I6,11H QUADRATICF17.5[ 00072800
|
|
932 FORMAT%I6,7H CUBICF21.5[ 00072900
|
|
933 FORMAT%I6,11H REMAINDERF17.5[ 00073000
|
|
934 FORMAT%I6,16H TOTAL F24.5[ 00073100
|
|
935 FORMAT%1H0[ 00073200
|
|
936 FORMAT%21H ORDERED VARIABLES...4I3[ 00073300
|
|
100 RETURN 00073400
|
|
END 00073500
|
|
CMATRIX SUBROUTINE MATRIX FOR BMD02V 9-3-6300073600
|
|
SUBROUTINE MATRIX 00073700
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00073800
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00073900
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[ 00074000
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00074100
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00074200
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00074300
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00074400
|
|
I4#I2-1 00074500
|
|
I5#I3-1 00074600
|
|
DO 5 I#1,9 00074700
|
|
DO 5 J#1,9 00074800
|
|
C%I,J[#0.0 00074900
|
|
5 D%I,J[#0.0 00075000
|
|
DO 10 I#1,I4 00075100
|
|
DO 10 J#1,I3 00075200
|
|
DO 10 K#1,I2 00075300
|
|
10 C%I,J[#C%I,J[&CO%I2,I,K[*ST%IN,K,J[ 00075400
|
|
DO 16 I#1,I4 00075500
|
|
DO 16 J#1,I5 00075600
|
|
DO 15 K#1,I3 00075700
|
|
15 D%I,J[#D%I,J[&C%I,K[*CO%I3,J,K[ 00075800
|
|
16 D%I,J[#D%I,J[*D%I,J[ 00075900
|
|
KL#NK/%I2*I3[ 00076000
|
|
DIV#KL 00076100
|
|
OR#NR 00076200
|
|
DIV#DIV/OR 00076300
|
|
DO 30 I#1,I4 00076400
|
|
DO 30 J#1,I5 00076500
|
|
30 D%I,J[#D%I,J[/DIV 00076600
|
|
WRITE (6, 908) 00076700
|
|
WRITE (6, 900) NO2 00076800
|
|
WRITE (6, 901 ) 00076900
|
|
WRITE (6, 902) NO1 00077000
|
|
TOTAL#0.0 00077100
|
|
DO 40 J#1,I5 00077200
|
|
40 TOTAL#TOTAL&D%1,J[ 00077300
|
|
C%1,9[#TOTAL-D%1,1[-D%1,2[-D%1,3[ 00077400
|
|
WRITE (6, 903) D%1,1[,D%1,2[,D%1,3[,C%1,9[,TOTAL 00077500
|
|
TOTAL#0.0 00077600
|
|
DO 41 J#1,I5 00077700
|
|
41 TOTAL#TOTAL&D%2,J[ 00077800
|
|
C%2,9[#TOTAL-D%2,1[-D%2,2[-D%2,3[ 00077900
|
|
WRITE (6, 904) D%2,1[,D%2,2[,D%2,3[,C%2,9[,TOTAL 00078000
|
|
TOTAL#0.0 00078100
|
|
DO 42 J#1,I5 00078200
|
|
42 TOTAL#TOTAL&D%3,J[ 00078300
|
|
C%3,9[#TOTAL-D%3,1[-D%3,2[-D%3,3[ 00078400
|
|
WRITE (6, 905) D%3,1[,D%3,2[,D%3,3[,C%3,9[,TOTAL 00078500
|
|
TOTAL#0.0 00078600
|
|
DO 46 J#1,3 00078700
|
|
C%1,J[#0.0 00078800
|
|
DO 45 I#1,I4 00078900
|
|
45 C%1,J[#C%1,J[&D%I,J[ 00079000
|
|
C%2,J[#C%1,J[-D%1,J[-D%2,J[-D%3,J[ 00079100
|
|
46 TOTAL#TOTAL&C%2,J[ 00079200
|
|
C%4,9[#0.0 00079300
|
|
IF%I5-3[ 48,48,464 00079400
|
|
464 IF%I4-3[ 48,48,465 00079500
|
|
465 DO 47 J#4,I5 00079600
|
|
DO 47 I#4,I4 00079700
|
|
47 C%4,9[#C%4,9[&D%I,J[ 00079800
|
|
48 TOTAL#TOTAL&C%4,9[ 00079900
|
|
WRITE (6,906)C%2,1[,C%2,2[,C%2,3[,C%4,9[,TOTAL 00080000
|
|
C%1,4[#C%1,9[&C%2,9[&C%3,9[&C%4,9[ 00080100
|
|
TOTAL#C%1,1[&C%1,2[&C%1,3[&C%1,4[ 00080200
|
|
WRITE (6, 907) C%1,1[,C%1,2[,C%1,3[,C%1,4[,TOTAL 00080300
|
|
900 FORMAT%1H047X,8HVARIABLEI3,9H ORDERED[ 00080400
|
|
901 FORMAT%28X,6HLINEAR7X,9HQUADRATIC9X,5HCUBIC7X,9HREMAINDER9X,5HTOTA00080500
|
|
1L[ 00080600
|
|
902 FORMAT%9H VARIABLEI3,9H ORDERED[ 00080700
|
|
903 FORMAT%8X,12HLINEAR 5F15.5[ 00080800
|
|
904 FORMAT%8X,12HQUADRATIC 5F15.5[ 00080900
|
|
905 FORMAT%8X,12HCUBIC 5F15.5[ 00081000
|
|
906 FORMAT%8X,12HREMAINDER 5F15.5[ 00081100
|
|
907 FORMAT%8X,12HTOTAL 5F15.5[ 00081200
|
|
908 FORMAT%1H0[ 00081300
|
|
RETURN 00081400
|
|
END 00081500
|
|
CPART SUBROUTINE PART FOR BMD02V 9-5-6300081600
|
|
SUBROUTINE PART %MZ,I8,Y,ICE,KL[ 00081700
|
|
DIMENSION MZ%8[ 00081800
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00081900
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00082000
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[,FMT%60[ 00082100
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00082200
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00082300
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00082400
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00082500
|
|
GO TO %71,72,73,74,75,76[,KL 00082600
|
|
71 IF%Y-1.0[ 61, 60, 61 00082700
|
|
60 IA#I8 00082800
|
|
GO TO 299 00082900
|
|
61 IF%IA-8[ 130, 200, 200 00083000
|
|
72 IF%Y-1.0[ 63, 62, 63 00083100
|
|
62 IB#I8 00083200
|
|
IA#I8&1 00083300
|
|
GO TO 299 00083400
|
|
63 IF%IA-8[ 130, 93, 93 00083500
|
|
93 IF%IB-7[ 129, 200, 200 00083600
|
|
73 IF%Y-1.0[ 65, 64, 65 00083700
|
|
64 IC#I8 00083800
|
|
IB#I8&1 00083900
|
|
IA#I8&2 00084000
|
|
GO TO 299 00084100
|
|
65 IF%IA-8[ 130, 97, 97 00084200
|
|
97 IF%IB-7[ 129, 98, 98 00084300
|
|
98 IF%IC-6[ 128, 200, 200 00084400
|
|
74 IF%Y-1.0[ 67, 66, 67 00084500
|
|
66 ID#I8 00084600
|
|
IC#I8&1 00084700
|
|
IB#I8&2 00084800
|
|
IA#I8&3 00084900
|
|
GO TO 299 00085000
|
|
67 IF%IA-8[ 130, 103, 103 00085100
|
|
103 IF%IB-7[ 129, 104, 104 00085200
|
|
104 IF%IC-6[ 128, 105, 105 00085300
|
|
105 IF%ID-5[ 127, 200, 200 00085400
|
|
75 IF%Y-1.0[ 69,68,69 00085500
|
|
68 IE#I8 00085600
|
|
ID#I8&1 00085700
|
|
IC#I8&2 00085800
|
|
IB#I8&3 00085900
|
|
IA#I8&4 00086000
|
|
GO TO 299 00086100
|
|
69 IF%IA-8[ 130, 110, 110 00086200
|
|
110 IF%IB-7[ 129, 111, 111 00086300
|
|
111 IF%IC-6[ 128, 112, 112 00086400
|
|
112 IF%ID-5[ 127, 113, 113 00086500
|
|
113 IF%IE-4[ 126, 200, 200 00086600
|
|
76 IF%Y-1.0[ 81, 80, 81 00086700
|
|
80 IG#I8 00086800
|
|
IE#I8&1 00086900
|
|
ID#I8&2 00087000
|
|
IC#I8&3 00087100
|
|
IB#I8&4 00087200
|
|
IA#I8&5 00087300
|
|
GO TO 299 00087400
|
|
81 IF%IA-8[ 130, 120, 120 00087500
|
|
120 IF%IB-7[ 129, 121, 121 00087600
|
|
121 IF%IC-6[ 128, 122, 122 00087700
|
|
122 IF%ID-5[ 127, 123, 123 00087800
|
|
123 IF%IE-4[ 126, 124, 124 00087900
|
|
124 IF%IG-3[ 125, 200, 200 00088000
|
|
125 DO 136 I#IG,8 00088100
|
|
136 MZ%I[#0 00088200
|
|
IG#IG&1 00088300
|
|
MZ%IG[#1 00088400
|
|
IE#IG&1 00088500
|
|
MZ%IE[#1 00088600
|
|
10 ID#IE&1 00088700
|
|
MZ%ID[#1 00088800
|
|
20 IC#ID&1 00088900
|
|
MZ%IC[#1 00089000
|
|
30 IB#IC&1 00089100
|
|
MZ%IB[#1 00089200
|
|
40 IA#IB&1 00089300
|
|
MZ%IA[#1 00089400
|
|
GO TO 300 00089500
|
|
126 DO 135 I#IE,8 00089600
|
|
135 MZ%I[#0 00089700
|
|
IE#IE&1 00089800
|
|
MZ%IE[#1 00089900
|
|
GO TO 10 00090000
|
|
127 DO 134 I#ID,8 00090100
|
|
134 MZ%I[#0 00090200
|
|
ID#ID&1 00090300
|
|
MZ%ID[#1 00090400
|
|
GO TO 20 00090500
|
|
128 DO 133 I#IC,8 00090600
|
|
133 MZ%I[#0 00090700
|
|
IC#IC&1 00090800
|
|
MZ%IC[#1 00090900
|
|
GO TO 30 00091000
|
|
129 DO 132 I#IB,8 00091100
|
|
132 MZ%I[#0 00091200
|
|
IB#IB&1 00091300
|
|
MZ%IB[#1 00091400
|
|
GO TO 40 00091500
|
|
130 DO 131 I#IA,8 00091600
|
|
131 MZ%I[#0 00091700
|
|
IA#IA&1 00091800
|
|
MZ%IA[#1 00091900
|
|
GO TO 300 00092000
|
|
200 ICE#1 00092100
|
|
GO TO 300 00092200
|
|
299 Y#0.0 00092300
|
|
300 RETURN 00092400
|
|
END 00092500
|
|
CMEANS SUBROUTINE MEANS FOR BMD02V JULY 22, 1965 00092600
|
|
SUBROUTINE MEANS%PROB,TMEAN,MEAN[ 00092700
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00092800
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00092900
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[,FT%12[ 00093000
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00093100
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00093200
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00093300
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00093400
|
|
END FILE 3 00093500
|
|
REWIND 3 00093600
|
|
WRITE (6, 900) PROB 00093700
|
|
IF%MEAN-1[300,4,3 00093800
|
|
3 WRITE (6,901 ) 00093900
|
|
4 TSUM%1[#%&3HI3,[ 00094000
|
|
TSUM%2[#%&4H2I3,[ 00094100
|
|
TSUM%3[#%&4H3I3,[ 00094200
|
|
TSUM%4[#%&4H4I3,[ 00094300
|
|
TSUM%5[#%&4H5I3,[ 00094400
|
|
TSUM%6[#%&4H6I3,[ 00094500
|
|
TSUM%7[#%&4H7I3,[ 00094600
|
|
TSUM%8[#%&4H8I3,[ 00094700
|
|
PQ%1[#%&4H22X,[ 00094800
|
|
PQ%2[#%&4H19X,[ 00094900
|
|
PQ%3[#%&4H16X,[ 00095000
|
|
PQ%4[#%&4H13X,[ 00095100
|
|
PQ%5[#%&4H10X,[ 00095200
|
|
PQ%6[#%&3H7X,[ 00095300
|
|
PQ%7[#%&3H4X,[ 00095400
|
|
PQ%8[#%&3H1X,[ 00095500
|
|
BLANKS=(+6H ) 00095600
|
|
DO 6 I=1,12 00095700
|
|
6 FT(I)=BLANKS 00095800
|
|
FT%1[#%&1H%[ 00095900
|
|
FT%2[#%&3H1H [ 00096000
|
|
FT%3[#TSUM%NV[ 00096100
|
|
FT%4[#PQ%NV[ 00096200
|
|
FT%5[#%&6HF19.5[[ 00096300
|
|
IN#9-NV 00096400
|
|
DO 5 I#1,NV 00096500
|
|
MM%IN[#LF%I[ 00096600
|
|
5 IN#IN&1 00096700
|
|
IN#9-NV 00096800
|
|
READ (3) %X%I[, I#1,NK[ 00096900
|
|
I3#0 00097000
|
|
10 N%1[#1 00097100
|
|
20 N%2[#1 00097200
|
|
30 N%3[#1 00097300
|
|
40 N%4[#1 00097400
|
|
50 N%5[#1 00097500
|
|
60 N%6[#1 00097600
|
|
70 N%7[#1 00097700
|
|
80 N%8[#1 00097800
|
|
90 I3#I3&1 00097900
|
|
IF%MEAN-1[300,93,92 00098000
|
|
92 WRITE (6, FT)%N%I[, I#IN,8[,X%I3[ 00098100
|
|
93 IF%N%8[-MM%8[[ 110, 120, 120 00098200
|
|
110 N%8[#N%8[&1 00098300
|
|
GO TO 90 00098400
|
|
120 IF%IN-8[ 121, 200, 200 00098500
|
|
121 IF%N%7[-MM%7[[ 122, 130, 130 00098600
|
|
122 N%7[#N%7[&1 00098700
|
|
GO TO 80 00098800
|
|
130 IF%IN-7[ 131, 200, 200 00098900
|
|
131 IF%N%6[-MM%6[[ 132, 140, 140 00099000
|
|
132 N%6[#N%6[&1 00099100
|
|
GO TO 70 00099200
|
|
140 IF%IN-6[ 141, 200, 200 00099300
|
|
141 IF%N%5[-MM%5[[ 142, 150, 150 00099400
|
|
142 N%5[#N%5[&1 00099500
|
|
GO TO 60 00099600
|
|
150 IF%IN-5[ 151, 200, 200 00099700
|
|
GO TO 40 00099800
|
|
162 N%3[#N%3[&1 00099900
|
|
161 IF%N%3[-MM%3[[ 162, 170, 170 00100000
|
|
151 IF%N%4[-MM%4[[ 152, 160, 160 00100100
|
|
152 N%4[#N%4[&1 00100200
|
|
GO TO 50 00100300
|
|
160 IF%IN-4[ 161, 200, 200 00100400
|
|
170 IF%IN-3[ 171, 200, 200 00100500
|
|
171 IF%N%2[-MM%2[[ 172, 180, 180 00100600
|
|
172 N%2[#N%2[&1 00100700
|
|
GO TO 30 00100800
|
|
180 IF%IN-2[ 181, 200, 200 00100900
|
|
181 IF%N%1[-MM%1[[ 182, 200, 200 00101000
|
|
182 N%1[#N%1[&1 00101100
|
|
GO TO 20 00101200
|
|
200 IF%NV-1[ 300, 300, 205 00101300
|
|
205 WRITE (6,902 ) 00101400
|
|
WRITE (6, 903) 00101500
|
|
I2#1 00101600
|
|
DO 220 I#1,NV 00101700
|
|
IN#LF%I[ 00101800
|
|
SUM#NK/IN 00101900
|
|
TRASH#SUM*TMEAN 00102000
|
|
READ (3) TRA 00102100
|
|
TRA#%TRA&TRASH[/SUM 00102200
|
|
WRITE (6, 904) I,I2,TRA 00102300
|
|
DO 210 J#2,IN 00102400
|
|
READ (3) TRA 00102500
|
|
TRA#%TRA&TRASH[/SUM 00102600
|
|
210 WRITE (6, 905) J, TRA 00102700
|
|
WRITE (6, 906 ) 00102800
|
|
220 CONTINUE 00102900
|
|
300 REWIND 3 00103000
|
|
900 FORMAT%13H1PROBLEM NO. A2//[ 00103100
|
|
901 FORMAT%25H0 C E L L N U M B E R S11X,9HM E A N S[ 00103200
|
|
902 FORMAT%1H0/29H0 M A R G I N A L M E A N S[ 00103300
|
|
903 FORMAT%23H VARIABLES CATEGORIES13X,9HM E A N S[ 00103400
|
|
904 FORMAT%1H I6,6X,I6,7X,F19.5[ 00103500
|
|
905 FORMAT%1H 12X,I6,7X,F19.5[ 00103600
|
|
906 FORMAT%1H [ 00103700
|
|
RETURN 00103800
|
|
END 00103900
|
|
CREPORT SUBROUTINE REPORT FOR BMD02V 00104000
|
|
SUBROUTINE REPORT 00104100
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00104200
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00104300
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[ 00104400
|
|
3,ZZ%9,9[,C1%9[,C2%9[ 00104500
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00104600
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00104700
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00104800
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00104900
|
|
M1#NK/I2 00105000
|
|
M2#NK/I3 00105100
|
|
M3#NK/%I2*I3[ 00105200
|
|
A1#M1 00105300
|
|
A2#M2 00105400
|
|
A3#M3 00105500
|
|
DO 10 I#1,I2 00105600
|
|
DO 10 J#1,I3 00105700
|
|
10 ST%IN,I,J[#ST%IN,I,J[/A3 00105800
|
|
DO 15 I#1,I2 00105900
|
|
15 C%1,I[#ST1%KB1,I[/A1 00106000
|
|
DO 16 I#1,I3 00106100
|
|
16 C%2,I[#ST1%KB2,I[/A2 00106200
|
|
TOTAL#TEMP%1[ 00106300
|
|
DO 20 J#1,I3 00106400
|
|
DO 20 I#1,I2 00106500
|
|
20 ST%IN,I,J[#ST%IN,I,J[-C%2,J[-C%1,I[ 00106600
|
|
WRITE (6, 900) 00106700
|
|
WRITE (6, 901) NO1 00106800
|
|
WRITE (6, 902) NO2 00106900
|
|
WRITE (6,920 ) 00107000
|
|
DO 29 K#1,I2 00107100
|
|
29 WRITE (6, 904) %ST%IN,K,J[, J#1,I3[, C%1,K[ 00107200
|
|
WRITE (6,920 ) 00107300
|
|
WRITE (6, 904) %C%2,J[, J#1,I3[, TOTAL 00107400
|
|
900 FORMAT%26H0TABLE OF INTERACTIONS FOR[ 00107500
|
|
901 FORMAT%13H VARIABLEI2,7H %ROW[[ 00107600
|
|
902 FORMAT%13H VARIABLEI2,10H %COLUMN[[ 00107700
|
|
904 FORMAT%F11.5,9F12.5/%10F12.5[[ 00107800
|
|
920 FORMAT%1H [ 00107900
|
|
DO 31 I#1,I2 00108000
|
|
DO 30 J#1,I3 00108100
|
|
30 ZZ%I,J[#ST%IN,I,J[&C%2,J[&C%1,I[&TOTAL 00108200
|
|
31 C1%I[#C%1,I[&TOTAL 00108300
|
|
DO 40 J#1,I3 00108400
|
|
40 C2%J[#C%2,J[&TOTAL 00108500
|
|
WRITE ( 6,905 ) 00108600
|
|
905 FORMAT%15H0TABLE OF MEANS[ 00108700
|
|
WRITE (6,920 ) 00108800
|
|
DO 50 I#1,I2 00108900
|
|
50 WRITE (6,904)%ZZ%I,J[,J#1,I3[,C1%I[ 00109000
|
|
WRITE (6,920 ) 00109100
|
|
WRITE (6,904)%C2%J[,J#1,I3[,TOTAL 00109200
|
|
RETURN 00109300
|
|
END 00109400
|
|
CTABLE SUBROUTINE TABLE FOR BMD02V 9-5-6300109500
|
|
SUBROUTINE TABLE 00109600
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00109700
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00109800
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[ 00109900
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00110000
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00110100
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00110200
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00110300
|
|
GO TO %71,72,73,74,75,76,77[,IJK 00110400
|
|
71 IF%Z-1.0[ 61, 60, 61 00110500
|
|
60 KA#I5 00110600
|
|
GO TO 299 00110700
|
|
61 IF%KA-8[ 143, 200, 200 00110800
|
|
72 IF%Z-1.0[ 63, 62, 63 00110900
|
|
62 KB#I5 00111000
|
|
KA#I5&1 00111100
|
|
GO TO 299 00111200
|
|
63 IF%KA-8[ 143, 93, 93 00111300
|
|
93 IF%KB-7[ 142, 200, 200 00111400
|
|
73 IF%Z-1.0[ 65, 64, 65 00111500
|
|
64 KC#I5 00111600
|
|
KB#I5&1 00111700
|
|
KA#I5&2 00111800
|
|
GO TO 299 00111900
|
|
65 IF%KA-8[ 143, 97, 97 00112000
|
|
97 IF%KB-7[ 142, 98, 98 00112100
|
|
98 IF%KC-6[ 141, 200, 200 00112200
|
|
74 IF%Z-1.0[ 67, 66, 67 00112300
|
|
66 KD#I5 00112400
|
|
KC#I5&1 00112500
|
|
KB#I5&2 00112600
|
|
KA#I5&3 00112700
|
|
GO TO 299 00112800
|
|
67 IF%KA-8[ 143, 103, 103 00112900
|
|
103 IF%KB-7[ 142, 104, 104 00113000
|
|
104 IF%KC-6[ 141, 105, 105 00113100
|
|
105 IF%KD-5[ 140, 200, 200 00113200
|
|
75 IF%Z-1.0[ 69,68,69 00113300
|
|
68 KE#I5 00113400
|
|
KD#I5&1 00113500
|
|
KC#I5&2 00113600
|
|
KB#I5&3 00113700
|
|
KA#I5&4 00113800
|
|
GO TO 299 00113900
|
|
69 IF%KA-8[ 143, 110, 110 00114000
|
|
110 IF%KB-7[ 142, 111, 111 00114100
|
|
111 IF%KC-6[ 141, 112, 112 00114200
|
|
112 IF%KD-5[ 140, 113, 113 00114300
|
|
113 IF%KE-4[ 139, 200, 200 00114400
|
|
76 IF%Z-1.0[ 81, 80, 81 00114500
|
|
80 KF#I5 00114600
|
|
KE#I5&1 00114700
|
|
KD#I5&2 00114800
|
|
KC#I5&3 00114900
|
|
KB#I5&4 00115000
|
|
KA#I5&5 00115100
|
|
GO TO 299 00115200
|
|
81 IF%KA-8[ 143, 120, 120 00115300
|
|
120 IF%KB-7[ 142, 121, 121 00115400
|
|
121 IF%KC-6[ 141, 122, 122 00115500
|
|
122 IF%KD-5[ 140, 123, 123 00115600
|
|
123 IF%KE-4[ 139, 124, 124 00115700
|
|
124 IF%KF-3[ 138, 200, 200 00115800
|
|
77 IF%Z-1.0[ 83, 82, 83 00115900
|
|
82 KG#I5 00116000
|
|
KF#I5&1 00116100
|
|
KE#I5&2 00116200
|
|
KD#I5&3 00116300
|
|
KC#I5&4 00116400
|
|
KB#I5&5 00116500
|
|
KA#I5&6 00116600
|
|
GO TO 299 00116700
|
|
83 IF%KA-8[ 143, 131, 131 00116800
|
|
131 IF%KB-7[ 142, 132, 132 00116900
|
|
132 IF%KC-6[ 141, 133, 133 00117000
|
|
133 IF%KD-5[ 140, 134, 134 00117100
|
|
134 IF%KE-4[ 139, 135, 135 00117200
|
|
135 IF%KF-3[ 138, 136, 136 00117300
|
|
136 IF%KG-2[ 137, 200, 200 00117400
|
|
137 DO 150 I#KG,8 00117500
|
|
150 M%I[#0 00117600
|
|
KG#KG&1 00117700
|
|
M%KG[#1 00117800
|
|
KF#KG&1 00117900
|
|
M%KF[#1 00118000
|
|
10 KE#KF&1 00118100
|
|
M%KE[#1 00118200
|
|
20 KD#KE&1 00118300
|
|
M%KD[#1 00118400
|
|
30 KC#KD&1 00118500
|
|
M%KC[#1 00118600
|
|
40 KB#KC&1 00118700
|
|
M%KB[#1 00118800
|
|
50 KA#KB&1 00118900
|
|
M%KA[#1 00119000
|
|
GO TO 300 00119100
|
|
138 DO 149 I#KF,8 00119200
|
|
149 M%I[#0 00119300
|
|
KF#KF&1 00119400
|
|
M%KF[#1 00119500
|
|
GO TO 10 00119600
|
|
139 DO 148 I#KE,8 00119700
|
|
148 M%I[#0 00119800
|
|
KE#KE&1 00119900
|
|
M%KE[#1 00120000
|
|
GO TO 20 00120100
|
|
140 DO 147 I#KD,8 00120200
|
|
147 M%I[#0 00120300
|
|
KD#KD&1 00120400
|
|
M%KD[#1 00120500
|
|
GO TO 30 00120600
|
|
141 DO 146 I#KC,8 00120700
|
|
146 M%I[#0 00120800
|
|
KC#KC&1 00120900
|
|
M%KC[#1 00121000
|
|
GO TO 40 00121100
|
|
142 DO 145 I#KB,8 00121200
|
|
145 M%I[#0 00121300
|
|
KB#KB&1 00121400
|
|
M%KB[#1 00121500
|
|
GO TO 50 00121600
|
|
143 DO 144 I#KA,8 00121700
|
|
144 M%I[#0 00121800
|
|
KA#KA&1 00121900
|
|
M%KA[#1 00122000
|
|
GO TO 300 00122100
|
|
200 MARK#1 00122200
|
|
GO TO 300 00122300
|
|
299 Z#0.0 00122400
|
|
300 RETURN 00122500
|
|
END 00122600
|
|
CTPWD SUBROUTINE TPWD FOR BMD02V VERSION OF SEPT. 26, 1963 00122700
|
|
SUBROUTINE TPWD%NT1,NT2[ 00122800
|
|
IF%NT1[40,10,12 00122900
|
|
10 NT1#5 00123000
|
|
12 IF%NT1-NT2[14,19,14 00123100
|
|
14 IF%NT2-5[15,19,17 00123200
|
|
15 REWIND NT2 00123300
|
|
GO TO 19 00123400
|
|
17 CONTINUE 00123500
|
|
19 IF%NT1-5[18,24,18 00123600
|
|
18 IF%NT1-6[22,40,22 00123700
|
|
22 REWIND NT1 00123800
|
|
24 NT2#NT1 00123900
|
|
28 RETURN 00124000
|
|
40 WRITE (6,49 ) 00124100
|
|
CALL EXIT 00124200
|
|
49 FORMAT%25H ERROR ON TAPE ASSIGNMENT[ 00124300
|
|
END 00124400
|
|
CTRANS SUBROUTINE TRANS FOR BMD02V 9-5-6300124500
|
|
SUBROUTINE TRANS %DF,NDF,CONST,ONNN,IERROR,I,II[ 00124600
|
|
ASNF%XX[#ATAN %XX/SQRT %1.0-XX**2[[ 00124700
|
|
IERROR#0 00124800
|
|
GO TO %10,20,30,40,50,60,70,80,90,100[,NDF 00124900
|
|
10 IF%DF[99,32,8 00125000
|
|
8 DF#SQRT %DF[ 00125100
|
|
GO TO 200 00125200
|
|
20 IF%DF[99,11,12 00125300
|
|
11 DF#1.0 00125400
|
|
GO TO 200 00125500
|
|
12 DF#SQRT %DF[&SQRT %DF&1.0[ 00125600
|
|
GO TO 200 00125700
|
|
30 IF%DF[99,99,14 00125800
|
|
14 DF#ALOG10%DF[ 00125900
|
|
GO TO 200 00126000
|
|
40 DF#EXP %DF[ 00126100
|
|
GO TO 200 00126200
|
|
50 IF%DF[99,32,17 00126300
|
|
17 IF%DF-1.0[18,19,99 00126400
|
|
19 DF#3.14159265/2.0 00126500
|
|
GO TO 200 00126600
|
|
18 A#SQRT %DF[ 00126700
|
|
DF#ASNF%A[ 00126800
|
|
GO TO 200 00126900
|
|
60 A#DF/%ONNN&1.0[ 00127000
|
|
B#A&1.0/%ONNN&1.0[ 00127100
|
|
IF%A[99,23,24 00127200
|
|
23 IF%B[99,32,27 00127300
|
|
27 DF#ASNF%SQRT %B[[ 00127400
|
|
GO TO 200 00127500
|
|
24 IF%B[99,28,29 00127600
|
|
28 DF#ASNF%SQRT%A[[ 00127700
|
|
GO TO 200 00127800
|
|
29 A#SQRT %A[ 00127900
|
|
B#SQRT %B[ 00128000
|
|
DF#ASNF%A[&ASNF%B[ 00128100
|
|
GO TO 200 00128200
|
|
70 IF%DF[31,99,31 00128300
|
|
31 DF#1.0/DF 00128400
|
|
GO TO 200 00128500
|
|
80 DF#DF&CONST 00128600
|
|
GO TO 200 00128700
|
|
90 DF#DF*CONST 00128800
|
|
GO TO 200 00128900
|
|
100 IF%DF[33,32,33 00129000
|
|
32 DF#0.0 00129100
|
|
GO TO 200 00129200
|
|
33 DF#DF**CONST 00129300
|
|
GO TO 200 00129400
|
|
99 WRITE (6,900)II,I,DF 00129500
|
|
IERROR#-99 00129600
|
|
900 FORMAT%51H0ERROR OCCURRED DURING TRANS-GENERATION PASS NUMBERI3, 00129700
|
|
120H, REPLICATION NUMBERI5,1H./26H VALUE OF THIS REPLICATE #F15.5,200129800
|
|
2H. 36HTHIS IS THE FIRST ERROR ENCOUNTERED./41H PROGRAM WILL GO TO 00129900
|
|
3NEXT PROBLEM, IF ANY.[ 00130000
|
|
200 RETURN 00130100
|
|
END 00130200
|
|
CVFCHCK SUBROUTINE TO CHECK FOR PROPER NUMBER OF VARIABLE FORMAT CRDS00130300
|
|
SUBROUTINE VFCHCK%NVF[ 00130400
|
|
IF%NVF[10,10,20 00130500
|
|
10 WRITE (6,4000) 00130600
|
|
NVF#1 00130700
|
|
50 RETURN 00130800
|
|
C 00130900
|
|
20 IF%NVF-5[50,50,10 00131000
|
|
C 00131100
|
|
4000 FORMAT%1H023X71HNUMBER OF VARIABLE FORMAT CARDS INCORRECTLY SPECIF00131200
|
|
XIED, ASSUMED TO BE 1.[ 00131300
|
|
END 00131400
|
|
CXLOCAT SUBROUTINE XLOCAT FOR BMD02V 9-5-6300131500
|
|
SUBROUTINE XLOCAT %MEAN[ 00131600
|
|
DIMENSION X%18000[,L%9[,TEMP%12[,LO%8[,LM%9[,LL%8[,K5%9[,N%8[ 00131700
|
|
1,KM%8[,MM%9[,MI%8[,ST%6,9,9[,ST1%4,9[,M%8[,TSUM%255[,PQ%255[,MD%8[00131800
|
|
2,NORD%4[,NSIZE%4[,LF%8[,C%9,9[,CO%9,8,9[,D%9,9[ 00131900
|
|
COMMON X,L,TEMP,LO,LM,LL,K5,N,KM,MM,MI,ST,ST1,M,TSUM,PQ,MD,NORD, 00132000
|
|
1NSIZE,LF,C,CO,D,SUM,NV,NR,NK,NX,NVM,I5,MT,MT1,MARK,JW,Z,TRASH,TRA,00132100
|
|
2IJK,NO,IN,I2,I3,NDF,TOTAL,NO1,NO2,KB1,KB2, KA,KB,KC,KD,KE, 00132200
|
|
3KF,KG,IA,IB,IC,ID,IE,IG 00132300
|
|
DO 300 I#1,8 00132400
|
|
300 N%I[#0 00132500
|
|
KI#0 00132600
|
|
K#1 00132700
|
|
DO 302 I#I5,8 00132800
|
|
KI#KI&1 00132900
|
|
IF%M%I[[ 301, 301, 302 00133000
|
|
301 N%K[#KI 00133100
|
|
K#K&1 00133200
|
|
302 CONTINUE 00133300
|
|
DO 303 I#I5,8 00133400
|
|
MM%I[#L%I[-1 00133500
|
|
303 LM%I[#LO%I[ 00133600
|
|
DO 305 I#1,NV 00133700
|
|
K#I5&%I-1[ 00133800
|
|
305 K5%K[#I 00133900
|
|
DO 310 I#1,JW 00134000
|
|
DO 309 J#I5,8 00134100
|
|
IF%N%I[-K5%J[[ 309, 306, 309 00134200
|
|
306 K5%9[#K5%J[ 00134300
|
|
LM%9[#LM%J[ 00134400
|
|
MM%9[#MM%J[ 00134500
|
|
DO 307 J2#J,8 00134600
|
|
K5%J2[#K5%J2&1[ 00134700
|
|
LM%J2[#LM%J2&1[ 00134800
|
|
307 MM%J2[#MM%J2&1[ 00134900
|
|
GO TO 310 00135000
|
|
309 CONTINUE 00135100
|
|
310 CONTINUE 00135200
|
|
SUM#0.0 00135300
|
|
TRASH#0.0 00135400
|
|
351 N%1[#0 00135500
|
|
352 N%2[#0 00135600
|
|
353 N%3[#0 00135700
|
|
354 N%4[#0 00135800
|
|
355 N%5[#0 00135900
|
|
356 N%6[#0 00136000
|
|
357 N%7[#0 00136100
|
|
358 N%8[#0 00136200
|
|
359 DO 360 J#I5,8 00136300
|
|
KK#K5%J[ 00136400
|
|
360 KM%KK[#N%J[ 00136500
|
|
LX#1 00136600
|
|
DO 361 J#1,NV 00136700
|
|
361 LX#LX&KM%J[*LL%J[ 00136800
|
|
TRASH#TRASH&X%LX[ 00136900
|
|
IF%N%8[-MM%8[[ 362, 363, 363 00137000
|
|
362 N%8[#N%8[&1 00137100
|
|
GO TO 359 00137200
|
|
363 IF%JW-1[ 364,364,370 00137300
|
|
364 SUM#SUM&TRASH*TRASH 00137400
|
|
TRA#TRASH 00137500
|
|
TRASH#0.0 00137600
|
|
IF%MEAN[ 705, 705, 701 00137700
|
|
701 IF%JW-NVM[ 705, 702, 705 00137800
|
|
702 WRITE (3) TRA 00137900
|
|
705 IF%NV-3[ 365, 367, 370 00138000
|
|
365 IF%LM%7[-100[ 370,370,366 00138100
|
|
366 N1#N%7[&1 00138200
|
|
ST1%MT1,N1[#TRA 00138300
|
|
MI%1[#1 00138400
|
|
GO TO 370 00138500
|
|
367 IF%LM%6[-100[ 370,370, 368 00138600
|
|
368 IF%LM%7[-100[ 370, 370, 369 00138700
|
|
369 N1#N%6[&1 00138800
|
|
N2#N%7[&1 00138900
|
|
ST%MT,N1,N2[#TRA 00139000
|
|
MI%1[#2 00139100
|
|
370 IF%N%7[-MM%7[[ 371,372,372 00139200
|
|
371 N%7[#N%7[&1 00139300
|
|
GO TO 358 00139400
|
|
372 IF%JW-2[ 373,375,381 00139500
|
|
373 SUM#SUM&TRASH*TRASH 00139600
|
|
IF%MI%1[-1[ 501, 374, 381 00139700
|
|
501 IF%NV-2[ 70, 70, 381 00139800
|
|
374 MI%1[#0 00139900
|
|
MT1#MT1&1 00140000
|
|
GO TO 70 00140100
|
|
375 SUM#SUM&TRASH*TRASH 00140200
|
|
TRA#TRASH 00140300
|
|
TRASH#0.0 00140400
|
|
IF%MEAN[ 715, 715, 711 00140500
|
|
711 IF%JW-NVM[ 715, 712, 715 00140600
|
|
712 WRITE (3) TRA 00140700
|
|
715 IF%NV-4[ 376,378,381 00140800
|
|
376 IF%LM%6[-100[ 381,381,377 00140900
|
|
377 N1#N%6[&1 00141000
|
|
ST1%MT1,N1[#TRA 00141100
|
|
MI%2[#1 00141200
|
|
GO TO 381 00141300
|
|
378 IF%LM%5[-100[ 381,381,379 00141400
|
|
379 IF%LM%6[-100[ 381,381,380 00141500
|
|
380 N1#N%5[&1 00141600
|
|
N2#N%6[&1 00141700
|
|
ST%MT,N1,N2[#TRA 00141800
|
|
MI%2[#2 00141900
|
|
381 IF%N%6[-MM%6[[ 382,383, 383 00142000
|
|
382 N%6[#N%6[&1 00142100
|
|
GO TO 357 00142200
|
|
383 IF%JW-3[ 384, 388, 394 00142300
|
|
384 SUM#SUM&TRASH*TRASH 00142400
|
|
IF%MI%1[-2[ 385, 387, 387 00142500
|
|
385 IF%MI%2[-1[ 502, 386, 394 00142600
|
|
502 IF%NV-3[ 70, 70, 394 00142700
|
|
386 MI%2[#0 00142800
|
|
MT1#MT1&1 00142900
|
|
GO TO 70 00143000
|
|
387 MI%1[#0 00143100
|
|
MT#MT&1 00143200
|
|
GO TO 70 00143300
|
|
388 SUM#SUM&TRASH*TRASH 00143400
|
|
TRA#TRASH 00143500
|
|
TRASH#0.0 00143600
|
|
IF%MEAN[ 725, 725, 721 00143700
|
|
721 IF%JW-NVM[ 725, 722, 725 00143800
|
|
722 WRITE (3) TRA 00143900
|
|
725 IF%NV-5[ 389, 391,394 00144000
|
|
389 IF%LM%5[-100[ 394, 394, 390 00144100
|
|
390 N1#N%5[&1 00144200
|
|
ST1%MT1,N1[#TRA 00144300
|
|
MI%3[#1 00144400
|
|
GO TO 394 00144500
|
|
391 IF%LM%4[-100[ 394, 394, 392 00144600
|
|
392 IF%LM%5[-100[ 394, 394, 393 00144700
|
|
393 N1#N%4[&1 00144800
|
|
N2#N%5[&1 00144900
|
|
ST%MT,N1,N2[#TRA 00145000
|
|
MI%3[#2 00145100
|
|
394 IF%N%5[-MM%5[[ 395,396,396 00145200
|
|
395 N%5[#N%5[&1 00145300
|
|
GO TO 356 00145400
|
|
396 IF%JW-4[ 397, 401, 407 00145500
|
|
397 SUM#SUM&TRASH*TRASH 00145600
|
|
IF%MI%2[-2[ 398,400, 400 00145700
|
|
398 IF%MI%3[-1[ 503, 399, 407 00145800
|
|
503 IF%NV-4[ 70, 70, 407 00145900
|
|
399 MI%3[#0 00146000
|
|
MT1#MT1&1 00146100
|
|
GO TO 70 00146200
|
|
400 MI%2[#0 00146300
|
|
MT#MT&1 00146400
|
|
GO TO 70 00146500
|
|
401 SUM#SUM&TRASH*TRASH 00146600
|
|
TRA#TRASH 00146700
|
|
TRASH#0.0 00146800
|
|
IF%MEAN[ 735, 735, 731 00146900
|
|
731 IF%JW-NVM[ 735, 732, 735 00147000
|
|
732 WRITE (3) TRA 00147100
|
|
735 IF%NV-6[ 402, 404, 407 00147200
|
|
402 IF%LM%4[-100[ 407,407,403 00147300
|
|
403 N1#N%4[&1 00147400
|
|
ST1%MT1,N1[#TRA 00147500
|
|
MI%4[#1 00147600
|
|
GO TO 407 00147700
|
|
404 IF%LM%3[-100[ 407,407,405 00147800
|
|
405 IF%LM%4[-100[ 407,407,406 00147900
|
|
406 N1#N%3[&1 00148000
|
|
N2#N%4[&1 00148100
|
|
ST%MT,N1,N2[#TRA 00148200
|
|
MI%4[#2 00148300
|
|
407 IF%N%4[-MM%4[[ 408,409,409 00148400
|
|
408 N%4[#N%4[&1 00148500
|
|
GO TO 355 00148600
|
|
409 IF%JW-5[ 410,414,420 00148700
|
|
410 SUM#SUM&TRASH*TRASH 00148800
|
|
IF%MI%3[-2[ 411,413,413 00148900
|
|
411 IF%MI%4[-1[ 504, 412, 420 00149000
|
|
504 IF%NV-5[ 70, 70, 420 00149100
|
|
412 MI%4[#0 00149200
|
|
MT1#MT1&1 00149300
|
|
GO TO 70 00149400
|
|
413 MI%3[#0 00149500
|
|
MT#MT&1 00149600
|
|
GO TO 70 00149700
|
|
414 SUM#SUM&TRASH*TRASH 00149800
|
|
TRA#TRASH 00149900
|
|
TRASH#0.0 00150000
|
|
IF%MEAN[ 745, 745, 741 00150100
|
|
741 IF%JW-NVM[ 745, 742, 745 00150200
|
|
742 WRITE (3) TRA 00150300
|
|
745 IF%NV-7[ 415, 417, 420 00150400
|
|
415 IF%LM%3[-100[ 420, 420, 416 00150500
|
|
416 N1#N%3[&1 00150600
|
|
ST1%MT1,N1[#TRA 00150700
|
|
MI%5[#1 00150800
|
|
GO TO 420 00150900
|
|
417 IF%LM%2[-100[ 420,420,418 00151000
|
|
418 IF%LM%3[-100[ 420,420,419 00151100
|
|
419 N1#N%2[&1 00151200
|
|
N2#N%3[&1 00151300
|
|
ST%MT,N1,N2[#TRA 00151400
|
|
MI%5[#2 00151500
|
|
420 IF%N%3[-MM%3[[ 421,422,422 00151600
|
|
421 N%3[#N%3[&1 00151700
|
|
GO TO 354 00151800
|
|
422 IF%JW-6[ 423,427,433 00151900
|
|
423 SUM#SUM&TRASH*TRASH 00152000
|
|
IF%MI%4[-2[ 424,426,426 00152100
|
|
424 IF%MI%5[-1[ 505, 425, 433 00152200
|
|
505 IF%NV-6[ 70, 70, 433 00152300
|
|
425 MI%5[#0 00152400
|
|
MT1#MT1&1 00152500
|
|
GO TO 70 00152600
|
|
426 MI%4[#0 00152700
|
|
MT#MT&1 00152800
|
|
GO TO 70 00152900
|
|
427 SUM#SUM&TRASH*TRASH 00153000
|
|
TRA#TRASH 00153100
|
|
TRASH#0.0 00153200
|
|
IF%MEAN[ 755, 755, 751 00153300
|
|
751 IF%JW-NVM[ 755, 752, 755 00153400
|
|
752 WRITE (3) TRA 00153500
|
|
755 IF%NV-8[ 428,430,430 00153600
|
|
428 IF%LM%2[-100[ 433,433,429 00153700
|
|
429 N1#N%2[&1 00153800
|
|
ST1%MT1,N1[#TRA 00153900
|
|
MI%6[#1 00154000
|
|
GO TO 433 00154100
|
|
430 IF%LM%1[-100[ 433,433,431 00154200
|
|
431 IF%LM%2[-100[ 433,433,432 00154300
|
|
432 N1#N%1[&1 00154400
|
|
N2#N%2[&1 00154500
|
|
ST%MT,N1,N2[#TRA 00154600
|
|
MI%6[#2 00154700
|
|
433 IF%N%2[-MM%2[[ 434,435,435 00154800
|
|
434 N%2[#N%2[&1 00154900
|
|
GO TO 353 00155000
|
|
435 IF%JW-7[ 436,440,440 00155100
|
|
436 SUM#SUM&TRASH*TRASH 00155200
|
|
IF%MI%5[-2[ 437,439,439 00155300
|
|
437 IF%MI%6[-1[ 506, 438, 443 00155400
|
|
506 IF%NV-7[ 70, 70, 443 00155500
|
|
438 MI%6[#0 00155600
|
|
MT1#MT1&1 00155700
|
|
GO TO 70 00155800
|
|
439 MI%5[#0 00155900
|
|
MT#MT&1 00156000
|
|
GO TO 70 00156100
|
|
440 SUM#SUM&TRASH*TRASH 00156200
|
|
TRA#TRASH 00156300
|
|
TRASH#0.0 00156400
|
|
IF%MEAN[ 765, 765, 761 00156500
|
|
761 IF%JW-NVM[ 765, 762, 765 00156600
|
|
762 WRITE (3) TRA 00156700
|
|
765 IF%NV-8[ 70,441,441 00156800
|
|
441 IF%LM%1[-100[ 443,443,442 00156900
|
|
442 N1#N%1[&1 00157000
|
|
ST1%MT1,N1[#TRA 00157100
|
|
MI%7[#1 00157200
|
|
443 IF%N%1[-MM%1[[ 444,445,445 00157300
|
|
444 N%1[#N%1[&1 00157400
|
|
GO TO 352 00157500
|
|
445 IF%MI%6[-2[ 446,448,448 00157600
|
|
446 IF%MI%7[-1[70,447,447 00157700
|
|
447 MI%7[#0 00157800
|
|
MT1#MT1&1 00157900
|
|
GO TO 70 00158000
|
|
448 MI%6[#0 00158100
|
|
MT#MT&1 00158200
|
|
70 RETURN 00158300
|
|
END 00158400
|