Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-01 - 43,50144/multry.bas
There are 2 other files named multry.bas in the archive. Click here to see a list.
10 DIMV(8,36),O(8,36)
11 DIMD(7),E(7)
12 DIMF(6),G(6),U(6),H(6),B(6),S(6,6),T(6,6)
15 PRINT
16 PRINT"MODIFIED 07/24/68"
17 PRINT
18 PRINT
20 LETR=0
21 LETR=R+1
22 LETF(R)=0
23 LETG(R)=0
24 LETU(R)=1
25 LETH(R)=0
26 LETB(R)=0
27 LETJ=0
28 LETJ=J+1
29 LETS(R,J)=0
30 LETT(R,J)=0
31 IFJ<6THEN28
32 IFR<6THEN21
33 LETI=0
34 LETJ=0
35 LETJ=J+1
36 LETV(I,J)=1E37
37 LETO(I,J)=1E37
38 IFJ<36THEN35
39 LETI=I+1
40 IFI<8THEN34
41 LETI=0
42 LETE(I)=0
43 LETD(I)=0
44 LETI=I+1
45 IFI>7THEN42
46 LETL=0
47 LETE(0)=99
48 LETI=0
50 LETI=I+1
60 READV(0,I)
70 IFV(0,I)<>1E37THEN50
80 LETP=I-1
100 GOTO2000
109 DATA 1E37
110 LETU(1)=0
120 LETU(2)=0
130 LETU(3)=0
140 LETU(4)=0
150 LETU(5)=0
160 LETU(6)=0
170 FORJ=6TO1STEP-1
171 IFU(J)=1THEN173
172 NEXTJ
173 LETQ=J
174 FORK=1TOQ
175 FORI=1TOP
176 READO(K,I)
177 NEXTI
178 NEXTK
179 LETI=0
180 LETI=I+1
181 LETX1=O(1,I)
182 LETX2=O(2,I)
183 LETX3=O(3,I)
184 LETX4=O(4,I)
185 LETX5=O(5,I)
186 LETX6=O(6,I)
187 GOSUB210
188 LETO(1,I)=W1
189 LETO(2,I)=W2
190 LETO(3,I)=W3
191 LETO(4,I)=W4
192 LETO(5,I)=W5
193 LETO(6,I)=W6
194 IFI>PTHEN200
195 GOTO180
200 LETL=1
201 READM
202 IFL<>1THEN205
203 LETN1=Q
204 GOTO435
205 LETN1=M
206 GOTO435
210 LETW1=X1
220 LETW2=X2
230 LETW3=X3
240 LETW4=X4
250 LETW5=X5
260 LETW6=X6
270 RETURN
435 LETN=N1
440 FORI=1TO6
441 IFI<=NTHEN444
442 LETD(I)=0
443 GOTO445
444 LETD(I)=I
445 NEXTI
446 LETQ1=1
447 LETY=7
448 GOSUB1200
449 GOSUB1300
450 GOTO495
451 LETQ1=0
495 FORI=0TON
500 LETF(I)=0
505 FORJ=0TON
510 LETS(I,J)=0
515 NEXTJ
520 NEXTI
525 FORK=1TOP
530 FORI=0TON
535 LETF(I)=F(I)+V(I,K)
540 FORJ=0TON
545 LETS(I,J)=S(I,J)+V(I,K)*V(J,K)
550 NEXTJ
555 NEXTI
560 NEXTK
561 FORI=0TON
562 LETF(I)=F(I)/P
563 NEXTI
565 FORI=0TON
570 FORJ=0TON
575 LETT(I,J)=S(I,J)-P*F(I)*F(J)
580 NEXTJ
585 NEXTI
590 FORJ=0TON
595 LETG(J)=SQR(T(J,J))/P
600 NEXTJ
610 GOSUB1000
620 FORI=1TON
625 LETB(I)=0
630 FORJ=1TON
635 LETB(I)=B(I)+T(0,J)*T(I,J)
640 NEXTJ
645 NEXTI
646 LETJ=0
647 LETJ=J+1
648 LETH(J)=B(J)/(G(J)+1E-36)
649 IFJ<NTHEN647
650 LETA1=F(0)
655 LETS2=0
660 FORI=1TON
665 LETA1=A1-B(I)*F(I)
670 NEXTI
675 FORI=1TOP
680 LETV(7,I)=A1
685 FORJ=1TON
690 LETV(7,I)=V(7,I)+B(J)*V(J,I)
693 NEXTJ
695 LETV(8,I)=V(7,I)-V(0,I)
700 LETS2=S2+V(8,I)^2
710 NEXTI
715 LETS2=S2/(P-N-1)
720 LETW4=(P-N-1)*S2
725 LETW5=(P*G(0))^2
730 LETW3=1-(W4/W5)
735 LETW1=((P*G(0))^2)/(P-1)
740 IFS2>=W1THEN755
745 LETS2=S2
750 GOTO760
755 LETS2=W1
760 LETR2=1-(S2/W1)
765 LETS2=SQR(S2)
770 FORJ=0TON
775 LETG(J)=G(J)*SQR(P)
780 NEXTJ
785 IFQ1<>0THEN800
790 GOSUB1400
795 GOTO835
800 PRINT
805 PRINT"INDEXES OF DETERMINATION:"
810 PRINT"(ADJUSTED)","(UNADJUSTED)","   VARIABLES IN ANALYSIS"
815 PRINT
820 GOSUB1400
825 LETN=N-1
830 LETI1=0
835 LETI1=I1+1
840 IFI1=7THEN870
845 IFD(I1)=0THEN835
850 LETY=I1
855 GOSUB1200
860 GOSUB1300
865 GOTO451
870 IFN=1THEN905
875 PRINT
880 PRINT"     WHICH VARIABLE IS TO BE ELIMINATED";
885 INPUTQ
890 IFQ=0THEN905
895 LETD(Q)=0
900 GOTO446
905 PRINT
910 PRINT
915 PRINT"FOR DETAILED OUTPUT, ENTER DATA INTO 'MULREG' PROGRAM,"
920 PRINT"AND RUN THE SELECTED COMBINATION OF VARIABLES."
925 GOTO5000
1000 FORK=1TON
1010 FORJ=1TON
1030 IFJ=KTHEN1050
1040 LETT(K,J)=T(K,J)/T(K,K)
1050 NEXTJ
1060 LETT(K,K)=1/T(K,K)
1070 FORI=1TON
1080 IFI=KTHEN1140
1090 FORJ=1TON
1100 IFJ=KTHEN1120
1110 LETT(I,J)=T(I,J)-T(I,K)*T(K,J)
1120 NEXTJ
1130 LETT(I,K)=-(T(I,K)*T(K,K))
1140 NEXTI
1150 NEXTK
1160 RETURN
1200 LETI=0
1201 LETI=I+1
1202 IFI>6THEN1220
1210 LETE(I)=D(I)
1215 GOTO1201
1220 LETE(Y)=0
1225 LETK1=0
1230 LETK2=0
1235 LETK1=K1+1
1240 IFE(K1)=0THEN1235
1245 LETK2=K2+1
1250 LETE(K2)=E(K1)
1255 IFK2<NTHEN1235
1260 RETURN
1300 FORI=1TOP
1310 FORJ=1TON
1320 LETV(J,I)=O(E(J),I)
1330 NEXTJ
1340 NEXTI
1350 RETURN
1400 PRINTR2,W3,"   ";
1410 FORI=1TON
1420 PRINTE(I);
1430 NEXTI
1440 PRINT
1450 RETURN
2000 PRINT
2100 PRINT"LIST 'MULTEX***' FOR INSTRUCTIONS ON USING THIS PROGRAM"
4999 DATA 0
5000 END