PROGRAM INVEST INTEGER SEL 43(150),SEL 44(150) DOUBLE PRECISION KA 1(6, 9), KB 1A(6, 8), KC 1A(6, 8), KD 1(6, 4), * KB 1B(6, 8), KC 1B(6, 8), * KA 2(6, 9), KB 2 (6, 7), KC 2 (6, 7), KD 2(6, 9), * KA 3(6, 8), KB 3 (6, 7), KC 3 (6, 7), KD 3(6, 8), * KA 4(6,11), KB 4 (6,10), KC 4 (6,10), KD 4(6,11), * KD 5(6, 7), * KA 6(6,18), KB 6 (6,17), KC 6 (6,17), KD 6(6,18), * KA 7(6, 9), KB 7 (6, 8), KC 7 (6, 8), KD 7(6, 9), * KA 8(6, 9), KB 8 (6, 8), KC 8 (6, 8), KD 8(6, 9), * KA 9(6,12), KB 9 (6,11), KC 9 (6,11), KD 9(6,12), * KA10(6, 4), KB10 (6, 3), KC10 (6, 3), KD10(6, 4), * KA11(6,17), KB11 (6,17), KC11 (6,17), KD11(6,17), * KA12(6, 7), KB12 (6, 7), KC12 (6, 7), KD12(6, 7), * KA13(6, 4), KB13 (6, 3), KC13 (6, 3), KD13(6, 4), * KA14(6, 4), KB14A(6, 3), KC14A(6, 3), KD14(6, 4), * KB14B(6, 3), KC14B(6, 3), * KA15(6, 3), KB15 (6, 3), KC15 (6, 3), KD15(6, 3), * KA16(6, 7), KB16 (6, 6), KC16 (6, 6), KD16(6, 7), * KA17(6, 7), KB17 (6, 7), KC17 (6, 7), KD17(6, 7) DOUBLE PRECISION KA18(6,18), KB18 (6,18), KC18 (6,18), KD18(6,18), * KA19(6,18), KB19 (6,18), KC19 (6,18), * KA20(6,12), KB20 (6,12), KC20 (6,12) DOUBLE PRECISION KA21(6, 6), KB21 (6, 5), KC21 (6, 5), KD21(6, 6), * KA22(6,10), KB22 (6, 9), KC22 (6, 9), KD22(6,10), * KA23(6, 6), KB23 (6, 5), KC23 (6, 5), KD23(6, 6), * KA24(6, 8), KB24 (6, 7), KC24 (6, 7), KD24(6, 8), * KA25(6,10), KB25 (6, 9), KC25 (6, 9), KA32(6, 6), * KA26(6, 6), KB26 (6, 6), KC26 (6, 6), KA33(6, 4), * KA27(6, 5), KB27 (6, 4), KC27 (6, 4), KA34(6, 6), * KA28(6, 5), KB28 (6, 4), KC28 (6, 4), KA35(6, 6), * KA29(6, 1), KB29 (6, 1), KC29 (6, 1), KA36(6, 6), * KA30(6, 6), KB30 (6, 6), KC30 (6, 6), KA37(6, 9), * KA31(6, 6), KB31 (6, 6), KC31 (6, 6), KA38(6,13) DOUBLE PRECISION KX (6,18), SM43 (5, 3), SM44 (5, 3) DOUBLE PRECISION TBM(20, 40),TRM(20, 40), SBK( 6, 25),SRK( 6, 25), * IVH(19, 21),IQH(19, 21), IHC( 6,250),IHD( 6,250), * IVR(19, 21),IQR(19, 21), IMC(20,250),IMD(20,250), * IVM(19, 21),IQM(19, 21), IXC(20,250),IXD(20,250), * IDM(19, 21),IDR(19, 21), ITC(20,250), * ITM(19, 21),IVT(19, 21) DOUBLE PRECISION INV(19, 66, 6), BLD(19, 66, 6),PLT(19, 66, 6), * GCF(19, 66, 10), DEP(19, 66, 10),KAP(19, 66, 10), * AUX(20,196) DOUBLE PRECISION KVC(20,250),KVD(20,250), MVC(20,250),MVD(20,250), * KPC(20,250),KPD(20,250), XVC(20,250),XVD(20,250), * PVC(20,250),TVC(20,250), SVC(20,250),SVD(20,250), * PCK( 6,250),PXK( 6,250) DOUBLE PRECISION MCA(20,208),MCB(20,208), RCA(20,208), * MFA(20, 71), RFA(20, 71), * MGA(20,215),MGB(20,215), RGA(20,215), * MSA(20,121),MSB(20,121), MSC(20,121),RSA(20,121), * MUA(20,324),MUB(20,324), RUA(20,324) DOUBLE PRECISION XCA(20,182),XCB(20,182), ZCA(20,182), * XFA(20, 71), ZFA(20, 71), * XGA(20,216),XGB(20,216), ZGA(20,216), * XSA(20,123),XSB(20,123), XSC(20,123),ZSA(20,123), * XUA(20,308),XUB(20,308), XUC(20,308),ZUA(20,308) DOUBLE PRECISION KA (6,18,38), KD (6,18,38), M(20,939,4), * KBA(6,18,38), KBB(6,18,38), X(20,900,4), * KCA(6,18,38), KCB(6,18,38), T(20,738 ), * PKA(6,18,38), PKB(6,18,38) COMMON KA, KBA, KBB, KCA, KCB, KD, PKA, PKB, X,M, XUC, T, * RCA, RFA, RGA, RSA, RUA, ZCA, ZFA, ZGA, ZSA, ZUA, TVC, * MCA, MFA, MGA, MSA, MUA, XCA, XFA, XGA, XSA, XUA, SVC, * MCB, MSC, MGB, MSB, MUB, XCB, XSC, XGB, XSB, XUB, SVD, * KVC, KVD, KPC, KPD, PCK, MVC, MVD, XVC, XVD, PXK, PVC COMMON /DATDEC/ KA 1 , KA 2, KA 3, KA 4, KA 6 , KA 7, KA 8, * KB 1A, KB 2, KB 3, KB 4, KB 6 , KB 7, KB 8, * KC 1A, KC 2, KC 3, KC 4, KC 6 , KC 7, KC 8, * KD 1 , KD 2, KD 3, KD 4, KD 5, KD 6 , KD 7, KD 8, * KA 9 , KA10, KA11, KA12, KA13, KA14 , KA15, KA16, * KB 9 , KB10, KB11, KB12, KB13, KB14A, KB15, KB16, * KC 9 , KC10, KC11, KC12, KC13, KC14A, KC15, KC16, * KD 9 , KD10, KD11, KD12, KD13, KD14 , KD15, KD16, * KA17 , KA18, KA19, KA20, KA21, KA22 , KA23, KA24, * KB17 , KB18, KB19, KB20, KB21, KB22 , KB23, KB24, * KC17 , KC18, KC19, KC20, KC21, KC22 , KC23, KC24, * KD17 , KD18, KD21, KD22 , KD23, KD24, * KA25 , KA26, KA27, KA28, KA29, KA30 , KA31, KA32, * KB25 , KB26, KB27, KB28, KB29, KB30 , KB31, KA33, * KC25 , KC26, KC27, KC28, KC29, KC30 , KC31, KA34, * KB 1B, KA35, KA36, KA37, KA38, KB14B, * KC 1B, SM43, SM44, KC14B c these calls inserted to assign values to mca and rca call zeroer(mca,208) call zeroer(rca,208) CALL ZEROER (MCB,208) call COPYM2 (M, MCA,MCB, RCA, 1,208,939) write(6,10) 10 format (1h ,'f029 passed') END SUBROUTINE ZEROER (MPORT,NJT) DOUBLE PRECISION MPORT(20,NJT) DOUBLE PRECISION KVC(20,250),KVD(20,250), MVC(20,250),MVD(20,250), * KPC(20,250),KPD(20,250), XVC(20,250),XVD(20,250), * PVC(20,250),TVC(20,250), SVC(20,250),SVD(20,250), * PCK( 6,250),PXK( 6,250) DOUBLE PRECISION MCA(20,208),MCB(20,208), RCA(20,208), * MFA(20, 71), RFA(20, 71), * MGA(20,215),MGB(20,215), RGA(20,215), * MSA(20,121),MSB(20,121), MSC(20,121),RSA(20,121), * MUA(20,324),MUB(20,324), RUA(20,324) DOUBLE PRECISION XCA(20,182),XCB(20,182), ZCA(20,182), * XFA(20, 71), ZFA(20, 71), * XGA(20,216),XGB(20,216), ZGA(20,216), * XSA(20,123),XSB(20,123), XSC(20,123),ZSA(20,123), * XUA(20,308),XUB(20,308), XUC(20,308),ZUA(20,308) DOUBLE PRECISION KA (6,18,38), KD (6,18,38), M(20,939,4), * KBA(6,18,38), KBB(6,18,38), X(20,900,4), * KCA(6,18,38), KCB(6,18,38), T(20,738 ), * PKA(6,18,38), PKB(6,18,38) COMMON KA, KBA, KBB, KCA, KCB, KD, PKA, PKB, X,M, XUC, T, * RCA, RFA, RGA, RSA, RUA, ZCA, ZFA, ZGA, ZSA, ZUA, TVC, * MCA, MFA, MGA, MSA, MUA, XCA, XFA, XGA, XSA, XUA, SVC, * MCB, MSC, MGB, MSB, MUB, XCB, XSC, XGB, XSB, XUB, SVD, * KVC, KVD, KPC, KPD, PCK, MVC, MVD, XVC, XVD, PXK, PVC DO 10, I=1, 20 DO 10, J=1,NJT 10 MPORT(I,J)=0.0D0 RETURN END C Procedure to collect double-quantity trade data into single arrays: SUBROUTINE COPYM2 (TRADES, QTRAD1, QTRAD2, VTRADE, N1, N2, NT) DOUBLE PRECISION QTRAD1(20,N2-N1+1), TRADES(20,NT, 4), * QTRAD2(20,N2-N1+1), VTRADE(20,N2-N1+1) INTEGER N1,N2,NT DO 100, I=1,20 DO 100, J=1,N2-N1+1 TRADES(I,J+N1-1,1)=VTRADE(I,J) TRADES(I,J+N1-1,2)=QTRAD1(I,J) 100 TRADES(I,J+N1-1,3)=QTRAD2(I,J) RETURN END