C 
Cc 
6l 
CONTINUE 
03(1) 
04 (1) 
PIH = 
-OS(I) 
06(T) 
= Z22(1) *VEL 
= E(1)#03(T) 
PIv2, 
= B(I)*(HwW(1)-B(I) #TAs20) 
= 05(1) *E(1)*.5 
CONTINUE 
RHOG=RHU#GRAVTY 
# #® # # ® SET UP THE FUNCTIONS FUR THE INTEGRALS (PAGE 5 UF NOTES)FUNCT 85 
PIH = 
KPI = 
PI7/2, 
KAR#PI 
EVALUATE INTEGRALS USING TRAP METHOD 
91 
93 
9% 
I = 1 
INDEX 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CONTI 
MASS 
QA = 
OMU 
EDMU 
E2DMu 
E30MuU 
BF = 
BMM = 
ZMA = 
ZwWMA 
EMAS 
ZZWMA 
ZWEMA 
Z2WMA 
E2MAZ 
CONTI 
IF ( 
INDEX 
I2zI 
GO TO 
= } 
TRAP (MACINDEX) »DIFF (1) »KTT (1) » TMASS) 
TRAP (QPART (INDEX) sOIFF (I) sKTT (1) 9QAl) 
TRAP (C1 (INDEX) oOIFF (1) oKTT (1) CEA) 
TRAP (C2 (INDEX) »DIFF (1) oKTT (1) »CE2A) 
TRAP (IPART (INDEX) sDIFF (1) sKTT(1I) 9 IAA) 
TRAP (D1 (INDEX) »OIFF (1) 9K TT (1) 9OMUA) 
TRAP (D2 (INDEX) »DOIFF (1) »KTT (1) 9EDMUA) 
TRAP (N3 (INDEX) »DIFF (1) »KTT( 1) 9E2DMUA) 
TRAP (N& (INDEX) oDIFF (1) »9KTT (1) s9E3DMUA) 
TRAP (DS (INDEX) »DIFF (1) oKTT(I) 9 BFA) 
TRAP (N6 (INDEX) oDIFF (1) »KTT (1) »BMMA) 
TRAP (Z1 (INDEX) oDIFF (1) sKTT (1) 9 ZMAA) 
TRAP (72 (INDEX) oDIFF (1) oKTT (1) 9 ZWMAA) 
TRAP (73 (INDEX) »OIFF (1) oKTT(1) 9EMASA) 
TRAP (74( INDEX) »DIFF (1) »KTT(1) 9 ZZWMAA) 
TRAP (75 (INDEX) »DIFF (1) »9KTT (1) 9 ZWEMAA) 
TRAP (76( INDEX) »DIFF (1) »KTT(1) »9Z2WwMAA) 
TRAP (77 (INDEX) »DIFF (I) »KTT (I) s9ECMAZA) 
NUE 
= MASS ¢« TMASS 
QA + QAl 
IA + TAA 
CE + CEA 
cE2 + CE2A 
OMU + DMUA 
= EDM! + EDMUA 
= E20MU «+ E20MUA 
= E30MU «+ E3DMUA 
BF + CHOG*BFA 
BMM + RHOG#®BMMA 
ZMA+ZMAA 
= ZwMA+ZWMAA 
= EMAS+EMASA 
ZZwMA+ZZWMAA 
ZWEMA+Z2WEMAA 
Z2wMA+Z2WMAA 
E2MAZ+E2MAZA 
NUE 
1.GE.11)GO 10 92 
2 INDEX+KTT (I) =} 
ol 
91 
92 CONTINUE 
61 
FUNCT 77 
FUNCT 78 
FUNCT 79 
FUNCT 80 
FUNCT 81 
FUNCT 82 
FUNCT 83 
FUNCT 84 
FUNCT 86 
FUNCT 87 
FUNCT 88 
FUNCT 89 
FUNCT 90 
FUNCT 91 
FUNCT 92 
FUNCT 93 
FUNCT 94 
FUNCT 95 
FUNCT 96 
FUNCT 97 
FUNCT 98 
FUNCT 99 
FUNCT100 
FUNCT101 
FUNCT102 
FUNCT103 
FUNCT104 
FUNCT105 
FUNCT106 
FUNCT107 
FUNCT108 
FUNCT109 
FUNCT110 
FUNCT111 
FUNCT11l2 
FUNCT113 
FUNCT114 
FUNCT115 
FUNCT116 
FUNCT117 
FUNCT118 
FUNCT119 
FUNCT120 
FUNCT121 
FUNCT122 
FUNCT123 
FUNCT124 
FUNCT12S 
FUNCT126 
FUNCT127 
FUNCT128 
FUNCT129 
FUNCT130 
FUNCT131 
FUNCT132 
FUNCT133 
FUNCT134 
FUNCT135 
