Télécharger jelc3d.eso

Retour à la liste

Numérotation des lignes :

jelc3d
  1. C JELC3D SOURCE PV090527 23/01/27 21:15:45 11574
  2. subroutine jelc3d(NBRINC,NINC,NDIMG,NDIMA,FRAC,MK,MG,MbG,JEA,
  3. # affiche,err1)
  4.  
  5. c Jacibienne d epse / d eps imp dans le cas d un seul type d inclusion
  6.  
  7. c organisation de la Jacobienne
  8. c ligne : [eie1, eie2, eie3, ere1, ere2, ere3, eoe1, eoe2, eoe3, eme1, eme2, eme3]
  9. c 1 2 3 4 5 6 7 8 9 10 11 12 = NDIMG
  10. c colonne : [eia1, eia2, eia3, era1, era2, era3, eoa1, eoa2, eoa3, e00a1, e00a2, e00a3, eca1, eca2, eca3, eh1, eh2, eh3, bwiPwi, bgiVgi, bgiVdi, bwmPwm, bgmVgm, bgmVdm]
  11. c 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
  12. implicit real*8 (a-h,o-z)
  13. implicit integer (i-n)
  14.  
  15. logical affiche
  16. integer NBRINC,NDIMG,NDIMA,NINC,err1
  17. c jacobienne depse/depsimp si un seul type dinclusion
  18. real*8 JEA(NDIMG,NDIMA)
  19. real*8 FRAC(0:NBRINC),MK(0:NBRINC),MG(0:NBRINC),MBG(0:NBRINC,0:1)
  20. real*8 Mgm,Mgi
  21.  
  22.  
  23. real*8 f,Km,Gm,Kgm,Ki,Gi,Kgi
  24.  
  25.  
  26. if((NINC.eq.1).and.(NBRINC.eq.1).and.
  27. # (NDIMG.ge.12).and.(NDIMA.ge.24)) then
  28.  
  29. f=FRAC(1)
  30. Km=MK(0)
  31. Gm=MG(0)
  32. Mgm=MbG(0,1)
  33. Ki=MK(1)
  34. Gi=MG(1)
  35. Mgi=MbG(1,1)
  36.  
  37. t1 = Gm ** 2
  38. t3 = 0.24D2 * t1 * Ki
  39. t4 = 0.48D2 * Gi
  40. t5 = 0.12D2 * Ki
  41. t8 = Ki * Gi
  42. t9 = 0.72D2 * t8
  43. t12 = 0.2D1 * Km
  44. t13 = 0.3D1 * Ki
  45. t14 = t12 - t13
  46. t19 = f ** 2
  47. t21 = 0.64D2 * Gi
  48. t22 = 0.16D2 * Ki
  49. t25 = 0.32D2 * Gi
  50. t26 = 0.18D2 * Ki
  51. t31 = Gi * Km
  52. t32 = Ki * t31
  53. t33 = 0.36D2 * t32
  54. t42 = 0.2D1 / 0.3D1 * t1
  55. t44 = -Gi - 0.3D1 / 0.4D1 * Km
  56. t46 = t31 / 0.2D1
  57. t48 = 0.1D1 / ((Gm + Km / 0.2D1) * (-Gm + Gi) * f - t42 + Gm * t44
  58. # - t46)
  59. t49 = t48 * (t19 * (-t3 + Gm * (Km * (-t4 - t5) + t9) - 0.12D2 * G
  60. #i * Km * t14) + f * (t1 * (-t21 - t22) + Gm * (Km * (-t25 - t26) -
  61. # t9) - t33))
  62. t54 = 0.1D1 / ((Ki - Km) * f - 0.4D1 / 0.3D1 * Gm - Ki)
  63. t55 = t54 * t49
  64. t56 = 0.8D1 / 0.9D1 * t55
  65. t57 = 0.24D2 * Gi
  66. t58 = t57 - t5
  67. t59 = Km * t58
  68. t61 = Km ** 2
  69. t62 = Gi * t61
  70. t63 = 0.12D2 * t62
  71. t66 = t25 - t22
  72. t68 = 0.16D2 * Gi
  73. t69 = t68 - t26
  74. t75 = t48 * (t19 * (Gm * t59 - t3 + t63) + f * (Gm * Km * t69 + t1
  75. # * t66))
  76. t76 = t54 * t75
  77. t77 = 0.4D1 / 0.9D1 * t76
  78. t79 = t1 * Gm
  79. t84 = 0.2D1 / 0.9D1 * t55
  80. t85 = -0.16D2 + t84 + t77
  81. t86 = Gi * t85
  82. t87 = 0.5D1 / 0.9D1 * t55
  83. t88 = 0.5D1 / 0.9D1 * t76
  84. t96 = t61 * (t55 / 0.8D1 + t76 / 0.4D1)
  85. t97 = 0.5D1 / 0.18D2 * t55
  86. t99 = Gi * (-0.8D1 + t97 + t88)
  87. t100 = 0.5D1 / 0.12D2 * t55
  88. t101 = 0.5D1 / 0.12D2 * t76
  89. t108 = t55 / 0.72D2
  90. t109 = t76 / 0.36D2
  91. t110 = -t108 - t109
  92. t112 = 0.6D1 * t110 * t62
  93. t116 = 0.1D1 / (0.4D1 * Gm + t13)
  94. t125 = 0.1D1 / (0.12D2 * Gi * Gm + 0.9D1 * Gm * Km + 0.8D1 * t1 +
  95. #0.6D1 * t31)
  96. t126 = t125 * t116 * (t79 * (t56 - t77 - 0.32D2) + t1 * (Km * (-0.
  97. #36D2 + 0.5D1 / 0.6D1 * t55) + t86 + Ki * (-0.16D2 + t87 - t88)) +
  98. #Gm * (t96 + Km * (t99 + Ki * (-0.18D2 + t100 - t101))) - t112)
  99. t127 = 0.2D1 / 0.3D1 * t76
  100. t129 = t79 * (t127 - t84)
  101. t130 = Km * t54
  102. t137 = t1 * (0.5D1 / 0.6D1 * t130 * t75 + t86 + Ki * (0.8D1 + 0.5D
  103. #1 / 0.18D2 * t76 - t97))
  104. t145 = Gm * (t96 + Km * (t99 + Ki * (0.9D1 + 0.5D1 / 0.24D2 * t76
  105. #- 0.5D1 / 0.24D2 * t55)))
  106. t148 = t125 * t116 * (t129 + t137 + t145 - t112)
  107. t149 = 0.16D2 * t79
  108. t150 = 0.24D2 * Km
  109. t153 = t1 * (t150 - t68 - 0.4D1 * Ki)
  110. t154 = 0.18D2 * t61
  111. t155 = 0.3D1 / 0.4D1 * Ki
  112. t157 = Km * (Gi - t155)
  113. t159 = 0.36D2 * t8
  114. t161 = Gm * (t154 + 0.16D2 * t157 - t159)
  115. t162 = Km / 0.3D1
  116. t169 = t19 * f
  117. t171 = 0.12D2 * Km
  118. t186 = t19 * (t149 + t1 * (t171 + t21 + 0.28D2 * Ki) + Gm * (-t154
  119. # + (0.8D1 * Gi + 0.39D2 * Ki) * Km + t9) - 0.12D2 * Gi * Km * (Km
  120. #- t13))
  121. t187 = 0.32D2 * t79
  122. t189 = 0.24D2 * Ki
  123. t199 = f * (-t187 + t1 * (-0.36D2 * Km - t4 - t189) + Gm * (Km * (
  124. #-t57 - 0.27D2 * Ki) - t159) - 0.18D2 * t32)
  125. t201 = t48 * (t169 * (t149 + t153 + t161 - 0.36D2 * Gi * Km * (-t1
  126. #62 + Ki / 0.2D1)) + t186 + t199)
  127. t202 = t54 * t201
  128. t205 = -t68 + 0.20D2 * Ki
  129. t207 = 0.9D1 * t61
  130. t211 = 0.6D1 * t62
  131. t224 = t48 * (t169 * (t149 + t1 * t205 + Gm * (-t207 - 0.20D2 * t1
  132. #57) - t211) + t19 * (-t149 - t1 * t205 + Gm * (t207 + (0.20D2 * Gi
  133. # - 0.15D2 * Ki) * Km) + t211))
  134. t225 = t54 * t224
  135. t226 = 0.4D1 / 0.9D1 * t225
  136. t228 = t79 * (0.8D1 / 0.9D1 * t202 - t226)
  137. t231 = 0.2D1 / 0.9D1 * t202
  138. t232 = t231 + t226
  139. t233 = Gi * t232
  140. t234 = t202 - t225
  141. t237 = t1 * (0.5D1 / 0.6D1 * t130 * t201 + t233 + 0.5D1 / 0.9D1 *
  142. #Ki * t234)
  143. t239 = t225 / 0.4D1
  144. t241 = t61 * (t202 / 0.8D1 + t239)
  145. t243 = 0.5D1 / 0.9D1 * t225
  146. t245 = Gi * (0.5D1 / 0.18D2 * t202 + t243)
  147. t250 = Gm * (t241 + Km * (t245 + 0.5D1 / 0.12D2 * Ki * t234))
  148. t251 = t202 / 0.72D2
  149. t252 = t225 / 0.36D2
  150. t253 = -t251 - t252
  151. t255 = 0.6D1 * t253 * t62
  152. t259 = 0.2D1 / 0.3D1 * t225
  153. t260 = 0.3D1 / 0.2D1 * Ki
  154. t261 = -Km + t260
  155. t267 = t48 * (t169 * (-0.12D2 * t261 * t31 + t149 + t153 + t161) +
  156. # t186 + t199)
  157. t268 = t54 * t267
  158. t269 = 0.2D1 / 0.9D1 * t268
  159. t271 = t79 * (t259 - t269)
  160. t273 = 0.5D1 / 0.6D1 * t130 * t224
  161. t274 = t226 + t269
  162. t275 = Gi * t274
  163. t276 = t225 - t268
  164. t279 = t1 * (t273 + t275 + 0.5D1 / 0.18D2 * Ki * t276)
  165. t282 = t61 * (t239 + t268 / 0.8D1)
  166. t285 = Gi * (t243 + 0.5D1 / 0.18D2 * t268)
  167. t290 = Gm * (t282 + Km * (t285 + 0.5D1 / 0.24D2 * Ki * t276))
  168. t291 = t268 / 0.72D2
  169. t292 = -t252 - t291
  170. t294 = 0.6D1 * t292 * t62
  171. t297 = t125 * t116 * (t271 + t279 + t290 - t294)
  172. t301 = t1 * (0.48D2 * Km - t25 - 0.8D1 * Ki)
  173. t302 = 0.36D2 * t61
  174. t305 = Gm * (t302 + 0.32D2 * t157 - t9)
  175. t310 = 0.36D2 * (-0.2D1 / 0.3D1 * Km + Ki) * Km * Gi
  176. t330 = t19 * (t1 * (-t150 + 0.160D3 * Gi + 0.88D2 * Ki) + Gm * (-0
  177. #.72D2 * t61 + Km * (t25 + 0.114D3 * Ki) + 0.144D3 * t8) - 0.12D2 *
  178. # Gi * Km * (t12 - 0.6D1 * Ki))
  179. t343 = f * (-0.128D3 * t79 + t1 * (-0.144D3 * Km - 0.128D3 * Gi -
  180. #0.80D2 * Ki) + Gm * (Km * (-t21 - 0.90D2 * Ki) - t9) - t33)
  181. t345 = t48 * (t169 * (t187 + t301 + t305 - t310) + t330 + t343)
  182. t346 = t54 * t345
  183. t347 = 0.8D1 / 0.9D1 * t346
  184. t349 = -t25 + 0.40D2 * Ki
  185. t350 = t1 * t349
  186. t356 = 0.64D2 * t79
  187. t361 = 0.48D2 * Ki
  188. t374 = t48 * (t169 * (t187 + t350 + Gm * (-t154 - 0.40D2 * t157) -
  189. # t63) + t19 * (-t356 + t1 * (t21 - 0.56D2 * Ki) + Gm * (t302 + Km
  190. #* (0.56D2 * Gi - t361)) + t63) + f * (-Gm * Km * t69 - t1 * t66))
  191. t375 = t54 * t374
  192. t376 = 0.4D1 / 0.9D1 * t375
  193. t382 = 0.2D1 / 0.9D1 * t346
  194. t383 = t382 + t376 + 0.16D2
  195. t384 = Gi * t383
  196. t385 = 0.5D1 / 0.9D1 * t346
  197. t386 = 0.5D1 / 0.9D1 * t375
  198. t392 = t375 / 0.4D1
  199. t394 = t61 * (t346 / 0.8D1 + t392)
  200. t395 = 0.5D1 / 0.18D2 * t346
  201. t397 = Gi * (t395 + t386 + 0.8D1)
  202. t398 = 0.5D1 / 0.12D2 * t346
  203. t399 = 0.5D1 / 0.12D2 * t375
  204. t406 = t346 / 0.72D2
  205. t407 = t375 / 0.36D2
  206. t408 = -t406 - t407
  207. t410 = 0.6D1 * t408 * t62
  208. t414 = 0.2D1 / 0.3D1 * t375
  209. t416 = -0.12D2 * t14 * t31
  210. t420 = t48 * (t169 * (t187 + t301 + t305 - t416) + t330 + t343)
  211. t421 = t54 * t420
  212. t422 = 0.2D1 / 0.9D1 * t421
  213. t424 = t79 * (t414 - t422)
  214. t426 = 0.5D1 / 0.6D1 * t130 * t374
  215. t427 = t376 + t422 + 0.16D2
  216. t428 = Gi * t427
  217. t429 = 0.5D1 / 0.18D2 * t375
  218. t430 = 0.5D1 / 0.18D2 * t421
  219. t434 = t1 * (t426 + t428 + Ki * (-0.8D1 + t429 - t430))
  220. t437 = t61 * (t392 + t421 / 0.8D1)
  221. t439 = Gi * (t386 + t430 + 0.8D1)
  222. t440 = 0.5D1 / 0.24D2 * t375
  223. t447 = Gm * (t437 + Km * (t439 + Ki * (-0.9D1 + t440 - 0.5D1 / 0.2
  224. #4D2 * t421)))
  225. t448 = t421 / 0.72D2
  226. t449 = -t407 - t448
  227. t451 = 0.6D1 * t449 * t62
  228. t454 = t125 * t116 * (t424 + t434 + t447 - t451)
  229. t455 = 0.48D2 * t79
  230. t456 = 0.72D2 * Km
  231. t458 = t1 * (-t456 + t4 + t5)
  232. t461 = 0.108D3 * t8
  233. t463 = Gm * (-0.54D2 * t61 - 0.48D2 * t157 + t461)
  234. t487 = t19 * (-t149 + t1 * (0.84D2 * Km - 0.224D3 * Gi - 0.164D3 *
  235. # Ki) + Gm * (0.126D3 * t61 + Km * (-t21 - 0.177D3 * Ki) - 0.216D3
  236. #* t8) - 0.12D2 * Gi * Km * (-t12 + 0.9D1 * Ki))
  237. t488 = 0.224D3 * t79
  238. t493 = t1 * (0.228D3 * Km + 0.176D3 * Gi + 0.152D3 * Ki)
  239. t499 = Gm * (t154 + (0.112D3 * Gi + 0.141D3 * Ki) * Km + t461)
  240. t500 = 0.9D1 / 0.2D1 * Ki
  241. t505 = f * (t488 + t493 + t499 - 0.12D2 * (-Km - t500) * t31)
  242. t508 = t54 * t48 * (t169 * (0.36D2 * Gi * Km * t261 - t455 + t458
  243. #+ t463) + t487 + t505)
  244. t538 = 0.40D2 * Gi - 0.30D2 * Ki
  245. t545 = t48 * (t169 * (-t455 + t1 * (t4 - 0.60D2 * Ki) + Gm * (0.27
  246. #D2 * t61 + 0.60D2 * t157) + 0.18D2 * t62) + t19 * (0.80D2 * t79 +
  247. #t1 * (-0.80D2 * Gi + 0.100D3 * Ki) + Gm * (-0.45D2 * t61 + (-0.100
  248. #D3 * Gi + 0.75D2 * Ki) * Km) - 0.30D2 * t62) + f * (-t187 - t1 * t
  249. #349 + Gm * (Km * t538 + t154) + t63))
  250. t546 = t54 * t545
  251. t547 = 0.4D1 / 0.9D1 * t546
  252. t549 = t79 * (0.8D1 / 0.9D1 * t508 - t547 - 0.64D2)
  253. t553 = 0.2D1 / 0.9D1 * t508
  254. t554 = t553 - 0.16D2 + t547
  255. t555 = Gi * t554
  256. t557 = 0.5D1 / 0.9D1 * t546
  257. t561 = t1 * (Km * (-0.60D2 + 0.5D1 / 0.6D1 * t508) + t555 + Ki * (
  258. #-0.40D2 + 0.5D1 / 0.9D1 * t508 - t557))
  259. t563 = t546 / 0.4D1
  260. t565 = t61 * (t508 / 0.8D1 - 0.9D1 + t563)
  261. t566 = 0.5D1 / 0.18D2 * t508
  262. t568 = Gi * (t566 - 0.20D2 + t557)
  263. t570 = 0.5D1 / 0.12D2 * t546
  264. t576 = Gm * (t565 + Km * (t568 + Ki * (-0.30D2 + 0.5D1 / 0.12D2 *
  265. #t508 - t570)))
  266. t577 = t508 / 0.72D2
  267. t578 = t546 / 0.36D2
  268. t579 = 0.1D1 - t577 - t578
  269. t581 = 0.6D1 * t579 * t62
  270. t585 = 0.2D1 / 0.3D1 * t546
  271. t591 = t169 * (-t455 + t458 + t463 - 0.12D2 * (0.3D1 * Km - t500)
  272. #* t31)
  273. t594 = t54 * t48 * (t591 + t487 + t505)
  274. t595 = 0.2D1 / 0.9D1 * t594
  275. t597 = t79 * (t585 - t595 + 0.16D2)
  276. t599 = 0.5D1 / 0.6D1 * t130 * t545
  277. t600 = t547 + t595 - 0.16D2
  278. t601 = Gi * t600
  279. t602 = 0.5D1 / 0.18D2 * t546
  280. t603 = 0.5D1 / 0.18D2 * t594
  281. t607 = t1 * (t599 + t601 + Ki * (0.20D2 + t602 - t603))
  282. t610 = t61 * (t563 + t594 / 0.8D1 - 0.9D1)
  283. t612 = Gi * (t557 + t603 - 0.20D2)
  284. t613 = 0.5D1 / 0.24D2 * t546
  285. t620 = Gm * (t610 + Km * (t612 + Ki * (0.15D2 + t613 - 0.5D1 / 0.2
  286. #4D2 * t594)))
  287. t621 = t594 / 0.72D2
  288. t622 = 0.1D1 - t578 - t621
  289. t624 = 0.6D1 * t622 * t62
  290. t627 = t125 * t116 * (t597 + t607 + t620 - t624)
  291. t636 = t54 * t48 * (t591 + t487 + f * (t488 + t493 + t499 - 0.36D2
  292. # * Gi * Km * (-t162 - t260)))
  293. t637 = 0.2D1 / 0.9D1 * t636
  294. t639 = t79 * (t585 - t637 + 0.16D2)
  295. t640 = t547 + t637 - 0.16D2
  296. t641 = Gi * t640
  297. t642 = 0.5D1 / 0.18D2 * t636
  298. t646 = t1 * (t599 + t641 + Ki * (0.20D2 + t602 - t642))
  299. t649 = t61 * (t563 + t636 / 0.8D1 - 0.9D1)
  300. t651 = Gi * (t557 + t642 - 0.20D2)
  301. t658 = Gm * (t649 + Km * (t651 + Ki * (0.15D2 + t613 - 0.5D1 / 0.2
  302. #4D2 * t636)))
  303. t659 = t636 / 0.72D2
  304. t660 = 0.1D1 - t578 - t659
  305. t662 = 0.6D1 * t660 * t62
  306. t665 = t125 * t116 * (t639 + t646 + t658 - t662)
  307. t668 = -Gi - Ki
  308. t683 = t54 * t48 * (t19 * (t1 * (t456 - t361) + Gm * (0.24D2 * Km
  309. #* t668 + t302) - t63) + f * (t356 + 0.48D2 * t1 * (Km + Ki) + Gm *
  310. # (-0.24D2 * Km * t668 + t154) + t63))
  311. t684 = 0.8D1 / 0.9D1 * t683
  312. t694 = t48 * (t19 * (-Gm * Km * t58 + t3 - t63) + f * (-t187 - t3
  313. #+ Gm * (t154 + t59) + t63))
  314. t695 = t54 * t694
  315. t696 = 0.4D1 / 0.9D1 * t695
  316. t699 = 0.5D1 / 0.6D1 * t683
  317. t702 = 0.2D1 / 0.9D1 * t683
  318. t703 = t702 + t696
  319. t704 = Gi * t703
  320. t705 = 0.5D1 / 0.9D1 * t683
  321. t706 = 0.5D1 / 0.9D1 * t695
  322. t714 = t61 * (t683 / 0.8D1 - 0.9D1 + t695 / 0.4D1)
  323. t715 = 0.5D1 / 0.18D2 * t683
  324. t717 = Gi * (t715 - 0.12D2 + t706)
  325. t718 = 0.5D1 / 0.12D2 * t683
  326. t719 = 0.5D1 / 0.12D2 * t695
  327. t726 = t683 / 0.72D2
  328. t727 = t695 / 0.36D2
  329. t728 = 0.1D1 - t726 - t727
  330. t730 = 0.6D1 * t728 * t62
  331. t733 = t125 * t116 * (t79 * (t684 - t696 - 0.32D2) + t1 * (Km * (-
  332. #0.24D2 + t699) + t704 + Ki * (-0.24D2 + t705 - t706)) + Gm * (t714
  333. # + Km * (t717 + Ki * (-0.12D2 + t718 - t719))) - t730)
  334. t734 = 0.2D1 / 0.3D1 * t695
  335. t736 = t79 * (t734 - t702 + 0.16D2)
  336. t743 = t1 * (0.5D1 / 0.6D1 * t130 * t694 + t704 + Ki * (0.12D2 + 0
  337. #.5D1 / 0.18D2 * t695 - t715))
  338. t751 = Gm * (t714 + Km * (t717 + Ki * (0.6D1 + 0.5D1 / 0.24D2 * t6
  339. #95 - 0.5D1 / 0.24D2 * t683)))
  340. t754 = t125 * t116 * (t736 + t743 + t751 - t730)
  341. t758 = Gm * (t302 + Km * (t25 - t189) - t9)
  342. t765 = 0.96D2 * (-Gm * t44 + t42 + t46) * (t155 + Gm)
  343. t767 = t48 * (f * (t187 + t301 + t758 - t416) + t765)
  344. t768 = t54 * t767
  345. t774 = f * (t187 + t350 + Gm * (-Km * t538 - t154) - t63)
  346. t775 = t54 * t48
  347. t776 = t775 * t774
  348. t777 = 0.4D1 / 0.9D1 * t776
  349. t779 = t79 * (0.8D1 / 0.9D1 * t768 - t777)
  350. t782 = 0.2D1 / 0.9D1 * t768
  351. t783 = t782 + t777
  352. t784 = Gi * t783
  353. t785 = t768 - t776
  354. t788 = t1 * (0.5D1 / 0.6D1 * t130 * t767 + t784 + 0.5D1 / 0.9D1 *
  355. #Ki * t785)
  356. t790 = t776 / 0.4D1
  357. t792 = t61 * (t768 / 0.8D1 + t790)
  358. t794 = 0.5D1 / 0.9D1 * t776
  359. t796 = Gi * (0.5D1 / 0.18D2 * t768 + t794)
  360. t801 = Gm * (t792 + Km * (t796 + 0.5D1 / 0.12D2 * Ki * t785))
  361. t802 = t768 / 0.72D2
  362. t803 = t776 / 0.36D2
  363. t804 = -t802 - t803
  364. t806 = 0.6D1 * t804 * t62
  365. t809 = t125 * t116 * (t779 + t788 + t801 - t806)
  366. t810 = 0.2D1 / 0.3D1 * t776
  367. t812 = t79 * (t810 - t782)
  368. t815 = 0.5D1 / 0.6D1 * Km * t775 * t774
  369. t818 = t1 * (t815 + t784 - 0.5D1 / 0.18D2 * Ki * t785)
  370. t823 = Gm * (t792 + Km * (t796 - 0.5D1 / 0.24D2 * Ki * t785))
  371. t826 = t125 * t116 * (t812 + t818 + t823 - t806)
  372. t830 = t48 * (f * (t187 + t301 + t758 - t310) + t765)
  373. t831 = t54 * t830
  374. t832 = 0.2D1 / 0.9D1 * t831
  375. t834 = t79 * (t810 - t832)
  376. t835 = t777 + t832
  377. t836 = Gi * t835
  378. t837 = t776 - t831
  379. t840 = t1 * (t815 + t836 + 0.5D1 / 0.18D2 * Ki * t837)
  380. t843 = t61 * (t790 + t831 / 0.8D1)
  381. t846 = Gi * (t794 + 0.5D1 / 0.18D2 * t831)
  382. t851 = Gm * (t843 + Km * (t846 + 0.5D1 / 0.24D2 * Ki * t837))
  383. t852 = t831 / 0.72D2
  384. t853 = -t803 - t852
  385. t855 = 0.6D1 * t853 * t62
  386. t858 = t125 * t116 * (t834 + t840 + t851 - t855)
  387. t859 = 0.24D2 * t1
  388. t860 = -t171 + t57
  389. t862 = 0.12D2 * t31
  390. t865 = 0.16D2 * t1
  391. t867 = -0.18D2 * Km - t57
  392. t872 = t48 * (t19 * (Gm * t860 - t859 + t862) + f * (Gm * t867 - t
  393. #862 - t865))
  394. t873 = t79 * t54
  395. t875 = 0.4D1 / 0.9D1 * t873 * t872
  396. t876 = t130 * t872
  397. t878 = Gi * t54
  398. t879 = t878 * t872
  399. t880 = 0.2D1 / 0.3D1 * t879
  400. t882 = t1 * (0.5D1 / 0.6D1 * t876 + t880 + 0.8D1)
  401. t883 = t61 * t54
  402. t884 = t883 * t872
  403. t889 = 0.12D2 * Gi
  404. t891 = Gm * (0.3D1 / 0.8D1 * t884 + Km * (0.5D1 / 0.6D1 * t879 + 0
  405. #.9D1) + t889)
  406. t892 = t876 / 0.24D2
  407. t898 = t125 * t116 * (t875 + t882 + t891 - 0.6D1 * (-t892 - 0.1D1)
  408. # * t31)
  409. t899 = t1 * Mgi
  410. t900 = 0.24D2 * t899
  411. t901 = Mgi * Gi
  412. t902 = 0.24D2 * t901
  413. t903 = Km * Mgi
  414. t905 = t902 - 0.12D2 * t903
  415. t908 = 0.12D2 * Gi * t903
  416. t911 = 0.16D2 * t899
  417. t913 = -t902 - 0.18D2 * t903
  418. t918 = t48 * (t19 * (Gm * t905 - t900 + t908) + f * (Gm * t913 - t
  419. #908 - t911))
  420. t920 = 0.4D1 / 0.9D1 * t873 * t918
  421. t921 = t130 * t918
  422. t923 = t878 * t918
  423. t924 = 0.2D1 / 0.3D1 * t923
  424. t925 = 0.8D1 * Mgi
  425. t927 = t1 * (0.5D1 / 0.6D1 * t921 + t924 + t925)
  426. t928 = t883 * t918
  427. t931 = 0.9D1 * Mgi
  428. t934 = 0.12D2 * t901
  429. t936 = Gm * (0.3D1 / 0.8D1 * t928 + Km * (0.5D1 / 0.6D1 * t923 + t
  430. #931) + t934)
  431. t937 = t921 / 0.24D2
  432. t943 = t125 * t116 * (t920 + t927 + t936 - 0.6D1 * (-t937 - Mgi) *
  433. # t31)
  434. t951 = t48 * (t19 * (-Gm * t905 + t900 - t908) + f * (-Gm * t913 +
  435. # t908 + t911))
  436. t953 = 0.4D1 / 0.9D1 * t873 * t951
  437. t954 = t130 * t951
  438. t956 = t878 * t951
  439. t957 = 0.2D1 / 0.3D1 * t956
  440. t959 = t1 * (0.5D1 / 0.6D1 * t954 + t957 - t925)
  441. t960 = t883 * t951
  442. t966 = Gm * (0.3D1 / 0.8D1 * t960 + Km * (0.5D1 / 0.6D1 * t956 - t
  443. #931) - t934)
  444. t967 = t954 / 0.24D2
  445. t973 = t125 * t116 * (t953 + t959 + t966 - 0.6D1 * (-t967 + Mgi) *
  446. # t31)
  447. t981 = t48 * (t19 * (-Gm * t860 + t859 - t862) + f * (-Gm * t867 +
  448. # t862 + t865))
  449. t983 = 0.4D1 / 0.9D1 * t873 * t981
  450. t984 = t130 * t981
  451. t986 = t878 * t981
  452. t987 = 0.2D1 / 0.3D1 * t986
  453. t989 = t1 * (0.5D1 / 0.6D1 * t984 + t987 - 0.8D1)
  454. t990 = t883 * t981
  455. t996 = Gm * (0.3D1 / 0.8D1 * t990 + Km * (0.5D1 / 0.6D1 * t986 - 0
  456. #.9D1) - t889)
  457. t997 = t984 / 0.24D2
  458. t1003 = t125 * t116 * (t983 + t989 + t996 - 0.6D1 * (-t997 + 0.1D1
  459. #) * t31)
  460. t1004 = t1 * Mgm
  461. t1005 = 0.24D2 * t1004
  462. t1006 = Mgm * Gi
  463. t1007 = 0.24D2 * t1006
  464. t1008 = Km * Mgm
  465. t1010 = -t1007 + 0.12D2 * t1008
  466. t1013 = 0.12D2 * Gi * t1008
  467. t1016 = 0.16D2 * t1004
  468. t1018 = t1007 + 0.18D2 * t1008
  469. t1023 = t48 * (t19 * (Gm * t1010 + t1005 - t1013) + f * (Gm * t101
  470. #8 + t1013 + t1016))
  471. t1025 = 0.4D1 / 0.9D1 * t873 * t1023
  472. t1026 = t130 * t1023
  473. t1028 = t878 * t1023
  474. t1029 = 0.2D1 / 0.3D1 * t1028
  475. t1030 = 0.8D1 * Mgm
  476. t1032 = t1 * (0.5D1 / 0.6D1 * t1026 + t1029 - t1030)
  477. t1033 = t883 * t1023
  478. t1036 = 0.9D1 * Mgm
  479. t1039 = 0.12D2 * t1006
  480. t1041 = Gm * (0.3D1 / 0.8D1 * t1033 + Km * (0.5D1 / 0.6D1 * t1028
  481. #- t1036) - t1039)
  482. t1042 = t1026 / 0.24D2
  483. t1048 = t125 * t116 * (t1025 + t1032 + t1041 - 0.6D1 * (-t1042 + M
  484. #gm) * t31)
  485. t1056 = t48 * (t19 * (-Gm * t1010 - t1005 + t1013) + f * (-Gm * t1
  486. #018 - t1013 - t1016))
  487. t1058 = 0.4D1 / 0.9D1 * t873 * t1056
  488. t1059 = t130 * t1056
  489. t1061 = t878 * t1056
  490. t1062 = 0.2D1 / 0.3D1 * t1061
  491. t1064 = t1 * (0.5D1 / 0.6D1 * t1059 + t1062 + t1030)
  492. t1065 = t883 * t1056
  493. t1071 = Gm * (0.3D1 / 0.8D1 * t1065 + Km * (0.5D1 / 0.6D1 * t1061
  494. #+ t1036) + t1039)
  495. t1072 = t1059 / 0.24D2
  496. t1078 = t125 * t116 * (t1058 + t1064 + t1071 - 0.6D1 * (-t1072 - M
  497. #gm) * t31)
  498. t1080 = t79 * (-t231 + t259)
  499. t1091 = t125 * t116 * (t1080 + t1 * (t273 - 0.5D1 / 0.18D2 * Ki *
  500. #t234 + t233) + Gm * (t241 + Km * (-0.5D1 / 0.24D2 * Ki * t234 + t2
  501. #45)) - t255)
  502. t1107 = t125 * t116 * (t79 * (-t226 + 0.8D1 / 0.9D1 * t268) + t1 *
  503. # (0.5D1 / 0.6D1 * t130 * t267 - 0.5D1 / 0.9D1 * Ki * t276 + t275)
  504. #+ Gm * (t282 + Km * (-0.5D1 / 0.12D2 * Ki * t276 + t285)) - t294)
  505. t1109 = t79 * (-t382 + t414)
  506. t1123 = t125 * t116 * (t1109 + t1 * (t426 + Ki * (-t395 + t429 - 0
  507. #.8D1) + t384) + Gm * (t394 + Km * (Ki * (-0.5D1 / 0.24D2 * t346 +
  508. #t440 - 0.9D1) + t397)) - t410)
  509. t1124 = 0.8D1 / 0.9D1 * t421
  510. t1130 = 0.5D1 / 0.9D1 * t421
  511. t1135 = 0.5D1 / 0.12D2 * t421
  512. t1144 = t125 * t116 * (t79 * (-t376 + t1124 + 0.32D2) + t1 * (Km *
  513. # (0.5D1 / 0.6D1 * t421 + 0.36D2) + Ki * (-t386 + t1130 + 0.16D2) +
  514. # t428) + Gm * (t437 + Km * (Ki * (-t399 + t1135 + 0.18D2) + t439))
  515. # - t451)
  516. t1146 = t79 * (-t553 + 0.16D2 + t585)
  517. t1160 = t125 * t116 * (t1146 + t1 * (t599 + Ki * (-t566 + 0.20D2 +
  518. # t602) + t555) + Gm * (t565 + Km * (Ki * (-0.5D1 / 0.24D2 * t508 +
  519. # 0.15D2 + t613) + t568)) - t581)
  520. t1181 = t125 * t116 * (t79 * (-t547 + 0.8D1 / 0.9D1 * t594 - 0.64D
  521. #2) + t1 * (Km * (0.5D1 / 0.6D1 * t594 - 0.60D2) + Ki * (-t557 + 0.
  522. #5D1 / 0.9D1 * t594 - 0.40D2) + t601) + Gm * (t610 + Km * (Ki * (-t
  523. #570 + 0.5D1 / 0.12D2 * t594 - 0.30D2) + t612)) - t624)
  524. t1202 = t125 * t116 * (t79 * (-t547 + 0.8D1 / 0.9D1 * t636 - 0.64D
  525. #2) + t1 * (Km * (-0.60D2 + 0.5D1 / 0.6D1 * t636) + t641 + Ki * (-0
  526. #.40D2 - t557 + 0.5D1 / 0.9D1 * t636)) + Gm * (t649 + Km * (t651 +
  527. #Ki * (-0.30D2 - t570 + 0.5D1 / 0.12D2 * t636))) - t662)
  528. t1218 = t125 * t116 * (t79 * (-t777 + 0.8D1 / 0.9D1 * t831) + t1 *
  529. # (0.5D1 / 0.6D1 * t130 * t830 + t836 - 0.5D1 / 0.9D1 * Ki * t837)
  530. #+ Gm * (t843 + Km * (t846 - 0.5D1 / 0.12D2 * Ki * t837)) - t855)
  531. t1230 = t1 * (Gi * (-0.48D2 + 0.4D1 / 0.3D1 * t55 - t127) + Km * (
  532. #t55 / 0.2D1 - t127) + Ki * (-0.16D2 + t55 / 0.18D2 + 0.17D2 / 0.18
  533. #D2 * t76))
  534. t1234 = t61 * (-t55 / 0.4D1 - t76 / 0.2D1)
  535. t1249 = Gm * (t1234 + Km * (-Gi * t85 + Ki * (-0.18D2 + 0.19D2 / 0
  536. #.24D2 * t55 + t76 / 0.3D1)) - 0.36D2 * Gi * Ki * (0.5D1 / 0.3D1 -
  537. #t55 / 0.27D2 - t76 / 0.216D3))
  538. t1253 = Km * (t55 / 0.108D3 + t76 / 0.54D2)
  539. t1262 = t55 / 0.3D1
  540. t1264 = t79 * (t76 / 0.9D1 + t1262)
  541. t1275 = t1 * (Gi * (0.24D2 + t76 - t1262) + Km * (t76 / 0.6D1 - t1
  542. #262) + Ki * (-0.16D2 + 0.19D2 / 0.36D2 * t76 + 0.17D2 / 0.36D2 * t
  543. #55))
  544. t1291 = Gm * (t1234 + Km * (Gi * (-0.8D1 - t77 - t84) + Ki * (-0.1
  545. #8D2 + 0.23D2 / 0.24D2 * t76 + t55 / 0.6D1)) - 0.36D2 * Gi * Ki * (
  546. #0.1D1 / 0.6D1 - 0.17D2 / 0.432D3 * t76 - t55 / 0.432D3))
  547. t1299 = t125 * t116 * (t1264 + t1275 + t1291 - 0.18D2 * (t1253 + K
  548. #i * (0.1D1 - t109 - t108)) * t31)
  549. t1307 = 0.17D2 / 0.18D2 * t225
  550. t1313 = t225 / 0.2D1
  551. t1315 = t61 * (-t202 / 0.4D1 - t1313)
  552. t1316 = -Gi * t232
  553. t1318 = t225 / 0.3D1
  554. t1324 = t225 / 0.216D3
  555. t1332 = t225 / 0.54D2
  556. t1342 = t225 / 0.9D1
  557. t1343 = t268 / 0.3D1
  558. t1345 = t79 * (t1342 + t1343)
  559. t1348 = t225 / 0.6D1
  560. t1351 = 0.19D2 / 0.36D2 * t225
  561. t1356 = t1 * (Gi * (t225 - t1343) + Km * (t1348 - t1343) + Ki * (t
  562. #1351 + 0.17D2 / 0.36D2 * t268))
  563. t1359 = t61 * (-t1313 - t268 / 0.4D1)
  564. t1360 = -Gi * t274
  565. t1361 = 0.23D2 / 0.24D2 * t225
  566. t1367 = 0.17D2 / 0.432D3 * t225
  567. t1374 = Gm * (t1359 + Km * (t1360 + Ki * (t1361 + t268 / 0.6D1)) -
  568. # 0.36D2 * Gi * Ki * (-t1367 - t268 / 0.432D3))
  569. t1384 = t125 * t116 * (t1345 + t1356 + t1374 - 0.18D2 * (Km * (t13
  570. #32 + t268 / 0.108D3) + Ki * t292) * t31)
  571. t1392 = 0.17D2 / 0.18D2 * t375
  572. t1398 = t375 / 0.2D1
  573. t1400 = t61 * (-t346 / 0.4D1 - t1398)
  574. t1403 = t375 / 0.3D1
  575. t1409 = t375 / 0.216D3
  576. t1417 = t375 / 0.54D2
  577. t1427 = t375 / 0.9D1
  578. t1428 = t421 / 0.3D1
  579. t1430 = t79 * (t1427 + t1428)
  580. t1433 = t375 / 0.6D1
  581. t1436 = 0.19D2 / 0.36D2 * t375
  582. t1441 = t1 * (Gi * (t375 - t1428 - 0.24D2) + Km * (t1433 - t1428)
  583. #+ Ki * (t1436 + 0.17D2 / 0.36D2 * t421 + 0.16D2))
  584. t1444 = t61 * (-t1398 - t421 / 0.4D1)
  585. t1447 = 0.23D2 / 0.24D2 * t375
  586. t1453 = 0.17D2 / 0.432D3 * t375
  587. t1460 = Gm * (t1444 + Km * (Gi * (-t376 - t422 + 0.8D1) + Ki * (0.
  588. #18D2 + t1447 + t421 / 0.6D1)) - 0.36D2 * Gi * Ki * (-0.1D1 / 0.6D1
  589. # - t1453 - t421 / 0.432D3))
  590. t1471 = t125 * t116 * (t1430 + t1441 + t1460 - 0.18D2 * (Km * (t14
  591. #17 + t421 / 0.108D3) + Ki * (-t407 - t448 - 0.1D1)) * t31)
  592. t1479 = 0.17D2 / 0.18D2 * t546
  593. t1485 = t546 / 0.2D1
  594. t1487 = t61 * (0.18D2 - t508 / 0.4D1 - t1485)
  595. t1488 = -Gi * t554
  596. t1490 = t546 / 0.3D1
  597. t1496 = t546 / 0.216D3
  598. t1504 = t546 / 0.54D2
  599. t1514 = t546 / 0.9D1
  600. t1515 = t594 / 0.3D1
  601. t1517 = t79 * (t1514 + t1515 - 0.24D2)
  602. t1520 = t546 / 0.6D1
  603. t1523 = 0.19D2 / 0.36D2 * t546
  604. t1528 = t1 * (Gi * (t546 - t1515 + 0.24D2) + Km * (t1520 - t1515 +
  605. # 0.24D2) + Ki * (t1523 + 0.17D2 / 0.36D2 * t594 - 0.34D2))
  606. t1531 = t61 * (0.18D2 - t1485 - t594 / 0.4D1)
  607. t1532 = -Gi * t600
  608. t1533 = 0.23D2 / 0.24D2 * t546
  609. t1539 = 0.17D2 / 0.432D3 * t546
  610. t1546 = Gm * (t1531 + Km * (t1532 + Ki * (-0.12D2 + t1533 + t594 /
  611. # 0.6D1)) - 0.36D2 * Gi * Ki * (0.1D1 / 0.6D1 - t1539 - t594 / 0.43
  612. #2D3))
  613. t1557 = t636 / 0.3D1
  614. t1559 = t79 * (t1514 + t1557 - 0.24D2)
  615. t1568 = t1 * (Gi * (t546 - t1557 + 0.24D2) + Km * (t1520 - t1557 +
  616. # 0.24D2) + Ki * (t1523 + 0.17D2 / 0.36D2 * t636 - 0.34D2))
  617. t1571 = t61 * (0.18D2 - t1485 - t636 / 0.4D1)
  618. t1572 = -Gi * t640
  619. t1584 = Gm * (t1571 + Km * (t1572 + Ki * (-0.12D2 + t1533 + t636 /
  620. # 0.6D1)) - 0.36D2 * Gi * Ki * (0.1D1 / 0.6D1 - t1539 - t636 / 0.43
  621. #2D3))
  622. t1606 = t1 * (Gi * (0.4D1 / 0.3D1 * t683 - 0.48D2 - t734) + Km * (
  623. #t683 / 0.2D1 - 0.36D2 - t734) + Ki * (t683 / 0.18D2 + 0.12D2 + 0.1
  624. #7D2 / 0.18D2 * t695))
  625. t1610 = t61 * (0.18D2 - t683 / 0.4D1 - t695 / 0.2D1)
  626. t1625 = Gm * (t1610 + Km * (-Gi * t703 + Ki * (-0.39D2 + 0.19D2 /
  627. #0.24D2 * t683 + t695 / 0.3D1)) - 0.36D2 * Gi * Ki * (0.1D1 - t683
  628. #/ 0.27D2 - t695 / 0.216D3))
  629. t1629 = Km * (t683 / 0.108D3 - 0.2D1 / 0.3D1 + t695 / 0.54D2)
  630. t1638 = t683 / 0.3D1
  631. t1640 = t79 * (t695 / 0.9D1 + t1638 - 0.24D2)
  632. t1651 = t1 * (Gi * (t695 - t1638) + Km * (t695 / 0.6D1 - t1638 + 0
  633. #.24D2) + Ki * (0.19D2 / 0.36D2 * t695 + 0.17D2 / 0.36D2 * t683 - 0
  634. #.18D2))
  635. t1667 = Gm * (t1610 + Km * (Gi * (-t696 - t702 + 0.24D2) + Ki * (0
  636. #.6D1 + 0.23D2 / 0.24D2 * t695 + t683 / 0.6D1)) - 0.36D2 * Gi * Ki
  637. #* (-0.17D2 / 0.432D3 * t695 - t683 / 0.432D3))
  638. t1675 = t125 * t116 * (t1640 + t1651 + t1667 - 0.18D2 * (t1629 + K
  639. #i * (-t727 - t726)) * t31)
  640. t1683 = 0.17D2 / 0.18D2 * t776
  641. t1687 = t1 * (Gi * (0.4D1 / 0.3D1 * t768 - t810) + Km * (t768 / 0.
  642. #2D1 - t810) + Ki * (t768 / 0.18D2 + t1683))
  643. t1689 = t776 / 0.2D1
  644. t1691 = t61 * (-t768 / 0.4D1 - t1689)
  645. t1692 = -Gi * t783
  646. t1694 = t776 / 0.3D1
  647. t1700 = t776 / 0.216D3
  648. t1706 = Gm * (t1691 + Km * (t1692 + Ki * (0.19D2 / 0.24D2 * t768 +
  649. # t1694)) - 0.36D2 * Gi * Ki * (-t768 / 0.27D2 - t1700))
  650. t1708 = t776 / 0.54D2
  651. t1714 = 0.18D2 * (Km * (t768 / 0.108D3 + t1708) + Ki * t804) * t31
  652. t1718 = t776 / 0.9D1
  653. t1719 = t768 / 0.3D1
  654. t1721 = t79 * (t1718 + t1719)
  655. t1724 = t776 / 0.6D1
  656. t1727 = 0.19D2 / 0.36D2 * t776
  657. t1732 = t1 * (Gi * (t776 - t1719) + Km * (t1724 - t1719) + Ki * (t
  658. #1727 + 0.17D2 / 0.36D2 * t768))
  659. t1733 = 0.23D2 / 0.24D2 * t776
  660. t1739 = 0.17D2 / 0.432D3 * t776
  661. t1746 = Gm * (t1691 + Km * (t1692 + Ki * (t1733 + t768 / 0.6D1)) -
  662. # 0.36D2 * Gi * Ki * (-t1739 - t768 / 0.432D3))
  663. t1750 = t831 / 0.3D1
  664. t1752 = t79 * (t1718 + t1750)
  665. t1761 = t1 * (Gi * (t776 - t1750) + Km * (t1724 - t1750) + Ki * (t
  666. #1727 + 0.17D2 / 0.36D2 * t831))
  667. t1764 = t61 * (-t1689 - t831 / 0.4D1)
  668. t1765 = -Gi * t835
  669. t1777 = Gm * (t1764 + Km * (t1765 + Ki * (t1733 + t831 / 0.6D1)) -
  670. # 0.36D2 * Gi * Ki * (-t1739 - t831 / 0.432D3))
  671. t1789 = Ki * t54
  672. t1790 = t1789 * t872
  673. t1792 = t1 * (t880 - t876 / 0.6D1 + t1790 - 0.16D2)
  674. t1797 = t1790 / 0.24D2
  675. t1802 = Gm * (-0.3D1 / 0.4D1 * t884 + Km * (-t880 + 0.9D1 / 0.8D1
  676. #* t1790 - 0.18D2) - 0.36D2 * Gi * (-t1797 + 0.2D1 / 0.3D1))
  677. t1811 = t1789 * t918
  678. t1812 = 0.16D2 * Mgi
  679. t1814 = t1 * (t924 - t921 / 0.6D1 + t1811 - t1812)
  680. t1817 = 0.18D2 * Mgi
  681. t1820 = t1811 / 0.24D2
  682. t1821 = 0.2D1 / 0.3D1 * Mgi
  683. t1826 = Gm * (-0.3D1 / 0.4D1 * t928 + Km * (-t924 + 0.9D1 / 0.8D1
  684. #* t1811 - t1817) - 0.36D2 * Gi * (-t1820 + t1821))
  685. t1835 = t1789 * t951
  686. t1837 = t1 * (t957 - t954 / 0.6D1 + t1835 + t1812)
  687. t1842 = t1835 / 0.24D2
  688. t1847 = Gm * (-0.3D1 / 0.4D1 * t960 + Km * (-t957 + 0.9D1 / 0.8D1
  689. #* t1835 + t1817) - 0.36D2 * Gi * (-t1842 - t1821))
  690. t1856 = t1789 * t981
  691. t1858 = t1 * (t987 - t984 / 0.6D1 + t1856 + 0.16D2)
  692. t1863 = t1856 / 0.24D2
  693. t1868 = Gm * (-0.3D1 / 0.4D1 * t990 + Km * (-t987 + 0.9D1 / 0.8D1
  694. #* t1856 + 0.18D2) - 0.36D2 * Gi * (-t1863 - 0.2D1 / 0.3D1))
  695. t1877 = t1789 * t1023
  696. t1878 = 0.16D2 * Mgm
  697. t1880 = t1 * (t1029 - t1026 / 0.6D1 + t1877 + t1878)
  698. t1883 = 0.18D2 * Mgm
  699. t1886 = t1877 / 0.24D2
  700. t1887 = 0.2D1 / 0.3D1 * Mgm
  701. t1892 = Gm * (-0.3D1 / 0.4D1 * t1033 + Km * (-t1029 + 0.9D1 / 0.8D
  702. #1 * t1877 + t1883) - 0.36D2 * Gi * (-t1886 - t1887))
  703. t1901 = t1789 * t1056
  704. t1903 = t1 * (t1062 - t1059 / 0.6D1 + t1901 - t1878)
  705. t1908 = t1901 / 0.24D2
  706. t1913 = Gm * (-0.3D1 / 0.4D1 * t1065 + Km * (-t1062 + 0.9D1 / 0.8D
  707. #1 * t1901 - t1883) - 0.36D2 * Gi * (-t1908 + t1887))
  708. t1921 = Km * t110
  709. t1922 = t55 / 0.48D2
  710. t1923 = t76 / 0.24D2
  711. t1932 = t125 * t116 * (t1264 + t1275 + t1291 + 0.12D2 * Gi * (t192
  712. #1 + Ki * (-0.3D1 / 0.2D1 + t1922 + t1923)) * Km)
  713. t1941 = t125 * t116 * (t129 + t1230 + t1249 + 0.12D2 * Gi * (t1921
  714. # + Ki * (t1923 + t1922)) * Km)
  715. t1942 = t202 / 0.3D1
  716. t1968 = t225 / 0.24D2
  717. t1977 = t125 * t116 * (t79 * (t1942 + t1342) + t1 * (Gi * (-t1942
  718. #+ t225) + Km * (-t1942 + t1348) + Ki * (0.17D2 / 0.36D2 * t202 + t
  719. #1351)) + Gm * (t1315 + Km * (t1316 + Ki * (t202 / 0.6D1 + t1361))
  720. #- 0.36D2 * Ki * (-t202 / 0.432D3 - t1367) * Gi) + 0.12D2 * Gi * (K
  721. #m * t253 + Ki * (t202 / 0.48D2 + t1968)) * Km)
  722. t2008 = 0.12D2 * Gi * (Km * t292 + Ki * (t1968 + t268 / 0.48D2)) *
  723. # Km
  724. t2011 = t125 * t116 * (t271 + t1 * (Gi * (-t259 + 0.4D1 / 0.3D1 *
  725. #t268) + Km * (-t259 + t268 / 0.2D1) + Ki * (t1307 + t268 / 0.18D2)
  726. #) + Gm * (t1359 + Km * (t1360 + Ki * (t1318 + 0.19D2 / 0.24D2 * t2
  727. #68)) - 0.36D2 * Ki * (-t1324 - t268 / 0.27D2) * Gi) + t2008)
  728. t2014 = t125 * t116 * (t1345 + t1356 + t1374 + t2008)
  729. t2015 = t346 / 0.3D1
  730. t2043 = t375 / 0.24D2
  731. t2052 = t125 * t116 * (t79 * (t2015 + t1427) + t1 * (Gi * (-t2015
  732. #+ t375 - 0.24D2) + Km * (-t2015 + t1433) + Ki * (0.17D2 / 0.36D2 *
  733. # t346 + t1436 + 0.16D2)) + Gm * (t1400 + Km * (Gi * (-t382 - t376
  734. #+ 0.8D1) + Ki * (0.18D2 + t346 / 0.6D1 + t1447)) - 0.36D2 * Ki * (
  735. #-0.1D1 / 0.6D1 - t346 / 0.432D3 - t1453) * Gi) + 0.12D2 * Gi * (Km
  736. # * t408 + Ki * (0.3D1 / 0.2D1 + t346 / 0.48D2 + t2043)) * Km)
  737. t2077 = Km * t449
  738. t2078 = t421 / 0.48D2
  739. t2087 = t125 * t116 * (t424 + t1 * (Gi * (-t414 + 0.4D1 / 0.3D1 *
  740. #t421 + 0.48D2) + Km * (-t414 + t421 / 0.2D1) + Ki * (t1392 + t421
  741. #/ 0.18D2 + 0.16D2)) + Gm * (t1444 + Km * (-Gi * t427 + Ki * (0.18D
  742. #2 + t1403 + 0.19D2 / 0.24D2 * t421)) - 0.36D2 * Ki * (-0.5D1 / 0.3
  743. #D1 - t1409 - t421 / 0.27D2) * Gi) + 0.12D2 * Gi * (t2077 + Ki * (t
  744. #2043 + t2078)) * Km)
  745. t2096 = t125 * t116 * (t1430 + t1441 + t1460 + 0.12D2 * Gi * (t207
  746. #7 + Ki * (0.3D1 / 0.2D1 + t2043 + t2078)) * Km)
  747. t2097 = t508 / 0.3D1
  748. t2123 = t546 / 0.24D2
  749. t2132 = t125 * t116 * (t79 * (t2097 - 0.24D2 + t1514) + t1 * (Gi *
  750. # (-t2097 + 0.24D2 + t546) + Km * (-t2097 + 0.24D2 + t1520) + Ki *
  751. #(0.17D2 / 0.36D2 * t508 - 0.34D2 + t1523)) + Gm * (t1487 + Km * (t
  752. #1488 + Ki * (-0.12D2 + t508 / 0.6D1 + t1533)) - 0.36D2 * Ki * (0.1
  753. #D1 / 0.6D1 - t508 / 0.432D3 - t1539) * Gi) + 0.12D2 * Gi * (Km * t
  754. #579 + Ki * (-0.3D1 / 0.2D1 + t508 / 0.48D2 + t2123)) * Km)
  755. t2163 = 0.12D2 * Gi * (Km * t622 + Ki * (-0.3D1 / 0.2D1 + t2123 +
  756. #t594 / 0.48D2)) * Km
  757. t2174 = 0.12D2 * Gi * (Km * t660 + Ki * (-0.3D1 / 0.2D1 + t2123 +
  758. #t636 / 0.48D2)) * Km
  759. t2178 = Km * t728
  760. t2179 = t683 / 0.48D2
  761. t2180 = t695 / 0.24D2
  762. t2189 = t125 * t116 * (t1640 + t1651 + t1667 + 0.12D2 * Gi * (t217
  763. #8 + Ki * (t2179 + t2180)) * Km)
  764. t2198 = t125 * t116 * (t736 + t1606 + t1625 + 0.12D2 * Gi * (t2178
  765. # + Ki * (-0.3D1 / 0.2D1 + t2180 + t2179)) * Km)
  766. t2201 = t776 / 0.24D2
  767. t2207 = 0.12D2 * Gi * (Km * t804 + Ki * (t768 / 0.48D2 + t2201)) *
  768. # Km
  769. t2210 = t125 * t116 * (t1721 + t1732 + t1746 + t2207)
  770. t2221 = 0.12D2 * Gi * (Km * t853 + Ki * (t2201 + t831 / 0.48D2)) *
  771. # Km
  772. t2232 = t125 * t116 * (t875 + t1792 + t1802 + 0.12D2 * Gi * (-t892
  773. # + t1790 / 0.16D2 - 0.1D1) * Km)
  774. t2240 = t125 * t116 * (t920 + t1814 + t1826 + 0.12D2 * Gi * (-t937
  775. # + t1811 / 0.16D2 - Mgi) * Km)
  776. t2248 = t125 * t116 * (t953 + t1837 + t1847 + 0.12D2 * Gi * (-t967
  777. # + t1835 / 0.16D2 + Mgi) * Km)
  778. t2256 = t125 * t116 * (t983 + t1858 + t1868 + 0.12D2 * Gi * (-t997
  779. # + t1856 / 0.16D2 + 0.1D1) * Km)
  780. t2264 = t125 * t116 * (t1025 + t1880 + t1892 + 0.12D2 * Gi * (-t10
  781. #42 + t1877 / 0.16D2 + Mgm) * Km)
  782. t2272 = t125 * t116 * (t1058 + t1903 + t1913 + 0.12D2 * Gi * (-t10
  783. #72 + t1901 / 0.16D2 - Mgm) * Km)
  784. t2329 = t79 * (t56 - t77)
  785. t2337 = t1 * (0.5D1 / 0.6D1 * t130 * t49 + Gi * (0.32D2 + t84 + t7
  786. #7) + Ki * (0.8D1 + t87 - t88))
  787. t2345 = Gm * (t96 + Km * (Gi * (0.16D2 + t97 + t88) + Ki * (0.9D1
  788. #+ t100 - t101)) + t159)
  789. t2348 = -t55 / 0.216D3 - t76 / 0.108D3
  790. t2362 = t125 * t116 * (-0.18D2 * Gi * t61 * t2348 + t129 + t137 +
  791. #t145)
  792. t2364 = t225 / 0.108D3
  793. t2379 = t125 * t116 * (t271 + t279 + t290 - 0.18D2 * Gi * t61 * (-
  794. #t2364 - t268 / 0.216D3))
  795. t2399 = t375 / 0.108D3
  796. t2416 = t125 * t116 * (t424 + t434 + t447 - 0.18D2 * Gi * t61 * (-
  797. #t2399 - t421 / 0.216D3))
  798. t2418 = t546 / 0.108D3
  799. t2443 = t79 * (t684 - 0.64D2 - t696)
  800. t2451 = t1 * (Km * (-0.60D2 + t699) + Gi * (t702 - 0.48D2 + t696)
  801. #+ Ki * (t705 - 0.48D2 - t706))
  802. t2459 = Gm * (t714 + Km * (Gi * (t715 - 0.36D2 + t706) + Ki * (t71
  803. #8 - 0.39D2 - t719)) - t159)
  804. t2462 = 0.1D1 / 0.3D1 - t683 / 0.216D3 - t695 / 0.108D3
  805. t2476 = t125 * t116 * (-0.18D2 * Gi * t61 * t2462 + t736 + t743 +
  806. #t751)
  807. t2478 = t776 / 0.108D3
  808. t2482 = 0.18D2 * Gi * t61 * (-t768 / 0.216D3 - t2478)
  809. t2506 = Mgi / 0.3D1
  810. t2531 = Mgm / 0.3D1
  811. t2553 = t125 * t116 * (t2329 + t2337 + t2345 - 0.6D1 * Gi * (t1921
  812. # - t13) * Km)
  813. t2578 = t125 * t116 * (t79 * (-t376 + t1124) + t1 * (0.5D1 / 0.6D1
  814. # * t130 * t420 + Gi * (t376 + t422 - 0.32D2) + Ki * (-t386 + t1130
  815. # - 0.8D1)) + Gm * (t437 + Km * (Gi * (t386 + t430 - 0.16D2) + Ki *
  816. # (-t399 + t1135 - 0.9D1)) - t159) - 0.6D1 * Gi * (t2077 + t13) * K
  817. #m)
  818. t2585 = t125 * t116 * (t2443 + t2451 + t2459 - 0.6D1 * Gi * (t2178
  819. # + t13) * Km)
  820. t2586 = t76 / 0.72D2
  821. t2587 = t225 / 0.72D2
  822. t2588 = t375 / 0.72D2
  823. t2590 = t546 / 0.72D2
  824. t2591 = t726 - 0.1D1
  825. t2592 = t695 / 0.72D2
  826. t2593 = t776 / 0.72D2
  827. t2595 = t54 * t872 / 0.72D2
  828. t2597 = t54 * t918 / 0.72D2
  829. t2599 = t54 * t951 / 0.72D2
  830. t2601 = t54 * t981 / 0.72D2
  831. t2603 = t54 * t1023 / 0.72D2
  832. t2605 = t54 * t1056 / 0.72D2
  833. JEA(1,1) = t126
  834. JEA(1,2) = t148
  835. JEA(1,3) = t148
  836. JEA(1,4) = t125 * t116 * (t228 + t237 + t250 - t255)
  837. JEA(1,5) = t297
  838. JEA(1,6) = t297
  839. JEA(1,7) = t125 * t116 * (t79 * (t347 - t376 + 0.32D2) + t1 * (Km
  840. #* (0.36D2 + 0.5D1 / 0.6D1 * t346) + t384 + Ki * (0.16D2 + t385 - t
  841. #386)) + Gm * (t394 + Km * (t397 + Ki * (0.18D2 + t398 - t399))) -
  842. #t410)
  843. JEA(1,8) = t454
  844. JEA(1,9) = t454
  845. JEA(1,10) = t125 * t116 * (t549 + t561 + t576 - t581)
  846. JEA(1,11) = t627
  847. JEA(1,12) = t665
  848. JEA(1,13) = t733
  849. JEA(1,14) = t754
  850. JEA(1,15) = t754
  851. JEA(1,16) = t809
  852. JEA(1,17) = t826
  853. JEA(1,18) = t858
  854. JEA(1,19) = t898
  855. JEA(1,20) = t943
  856. JEA(1,21) = t973
  857. JEA(1,22) = t1003
  858. JEA(1,23) = t1048
  859. JEA(1,24) = t1078
  860. JEA(2,1) = t148
  861. JEA(2,2) = t126
  862. JEA(2,3) = t148
  863. JEA(2,4) = t1091
  864. JEA(2,5) = t1107
  865. JEA(2,6) = t297
  866. JEA(2,7) = t1123
  867. JEA(2,8) = t1144
  868. JEA(2,9) = t454
  869. JEA(2,10) = t1160
  870. JEA(2,11) = t1181
  871. JEA(2,12) = t665
  872. JEA(2,13) = t754
  873. JEA(2,14) = t733
  874. JEA(2,15) = t754
  875. JEA(2,16) = t826
  876. JEA(2,17) = t809
  877. JEA(2,18) = t858
  878. JEA(2,19) = t898
  879. JEA(2,20) = t943
  880. JEA(2,21) = t973
  881. JEA(2,22) = t1003
  882. JEA(2,23) = t1048
  883. JEA(2,24) = t1078
  884. JEA(3,1) = t148
  885. JEA(3,2) = t148
  886. JEA(3,3) = t126
  887. JEA(3,4) = t1091
  888. JEA(3,5) = t297
  889. JEA(3,6) = t1107
  890. JEA(3,7) = t1123
  891. JEA(3,8) = t454
  892. JEA(3,9) = t1144
  893. JEA(3,10) = t1160
  894. JEA(3,11) = t627
  895. JEA(3,12) = t1202
  896. JEA(3,13) = t754
  897. JEA(3,14) = t754
  898. JEA(3,15) = t733
  899. JEA(3,16) = t826
  900. JEA(3,17) = t826
  901. JEA(3,18) = t1218
  902. JEA(3,19) = t898
  903. JEA(3,20) = t943
  904. JEA(3,21) = t973
  905. JEA(3,22) = t1003
  906. JEA(3,23) = t1048
  907. JEA(3,24) = t1078
  908. JEA(4,1) = t125 * t116 * (t129 + t1230 + t1249 - 0.18D2 * (Ki * t1
  909. #10 + t1253) * t31)
  910. JEA(4,2) = t1299
  911. JEA(4,3) = t1299
  912. JEA(4,4) = t125 * t116 * (t1080 + t1 * (Gi * (0.4D1 / 0.3D1 * t202
  913. # - t259) + Km * (t202 / 0.2D1 - t259) + Ki * (t202 / 0.18D2 + t130
  914. #7)) + Gm * (t1315 + Km * (t1316 + Ki * (0.19D2 / 0.24D2 * t202 + t
  915. #1318)) - 0.36D2 * Gi * Ki * (-t202 / 0.27D2 - t1324)) - 0.18D2 * (
  916. #Km * (t202 / 0.108D3 + t1332) + Ki * t253) * t31)
  917. JEA(4,5) = t1384
  918. JEA(4,6) = t1384
  919. JEA(4,7) = t125 * t116 * (t1109 + t1 * (Gi * (0.4D1 / 0.3D1 * t346
  920. # - t414 + 0.48D2) + Km * (t346 / 0.2D1 - t414) + Ki * (t346 / 0.18
  921. #D2 + t1392 + 0.16D2)) + Gm * (t1400 + Km * (-Gi * t383 + Ki * (0.1
  922. #8D2 + 0.19D2 / 0.24D2 * t346 + t1403)) - 0.36D2 * Gi * Ki * (-0.5D
  923. #1 / 0.3D1 - t346 / 0.27D2 - t1409)) - 0.18D2 * (Km * (t346 / 0.108
  924. #D3 + t1417) + Ki * t408) * t31)
  925. JEA(4,8) = t1471
  926. JEA(4,9) = t1471
  927. JEA(4,10) = t125 * t116 * (t1146 + t1 * (Gi * (0.4D1 / 0.3D1 * t50
  928. #8 - 0.96D2 - t585) + Km * (t508 / 0.2D1 - 0.36D2 - t585) + Ki * (t
  929. #508 / 0.18D2 - 0.4D1 + t1479)) + Gm * (t1487 + Km * (t1488 + Ki *
  930. #(-0.57D2 + 0.19D2 / 0.24D2 * t508 + t1490)) - 0.36D2 * Gi * Ki * (
  931. #0.8D1 / 0.3D1 - t508 / 0.27D2 - t1496)) - 0.18D2 * (Km * (t508 / 0
  932. #.108D3 - 0.2D1 / 0.3D1 + t1504) + Ki * t579) * t31)
  933. JEA(4,11) = t125 * t116 * (t1517 + t1528 + t1546 - 0.18D2 * (Km *
  934. #(t1504 + t594 / 0.108D3 - 0.2D1 / 0.3D1) + Ki * t622) * t31)
  935. JEA(4,12) = t125 * t116 * (t1559 + t1568 + t1584 - 0.18D2 * (Km *
  936. #(t1504 + t636 / 0.108D3 - 0.2D1 / 0.3D1) + Ki * t660) * t31)
  937. JEA(4,13) = t125 * t116 * (t736 + t1606 + t1625 - 0.18D2 * (Ki * t
  938. #728 + t1629) * t31)
  939. JEA(4,14) = t1675
  940. JEA(4,15) = t1675
  941. JEA(4,16) = t125 * t116 * (t812 + t1687 + t1706 - t1714)
  942. JEA(4,17) = t125 * t116 * (t1721 + t1732 + t1746 - t1714)
  943. JEA(4,18) = t125 * t116 * (t1752 + t1761 + t1777 - 0.18D2 * (Km *
  944. #(t1708 + t831 / 0.108D3) + Ki * t853) * t31)
  945. JEA(4,19) = t125 * t116 * (t875 + t1792 + t1802 - 0.18D2 * (t876 /
  946. # 0.36D2 - t1797 + 0.2D1 / 0.3D1) * t31)
  947. JEA(4,20) = t125 * t116 * (t920 + t1814 + t1826 - 0.18D2 * (t921 /
  948. # 0.36D2 - t1820 + t1821) * t31)
  949. JEA(4,21) = t125 * t116 * (t953 + t1837 + t1847 - 0.18D2 * (t954 /
  950. # 0.36D2 - t1842 - t1821) * t31)
  951. JEA(4,22) = t125 * t116 * (t983 + t1858 + t1868 - 0.18D2 * (t984 /
  952. # 0.36D2 - t1863 - 0.2D1 / 0.3D1) * t31)
  953. JEA(4,23) = t125 * t116 * (t1025 + t1880 + t1892 - 0.18D2 * (t1026
  954. # / 0.36D2 - t1886 - t1887) * t31)
  955. JEA(4,24) = t125 * t116 * (t1058 + t1903 + t1913 - 0.18D2 * (t1059
  956. # / 0.36D2 - t1908 + t1887) * t31)
  957. JEA(5,1) = t1932
  958. JEA(5,2) = t1941
  959. JEA(5,3) = t1932
  960. JEA(5,4) = t1977
  961. JEA(5,5) = t2011
  962. JEA(5,6) = t2014
  963. JEA(5,7) = t2052
  964. JEA(5,8) = t2087
  965. JEA(5,9) = t2096
  966. JEA(5,10) = t2132
  967. JEA(5,11) = t125 * t116 * (t597 + t1 * (Gi * (-t585 + 0.4D1 / 0.3D
  968. #1 * t594 - 0.96D2) + Km * (-t585 + t594 / 0.2D1 - 0.36D2) + Ki * (
  969. #t1479 + t594 / 0.18D2 - 0.4D1)) + Gm * (t1531 + Km * (t1532 + Ki *
  970. # (-0.57D2 + t1490 + 0.19D2 / 0.24D2 * t594)) - 0.36D2 * Ki * (0.8D
  971. #1 / 0.3D1 - t1496 - t594 / 0.27D2) * Gi) + t2163)
  972. JEA(5,12) = t125 * t116 * (t1559 + t1568 + t1584 + t2174)
  973. JEA(5,13) = t2189
  974. JEA(5,14) = t2198
  975. JEA(5,15) = t2189
  976. JEA(5,16) = t2210
  977. JEA(5,17) = t125 * t116 * (t812 + t1687 + t1706 + t2207)
  978. JEA(5,18) = t125 * t116 * (t1752 + t1761 + t1777 + t2221)
  979. JEA(5,19) = t2232
  980. JEA(5,20) = t2240
  981. JEA(5,21) = t2248
  982. JEA(5,22) = t2256
  983. JEA(5,23) = t2264
  984. JEA(5,24) = t2272
  985. JEA(6,1) = t1932
  986. JEA(6,2) = t1932
  987. JEA(6,3) = t1941
  988. JEA(6,4) = t1977
  989. JEA(6,5) = t2014
  990. JEA(6,6) = t2011
  991. JEA(6,7) = t2052
  992. JEA(6,8) = t2096
  993. JEA(6,9) = t2087
  994. JEA(6,10) = t2132
  995. JEA(6,11) = t125 * t116 * (t1517 + t1528 + t1546 + t2163)
  996. JEA(6,12) = t125 * t116 * (t639 + t1 * (Gi * (-t585 + 0.4D1 / 0.3D
  997. #1 * t636 - 0.96D2) + Km * (-t585 + t636 / 0.2D1 - 0.36D2) + Ki * (
  998. #t1479 + t636 / 0.18D2 - 0.4D1)) + Gm * (t1571 + Km * (t1572 + Ki *
  999. # (-0.57D2 + t1490 + 0.19D2 / 0.24D2 * t636)) - 0.36D2 * Ki * (0.8D
  1000. #1 / 0.3D1 - t1496 - t636 / 0.27D2) * Gi) + t2174)
  1001. JEA(6,13) = t2189
  1002. JEA(6,14) = t2189
  1003. JEA(6,15) = t2198
  1004. JEA(6,16) = t2210
  1005. JEA(6,17) = t2210
  1006. JEA(6,18) = t125 * t116 * (t834 + t1 * (Gi * (-t810 + 0.4D1 / 0.3D
  1007. #1 * t831) + Km * (-t810 + t831 / 0.2D1) + Ki * (t1683 + t831 / 0.1
  1008. #8D2)) + Gm * (t1764 + Km * (t1765 + Ki * (t1694 + 0.19D2 / 0.24D2
  1009. #* t831)) - 0.36D2 * Ki * (-t1700 - t831 / 0.27D2) * Gi) + t2221)
  1010. JEA(6,19) = t2232
  1011. JEA(6,20) = t2240
  1012. JEA(6,21) = t2248
  1013. JEA(6,22) = t2256
  1014. JEA(6,23) = t2264
  1015. JEA(6,24) = t2272
  1016. JEA(7,1) = t125 * t116 * (t2329 + t2337 + t2345 - 0.18D2 * Gi * Km
  1017. # * (Km * t2348 - Ki))
  1018. JEA(7,2) = t2362
  1019. JEA(7,3) = t2362
  1020. JEA(7,4) = t125 * t116 * (t228 + t237 + t250 - 0.18D2 * Gi * t61 *
  1021. # (-t202 / 0.216D3 - t2364))
  1022. JEA(7,5) = t2379
  1023. JEA(7,6) = t2379
  1024. JEA(7,7) = t125 * t116 * (t79 * (t347 - t376) + t1 * (0.5D1 / 0.6D
  1025. #1 * t130 * t345 + Gi * (t382 + t376 - 0.32D2) + Ki * (t385 - t386
  1026. #- 0.8D1)) + Gm * (t394 + Km * (Gi * (t395 + t386 - 0.16D2) + Ki *
  1027. #(t398 - t399 - 0.9D1)) - t159) - 0.18D2 * Gi * Km * (Km * (-t346 /
  1028. # 0.216D3 - t2399) + Ki))
  1029. JEA(7,8) = t2416
  1030. JEA(7,9) = t2416
  1031. JEA(7,10) = t125 * t116 * (t549 + t561 + t576 - 0.18D2 * Gi * t61
  1032. #* (0.1D1 / 0.3D1 - t508 / 0.216D3 - t2418))
  1033. JEA(7,11) = t125 * t116 * (t597 + t607 + t620 - 0.18D2 * Gi * t61
  1034. #* (0.1D1 / 0.3D1 - t2418 - t594 / 0.216D3))
  1035. JEA(7,12) = t125 * t116 * (t639 + t646 + t658 - 0.18D2 * Gi * t61
  1036. #* (0.1D1 / 0.3D1 - t2418 - t636 / 0.216D3))
  1037. JEA(7,13) = t125 * t116 * (t2443 + t2451 + t2459 - 0.18D2 * Gi * K
  1038. #m * (Km * t2462 + Ki))
  1039. JEA(7,14) = t2476
  1040. JEA(7,15) = t2476
  1041. JEA(7,16) = t125 * t116 * (t779 + t788 + t801 - t2482)
  1042. JEA(7,17) = t125 * t116 * (t812 + t818 + t823 - t2482)
  1043. JEA(7,18) = t125 * t116 * (t834 + t840 + t851 - 0.18D2 * Gi * t61
  1044. #* (-t2478 - t831 / 0.216D3))
  1045. JEA(7,19) = t125 * t116 * (t875 + t882 + t891 - 0.18D2 * Gi * Km *
  1046. # (-t876 / 0.72D2 - 0.1D1 / 0.3D1))
  1047. JEA(7,20) = t125 * t116 * (t920 + t927 + t936 - 0.18D2 * Gi * Km *
  1048. # (-t921 / 0.72D2 - t2506))
  1049. JEA(7,21) = t125 * t116 * (t953 + t959 + t966 - 0.18D2 * Gi * Km *
  1050. # (-t954 / 0.72D2 + t2506))
  1051. JEA(7,22) = t125 * t116 * (t983 + t989 + t996 - 0.18D2 * Gi * Km *
  1052. # (-t984 / 0.72D2 + 0.1D1 / 0.3D1))
  1053. JEA(7,23) = t125 * t116 * (t1025 + t1032 + t1041 - 0.18D2 * Gi * K
  1054. #m * (-t1026 / 0.72D2 + t2531))
  1055. JEA(7,24) = t125 * t116 * (t1058 + t1064 + t1071 - 0.18D2 * Gi * K
  1056. #m * (-t1059 / 0.72D2 - t2531))
  1057. JEA(8,1) = t148
  1058. JEA(8,2) = t2553
  1059. JEA(8,3) = t148
  1060. JEA(8,4) = t1091
  1061. JEA(8,5) = t1107
  1062. JEA(8,6) = t297
  1063. JEA(8,7) = t1123
  1064. JEA(8,8) = t2578
  1065. JEA(8,9) = t454
  1066. JEA(8,10) = t1160
  1067. JEA(8,11) = t1181
  1068. JEA(8,12) = t665
  1069. JEA(8,13) = t754
  1070. JEA(8,14) = t2585
  1071. JEA(8,15) = t754
  1072. JEA(8,16) = t826
  1073. JEA(8,17) = t809
  1074. JEA(8,18) = t858
  1075. JEA(8,19) = t898
  1076. JEA(8,20) = t943
  1077. JEA(8,21) = t973
  1078. JEA(8,22) = t1003
  1079. JEA(8,23) = t1048
  1080. JEA(8,24) = t1078
  1081. JEA(9,1) = t148
  1082. JEA(9,2) = t148
  1083. JEA(9,3) = t2553
  1084. JEA(9,4) = t1091
  1085. JEA(9,5) = t297
  1086. JEA(9,6) = t1107
  1087. JEA(9,7) = t1123
  1088. JEA(9,8) = t454
  1089. JEA(9,9) = t2578
  1090. JEA(9,10) = t1160
  1091. JEA(9,11) = t627
  1092. JEA(9,12) = t1202
  1093. JEA(9,13) = t754
  1094. JEA(9,14) = t754
  1095. JEA(9,15) = t2585
  1096. JEA(9,16) = t826
  1097. JEA(9,17) = t826
  1098. JEA(9,18) = t1218
  1099. JEA(9,19) = t898
  1100. JEA(9,20) = t943
  1101. JEA(9,21) = t973
  1102. JEA(9,22) = t1003
  1103. JEA(9,23) = t1048
  1104. JEA(9,24) = t1078
  1105. JEA(10,1) = t108
  1106. JEA(10,2) = t2586
  1107. JEA(10,3) = t2586
  1108. JEA(10,4) = t251
  1109. JEA(10,5) = t2587
  1110. JEA(10,6) = t2587
  1111. JEA(10,7) = t406
  1112. JEA(10,8) = t2588
  1113. JEA(10,9) = t2588
  1114. JEA(10,10) = t577 - 0.1D1
  1115. JEA(10,11) = t2590
  1116. JEA(10,12) = t2590
  1117. JEA(10,13) = t2591
  1118. JEA(10,14) = t2592
  1119. JEA(10,15) = t2592
  1120. JEA(10,16) = t802
  1121. JEA(10,17) = t2593
  1122. JEA(10,18) = t2593
  1123. JEA(10,19) = t2595
  1124. JEA(10,20) = t2597
  1125. JEA(10,21) = t2599
  1126. JEA(10,22) = t2601
  1127. JEA(10,23) = t2603
  1128. JEA(10,24) = t2605
  1129. JEA(11,1) = t2586
  1130. JEA(11,2) = t108
  1131. JEA(11,3) = t2586
  1132. JEA(11,4) = t2587
  1133. JEA(11,5) = t291
  1134. JEA(11,6) = t2587
  1135. JEA(11,7) = t2588
  1136. JEA(11,8) = t448
  1137. JEA(11,9) = t2588
  1138. JEA(11,10) = t2590
  1139. JEA(11,11) = t621 - 0.1D1
  1140. JEA(11,12) = t2590
  1141. JEA(11,13) = t2592
  1142. JEA(11,14) = t2591
  1143. JEA(11,15) = t2592
  1144. JEA(11,16) = t2593
  1145. JEA(11,17) = t802
  1146. JEA(11,18) = t2593
  1147. JEA(11,19) = t2595
  1148. JEA(11,20) = t2597
  1149. JEA(11,21) = t2599
  1150. JEA(11,22) = t2601
  1151. JEA(11,23) = t2603
  1152. JEA(11,24) = t2605
  1153. JEA(12,1) = t2586
  1154. JEA(12,2) = t2586
  1155. JEA(12,3) = t108
  1156. JEA(12,4) = t2587
  1157. JEA(12,5) = t2587
  1158. JEA(12,6) = t291
  1159. JEA(12,7) = t2588
  1160. JEA(12,8) = t2588
  1161. JEA(12,9) = t448
  1162. JEA(12,10) = t2590
  1163. JEA(12,11) = t2590
  1164. JEA(12,12) = t659 - 0.1D1
  1165. JEA(12,13) = t2592
  1166. JEA(12,14) = t2592
  1167. JEA(12,15) = t2591
  1168. JEA(12,16) = t2593
  1169. JEA(12,17) = t2593
  1170. JEA(12,18) = t852
  1171. JEA(12,19) = t2595
  1172. JEA(12,20) = t2597
  1173. JEA(12,21) = t2599
  1174. JEA(12,22) = t2601
  1175. JEA(12,23) = t2603
  1176. JEA(12,24) = t2605
  1177.  
  1178.  
  1179.  
  1180. if(affiche) then
  1181. print*,'Dans JEA1inc3D matrice JEA'
  1182. do i=1,12
  1183. write(*,'(24(E10.3))') (JEA(i,j),j=1,24)
  1184. end do
  1185. end if
  1186.  
  1187. else
  1188. print*,'Pb de dimansion dans jelc3d'
  1189. err1=1
  1190. return
  1191. end if
  1192.  
  1193. return
  1194. end
  1195.  
  1196.  
  1197.  
  1198.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales