SUBROUTINE VECTW2(INDCAS,W,R1,R2) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION W(*) GOTO (10,20,30,40,50),INDCAS C ------------------------------------------- C CASE R1=R2=1/3 (H1=H2=H3) C ---------------------------- 10 CONTINUE W(1)= 0.005735253586193482089659 W(2)=-0.248356781363971259867437 W(3)= 1.725694444444444444444446 W(4)=-1.991056224423800892964037 W(5)= 1.265361779979356448519593 W(6)=-2.740451388888888888888889 W(7)= 2.361724033538074685924453 W(8)=-0.504519172426963574813342 W(9)= 1.1258680555555555555555 RETURN C ------------------------------- C CASE R1=R2=R (H1=H2 AND H3#H1) C ------------------------------- 20 CONTINUE R=R1 T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R-1 T5 = R**2 T6 = T5**2 T7 = T6*T5 T8 = T6*R T12 = T5*R T20 = 4*T5-4*R+1 T22 = 1/(66*T8-291*T6+429*T12-293*T5+93*R-12) T24 = 1/T12 T25 = T20*T22*T24 T27 = 16+T1 T45 = T6**2 T46 = T6*T12 T62 = 1/T3 T64 = T24*T62*T22 T97 = T62*T22 W(1) = -T2*T3*(165*T7-1140*T8+(438*T1+1758)*T6+(-978*T1-968)*T12+( #843*T1+23)*T5+(-318*T1+132)*R+45*T1-30)*T25/360 W(2) = T27*(165*T7-1140*T8+(-438*T1+1758)*T6+(978*T1-968)*T12+(-84 #3*T1+23)*T5+(318*T1+132)*R-45*T1-30)*T3*T25/360 W(3) = T3*T20*(33*T7-228*T8-174*T6+980*T12-1007*T5+408*R-60)*T22*T #24/18 W(4) = -T2*(165*T45-1470*T46+(42*T1+3387)*T7+(-3656-306*T1)*T8+(84 #3*T1+3543)*T6+(-1158*T1-3738)*T12+(2641+846*T1)*T5+(-912-312*T1)*R #+45*T1+120)*T20*T64/360 W(5) = T27*(165*T45-1470*T46+(-42*T1+3387)*T7+(306*T1-3656)*T8+(35 #43-843*T1)*T6+(1158*T1-3738)*T12+(2641-846*T1)*T5+(-912+312*T1)*R+ #120-45*T1)*T20*T64/360 W(6) = T20*(33*T45-294*T46+627*T7-364*T8-303*T6+642*T12-487*T5-30+ #192*R)*T64/18 W(7) = T2*(132*T46-1374*T7+(4422+12*T1)*T8+(-6748-42*T1)*T6+(5610+ #57*T1)*T12+(-36*T1-2598)*T5+(9*T1+624)*R-60)*T97/36 W(8) = -T27*(132*T46-1374*T7+(-12*T1+4422)*T8+(-6748+42*T1)*T6+(56 #10-57*T1)*T12+(36*T1-2598)*T5+(-9*T1+624)*R-60)*T97/36 W(9) = -2.D0/9.D0*(66*T46-687*T7+2175*T8-3248*T6+2634*T12-1191*T5+ #285*R-30)*T62*T22 RETURN C ---------------------------- C CASE R1=R AND R2=(1-R)/2 C (H2=H3 AND H1#H2) C ---------------------------- 30 CONTINUE R=R1 T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R**2 T4 = T3**2 T5 = T3*R T12 = T4-4*T5+6*T3-4*R+1 T15 = T4*R T17 = 1/(33*T15+126*T4+24*T5+14*T3-25*R+20) T20 = 1/(R+1) T21 = T20*T17/R T23 = 16+T1 T35 = T4*T5 T36 = T4*T3 T49 = T20*T17 W(1) = -T21*T12*(33*T4-6*T5+T3*(12*T1-33)+R*(-52+18*T1)+6+6*T1)*T2 #/36 W(2) = T21*T12*(33*T4-6*T5+T3*(-12*T1-33)+R*(-52-18*T1)+6-6*T1)*T2 #3/36 W(3) = T21*(33*T4-6*T5-105*T3-160*R-30)*T12/9 W(4) = T49*(165*T35-855*T36+T15*(168*T1+1098)+T4*(792*T1-1568)+T5* #(264*T1+49)+T3*(-2375+120*T1)+R*(240*T1+3040)-1090-240*T1)*T2/360 W(5) = -T49*(165*T35-855*T36+T15*(-168*T1+1098)+T4*(-1568-792*T1)+ #T5*(-264*T1+49)+T3*(-120*T1-2375)+R*(3040-240*T1)+240*T1-1090)*T23 #/360 W(6) = -T20*T17*(70-1264*T4+18*T15-171*T36+33*T35-307*T5-619*T3+32 #0*R)/18 W(7) = T17*(165*T36-1020*T15+T4*(2022+1752*T1)+T5*(346+336*T1)+T3* #(855+360*T1)+R*(-830-480*T1)+240*T1+190)*T2/360 W(8) = -T17*(165*T36-1020*T15+T4*(2022-1752*T1)+T5*(346-336*T1)+T3 #*(855-360*T1)+R*(-830+480*T1)-240*T1+190)*T23/360 W(9) = -T17*(33*T36-204*T15-1698*T4-334*T5-261*T3+410*R-250)/18 RETURN C ---------------------------- C CASE R2=R AND R1=(1-R)/2 C (H1=H3 AND H2#H1) C ---------------------------- 40 CONTINUE R=R2 T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R**2 T4 = T3*R T5 = T3**2 T6 = T4*T5 T7 = T5*T3 T9 = T5*R T21 = R-1 T23 = 1/(33*T9+93*T5+30*T4+14*T3-15*R+5) T25 = R+1 T26 = T25**2 T27 = 1/T26 T30 = 16+T1 T48 = T27*T23 T62 = T3-2*R+1 T66 = T27*T23/R T85 = T5**2 W(1) = T27*T23*T21*(165*T6+795*T7+T9*(240*T1+945)+T5*(564*T1-641)+ #T4*(168*T1-217)+T3*(-327+48*T1)+R*(-120*T1+355)-115+60*T1)*T2/360 W(2) = -T27*T23*(165*T6+795*T7+T9*(-240*T1+945)+T5*(-641-564*T1)+T #4*(-168*T1-217)+T3*(-48*T1-327)+R*(355+120*T1)-60*T1-115)*T21*T30/ #360 W(3) = -T48*(33*T6+159*T7-99*T9-805*T5-245*T4-123*T3+215*R-95)*T21 #/18 W(4) = -T66*T62*(33*T6+225*T7+T9*(3*T1+405)+T5*(3*T1+137)+T4*(-6*T #1+35)+T3*(-61-6*T1)+R*(3*T1+23)+3+3*T1)*T2/36 W(5) = T66*T62*(33*T6+225*T7+T9*(-3*T1+405)+T5*(137-3*T1)+T4*(6*T1 #+35)+T3*(-61+6*T1)+R*(23-3*T1)-3*T1+3)*T30/36 W(6) = T66*(33*T6+225*T7+387*T9+119*T5+71*T4-25*T3+5*R-15)*T62/9 W(7) = T48*(165*T85+630*T6+T7*(240*T1-1410)+T9*(1116*T1-3914)+T5*( #1096+1836*T1)+T4*(600*T1+130)+T3*(1138+168*T1)+R*(-430-180*T1)+60* #T1+35)*T2/360 W(8) = -T48*(165*T85+630*T6+T7*(-240*T1-1410)+T9*(-1116*T1-3914)+T #5*(-1836*T1+1096)+T4*(-600*T1+130)+T3*(1138-168*T1)+R*(-430+180*T1 #)-60*T1+35)*T30/360 W(9) = -T23/T25*(33*T6+93*T7-663*T9-1459*T5-525*T4-169*T3+195*R-65 #)/18 RETURN C ---------------------------- C GENERAL CASE C ---------------------------- 50 CONTINUE T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R1**2 T4 = T3**2 T5 = T4**2 T7 = T3*R1 T8 = T4*T7 T9 = T8*(-70+30*R2) T10 = R2**2 T12 = T4*T3 T14 = T10*R2 T15 = R2*T1 T17 = T4*R1 T19 = T10*T1 T20 = T10**2 T23 = T14*T1 T24 = T20*R2 T27 = T20*T10 T28 = T20*T1 T29 = -45+603*T19-654*T14+10*T27+30*T1+108*T10-270*T15+549*T20+180 #*R2+144*T28-504*T23-150*T24 T31 = T24*T1 T32 = -186*T28+84*T24+54*T31-18*T10-276*T20+45*T15-40*T27+303*T23- #30*R2-198*T19+198*T14 T34 = T27*T1 T35 = 10*T5+T9+T12*(75*T10+135-270*R2+30*T1)+T17*(-40-120*T1+540*R #2+100*T14-600*T10+135*T15)+T4*(252*T19+180*T1-270*R2+1092*T10+75*T #20-140-700*T14-450*T15)+T7*(-672*T19+540*T15-450*T20-120*T1-180*R2 #+150-672*T10+249*T23+30*T24+1104*T14)+T3*T29+R1*T32+12*T27-48*T14- #24*T31+60*T20-48*T23+15*T19+60*T28+15*T10+12*T34-24*T24 T37 = R1-1 T42 = (-1+R2)**2 T43 = T3+(-2+2*R2)*R1+T42 T56 = T20*T14 T59 = 2*T5+T8*(8*R2-8)+T12*(15+21*T10-46*R2)+T17*(99*R2-17-117*T10 #+35*T14)+T4*(35*T20-170*T14-101*R2+225*T10+11)+T7*(-3+49*R2+21*T24 #-192*T10-145*T20+270*T14)+T3*(72*T10-9*R2-177*T14+8*T27-72*T24+180 #*T20)+R1*(-9*T10+63*T24+45*T14-81*T20+2*T56-20*T27)-3*T14-18*T24+6 #*T27-4*T56+9*T20 T60 = 1/T59 T63 = 1/(R1+R2) T64 = 1/R1 T66 = T64*T63*T60*T43 T68 = 16+T1 T77 = -45-603*T19-654*T14+10*T27-30*T1+108*T10+270*T15+549*T20+180 #*R2-144*T28+504*T23-150*T24 T79 = 186*T28+84*T24-54*T31-18*T10-276*T20-45*T15-40*T27-303*T23-3 #0*R2+198*T19+198*T14 T81 = 10*T5+T9+T12*(75*T10+135-270*R2-30*T1)+T17*(-40+120*T1+540*R #2+100*T14-600*T10-135*T15)+T4*(-252*T19-180*T1-270*R2+1092*T10+75* #T20-140-700*T14+450*T15)+T7*(672*T19-540*T15-450*T20+120*T1-180*R2 #+150-672*T10-249*T23+30*T24+1104*T14)+T3*T77+R1*T79+12*T27-48*T14+ #24*T31+60*T20+48*T23-15*T19-60*T28+15*T10-12*T34-24*T24 T100 = 2*T5+T8*(6*R2-14)+T12*(-54*R2+15*T10-9)+T17*(20*T14-120*T10 #+136-54*R2)+T4*(15*T20+486*R2-140*T14-84*T10-244)+T7*(6*T24-90*T20 #+174-78*T14-684*R2+672*T10)+T3*(-45+474*T14+2*T27+360*R2-63*T20-70 #2*T10-30*T24)+R1*(168*T20-48*T24+234*T10-324*T14-60*R2-8*T27)-12*T #27+48*T14+24*T24-15*T10-60*T20 T102 = T63*T60 T113 = -168*T23+138+2478*T10-1698*R2+75*T27+663*T19+138*T1+9*T28-1 #668*T14+1869*T20-900*T24-648*T15 T115 = -72*T1-930*T19-1808*T20-72-600*T27-48*T28+429*T23+1950*T24- #3240*T10+561*T15+30*T56+1496*R2+2064*T14 T117 = T20**2 T118 = -210*T56+15*T1-672*R2+2373*T10+10*T117-510*T23-2040*T14+678 #*T19+99*T28+879*T20+15+1095*T27-1470*T24-252*T15 T120 = -90*T28-60*T117+120*R2+45*T15-240*T19+285*T23-840*T10-390*T #20-690*T27+1140*T14+330*T24+240*T56 T122 = T5*(6+10*T10+6*T1)+T8*(78*R2-30-90*T10+18*T15-30*T1+30*T14) #+T12*(432*T10-138*T15-330*T14+81+81*T1+75*T20+27*T19-398*R2)+T17*( #-750*T20+1074*R2+414*T15-1228*T10-138*T1-228*T19+24*T23+100*T24+10 #44*T14-138)+T4*T113+T7*T115+T3*T118+R1*T120+135*T20+30*T24+60*T117 #+30*T28-240*T14+30*T19+105*T10+180*T27-20*T56-60*T23 T125 = 1/T37 T126 = T60*T125 T127 = 1/R2 T129 = T127*T63*T126 T139 = -168*T23-138-2478*T10+1698*R2-75*T27+663*T19+138*T1+9*T28+1 #668*T14-1869*T20+900*T24-648*T15 T141 = -72*T1-930*T19+1808*T20+72+600*T27-48*T28+429*T23-1950*T24+ #3240*T10+561*T15-30*T56-1496*R2-2064*T14 T143 = 210*T56+15*T1+672*R2-2373*T10-10*T117-510*T23+2040*T14+678* #T19+99*T28-879*T20-15-1095*T27+1470*T24-252*T15 T145 = -90*T28+60*T117-120*R2+45*T15-240*T19+285*T23+840*T10+390*T #20+690*T27-1140*T14-330*T24-240*T56 T147 = T5*(-6-10*T10+6*T1)+T8*(-78*R2+30+90*T10+18*T15-30*T1-30*T1 #4)+T12*(-432*T10-138*T15+330*T14-81+81*T1-75*T20+27*T19+398*R2)+T1 #7*(750*T20-1074*R2+414*T15+1228*T10-138*T1-228*T19+24*T23-100*T24- #1044*T14+138)+T4*T139+T7*T141+T3*T143+R1*T145-135*T20-30*T24-60*T1 #17+30*T28+240*T14+30*T19-105*T10-180*T27+20*T56-60*T23 T167 = T5*(-6+2*T10)+T8*(30+6*T14-18*T10-6*R2)+T12*(-66*T14+86*R2+ #54*T10-81+15*T20)+T17*(-282*R2-150*T20+28*T10+138+180*T14+20*T24)+ #T4*(-180*T24+15*T27+438*R2-300*T10-132*T14-138+363*T20)+T7*(-120*T #27+6*T56-102*T14+468*T10+390*T24-374*R2-304*T20+72)+T3*(2*T117+219 #*T27+204*T14-15-339*T10-42*T56+57*T20+168*R2-294*T24)+R1*(66*T24+3 #0*T20-12*T117-30*R2-138*T27+120*T10-114*T14+48*T56)+24*T14-9*T20-4 #*T56+12*T117+6*T24-15*T10+36*T27 T172 = T5*T3 T174 = T5*R1 T175 = T174*(-120+50*R2) T184 = 145*T27+14898*T10-2625*T24-4860*R2+30*T1-19660*T14+11529*T2 #0+505-90*T15-30*T23+9*T28+63*T19 T186 = -75+16878*T14+48*T23+1640*R2+50*T56-16073*T20+15*T15+6900*T #24-1200*T27-8238*T10-48*T28-18*T19 T192 = 10*T172+T175+T5*(-630*R2+145*T10+585+30*T1)+T8*(-1815*T10-1 #20*T1+2880*R2+45*T15-1415+280*T14)+T12*(-3255*T14-150*T15+7467*T10 #+1860-6325*R2+27*T19+180*T1+350*T20)+T17*(-14482*T10-72*T19+7470*R #2+11784*T14+9*T23-1350+180*T15-3675*T20-120*T1+280*T24)+T4*T184+T7 #*T186+T3*(2475*T27-225*R2-42*T23+10*T117+2250*T10-360*T56-7392*T14 #-7515*T24+10974*T20+99*T28)+R1*(3600*T24-60*T117-3510*T20-90*T28+1 #5*T23-225*T10+1440*T14+540*T56-1810*T27)-200*T56-75*T14-630*T24+60 #*T117+30*T28+405*T20+420*T27 T203 = 145*T27+14898*T10-2625*T24-4860*R2-30*T1-19660*T14+11529*T2 #0+505+90*T15+30*T23-9*T28-63*T19 T205 = -75+16878*T14-48*T23+1640*R2+50*T56-16073*T20-15*T15+6900*T #24-1200*T27-8238*T10+48*T28+18*T19 T211 = 10*T172+T175+T5*(-630*R2+145*T10+585-30*T1)+T8*(-1815*T10+1 #20*T1+2880*R2-45*T15-1415+280*T14)+T12*(-3255*T14+150*T15+7467*T10 #+1860-6325*R2-27*T19-180*T1+350*T20)+T17*(-14482*T10+72*T19+7470*R #2+11784*T14-9*T23-1350-180*T15-3675*T20+120*T1+280*T24)+T4*T203+T7 #*T205+T3*(2475*T27-225*R2+42*T23+10*T117+2250*T10-360*T56-7392*T14 #-7515*T24+10974*T20-99*T28)+R1*(3600*T24-60*T117-3510*T20+90*T28-1 #5*T23-225*T10+1440*T14+540*T56-1810*T27)-200*T56-75*T14-630*T24+60 #*T117-30*T28+405*T20+420*T27 T232 = 2*T172+T174*(10*R2-24)+T5*(81-126*R2+29*T10)+T8*(56*T14-363 #*T10+522*R2-139)+T12*(70*T20+156-651*T14-1085*R2+1461*T10)+T17*(-7 #35*T20+1278*R2+56*T24+2346*T14-2810*T10-126)+T4*(2904*T10+65-864*R #2+29*T27-525*T24-3896*T14+2295*T20)+T7*(-1626*T10+1380*T24-240*T27 #-3157*T20+310*R2+10*T56-15+3318*T14)+T3*(-45*R2-1428*T14+2*T117-72 #*T56+495*T27+450*T10-1503*T24+2076*T20)+R1*(720*T24-12*T117-594*T2 #0+108*T56+270*T14-362*T27-45*T10)-15*T14+45*T20-40*T56+12*T117-126 #*T24+84*T27 W(1) = -T66*T37*T35*T2/180 W(2) = T66*T37*T81*T68/180 W(3) = T64*T102*T100*T43*T37/9 W(4) = -T129*T43*T122*T2/180 W(5) = -T129*T43*T147*T68/180 W(6) = T127*T102*T125*T167*T43/9 W(7) = T126*T192*T2/180 W(8) = -T126*T211*T68/180 W(9) = -T125*T60*T232/9 RETURN END SUBROUTINE VECTWY(WY,R,T) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION WY(*) IF (R.EQ.0.5D0) THEN T1 = DSQRT(6.D0) T3 = -16+T1 T5 = T**2 T6 = T5**2 T10 = T5*T T22 = T1+16 WY(1) = 8.D0/45.D0*T3*(-15.D0/2.D0*T1+45.D0/4.D0)*T6+8.D0/45.D0*(1 #8*T1-103.D0/4.D0)*T3*T10+8.D0/45.D0*(153.D0/8.D0-117.D0/8.D0*T1)*T #3*T5+8.D0/45.D0*(-151.D0/32.D0+9.D0/2.D0*T1)*T3*T+8.D0/45.D0*T3*(- #3.D0/8.D0*T1+3.D0/32.D0) WY(2) = -8.D0/45.D0*(15.D0/2.D0*T1+45.D0/4.D0)*T22*T6-8.D0/45.D0*( #-18*T1-103.D0/4.D0)*T22*T10-8.D0/45.D0*(153.D0/8.D0+117.D0/8.D0*T1 #)*T22*T5-8.D0/45.D0*(-151.D0/32.D0-9.D0/2.D0*T1)*T22*T-8.D0/45.D0* #(3.D0/8.D0*T1+3.D0/32.D0)*T22 WY(3) = -805.D0/9.D0*T+19.D0/3.D0-4040.D0/9.D0*T10+200*T6+332*T5 WY(4) = 16.D0/45.D0*(15.D0/4.D0*T1+105.D0/8.D0)*T3*T6+16.D0/45.D0* #(-9*T1-227.D0/8.D0)*T3*T10+16.D0/45.D0*(327.D0/16.D0+117.D0/16.D0* #T1)*T3*T5+16.D0/45.D0*(-359.D0/64.D0-9.D0/4.D0*T1)*T3*T+16.D0/45.D #0*(3.D0/16.D0*T1+27.D0/64.D0)*T3 WY(5) = -16.D0/45.D0*(-15.D0/4.D0*T1+105.D0/8.D0)*T22*T6-16.D0/45. #D0*(9*T1-227.D0/8.D0)*T22*T10-16.D0/45.D0*(327.D0/16.D0-117.D0/16. #D0*T1)*T22*T5-16.D0/45.D0*(-359.D0/64.D0+9.D0/4.D0*T1)*T22*T-16.D0 #/45.D0*(-3.D0/16.D0*T1+27.D0/64.D0)*T22 WY(6) = -11.D0/9.D0*T-184.D0/9.D0*T10+40.D0/3.D0*T6+28.D0/3.D0*T5 ELSE T1 = -2*T+2 T2 = R**2 T3 = T2**2 T4 = T3*T2 T5 = T1*T4 T6 = T**2 T8 = T3*R T9 = (12*T6-12*T)*T8 T10 = T6*T T12 = (21*T-1-20*T10)*T3 T13 = DSQRT(6.D0) T14 = T13*T T15 = T13*T6 T16 = T6**2 T18 = T2*R T20 = T10*T13 T23 = T16*T13 T26 = T5+T9+T12+(-42*T6-15*T+40*T10-12*T14+3*T13+9*T15+10*T16+7)*T #18+(36*T6-30*T16-3-9*T+27*T15+6*T10-24*T20-3*T13)*T2+(15*T16+15*T2 #3+12*T14+18*T6-30*T10-27*T15-3*T)*R+10*T16+6*T6-9*T15+24*T20-16*T1 #0-15*T23 T27 = -16+T13 T32 = 1/T2/(2*T18-4*T2+9*R-5) T40 = T5+T9+T12+(-42*T6-15*T+40*T10+12*T14-3*T13-9*T15+10*T16+7)*T #18+(36*T6-30*T16-3-9*T-27*T15+6*T10+24*T20+3*T13)*T2+(15*T16-15*T2 #3-12*T14+18*T6-30*T10+27*T15-3*T)*R+10*T16+6*T6+9*T15-24*T20-16*T1 #0+15*T23 T41 = T13+16 T44 = -4*T+4 T45 = T3**2 T51 = T3*T18 T65 = T44*T45*R+(-12-12*T+24*T6)*T45+(10-40*T10+102*T-72*T6)*T51+( #-224*T+16+200*T10-12*T6+20*T16)*T4+(18-120*T16-102*T+624*T6-420*T1 #0)*T8+(26+360*T16+54*T6-320*T10-120*T)*T3+(-220*T10-96*T-180*T16-3 #2+528*T6)*T18+(-486*T6+15+495*T16-240*T10+216*T)*T2+(-600*T16+720* #T10-120*T)*R+150*T6-400*T10+250*T16 T69 = T1*T51 T71 = (-4*T-8+12*T6)*T4 T73 = (-20*T10-48*T6+57*T+11)*T8 T81 = (50*T10+90*T6-140*T16)*R T85 = 1/(2*T4-10*T8+27*T3-46*T18+46*T2-24*R+5) WY(1) = -T26*T27*T32/18 WY(2) = T40*T41*T32/18 WY(3) = T65/(2*T45-10*T51+27*T4-46*T8+46*T3-24*T18+5*T2)/9 WY(4) = (T69+T71+T73+(-123*T+3*T13+9*T15-7-12*T14+120*T10+10*T16)* #T3+(-214*T10+93*T+186*T6-70*T16-24*T20+36*T15-12*T14+5)*T18+(165*T #16+66*T10-186*T6+15*T23-24*T20-45*T+9*T15)*T2+T81-50*T10+50*T16)*T #27*T85/18 WY(5) = -(T69+T71+T73+(-123*T-3*T13-9*T15-7+12*T14+120*T10+10*T16) #*T3+(-214*T10+93*T+186*T6-70*T16+24*T20-36*T15+12*T14+5)*T18+(165* #T16+66*T10-186*T6-15*T23+24*T20-45*T-9*T15)*T2+T81-50*T10+50*T16)* #T41*T85/18 WY(6) = -T85*(T44*T51+(-44*T+24*T6+2)*T4+(4+6*T+120*T6-40*T10)*T8+ #(-48*T+20*T16+13-108*T6-120*T10)*T3+(264*T6-8+220*T10-102*T+40*T16 #)*T18+(-120*T16-480*T10+72*T+114*T6)*T2+(100*T10-144*T6+260*T16)*R #+80*T10-125*T16)/9 ENDIF RETURN END SUBROUTINE VECTWZ(INDCAS,W,R1,R2,T) IMPLICIT DOUBLE PRECISION (A-H,O-Z) DIMENSION W(*) GOTO (10,20,30,40,50),INDCAS C ------------------------------------------- C CASE R1=R2=1/3 (H1=H2=H3) C ---------------------------- 10 CONTINUE T1 = DSQRT(6.D0) T2 = T1-16 T5 = T**2 T6 = T5**2 T10 = T5*T T21 = 16+T1 W(1) = -T2*(3645*T1-2430)*T6/4320-T2*(4212-7938*T1)*T10/4320-T2*(- #1728+5832*T1)*T5/4320-T2*(-306-1611*T1)*T/4320-T2*(117*T1+202)/432 #0 W(2) = -T21*(3645*T1+2430)*T6/4320-T21*(-4212-7938*T1)*T10/4320-T2 #1*(1728+5832*T1)*T5/4320-T21*(306-1611*T1)*T/4320-T21*(117*T1-202) #/4320 W(3) = 26.D0/3.D0*T+48*T10-45.D0/2.D0*T6-34*T5-25.D0/54.D0 W(4) = (-45.D0/16.D0*T1+45)*T6-T2*(648*T1-11502)*T10/2160-T2*(6588 #-972*T1)*T5/2160-T2*(-1089+396*T1)*T/2160-T2*(-36*T1+29)/2160 W(5) = (45+45.D0/16.D0*T1)*T6+T21*(-648*T1-11502)*T10/2160+T21*(65 #88+972*T1)*T5/2160+T21*(-1089-396*T1)*T/2160+T21*(36*T1+29)/2160 W(6) = -77.D0/12.D0*T+49.D0/108.D0-57.D0/2.D0*T10+45.D0/4.D0*T6+23 #*T5 W(7) = T2*(3645*T1+9720)*T6/4320+T2*(-18792-6642*T1)*T10/4320+T2*( #11448+3888*T1)*T5/4320+T2*(-2484-819*T1)*T/4320+T2*(45*T1+140)/432 #0 W(8) = T21*(3645*T1-9720)*T6/4320+T21*(18792-6642*T1)*T10/4320+T21 #*(-11448+3888*T1)*T5/4320+T21*(2484-819*T1)*T/4320+T21*(45*T1-140) #/4320 W(9) = -9.D0/4.D0*T-39.D0/2.D0*T10+45.D0/4.D0*T6+11*T5+13.D0/108.D #0 RETURN C ------------------------------- C CASE R1=R2=R (H1=H2 AND H3#H1) C ------------------------------- 20 CONTINUE R=R1 T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R**2 T4 = T3**2 T5 = T4**2 T6 = T5*R T8 = (-5940*T+60)*T5 T9 = T**2 T10 = T1*T T12 = T3*R T13 = T4*T12 T15 = T1*T9 T16 = T9*T T18 = T4*T3 T20 = T9**2 T21 = T1*T16 T23 = T4*R T25 = T1*T20 T34 = 660*T6+T8+(12870*T9+8214*T-1752*T1-4107+3504*T10)*T13+(3672* #T1-11232*T15+144*T10-9900*T16+6577-19242*T9-1236*T)*T18+(-3312*T+1 #4400*T16-9792*T10-2970*T1+8568*T9-4365+15948*T15+2475*T20+10200*T2 #1)*T23+(-984*T9-3024*T15-2790*T25+1327-3465*T20+10698*T10+1032*T1- #4560*T16+1428*T-19260*T21)*T4+(-138*T1+2106*T9-168-6*T+12780*T21-2 #520*T16-4176*T10-6444*T15+495*T20+5895*T25)*T12+(820*T16-2280*T21- #4815*T25+1635*T20-1542*T9-36*T+594*T10+3708*T15)*T3+(1695*T25-930* #T20+288*T9+360*T16-540*T21-612*T15)*R-225*T25+150*T20+180*T21-120* #T16 T37 = 1/(66*T23-291*T4+429*T12-293*T3+93*R-12) T38 = 1/T12 T39 = T37*T38 T41 = T1+16 T56 = 660*T6+T8+(12870*T9+8214*T+1752*T1-4107-3504*T10)*T13+(-3672 #*T1+11232*T15-144*T10-9900*T16+6577-19242*T9-1236*T)*T18+(-3312*T+ #14400*T16+9792*T10+2970*T1+8568*T9-4365-15948*T15+2475*T20-10200*T #21)*T23+(-984*T9+3024*T15+2790*T25+1327-3465*T20-10698*T10-1032*T1 #-4560*T16+1428*T+19260*T21)*T4+(138*T1+2106*T9-168-6*T-12780*T21-2 #520*T16+4176*T10+6444*T15+495*T20-5895*T25)*T12+(820*T16+2280*T21+ #4815*T25+1635*T20-1542*T9-36*T-594*T10-3708*T15)*T3+(-1695*T25-930 #*T20+288*T9+360*T16+540*T21+612*T15)*R+225*T25+150*T20-180*T21-120 #*T16 T75 = 132*T6+(-1188*T+12)*T5+(1281+2574*T9-2562*T)*T13+(-420*T-309 #1-1980*T16+9630*T9)*T18+(2691+11088*T+495*T20-9360*T16-17424*T9)*T #23+(-12552*T+3432*T9+22200*T16-973+2655*T20)*T4+(5010*T+132+8154*T #9-15840*T16-6975*T20)*T12+(-720*T-4758*T9+6105*T20+2900*T16)*T3+(7 #20*T16+792*T9-2220*T20)*R-240*T16+300*T20 T78 = T5*T3 T80 = (-5940*T-600)*T6 T97 = 660*T78+T80+(-903+7626*T-168*T1+336*T10+12870*T9)*T5+(-1770* #T+816*T1-480*T10+3196-21168*T9-1728*T15-9900*T16)*T13+(-1278*T1-30 #24*T10+6660*T15+26370*T9+2475*T20-12444*T+2280*T21-4198+19980*T16) #*T18+(3088-10380*T21-14088*T9+7110*T10-810*T25+22860*T-38880*T16+9 #18*T1-5760*T20-5148*T15)*T23+(-5886*T10+45960*T16-3600*T15+14760*T #21-10110*T9+13590*T20-1055-330*T1-19326*T+3915*T25)*T4+(-18840*T20 #+2250*T10-8100*T21+7050*T-6390*T25+42*T1+132+20208*T9+6408*T15-265 #00*T16)*T12+(900*T21+13815*T20-306*T10-9210*T9+3100*T16+4740*T25-3 #060*T15-936*T)*T3+(720*T21+1368*T9+468*T15-4680*T20-1680*T25+1920* #T16)*R-180*T21+225*T25+600*T20-480*T16 T100 = 1/(R-1) T102 = T100*T37*T38 T120 = 660*T78+T80+(-903+7626*T+168*T1-336*T10+12870*T9)*T5+(-1770 #*T-816*T1+480*T10+3196-21168*T9+1728*T15-9900*T16)*T13+(1278*T1+30 #24*T10-6660*T15+26370*T9+2475*T20-12444*T-2280*T21-4198+19980*T16) #*T18+(3088+10380*T21-14088*T9-7110*T10+810*T25+22860*T-38880*T16-9 #18*T1-5760*T20+5148*T15)*T23+(5886*T10+45960*T16+3600*T15-14760*T2 #1-10110*T9+13590*T20-1055+330*T1-19326*T-3915*T25)*T4+(-18840*T20- #2250*T10+8100*T21+7050*T+6390*T25-42*T1+132+20208*T9-6408*T15-2650 #0*T16)*T12+(-900*T21+13815*T20+306*T10-9210*T9+3100*T16-4740*T25+3 #060*T15-936*T)*T3+(-720*T21+1368*T9-468*T15-4680*T20+1680*T25+1920 #*T16)*R+180*T21-225*T25+600*T20-480*T16 T141 = 132*T78+(-1188*T-120)*T6+(1122*T+2574*T9+21)*T5+(-2160*T9-1 #980*T16-340+222*T)*T13+(694+495*T20+1260*T16-2718*T9+1140*T)*T18+( #4680*T16-3960*T+3360*T9-484-180*T20)*T23+(2298*T9+3198*T-1980*T20+ #185-8520*T16)*T4+(-24+4420*T16-3648*T9-1290*T+3900*T20)*T12+(180*T #+1830*T9-460*T16-2925*T20)*T3+(-480*T16+1080*T20-288*T9)*R-150*T20 #+120*T16 T145 = (-504-2376*T)*T13 T161 = 1/(-1+2*R) T163 = T100*T161*T37 W(1) = -T2*T34*T39/360 W(2) = T41*T56*T39/360 W(3) = T75*T37*T38/18 W(4) = -T2*T97*T102/360 W(5) = T41*T120*T102/360 W(6) = T141*T100*T39/18 W(7) = T2*(264*T5+T145+(5148*T9-24*T1+246+4824*T+48*T10)*T18+(88+4 #2*T1-162*T15-3246*T-3960*T16-10800*T9+24*T10)*T23+(8280*T9+264*T-1 #62*T10+180*T15-62+156*T21-18*T1+8496*T16+990*T20)*T4+(-2160*T20-26 #4*T21-45*T25-1848*T9+90*T15+90*T10+246*T-7044*T16)*T12+(72*T21-156 #*T9+90*T25+2112*T16+1890*T20-108*T15)*T3+(-84*T16-45*T25+36*T21-66 #0*T20)*R+60*T20)*T163/36 W(8) = -T41*(264*T5+T145+(5148*T9+24*T1+246+4824*T-48*T10)*T18+(88 #-42*T1+162*T15-3246*T-3960*T16-10800*T9-24*T10)*T23+(8280*T9+264*T #+162*T10-180*T15-62-156*T21+18*T1+8496*T16+990*T20)*T4+(-2160*T20+ #264*T21+45*T25-1848*T9-90*T15-90*T10+246*T-7044*T16)*T12+(-72*T21- #156*T9-90*T25+2112*T16+1890*T20+108*T15)*T3+(-84*T16+45*T25-36*T21 #-660*T20)*R+60*T20)*T163/36 W(9) = -2.D0/9.D0*(132*T5+(-252-1188*T)*T13+(195+2268*T+2574*T9)*T #18+(-82-1695*T-4914*T9-1980*T16)*T23+(3600*T9+618*T+23+495*T20+378 #0*T16)*T4+(-2730*T16-945*T20-1194*T9-147*T)*T12+(840*T16+246*T9+67 #5*T20)*T3+(-150*T16-195*T20)*R+30*T20)*T100*T161*T37 RETURN C ---------------------------- C CASE R1=R AND R2=(1-R)/2 C (H2=H3 AND H1#H2) C ---------------------------- 30 CONTINUE R=R1 T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R**2 T4 = T3**2 T5 = T4**2 T7 = T3*R T8 = T4*T7 T9 = T8*(456-594*T) T10 = T**2 T11 = T*T1 T13 = T4*T3 T15 = T10*T T16 = T10*T1 T18 = T4*R T20 = T10**2 T21 = T15*T1 T24 = T1*T20 T31 = 33*T5+T9+T13*(2574*T10+24*T11-3354*T-12*T1+969)+T18*(-24*T1+ #488-3528*T-3960*T15+156*T11+6912*T10-162*T16)+T4*(-1140*T+1980*T20 #-486*T16+312*T21+3654*T10-12*T1+192*T11-4488*T15+67)+T7*(360*T20+7 #68*T10-144+12*T1-324*T16+438*T+528*T21-180*T24-1632*T15)+T3*(960*T #15-192*T11-180*T24+324*T16+678*T+19+24*T1-2166*T10+720*T20)+R*(32- #1320*T20+486*T16-384*T10+1752*T15-156*T11-156*T+180*T24-528*T21+12 #*T1)-24*T-24*T11-312*T15+180*T20+162*T16-312*T21+180*T24+162*T10 T33 = 1/R T35 = 1/(R+1) T38 = 1/(33*T18+126*T4+24*T7+14*T3-25*R+20) T39 = T38*T35*T33 T41 = T1+16 T54 = 33*T5+T9+T13*(2574*T10-24*T11-3354*T+12*T1+969)+T18*(24*T1+4 #88-3528*T-3960*T15-156*T11+6912*T10+162*T16)+T4*(-1140*T+1980*T20+ #486*T16-312*T21+3654*T10+12*T1-192*T11-4488*T15+67)+T7*(360*T20+76 #8*T10-144-12*T1+324*T16+438*T-528*T21+180*T24-1632*T15)+T3*(960*T1 #5+192*T11+180*T24-324*T16+678*T+19-24*T1-2166*T10+720*T20)+R*(32-1 #320*T20-486*T16-384*T10+1752*T15+156*T11-156*T-180*T24+528*T21-12* #T1)-24*T+24*T11-312*T15+180*T20-162*T16+312*T21-180*T24+162*T10 T69 = 33*T5+T9+T13*(-3498*T+1041+2574*T10)+T18*(-4464*T+632-3960*T #15+7884*T10)+T4*(-2292*T+6570*T10+1980*T20-6360*T15+139)+T7*(438*T #+1440*T20+2712*T10-4800*T15-216)+T3*(960*T15-125-4110*T10+1830*T+1 #800*T20)+R*(780*T-40+4920*T15-3300*T10-2400*T20)+120*T-900*T20-810 #*T10+1560*T15 T71 = T38*T35 T73 = T5*T3 T75 = T5*R T76 = T75*(1620-2970*T) T93 = 165*T73+T76+T5*(-2748-168*T1-5964*T+336*T11+12870*T10)+T8*(- #2420+22764*T-3456*T16-1440*T1-8136*T10-19800*T15+5184*T11)+T13*(-3 #768*T1+14782+22848*T11-6768*T10+9900*T20-29808*T16-46332*T+37320*T #15+9120*T21)+T18*(-4224*T1-85560*T15+154680*T10+54240*T21-50760*T- #26280*T20-404-77760*T16-6480*T24+35760*T11)+T4*(75060*T20-70680*T1 #5+36-36*T-30240*T24+7500*T10+78720*T21-1944*T1+23184*T11-69840*T16 #)+T7*(-96*T1-5196+5856*T11+37440*T21-41520*T20-25344*T16+11316*T-1 #8720*T24+38376*T10-20360*T15)+T3*(120*T1+21612*T-288*T11+7200*T21- #78432*T10-4752*T16+3800*T15+65940*T20-960*T24+1845)+R*(7080*T10-36 #00*T24+64280*T15+480*T21+2880*T16-66600*T20-720*T11-11070*T)+720*T #16+11070*T10-2880*T21+2400*T24+21900*T20-31880*T15 T95 = R-1 T96 = T95**2 T98 = 1/T96/T95 T99 = T38*T98 T100 = T35*T99 T118 = 165*T73+T76+T5*(-2748+168*T1-5964*T-336*T11+12870*T10)+T8*( #-2420+22764*T+3456*T16+1440*T1-8136*T10-19800*T15-5184*T11)+T13*(3 #768*T1+14782-22848*T11-6768*T10+9900*T20+29808*T16-46332*T+37320*T #15-9120*T21)+T18*(4224*T1-85560*T15+154680*T10-54240*T21-50760*T-2 #6280*T20-404+77760*T16+6480*T24-35760*T11)+T4*(75060*T20-70680*T15 #+36-36*T+30240*T24+7500*T10-78720*T21+1944*T1-23184*T11+69840*T16) #+T7*(96*T1-5196-5856*T11-37440*T21-41520*T20+25344*T16+11316*T+187 #20*T24+38376*T10-20360*T15)+T3*(-120*T1+21612*T+288*T11-7200*T21-7 #8432*T10+4752*T16+3800*T15+65940*T20+960*T24+1845)+R*(7080*T10+360 #0*T24+64280*T15-480*T21-2880*T16-66600*T20+720*T11-11070*T)-720*T1 #6+11070*T10+2880*T21-2400*T24+21900*T20-31880*T15 T139 = 33*T73+T75*(324-594*T)+T5*(-1596*T+2574*T10-348)+T8*(1244-1 #668*T+2520*T10-3960*T15)+T13*(7478-3480*T15+34416*T10-36684*T+1980 #*T20)+T18*(-82200*T15+4988+2520*T20+124248*T10-53064*T)+T4*(51300* #T20-27828*T-108600*T15+2340+85308*T10)+T7*(38088*T10-924-49000*T15 #+14160*T20-4764*T)+T3*(14340*T20+225-9984*T10-7880*T15+4668*T)+R*( #-9000*T20-2040*T10-1350*T+12280*T15)-2920*T15+1500*T20+1350*T10 T143 = T5*(1455-2970*T) T158 = 165*T75+T143+T8*(-3186*T+3504*T11-4107-1752*T1+12870*T10)+T #13*(-22464*T16+55902*T-7913+20112*T11-2568*T1-19800*T15-37134*T10) #+T18*(-4569+51750*T+40800*T21-1344*T1-153378*T10+16560*T11+9900*T2 #0-51408*T16+103200*T15)+T4*(-192*T1+39360*T21+26130*T+120360*T15-2 #2320*T24-26352*T16-83622*T10-66420*T20+7200*T11-827)+T7*(-4854*T-1 #584*T11+216*T1-5040*T24+831-31758*T10-7488*T16-14040*T20+16320*T21 #+54720*T15)+T3*(-3840*T21-18840*T20-5040*T24+24822*T10-395-7240*T1 #5-2262*T-432*T11-120*T1+6912*T16)+R*(-5910*T10+2370*T-3360*T21+600 #0*T24-2160*T16-14560*T15+19500*T20+720*T11)-720*T16+2880*T21+9080* #T15-6900*T20-2400*T24-2370*T10 T175 = 165*T75+T143+T8*(-3186*T-3504*T11-4107+1752*T1+12870*T10)+T #13*(22464*T16+55902*T-7913-20112*T11+2568*T1-19800*T15-37134*T10)+ #T18*(-4569+51750*T-40800*T21+1344*T1-153378*T10-16560*T11+9900*T20 #+51408*T16+103200*T15)+T4*(192*T1-39360*T21+26130*T+120360*T15+223 #20*T24+26352*T16-83622*T10-66420*T20-7200*T11-827)+T7*(-4854*T+158 #4*T11-216*T1+5040*T24+831-31758*T10+7488*T16-14040*T20-16320*T21+5 #4720*T15)+T3*(3840*T21-18840*T20+5040*T24+24822*T10-395-7240*T15-2 #262*T+432*T11+120*T1-6912*T16)+R*(-5910*T10+2370*T+3360*T21-6000*T #24+2160*T16-14560*T15+19500*T20-720*T11)+720*T16-2880*T21+9080*T15 #-6900*T20+2400*T24-2370*T10 T194 = 33*T75+T5*(291-594*T)+T8*(2574*T10-4842*T+1281)+T13*(-12954 #*T+1499+19530*T10-3960*T15)+T18*(699-9522*T+31014*T10-28320*T15+19 #80*T20)+T4*(-23160*T15-3414*T+14898*T10+65+13500*T20)+T7*(-93+930* #T+3240*T20+2634*T10-8640*T15)+T3*(65-3330*T10+3160*T15+2280*T20+66 #*T)+R*(-390*T+1410*T10+1120*T15-3300*T20)+1500*T20+390*T10-1640*T1 #5 W(1) = -T39*T31*T2/36 W(2) = T39*T54*T41/36 W(3) = T71*T33*T69/9 W(4) = T100*T93*T2/360 W(5) = -T100*T118*T41/360 W(6) = -T71*T98*T139/18 W(7) = T99*T158*T2/360 W(8) = -T99*T175*T41/360 W(9) = -T98*T38*T194/18 RETURN C ---------------------------- C CASE R2=R AND R1=(1-R)/2 C (H1=H3 AND H2#H1) C ---------------------------- 40 CONTINUE R=R2 T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R**2 T4 = T3**2 T5 = T4**2 T6 = T5*T3 T7 = T5*R T8 = T*T1 T9 = T**2 T12 = T1*T9 T14 = T3*R T15 = T4*T14 T17 = T9*T T18 = T17*T1 T19 = T9**2 T21 = T4*T3 T23 = T19*T1 T25 = T4*R T35 = 165*T6+630*T7+T5*(-2055-240*T1+9960*T+480*T8-7920*T9)+T15*(- #2376*T12-20736*T9+25296*T+3816*T8-6776-1116*T1)+T21*(69264*T9+5736 #*T8-936*T12-3840*T18-1596*T1-686+39600*T19-90240*T17-17664*T)+T25* #(15840*T23+5076+516*T1-192000*T17+36504*T12-10104*T8-43200*T18-481 #44*T+148704*T9+87840*T19)+T4*(57600*T17-30960*T19+44640*T23-23664* #T9-6096*T+1668*T1-21816*T8+3378-100800*T18+76536*T12)+T14*(-4800*T #19+780*T1+1000-8520*T8-32640*T18-720*T-4800*T9+8960*T17+14400*T23+ #26280*T12)+T3*(8712*T12-13440*T18+56960*T17+6720*T23+108*T1-887-38 #928*T9-1992*T8-28080*T19+10368*T)+R*(14400*T18-180*T1-9720*T12-320 #00*T17-7200*T23+70-4080*T+19680*T9+16800*T19+2520*T8)+85+360*T-840 #*T8+2400*T23-3360*T9-3600*T19-4800*T18+3240*T12+6400*T17+60*T1 T38 = 1/(R-1) T40 = 1/(5-15*R+14*T3+30*T14+93*T4+33*T25) T41 = T40*T38 T42 = R+1 T43 = T42**2 T45 = 1/T43/T42 T46 = T45*T41 T48 = T1+16 T65 = 165*T6+630*T7+T5*(-2055+240*T1+9960*T-480*T8-7920*T9)+T15*(2 #376*T12-20736*T9+25296*T-3816*T8-6776+1116*T1)+T21*(69264*T9-5736* #T8+936*T12+3840*T18+1596*T1-686+39600*T19-90240*T17-17664*T)+T25*( #-15840*T23+5076-516*T1-192000*T17-36504*T12+10104*T8+43200*T18-481 #44*T+148704*T9+87840*T19)+T4*(57600*T17-30960*T19-44640*T23-23664* #T9-6096*T-1668*T1+21816*T8+3378+100800*T18-76536*T12)+T14*(-4800*T #19-780*T1+1000+8520*T8+32640*T18-720*T-4800*T9+8960*T17-14400*T23- #26280*T12)+T3*(-8712*T12+13440*T18+56960*T17-6720*T23-108*T1-887-3 #8928*T9+1992*T8-28080*T19+10368*T)+R*(-14400*T18+180*T1+9720*T12-3 #2000*T17+7200*T23+70-4080*T+19680*T9+16800*T19-2520*T8)+85+360*T+8 #40*T8-2400*T23-3360*T9-3600*T19+4800*T18-3240*T12+6400*T17-60*T1 T84 = 33*T6+126*T7+T5*(-123+1416*T-1584*T9)+T15*(480*T-1296*T9-16) #+T21*(1778-13440*T17-10416*T+14976*T9+7920*T19)+T25*(13440*T17+396 #+2496*T-14064*T9-1440*T19)+T4*(-96576*T9-59760*T19+24960*T-1326+13 #2480*T17)+T14*(40960*T17-18240*T19-32496*T9-736+10080*T)+T3*(-307+ #4464*T+27520*T17-18240*T9-13680*T19)+R*(230-3840*T+15600*T9-23680* #T17+12000*T19)-55+1080*T-3600*T19-4560*T9+7040*T17 T86 = T45*T40 T91 = T15*(-5256*T9-460+5256*T) T95 = T25*(-5304*T+27624*T9+172+22320*T19-44640*T17) T99 = T14*(-6720*T17+76-1128*T+3360*T19+4488*T9) T103 = R*(20-360*T+1200*T19-2400*T17+1560*T9) T104 = 33*T6+192*T7+T5*(141+1590*T-1584*T9-3*T1+6*T8)+T91+T21*(-15 #888*T17+1248*T+72*T12-48*T18+7920*T19+6696*T9+12*T1-48*T8-258)+T95 #+T4*(-14256*T17-2436*T+108*T8+7200*T19-18*T1+144*T18-216*T12+9528* #T9+150)+T99+T3*(-4584*T9-3600*T19+1104*T-96*T8-63+7056*T17+216*T12 #-144*T18+12*T1)+T103-3+30*T-72*T12-3*T1+48*T18+48*T17+30*T8-72*T9 T106 = 1/R T107 = T106*T86 T117 = 33*T6+192*T7+T5*(141+1590*T-1584*T9+3*T1-6*T8)+T91+T21*(-15 #888*T17+1248*T-72*T12+48*T18+7920*T19+6696*T9-12*T1+48*T8-258)+T95 #+T4*(-14256*T17-2436*T-108*T8+7200*T19+18*T1-144*T18+216*T12+9528* #T9+150)+T99+T3*(-4584*T9-3600*T19+1104*T+96*T8-63+7056*T17-216*T12 #+144*T18-12*T1)+T103-3+30*T+72*T12+3*T1-48*T18+48*T17-30*T8-72*T9 T128 = 33*T6+192*T7+T5*(159+1554*T-1584*T9)+T91+T21*(-330+7920*T19 #-15600*T17+1536*T+6264*T9)+T95+T4*(10824*T9+7200*T19-3084*T-15120* #T17+258)+T99+T3*(-135-3600*T19+7920*T17+1680*T-5880*T9)+T103+15-15 #0*T+360*T9-240*T17 T148 = 165*T6+630*T7+T5*(-495-240*T1+480*T8-7920*T9+6840*T)+T15*(2 #376*T12-15984*T9+14304*T-936*T8-324*T1-2864)+T21*(636*T1+12456*T12 #-3840*T18+61056*T9-7656*T8+39600*T19-75840*T17+1546-23856*T)+T25*( #2364-119040*T17-15840*T23+444*T1-31296*T+20160*T18-3336*T8+56160*T #19+95856*T9-1944*T12)+T4*(-120576*T9+7416*T8+77760*T18+221760*T17- #198-120240*T19-41976*T12+19296*T-228*T1-44640*T23)+T14*(-960-33600 #*T19+3720*T8-14400*T23-36240*T9-300*T1+60160*T17-14760*T12+24960*T #18+10080*T)+T3*(13440*T18-108*T1+82240*T17-6720*T23+12912*T-53952* #T9-8712*T12+1992*T8-41520*T19-783)+R*(-14400*T18+180*T1-64000*T17+ #31200*T19+7200*T23+9720*T12+43920*T9+830-11520*T-2520*T8)+3240*T+8 #40*T8-235-12240*T9-8400*T19-60*T1+4800*T18-2400*T23+17600*T17-3240 #*T12 T167 = 165*T6+630*T7+T5*(-495+240*T1-480*T8-7920*T9+6840*T)+T15*(- #2376*T12-15984*T9+14304*T+936*T8+324*T1-2864)+T21*(-636*T1-12456*T #12+3840*T18+61056*T9+7656*T8+39600*T19-75840*T17+1546-23856*T)+T25 #*(2364-119040*T17+15840*T23-444*T1-31296*T-20160*T18+3336*T8+56160 #*T19+95856*T9+1944*T12)+T4*(-120576*T9-7416*T8-77760*T18+221760*T1 #7-198-120240*T19+41976*T12+19296*T+228*T1+44640*T23)+T14*(-960-336 #00*T19-3720*T8+14400*T23-36240*T9+300*T1+60160*T17+14760*T12-24960 #*T18+10080*T)+T3*(-13440*T18+108*T1+82240*T17+6720*T23+12912*T-539 #52*T9+8712*T12-1992*T8-41520*T19-783)+R*(14400*T18-180*T1-64000*T1 #7+31200*T19-7200*T23-9720*T12+43920*T9+830-11520*T+2520*T8)+3240*T #-840*T8-235-12240*T9-8400*T19+60*T1-4800*T18+2400*T23+17600*T17+32 #40*T12 T184 = 33*T7+93*T5+T15*(96+792*T-1584*T9)+T21*(-280+3192*T-4464*T9 #)+T25*(-174+1728*T9-10560*T17+1224*T+7920*T19)+T4*(19776*T9-3480*T #+22320*T19+114-37440*T17)+T14*(-1560*T+6480*T9-11520*T17+7200*T19+ #120)+T3*(48+3360*T19-888*T-6400*T17+3984*T9)+R*(1080*T-75+6720*T17 #-4320*T9-3600*T19)+25-360*T+1440*T9-2240*T17+1200*T19 W(1) = T46*T35*T2/360 W(2) = -T46*T65*T48/360 W(3) = -T86*T38*T84/18 W(4) = -T107*T104*T2/36 W(5) = T107*T117*T48/36 W(6) = T106*T45*T40*T128/9 W(7) = T46*T148*T2/360 W(8) = -T46*T167*T48/360 W(9) = -T41/T43*T184/18 RETURN C ---------------------------- C GENERAL CASE C ---------------------------- 50 CONTINUE T1 = DSQRT(6.D0) T2 = -16+T1 T3 = R2**2 T4 = T3*R2 T5 = T3**2 T6 = T5*T4 T7 = T6*T T8 = T*T5 T9 = T5*R2 T10 = T*T9 T11 = T**2 T12 = T11**2 T13 = T4*T12 T14 = T3*T12 T15 = T11*T T16 = T9*T15 T17 = T15*T5 T18 = T1*T12 T19 = T1*T9 T20 = T1*T15 T21 = T4*T11 T22 = T1*T5 T23 = 7500*T17+600*T18-360*T7-2988*T8-132*T10-2400*T13-240*T14-270 #0*T16+819*T19-120*T20+990*T4-1800*T21-1053*T22+360*T11 T24 = T9*T11 T25 = T5*T3 T26 = T20*T5 T27 = T1*T25 T28 = T1*T3 T29 = T28*T12 T30 = T3*T15 T31 = R2*T T32 = T1*T T33 = T32*T5 T34 = T28*T15 T35 = T1*T4 T36 = T35*T15 T37 = T1*T11 T38 = 2034*T9-5364*T24-669*T25+4500*T26-2178*T5+180*T15-264*T27-14 #4*T3+4035*T29-2040*T30-180*T31-4428*T33+8760*T34-11520*T36-270*T37 T40 = T25*T11 T41 = T3*T T42 = T35*T12 T43 = T32*T25 T44 = T35*T11 T45 = T22*T11 T46 = T32*T9 T47 = T19*T11 T48 = T11*T3 T49 = T15*R2 T50 = T5**2 T51 = T32*R2 T52 = T37*R2 T53 = 1890*T40-96*T41-1650*T42+528*T43-1170*T44+6426*T45+498*T46-3 #204*T47+2250*T48+1200*T49+10*T50+60*T6+450*T51+2430*T52 T54 = T37*T3 T55 = T32*T4 T56 = T4*T15 T57 = T5*T12 T58 = T32*T3 T59 = T5*T11 T60 = T*T25 T61 = R2*T11 T62 = T4*T T63 = T12*R2 T64 = T18*R2 T65 = T20*R2 T66 = -4230*T54+6120*T55-4020*T56-825*T12+1125*T57-2976*T58+5526*T #59+768*T60-1980*T61+2070*T62-99*T28+2250*T63+540*T35-2925*T64-1800 #*T65 T69 = R1**2 T70 = T69*R1 T72 = T6*T15 T73 = T50*T T74 = T1*T6 T75 = T9*T12 T76 = T6*T11 T77 = -4560*T17-150*T18+96*T7+2232*T8-1872*T10+1170*T13+990*T14+28 #20*T16-60*T73-351*T19+120*T20-78*T74+450*T75+540*T76-108*T4 T78 = 360*T21+198*T22-846*T9+3096*T24+549*T25-6360*T26+468*T5-180* #T15+204*T27-3510*T29-60*T30+3222*T33-1560*T34+6240*T36-1338*T40 T80 = T27*T11 T81 = 126*T41+3270*T42+444*T43-3690*T44-54*T45-1872*T46+2196*T47-7 #38*T48+480*T49+10*T50-198*T6-540*T52+2682*T54-1278*T80-1872*T55 T82 = T19*T15 T83 = T22*T12 T84 = T25*T15 T85 = T74*T T86 = 1140*T56+225*T12-1605*T57+396*T58-414*T59+2220*T82-1005*T83+ #114*T60-1000*T84+360*T61-1062*T62-1200*T63-48*T35+1425*T64-420*T65 #+156*T85 T90 = T50*T11 T91 = -36*T7+18*T8-54*T10+270*T13-75*T14+420*T16+(T23+T38+T53+T66) #*T70+120*T72+24*T73+(T77+T78+T81+T86)*T69-36*T90+150*T75-36*T76 T92 = T25*T12 T94 = T69**2 T95 = T94**2 T96 = T95*T69 T97 = (-120*T+50*R2+20)*T96 T99 = 510*T11+70*T9-30*T1+145*T25-2274*T5-1820*T15-2023*T3-660*T32 #+5820*T30-1500*T31+1788*T33 T101 = 4620*T34-540*T37+132*T41-6084*T44+2646*T48+1440*T49+3660*T5 #1-1890*T52+8496*T54-264*T55-5500*T56 T102 = T1*R2 T103 = -225*T12-4608*T58+6750*T59+540*T61-1164*T62-1683*T28-600*T6 #3+2160*T35-750*T64-6060*T65+465*T102 T106 = T94*R1 T108 = 1140*T17+180*T7-540*T8+918*T10-540*T13-105*T14-1980*T16-200 #*T72+24*T73+27*T19+18*T74+60*T90-240*T75-264*T76 T109 = T1*T50 T110 = 150*T92-12*T109+468*T21-9*T22+72*T9+36*T24-144*T25+2340*T26 #-24*T5-54*T27+1095*T29+180*T30-540*T33-420*T34-480*T36 T112 = 546*T40-1890*T42-396*T43+1278*T44-1674*T45+918*T46+36*T47-7 #2*T48-120*T49-32*T50+48*T6-342*T54+216*T80+144*T55-480*T56 T113 = T74*T11 T114 = T19*T12 T115 = T109*T T116 = T27*T15 T117 = 1185*T57-684*T59-324*T113-1680*T82+1260*T83-390*T114+24*T11 #5-396*T60+720*T116+520*T84+144*T62+150*T63-225*T64+180*T65+180*T85 T121 = -5+350*T5+150*T12+3840*T48-70*T4-1400*T49-2730*T62+810*T37+ #1104*T58+760*T15-360*T32+360*T20 T122 = 730*R2+675*T102-60*T-270*T51-1022*T3+2724*T41-180*T1-2640*T #61-180*T31-390*T11-1620*T52-552*T28 T124 = T94*T70 T126 = -60*T32+120*T1-65*T3+240*T+95+960*T31-400*T15+1440*T61+390* #T51-1500*T41-600*T11+280*T4-195*T102-270*T37-390*R2 T128 = T27*T12 T130 = T95*R1 T132 = T74*T15 T133 = -90*T92+T97+(600*T18+4548*T8-2310*T10+1125*T14+1680*T20+378 #0*T4-11004*T21-894*T22+660*T+270*R2+T99+T101+T103)*T106+(T108+T110 #+T112+T117)*R1+(T121+T122)*T124+T126*T95-72*T21-90*T128-324*T24+(- #30*T1-105-540*T31+120*T+60*T32+360*T11+145*T3)*T130+120*T132-75*T2 #9+60*T30 T135 = 18*T33+60*T34-120*T36-18*T40+270*T42+108*T43-72*T44+216*T45 #-54*T46-324*T47-18*T80-120*T56-405*T57 T136 = T95*T70 T137 = -15-150*T18-3150*T8-1320*T20-1938*T4+6480*T21-660*T-600*R2+ #810*T11+280*T9+120*T1+480*T15+2240*T3+840*T32-3600*T30+360*T31 T138 = -540*T37-1074*T41-6804*T48+2280*T49-1980*T51+4050*T52-4194* #T54+1776*T55-300*T12-444*T58+180*T61+4506*T62+1620*T28+450*T63-888 #*T35+1980*T65-855*T102 T140 = T94*T69 T142 = -5000*T17-900*T18-480*T8+2682*T10+1500*T13-1695*T14-591*T19 #-720*T20-2934*T4+5568*T21+1728*T22-180*T-60*R2 T143 = -1050*T11-1586*T9+4380*T24+95*T25+3678*T5+980*T15+869*T3-15 #45*T29+180*T32-360*T30+1080*T31+180*T33-11460*T34 T145 = 5880*T36+810*T37-312*T41+9738*T44-5454*T45+1182*T46-1122*T4 #8-3880*T49+50*T6-2250*T51-2430*T52-2412*T54-5904*T55 T146 = 8980*T56+975*T12+6528*T58-10434*T59-1110*T60+2100*T61-1764* #T62+714*T28-1050*T63-1764*T35+2475*T64+6120*T65-90*T102 T150 = T109*T11 T151 = 216*T59-36*T113+420*T82-405*T83+150*T114+24*T115+10*T136+(T #137+T138)*T140+108*T60-120*T116-120*T84-36*T85+(T142+T143+T145+T14 #6)*T94-36*T150 T169 = 2*T95+(8*R2-8)*T124+(-46*R2+21*T3+15)*T140+(99*R2-17+35*T4- #117*T3)*T106+(-101*R2+225*T3+35*T5+11-170*T4)*T94+(270*T4+49*R2-3- #192*T3-145*T5+21*T9)*T70+(-72*T9+180*T5-9*R2-177*T4+8*T25+72*T3)*T #69+(45*T4+2*T6-20*T25-81*T5+63*T9-9*T3)*R1-18*T9+9*T5+6*T25-4*T6-3 #*T4 T170 = 1/T169 T172 = 1/(R1+R2) T174 = 1/R1 T175 = T170*T172*T174 T177 = 16+T1 T178 = 60*T32-120*T1-65*T3+240*T+95+960*T31-400*T15+1440*T61-390*T #51-1500*T41-600*T11+280*T4+195*T102+270*T37-390*R2 T180 = T178*T95-36*T7+18*T8-54*T10+270*T13-75*T14+420*T16+120*T72+ #24*T73-36*T90+150*T75-36*T76-90*T92 T181 = 7500*T17-600*T18-360*T7-2988*T8-132*T10-2400*T13-240*T14-27 #00*T16-819*T19+120*T20+990*T4-1800*T21+1053*T22+360*T11 T182 = 2034*T9-5364*T24-669*T25-4500*T26-2178*T5+180*T15+264*T27-1 #44*T3-4035*T29-2040*T30-180*T31+4428*T33-8760*T34+11520*T36+270*T3 #7 T184 = 1890*T40-96*T41+1650*T42-528*T43+1170*T44-6426*T45-498*T46+ #3204*T47+2250*T48+1200*T49+10*T50+60*T6-450*T51-2430*T52 T185 = 4230*T54-6120*T55-4020*T56-825*T12+1125*T57+2976*T58+5526*T #59+768*T60-1980*T61+2070*T62+99*T28+2250*T63-540*T35+2925*T64+1800 #*T65 T189 = -5+350*T5+150*T12+3840*T48-70*T4-1400*T49-2730*T62-810*T37- #1104*T58+760*T15+360*T32-360*T20 T190 = 730*R2-675*T102-60*T+270*T51-1022*T3+2724*T41+180*T1-2640*T #61-180*T31-390*T11+1620*T52+552*T28 T193 = T97-72*T21+(T181+T182+T184+T185)*T70+(T189+T190)*T124+90*T1 #28-324*T24-120*T132+75*T29+60*T30-18*T33-60*T34+120*T36-18*T40 T195 = -270*T42-108*T43+72*T44-216*T45+54*T46+324*T47+18*T80-120*T #56-405*T57+216*T59+36*T113-420*T82+405*T83 T197 = 510*T11+70*T9+30*T1+145*T25-2274*T5-1820*T15-2023*T3+660*T3 #2+5820*T30-1500*T31-1788*T33 T199 = -4620*T34+540*T37+132*T41+6084*T44+2646*T48+1440*T49-3660*T #51+1890*T52-8496*T54+264*T55-5500*T56 T200 = -225*T12+4608*T58+6750*T59+540*T61-1164*T62+1683*T28-600*T6 #3-2160*T35+750*T64+6060*T65-465*T102 T204 = -5000*T17+900*T18-480*T8+2682*T10+1500*T13-1695*T14+591*T19 #+720*T20-2934*T4+5568*T21-1728*T22-180*T-60*R2 T205 = -1050*T11-1586*T9+4380*T24+95*T25+3678*T5+980*T15+869*T3+15 #45*T29-180*T32-360*T30+1080*T31-180*T33+11460*T34 T207 = -5880*T36-810*T37-312*T41-9738*T44+5454*T45-1182*T46-1122*T #48-3880*T49+50*T6+2250*T51+2430*T52+2412*T54+5904*T55 T208 = 8980*T56+975*T12-6528*T58-10434*T59-1110*T60+2100*T61-1764* #T62-714*T28-1050*T63+1764*T35-2475*T64-6120*T65+90*T102 T214 = -15+150*T18-3150*T8+1320*T20-1938*T4+6480*T21-660*T-600*R2+ #810*T11+280*T9-120*T1+480*T15+2240*T3-840*T32-3600*T30+360*T31 T215 = 540*T37-1074*T41-6804*T48+2280*T49+1980*T51-4050*T52+4194*T #54-1776*T55-300*T12+444*T58+180*T61+4506*T62-1620*T28+450*T63+888* #T35-1980*T65+855*T102 T218 = 1140*T17+180*T7-540*T8+918*T10-540*T13-105*T14-1980*T16-200 #*T72+24*T73-27*T19-18*T74+60*T90-240*T75-264*T76 T219 = 150*T92+12*T109+468*T21+9*T22+72*T9+36*T24-144*T25-2340*T26 #-24*T5+54*T27-1095*T29+180*T30+540*T33+420*T34+480*T36 T221 = 546*T40+1890*T42+396*T43-1278*T44+1674*T45-918*T46-36*T47-7 #2*T48-120*T49-32*T50+48*T6+342*T54-216*T80-144*T55-480*T56 T222 = 1185*T57-684*T59+324*T113+1680*T82-1260*T83+390*T114-24*T11 #5-396*T60-720*T116+520*T84+144*T62+150*T63+225*T64-180*T65-180*T85 T226 = -4560*T17+150*T18+96*T7+2232*T8-1872*T10+1170*T13+990*T14+2 #820*T16-60*T73+351*T19-120*T20+78*T74+450*T75+540*T76-108*T4 T227 = 360*T21-198*T22-846*T9+3096*T24+549*T25+6360*T26+468*T5-180 #*T15-204*T27+3510*T29-60*T30-3222*T33+1560*T34-6240*T36-1338*T40 T229 = 126*T41-3270*T42-444*T43+3690*T44+54*T45+1872*T46-2196*T47- #738*T48+480*T49+10*T50-198*T6+540*T52-2682*T54+1278*T80+1872*T55 T230 = 1140*T56+225*T12-1605*T57-396*T58-414*T59-2220*T82+1005*T83 #+114*T60-1000*T84+360*T61-1062*T62-1200*T63+48*T35-1425*T64+420*T6 #5-156*T85 T234 = -150*T114-24*T115+10*T136+108*T60+120*T116-120*T84+36*T85+( #-600*T18+4548*T8-2310*T10+1125*T14-1680*T20+3780*T4-11004*T21+894* #T22+660*T+270*R2+T197+T199+T200)*T106+(T204+T205+T207+T208)*T94+36 #*T150+(30*T1-105-540*T31+120*T-60*T32+360*T11+145*T3)*T130+(T214+T #215)*T140+(T218+T219+T221+T222)*R1+(T226+T227+T229+T230)*T69 T239 = 6720*T17-168*T7-3420*T8+1872*T10-3690*T13+4410*T14-2100*T16 #-12*T73+90*T75+108*T76+36*T4+4500*T21+252*T9-2016*T24-135*T25-144* #T5 T240 = -180*T15+1860*T30+1266*T40-450*T41-3366*T48+600*T49+2*T50+5 #4*T6-7260*T56+225*T12+885*T57-18*T59-510*T60-200*T84+720*T61+2034* #T62-1950*T63 T243 = -4824*T61+906*R2+1680*T15-1230*T62+1296*T21+3672*T48+678*T4 #+810*T11+120*T12+56*T9-1140*T-147+318*T41-630*T8+2448*T31-1920*T49 #-720*T30+90*T63-1496*T3 T245 = -280*T15+30*T12-546*T62-1050*T11+215+458*T3+1416*T61-14*T4+ #420*T-280*T49+288*T31+70*T5-664*R2-780*T41+768*T48 T247 = 36-1236*T8-462*T10+225*T14-1836*T4+5100*T21+924*T-504*R2+75 #0*T11+14*T9+29*T25+618*T5 T248 = -2380*T15+1615*T3-4380*T30-4692*T31+5556*T41-9666*T48+7560* #T49-1100*T56-765*T12+1350*T59+2376*T61+84*T62+780*T63 T251 = -2580*T17-180*T7+540*T8-918*T10+2160*T13-1335*T14+1620*T16- #40*T72-24*T73+12*T90+420*T75+336*T76+30*T92-1440*T21-18*T9 T252 = -36*T24+36*T25+6*T5+540*T30-150*T40+396*T48-240*T49+8*T50-1 #2*T6+480*T56-1275*T57+1872*T59+396*T60-760*T84-144*T62+300*T63 T255 = -1000*T17-312*T8-882*T10+300*T13+1515*T14+1530*T4-10572*T21 #-252*T+96*R2-1182*T11+392*T9+876*T24+19*T25-1338*T5 T256 = 1060*T15-683*T3+13680*T30+2916*T31-7896*T41+2670*T48-8120*T #49+10*T6-5260*T56+1275*T12+4458*T59-222*T60+3336*T61+6732*T62-3180 #*T63 T261 = (T239+T240)*T69+T243*T140+T245*T124+(T247+T248)*T106+(T251+ #T252)*R1+36*T7-18*T8+54*T10-270*T13+75*T14-420*T16+(T255+T256)*T94 #-120*T72-24*T73+36*T90+(56*T4-276*T31+120*T-13*T3-80*T15-300*T41+2 #04*T11-125+156*R2+288*T61)*T95 T266 = -3900*T17-72*T7+4716*T8-624*T10+1500*T13-4890*T14-540*T16-4 #50*T4+1044*T21+396*T11-576*T9+2772*T24+183*T25+828*T5+180*T15+90*T #3 T267 = -10920*T30-576*T31+378*T40+3552*T41+5526*T48+2400*T49+2*T50 #+12*T6+13020*T56-885*T12+225*T57-6606*T59-480*T60-3312*T61-6930*T6 #2+3960*T63 T270 = -150*T75+36*T76+90*T92+(4+10*R2-24*T)*T96+(-108*T31-48*T+72 #*T11+29*T3+15)*T130+(T266+T267)*T70+72*T21+324*T24-60*T30+18*T40+1 #20*T56+405*T57-216*T59+2*T136-108*T60+120*T84 T276 = T109*T12 T277 = T50*R2 T278 = T1*T277 T279 = T278*T15 T280 = -1980*T17-102*T7+648*T8-360*T10+1920*T13-2205*T14+180*T16-5 #40*T72+48*T73-33*T19-20*T74+378*T90-330*T75-528*T76+225*T92-11*T10 #9+2*T277-2*T278 T281 = T277*T T282 = T50*T3 T283 = T1*T282 T284 = -72*T281-6*T4-1908*T21+2*T282-2*T283+27*T22-225*T128-12*T9+ #252*T24+15*T25-2760*T26+18*T5+60*T15+540*T132-1560*T29-240*T30+972 #*T33-420*T34+2460*T36 T286 = 252*T40+126*T41+1845*T42+252*T43-1980*T44+936*T45-486*T46+7 #2*T47+1314*T48-360*T49+11*T50+20*T6-126*T52+1062*T54-414*T80-630*T #55+2040*T56-75*T12-465*T57 T287 = T278*T T288 = 126*T58+792*T59+528*T113+1200*T82-765*T83+405*T114-48*T115- #144*T60-1020*T116+960*T84-252*T61-504*T62+870*T63-9*T35+360*T64-12 #0*T65+102*T85-378*T150+72*T287 T292 = T283*T11 T293 = T50*T12 T294 = 1400*T17+324*T7-1206*T8+168*T10-1170*T13+2820*T14+1160*T16- #222*T73+93*T19-6*T74+300*T75+876*T76+10*T109+10*T277-10*T278+78*T4 #+1392*T21 T295 = -162*T22-108*T11+60*T9-24*T24+12*T25+2800*T26-108*T5-120*T1 #5-36*T27-18*T3+1800*T29+3160*T30+144*T31-1638*T33+2600*T34-4660*T3 #6-1320*T40-936*T41 T297 = -1320*T42+60*T43+1038*T44+972*T45+456*T46-1056*T47-2604*T48 #-440*T49-10*T50+6*T6+72*T51+498*T52-1788*T54+1392*T80+1944*T55-464 #0*T56+330*T12-150*T57 T298 = -792*T58+756*T59-876*T113-1580*T82+480*T83-300*T114+222*T11 #5-60*T60+1000*T116-1000*T84+1212*T61+1728*T62-18*T28-1830*T63+99*T #35-660*T64-280*T65-324*T85 T302 = 210*T13-75*T14-220*T72-30*T276+40*T279+66*T90+120*T75+120*T #76+165*T92+(T280+T284+T286+T288)*T69-12*T292+30*T293+(T294+T295+T2 #97+T298)*T70 T303 = -36*T21-165*T128-72*T24+220*T132-75*T29+60*T30+60*T34-120*T #36+90*T40+240*T42-54*T44+162*T45-198*T47 T305 = T50*T15 T306 = -108*T62-32*T102-72*T54-24*T51+29*T5+39+192*T61+6*T41+84*T3 #1+72*T48-106*R2 T307 = 132*T52-4*T3-144*T58+73*T28-120*T15-10*T4+78*T+10*T35-29*T2 #2-144*T11+108*T55 T310 = T282*T11 T313 = -40*T305-120*T56+(T306+T307)*T95-180*T57+108*T59-120*T113+2 #40*T82-270*T83+30*T114+120*T116-240*T84+12*T310+(-12*T-2*T28+2*T3+ #6)*T96 T314 = 600*T17-120*T7-108*T8+72*T10-1110*T13+750*T14-600*T16+200*T #72-66*T73-12*T90+30*T75+300*T76-210*T92-12*T281+612*T21 T315 = 210*T128+360*T24+1020*T26-200*T132+600*T29-360*T30-162*T33- #240*T34-180*T36+360*T40-1140*T42+738*T44-972*T45+198*T46+270*T47 T317 = T282*T T318 = -180*T48+120*T49-180*T54-200*T305-360*T80+54*T55+600*T57-12 #*T317-648*T59-300*T113-600*T82+750*T83-255*T114+66*T115-90*T60 T319 = T6*T12 T320 = T277*T11 T321 = T109*T15 T322 = T278*T11 T323 = T283*T T324 = T74*T12 T325 = 540*T116-360*T84+36*T62-150*T63-75*T64+60*T65+120*T85+90*T3 #19+108*T320+200*T321-108*T322+12*T323+12*T150+12*T287-90*T324 T329 = 10*T4-36*T31-24*T51+2*T3+24*T58-18-12*T+18*R2-2*T28-10*T35- #24*T41+72*T11+12*T102 T331 = -54-300*T8-56*T19-24*T4+288*T21+42*T22-216*T+198*R2+180*T11 #+56*T9-42*T5+360*T15-212*T3-80*T30+432*T31+300*T33 T332 = 80*T34+384*T41-288*T44+24*T48-280*T49+108*T51-228*T52+696*T #54-384*T55+60*T12-96*T58-1032*T61+192*T62-172*T28+120*T35-200*T65+ #72*T102 T335 = -6-720*T17-6*T8+660*T10+90*T13+60*T14+22*T19-56*T74+430*T4- #1908*T21-264*T22-252*T+106*R2 T336 = -396*T11+20*T9+1296*T24-70*T25+720*T26-126*T5+720*T15+70*T2 #7-326*T3+360*T29+3240*T30+1380*T31+210*T33+2400*T34 T338 = -1500*T36-2460*T41-90*T42+630*T43-2412*T44+1788*T45-744*T46 #-1296*T47+2160*T48-3240*T49+56*T6+564*T51-60*T52+1080*T54 T339 = 1188*T55+120*T56+420*T12-1500*T58-924*T59-630*T60+600*T61+6 #72*T62-268*T28-930*T63+425*T35-270*T64-1020*T65+47*T102 T343 = T277*T15 T344 = 620*T17-462*T7+498*T8-126*T10+60*T13-1380*T14-1100*T16-108* #T19+42*T74-29*T109-270*T4+1944*T21+306*T22+54*T-18*R2+414*T11 T345 = -42*T9-1452*T24+3*T25-1700*T26+225*T5-240*T15+6*T27+135*T3- #1095*T29-5680*T30-780*T31+960*T33-3560*T34+3720*T36+1350*T40+2208* #T41 T347 = 465*T42-642*T43+1566*T44-2190*T45+162*T46+1848*T47+1380*T48 #+2320*T49+29*T50-42*T6-336*T51-582*T52+420*T54-1350*T80-2316*T55+2 #880*T56 T348 = -525*T12+225*T57+1698*T58-780*T59+1100*T82-225*T83+624*T60- #1668*T61-2364*T62+117*T28+1920*T63-315*T35+555*T64+1040*T65+462*T8 #5-9*T102 T352 = 33+450*T8-546*T10+30*T14+70*T19-208*T4-384*T21+78*T22+360*T #-198*R2-18*T11 T353 = -70*T9+70*T25+6*T5-660*T15-70*T27+421*T3-30*T29-100*T30-122 #4*T31-618*T33-920*T34 T355 = 280*T36+696*T41+1392*T44-768*T45+546*T46-2166*T48+1880*T49- #360*T51+366*T52-1218*T54+36*T55 T356 = -280*T56-210*T12+684*T58+768*T59+948*T61+348*T62+272*T28+12 #0*T63-320*T35+90*T64+520*T65-90*T102 T360 = (T314+T315+T318+T325)*R1+T329*T130+30*T319+12*T320+40*T321- #12*T322+(T331+T332)*T124-66*T150+(T335+T336+T338+T339)*T106-30*T32 #4-40*T343+(T344+T345+T347+T348)*T94+(T352+T353+T355+T356)*T140 T364 = 1/R2 T367 = 1/(R1-1) T369 = T367*T170*T172 T372 = 33+450*T8-546*T10+30*T14-70*T19-208*T4-384*T21-78*T22+360*T #-198*R2-18*T11 T373 = -70*T9+70*T25+6*T5-660*T15+70*T27+421*T3+30*T29-100*T30-122 #4*T31+618*T33+920*T34 T375 = -280*T36+696*T41-1392*T44+768*T45-546*T46-2166*T48+1880*T49 #+360*T51-366*T52+1218*T54-36*T55 T376 = -280*T56-210*T12-684*T58+768*T59+948*T61+348*T62-272*T28+12 #0*T63+320*T35-90*T64-520*T65+90*T102 T380 = -54-300*T8+56*T19-24*T4+288*T21-42*T22-216*T+198*R2+180*T11 #+56*T9-42*T5+360*T15-212*T3-80*T30+432*T31-300*T33 T381 = -80*T34+384*T41+288*T44+24*T48-280*T49-108*T51+228*T52-696* #T54+384*T55+60*T12+96*T58-1032*T61+192*T62+172*T28-120*T35+200*T65 #-72*T102 T384 = 210*T13-75*T14-220*T72+30*T276-40*T279+66*T90+120*T75+120*T #76+165*T92+12*T292+30*T293+(T372+T373+T375+T376)*T140+(T380+T381)* #T124 T385 = -108*T62+32*T102+72*T54+24*T51+29*T5+39+192*T61+6*T41+84*T3 #1+72*T48-106*R2 T386 = -132*T52-4*T3+144*T58-73*T28-120*T15-10*T4+78*T-10*T35+29*T #22-144*T11-108*T55 T391 = 10*T4-36*T31+24*T51+2*T3-24*T58-18-12*T+18*R2+2*T28+10*T35- #24*T41+72*T11-12*T102 T393 = -1980*T17-102*T7+648*T8-360*T10+1920*T13-2205*T14+180*T16-5 #40*T72+48*T73+33*T19+20*T74+378*T90-330*T75-528*T76+225*T92+11*T10 #9+2*T277+2*T278 T394 = -72*T281-6*T4-1908*T21+2*T282+2*T283-27*T22+225*T128-12*T9+ #252*T24+15*T25+2760*T26+18*T5+60*T15-540*T132+1560*T29-240*T30-972 #*T33+420*T34-2460*T36 T396 = 252*T40+126*T41-1845*T42-252*T43+1980*T44-936*T45+486*T46-7 #2*T47+1314*T48-360*T49+11*T50+20*T6+126*T52-1062*T54+414*T80+630*T #55+2040*T56-75*T12-465*T57 T397 = -126*T58+792*T59-528*T113-1200*T82+765*T83-405*T114+48*T115 #-144*T60+1020*T116+960*T84-252*T61-504*T62+870*T63+9*T35-360*T64+1 #20*T65-102*T85+378*T150-72*T287 T401 = 1400*T17+324*T7-1206*T8+168*T10-1170*T13+2820*T14+1160*T16- #222*T73-93*T19+6*T74+300*T75+876*T76-10*T109+10*T277+10*T278+78*T4 #+1392*T21 T402 = 162*T22-108*T11+60*T9-24*T24+12*T25-2800*T26-108*T5-120*T15 #+36*T27-18*T3-1800*T29+3160*T30+144*T31+1638*T33-2600*T34+4660*T36 #-1320*T40-936*T41 T404 = 1320*T42-60*T43-1038*T44-972*T45-456*T46+1056*T47-2604*T48- #440*T49-10*T50+6*T6-72*T51-498*T52+1788*T54-1392*T80-1944*T55-4640 #*T56+330*T12-150*T57 T405 = 792*T58+756*T59+876*T113+1580*T82-480*T83+300*T114-222*T115 #-60*T60-1000*T116-1000*T84+1212*T61+1728*T62+18*T28-1830*T63-99*T3 #5+660*T64+280*T65+324*T85 T409 = 620*T17-462*T7+498*T8-126*T10+60*T13-1380*T14-1100*T16+108* #T19-42*T74+29*T109-270*T4+1944*T21-306*T22+54*T-18*R2+414*T11 T410 = -42*T9-1452*T24+3*T25+1700*T26+225*T5-240*T15-6*T27+135*T3+ #1095*T29-5680*T30-780*T31-960*T33+3560*T34-3720*T36+1350*T40+2208* #T41 T412 = -465*T42+642*T43-1566*T44+2190*T45-162*T46-1848*T47+1380*T4 #8+2320*T49+29*T50-42*T6+336*T51+582*T52-420*T54+1350*T80+2316*T55+ #2880*T56 T413 = -525*T12+225*T57-1698*T58-780*T59-1100*T82+225*T83+624*T60- #1668*T61-2364*T62-117*T28+1920*T63+315*T35-555*T64-1040*T65-462*T8 #5+9*T102 T417 = -6-720*T17-6*T8+660*T10+90*T13+60*T14-22*T19+56*T74+430*T4- #1908*T21+264*T22-252*T+106*R2 T418 = -396*T11+20*T9+1296*T24-70*T25-720*T26-126*T5+720*T15-70*T2 #7-326*T3-360*T29+3240*T30+1380*T31-210*T33-2400*T34 T420 = 1500*T36-2460*T41+90*T42-630*T43+2412*T44-1788*T45+744*T46+ #1296*T47+2160*T48-3240*T49+56*T6-564*T51+60*T52-1080*T54 T421 = -1188*T55+120*T56+420*T12+1500*T58-924*T59-630*T60+600*T61+ #672*T62+268*T28-930*T63-425*T35+270*T64+1020*T65-47*T102 T425 = (T385+T386)*T95+(-12*T+2*T28+2*T3+6)*T96+T391*T130+(T393+T3 #94+T396+T397)*T69+(T401+T402+T404+T405)*T70-36*T21+(T409+T410+T412 #+T413)*T94+(T417+T418+T420+T421)*T106+165*T128-72*T24-220*T132+75* #T29+60*T30 T427 = -210*T128+360*T24-1020*T26+200*T132-600*T29-360*T30+162*T33 #+240*T34+180*T36+360*T40+1140*T42-738*T44+972*T45-198*T46-270*T47 T429 = -180*T48+120*T49+180*T54-200*T305+360*T80-54*T55+600*T57-12 #*T317-648*T59+300*T113+600*T82-750*T83+255*T114-66*T115-90*T60 T430 = -540*T116-360*T84+36*T62-150*T63+75*T64-60*T65-120*T85+90*T #319+108*T320-200*T321+108*T322-12*T323-12*T150-12*T287+90*T324 T434 = (T314+T427+T429+T430)*R1-60*T34+120*T36+90*T40-240*T42+54*T #44-162*T45+198*T47-40*T305-120*T56-180*T57+108*T59+120*T113 T435 = -240*T82+270*T83-30*T114-120*T116-240*T84+12*T310+30*T319+1 #2*T320-40*T321+12*T322+66*T150+30*T324-40*T343 T441 = 210*T13-75*T14-220*T72+66*T90+120*T75+120*T76+165*T92+30*T2 #93-36*T21-72*T24+60*T30+90*T40-40*T305-120*T56-180*T57 T442 = 620*T17-462*T7+498*T8-126*T10+60*T13-1380*T14-1100*T16-270* #T4+1944*T21+54*T-18*R2+414*T11-42*T9-1452*T24+3*T25+225*T5-240*T15 T443 = 135*T3-5680*T30-780*T31+1350*T40+2208*T41+1380*T48+2320*T49 #+29*T50-42*T6+2880*T56-525*T12+225*T57-780*T59+624*T60-1668*T61-23 #64*T62+1920*T63 T446 = -1032*T61+60*T12-42*T5+192*T62+24*T48-54+288*T21+360*T15+19 #8*R2-80*T30+56*T9+384*T41+432*T31-280*T49-216*T-24*T4-212*T3+180*T #11-300*T8 T450 = -1980*T17-102*T7+648*T8-360*T10+1920*T13-2205*T14+180*T16-5 #40*T72+48*T73+378*T90-330*T75-528*T76+225*T92+2*T277-72*T281-6*T4- #1908*T21+2*T282-12*T9 T451 = 252*T24+15*T25+18*T5+60*T15-240*T30+252*T40+126*T41+1314*T4 #8-360*T49+11*T50+20*T6+2040*T56-75*T12-465*T57+792*T59-144*T60+960 #*T84-252*T61-504*T62+870*T63 T454 = 1400*T17+324*T7-1206*T8+168*T10-1170*T13+2820*T14+1160*T16- #222*T73+300*T75+876*T76+10*T277+78*T4+1392*T21-108*T11+60*T9-24*T2 #4+12*T25-108*T5 T455 = -120*T15-18*T3+3160*T30+144*T31-1320*T40-936*T41-2604*T48-4 #40*T49-10*T50+6*T6-4640*T56+330*T12-150*T57+756*T59-60*T60-1000*T8 #4+1212*T61+1728*T62-1830*T63 T458 = -6-720*T17-6*T8+660*T10+90*T13+60*T14+430*T4-1908*T21-252*T #+106*R2-396*T11+20*T9+1296*T24-70*T25-126*T5 T459 = 720*T15-326*T3+3240*T30+1380*T31-2460*T41+2160*T48-3240*T49 #+56*T6+120*T56+420*T12-924*T59-630*T60+600*T61+672*T62-930*T63 T462 = 33+450*T8-546*T10+30*T14-208*T4-384*T21+360*T-198*R2-18*T11 #-70*T9+70*T25+6*T5 T463 = -660*T15+421*T3-100*T30-1224*T31+696*T41-2166*T48+1880*T49- #280*T56-210*T12+768*T59+948*T61+348*T62+120*T63 T466 = -4*T3-120*T15-106*R2-144*T11+29*T5+78*T+72*T48+84*T31+6*T41 #+39-10*T4+192*T61-108*T62 T470 = 360*T24-360*T30+360*T40-180*T48+120*T49-200*T305+600*T57-12 #*T317-648*T59-90*T60-360*T84+36*T62-150*T63+90*T319+108*T320 T473 = 108*T59+(T442+T443)*T94-240*T84+T446*T124+12*T310+(-36*T31+ #10*T4+2*T3-24*T41-18-12*T+18*R2+72*T11)*T130+(T450+T451)*T69+(T454 #+T455)*T70+(T458+T459)*T106+30*T319+12*T320+(T462+T463)*T140+T466* #T95+(-12*T+6+2*T3)*T96+(T314+T470)*R1-40*T343 T478 = 39*T22+576*T9+11130*T24-550*T25-96*T5+320*T15-45*T29-16860* #T30+18*T33-360*T34+120*T36 T481 = 22720*T56-900*T12+126*T58-27864*T59-3420*T60-450*T61+2688*T #62+3375*T63-15*T35+225*T64+180*T65 T486 = (60*R2-120*T-10)*T96 T487 = 300*T18+12816*T8-5460*T10+1575*T14+120*T20-124*T4-27000*T21 #-18*T22-990*T11-1120*T9+425*T25+942*T5+1120*T15-119*T3+16200*T30+9 #00*T31+36*T33+300*T34-270*T37 T488 = 1176*T41-180*T44+18018*T48-10080*T49+270*T51+270*T52+378*T5 #4+18*T55-9100*T56+1725*T12-324*T58+13230*T59-1800*T61-9702*T62-39* #T28-2700*T63+51*T35-150*T64-780*T65 T491 = 21080*T17-1470*T7+1224*T8-4860*T10-6180*T13+4770*T14-7700*T #16+24*T19-7080*T21-21*T22+24*T9-17952*T24+225*T25+180*T26-66*T5+45 #*T29+6560*T30-216*T33+60*T34-240*T36 T492 = 6270*T40-30*T42+90*T44+252*T45-72*T47+798*T48-440*T49+60*T5 #0-140*T6-72*T54+72*T55-19140*T56+75*T12+2625*T57+19062*T59+4080*T6 #0-108*T62-1500*T63-75*T64+60*T65 T495 = 420*T17+75*T13+720*T16+200*T72+(-10500*T17-150*T18-8802*T8+ #9348*T10+2625*T13-4770*T14-9*T19+120*T20-15*T4+23418*T21+T478+66*T #41+198*T44-144*T45+18*T46-6432*T48+4080*T49+195*T6-270*T52+108*T54 #-162*T55+T481)*T94-60*T90+675*T75+270*T76-150*T92+T486+(T487+T488) #*T106+(T491+T492)*T70-360*T24+120*T26 T496 = 2880*T17-270*T7+360*T10-1500*T13+225*T14-5580*T16-1200*T72+ #60*T73+600*T90-1800*T75-1080*T76+600*T92-60*T281 T497 = -720*T24+180*T26+60*T36+2430*T40-75*T42-270*T45+90*T46-240* #T56+3375*T57-900*T59-180*T82+225*T83+2800*T84 T500 = -550*T3+120*T-400*T15-1980*T11+260+425*T4+2520*T31+150*T51- #75*T102+60*T32-2040*T41-225*R2-270*T37+1800*T61+60*T1 T502 = 540*R2-72*T28-1120*T4-1440*T31-30*T1-1200*T+630*T5-540*T52+ #144*T58+1890*T11-300*T32+150*T12 T503 = 270*T37+5280*T48-1800*T49-4230*T62+360*T20-125+90*T51+135*T #102+2560*T15-7920*T61+198*T3+6684*T41 T508 = -14280*T17+900*T7+522*T8+72*T10+4980*T13-1800*T14+11160*T16 #-420*T73-15*T19+1575*T75+2430*T76+10*T277+342*T21-60*T9+9342*T24-4 #80*T26-900*T30 T509 = 162*T33+60*T36-6480*T40+105*T42-108*T44+162*T45-108*T46+162 #*T47-10*T50+45*T6+5760*T56-4950*T57-3168*T59+60*T82-75*T83-1620*T6 #0-3700*T84+225*T63 T512 = -150*T18-5880*T8-600*T20+834*T4+10320*T21+720*T-235*R2+720* #T11+630*T9-1400*T5-3600*T15+276*T3+180*T32-5000*T30-1320*T31+270*T #37 T513 = -5886*T41-17664*T48+8240*T49-510*T51+540*T52-414*T54+72*T55 #-1050*T12+54*T58+8370*T61+11352*T62+111*T28+600*T63-36*T35+540*T65 #-60*T102 T516 = -90*T47-200*T305-900*T57+120*T82-150*T83+10*T136-900*T84+15 #0*T319+60*T320+(T496+T497)*R1+T500*T95+(T502+T503)*T124+(480*T-135 #-30*T1+360*T11-140*R2+60*T32+195*T3-660*T31)*T130+(T508+T509)*T69+ #(T512+T513)*T140 T520 = 1/(-1+R1+R2) T522 = T367*T520*T170 T526 = 21080*T17-1470*T7+1224*T8-4860*T10-6180*T13+4770*T14-7700*T #16-24*T19-7080*T21+21*T22+24*T9-17952*T24+225*T25-180*T26-66*T5-45 #*T29+6560*T30+216*T33-60*T34+240*T36 T527 = 6270*T40+30*T42-90*T44-252*T45+72*T47+798*T48-440*T49+60*T5 #0-140*T6+72*T54-72*T55-19140*T56+75*T12+2625*T57+19062*T59+4080*T6 #0-108*T62-1500*T63+75*T64-60*T65 T530 = -720*T24-180*T26-60*T36+2430*T40+75*T42+270*T45-90*T46-240* #T56+3375*T57-900*T59+180*T82-225*T83+2800*T84 T533 = -14280*T17+900*T7+522*T8+72*T10+4980*T13-1800*T14+11160*T16 #-420*T73+15*T19+1575*T75+2430*T76+10*T277+342*T21-60*T9+9342*T24+4 #80*T26-900*T30 T534 = -162*T33-60*T36-6480*T40-105*T42+108*T44-162*T45+108*T46-16 #2*T47-10*T50+45*T6+5760*T56-4950*T57-3168*T59-60*T82+75*T83-1620*T #60-3700*T84+225*T63 T537 = -300*T18+12816*T8-5460*T10+1575*T14-120*T20-124*T4-27000*T2 #1+18*T22-990*T11-1120*T9+425*T25+942*T5+1120*T15-119*T3+16200*T30+ #900*T31-36*T33-300*T34+270*T37 T538 = 1176*T41+180*T44+18018*T48-10080*T49-270*T51-270*T52-378*T5 #4-18*T55-9100*T56+1725*T12+324*T58+13230*T59-1800*T61-9702*T62+39* #T28-2700*T63-51*T35+150*T64+780*T65 T542 = -39*T22+576*T9+11130*T24-550*T25-96*T5+320*T15+45*T29-16860 #*T30-18*T33+360*T34-120*T36 T545 = 22720*T56-900*T12-126*T58-27864*T59-3420*T60-450*T61+2688*T #62+3375*T63+15*T35-225*T64-180*T65 T549 = 420*T17+75*T13+720*T16+200*T72-60*T90+675*T75+270*T76-150*T #92+(480*T-135+30*T1+360*T11-140*R2-60*T32+195*T3-660*T31)*T130+(T5 #26+T527)*T70+(T496+T530)*R1+(T533+T534)*T69+(T537+T538)*T106+(-105 #00*T17+150*T18-8802*T8+9348*T10+2625*T13-4770*T14+9*T19-120*T20-15 #*T4+23418*T21+T542+66*T41-198*T44+144*T45-18*T46-6432*T48+4080*T49 #+195*T6+270*T52-108*T54+162*T55+T545)*T94 T550 = 150*T18-5880*T8+600*T20+834*T4+10320*T21+720*T-235*R2+720*T #11+630*T9-1400*T5-3600*T15+276*T3-180*T32-5000*T30-1320*T31-270*T3 #7 T551 = -5886*T41-17664*T48+8240*T49+510*T51-540*T52+414*T54-72*T55 #-1050*T12-54*T58+8370*T61+11352*T62-111*T28+600*T63+36*T35-540*T65 #+60*T102 T554 = -550*T3+120*T-400*T15-1980*T11+260+425*T4+2520*T31-150*T51+ #75*T102-60*T32-2040*T41-225*R2+270*T37+1800*T61-60*T1 T556 = 540*R2+72*T28-1120*T4-1440*T31+30*T1-1200*T+630*T5+540*T52- #144*T58+1890*T11+300*T32+150*T12 T557 = -270*T37+5280*T48-1800*T49-4230*T62-360*T20-125-90*T51-135* #T102+2560*T15-7920*T61+198*T3+6684*T41 T560 = (T550+T551)*T140+T554*T95+T486+(T556+T557)*T124-360*T24-120 #*T26+90*T47-200*T305-900*T57-120*T82+150*T83+10*T136-900*T84+150*T #319+60*T320 T570 = 54*T11-936*T61-360*T49+80*T15-54*R2+30*T12-224*T4+11-846*T6 #2-396*T31+126*T5+126*T3+120*T+1164*T41+1056*T48 T572 = -1242*T41+25*R2-280*T5+348*T31-1000*T30+210*T4+120*T63+2184 #*T62-180*T11+126*T9-30*T12-72*T-78*T3+1000*T49+2064*T21-3036*T48+1 #026*T61-1176*T8 T574 = 210*T5-144*T31-1080*T49-86*T4+23*T3-1092*T10+315*T14-224*T9 #+2646*T59+85*T25-5184*T21 T575 = -360*T63+2520*T8-15*T12+3150*T48+126*T11-684*T61-1820*T56+6 #24*T41-1962*T62+80*T15+2880*T30 T578 = -2100*T17-1782*T8-2940*T30+126*T9+15*T4-1416*T48+405*T63+18 #48*T10+2226*T24-80*T15+4400*T56 T579 = -684*T60+600*T49-110*T25-138*T41-5400*T59+39*T6+732*T62+234 #*T61+525*T13-66*T5+4446*T21-900*T14 T582 = 4000*T17-294*T7+504*T8-972*T10-1200*T13+900*T14-1540*T16-15 #24*T21-24*T9-3504*T24+45*T25+12*T5 T583 = 1240*T30+1254*T40+246*T48-160*T49+12*T50-28*T6-3540*T56+15* #T12+525*T57+3510*T59+816*T60-108*T62-210*T63 T586 = 2*T277-2*T50+870*T13+198*T21+1080*T56+486*T76-900*T57+6*T9+ #45*T63+315*T75-828*T59+144*T10 T587 = 9*T6-360*T14-324*T60-180*T30-84*T73-90*T8-740*T84+2160*T16+ #180*T7+1674*T24-1296*T40-2280*T17 T590 = -36*T10-12*T281-216*T76+360*T17-120*T56+45*T14-360*T75-900* #T16+144*T59+120*T90-144*T24-240*T72+12*T73+120*T92+405*T57+560*T84 #+486*T40-210*T13-54*T7 T592 = 2*T136+(-24*T-2+12*R2)*T96+(72*T11+24*T+9+39*T3-28*R2-132*T #31)*T130+(85*T4+360*T61-48*T+45*R2-408*T41-80*T15-72*T11-20+324*T3 #1-110*T3)*T95+T570*T124+T572*T140+(T574+T575)*T106+(T578+T579)*T94 #+(T582+T583)*T70+(T586+T587)*T69+T590*R1 T593 = 30*T319-40*T305-12*T90-180*T84+36*T24+15*T13+40*T72+54*T76+ #135*T75+12*T320-30*T92-60*T17 W(1) = -T2*(T91+T133+T135+T151)*T175/180 W(2) = T177*(T180+T193+T195+T234)*T175/180 W(3) = (T261+T270)*T170*T172*T174/9 W(4) = -(2+3*T1)*(T302+T303+T313+T360)*T364*T369/36 W(5) = (3*T1-2)*(T384+T425+T434+T435)*T364*T369/36 W(6) = (T441+T473)*T364*T369/9 W(7) = T2*(T495+T516)*T522/180 W(8) = -T177*(T549+T560)*T522/180 W(9) = -(T592+T593)*T367*T520*T170/9 RETURN END