C JELC3D SOURCE PV090527 23/01/27 21:15:45 11574 subroutine jelc3d(NBRINC,NINC,NDIMG,NDIMA,FRAC,MK,MG,MbG,JEA, # affiche,err1) c Jacibienne d epse / d eps imp dans le cas d un seul type d inclusion c organisation de la Jacobienne c ligne : [eie1, eie2, eie3, ere1, ere2, ere3, eoe1, eoe2, eoe3, eme1, eme2, eme3] c 1 2 3 4 5 6 7 8 9 10 11 12 = NDIMG 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] 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 implicit real*8 (a-h,o-z) implicit integer (i-n) logical affiche integer NBRINC,NDIMG,NDIMA,NINC,err1 c jacobienne depse/depsimp si un seul type dinclusion real*8 JEA(NDIMG,NDIMA) real*8 FRAC(0:NBRINC),MK(0:NBRINC),MG(0:NBRINC),MBG(0:NBRINC,0:1) real*8 Mgm,Mgi real*8 f,Km,Gm,Kgm,Ki,Gi,Kgi if((NINC.eq.1).and.(NBRINC.eq.1).and. # (NDIMG.ge.12).and.(NDIMA.ge.24)) then f=FRAC(1) Km=MK(0) Gm=MG(0) Mgm=MbG(0,1) Ki=MK(1) Gi=MG(1) Mgi=MbG(1,1) t1 = Gm ** 2 t3 = 0.24D2 * t1 * Ki t4 = 0.48D2 * Gi t5 = 0.12D2 * Ki t8 = Ki * Gi t9 = 0.72D2 * t8 t12 = 0.2D1 * Km t13 = 0.3D1 * Ki t14 = t12 - t13 t19 = f ** 2 t21 = 0.64D2 * Gi t22 = 0.16D2 * Ki t25 = 0.32D2 * Gi t26 = 0.18D2 * Ki t31 = Gi * Km t32 = Ki * t31 t33 = 0.36D2 * t32 t42 = 0.2D1 / 0.3D1 * t1 t44 = -Gi - 0.3D1 / 0.4D1 * Km t46 = t31 / 0.2D1 t48 = 0.1D1 / ((Gm + Km / 0.2D1) * (-Gm + Gi) * f - t42 + Gm * t44 # - t46) t49 = t48 * (t19 * (-t3 + Gm * (Km * (-t4 - t5) + t9) - 0.12D2 * G #i * Km * t14) + f * (t1 * (-t21 - t22) + Gm * (Km * (-t25 - t26) - # t9) - t33)) t54 = 0.1D1 / ((Ki - Km) * f - 0.4D1 / 0.3D1 * Gm - Ki) t55 = t54 * t49 t56 = 0.8D1 / 0.9D1 * t55 t57 = 0.24D2 * Gi t58 = t57 - t5 t59 = Km * t58 t61 = Km ** 2 t62 = Gi * t61 t63 = 0.12D2 * t62 t66 = t25 - t22 t68 = 0.16D2 * Gi t69 = t68 - t26 t75 = t48 * (t19 * (Gm * t59 - t3 + t63) + f * (Gm * Km * t69 + t1 # * t66)) t76 = t54 * t75 t77 = 0.4D1 / 0.9D1 * t76 t79 = t1 * Gm t84 = 0.2D1 / 0.9D1 * t55 t85 = -0.16D2 + t84 + t77 t86 = Gi * t85 t87 = 0.5D1 / 0.9D1 * t55 t88 = 0.5D1 / 0.9D1 * t76 t96 = t61 * (t55 / 0.8D1 + t76 / 0.4D1) t97 = 0.5D1 / 0.18D2 * t55 t99 = Gi * (-0.8D1 + t97 + t88) t100 = 0.5D1 / 0.12D2 * t55 t101 = 0.5D1 / 0.12D2 * t76 t108 = t55 / 0.72D2 t109 = t76 / 0.36D2 t110 = -t108 - t109 t112 = 0.6D1 * t110 * t62 t116 = 0.1D1 / (0.4D1 * Gm + t13) t125 = 0.1D1 / (0.12D2 * Gi * Gm + 0.9D1 * Gm * Km + 0.8D1 * t1 + #0.6D1 * t31) t126 = t125 * t116 * (t79 * (t56 - t77 - 0.32D2) + t1 * (Km * (-0. #36D2 + 0.5D1 / 0.6D1 * t55) + t86 + Ki * (-0.16D2 + t87 - t88)) + #Gm * (t96 + Km * (t99 + Ki * (-0.18D2 + t100 - t101))) - t112) t127 = 0.2D1 / 0.3D1 * t76 t129 = t79 * (t127 - t84) t130 = Km * t54 t137 = t1 * (0.5D1 / 0.6D1 * t130 * t75 + t86 + Ki * (0.8D1 + 0.5D #1 / 0.18D2 * t76 - t97)) t145 = Gm * (t96 + Km * (t99 + Ki * (0.9D1 + 0.5D1 / 0.24D2 * t76 #- 0.5D1 / 0.24D2 * t55))) t148 = t125 * t116 * (t129 + t137 + t145 - t112) t149 = 0.16D2 * t79 t150 = 0.24D2 * Km t153 = t1 * (t150 - t68 - 0.4D1 * Ki) t154 = 0.18D2 * t61 t155 = 0.3D1 / 0.4D1 * Ki t157 = Km * (Gi - t155) t159 = 0.36D2 * t8 t161 = Gm * (t154 + 0.16D2 * t157 - t159) t162 = Km / 0.3D1 t169 = t19 * f t171 = 0.12D2 * Km t186 = t19 * (t149 + t1 * (t171 + t21 + 0.28D2 * Ki) + Gm * (-t154 # + (0.8D1 * Gi + 0.39D2 * Ki) * Km + t9) - 0.12D2 * Gi * Km * (Km #- t13)) t187 = 0.32D2 * t79 t189 = 0.24D2 * Ki t199 = f * (-t187 + t1 * (-0.36D2 * Km - t4 - t189) + Gm * (Km * ( #-t57 - 0.27D2 * Ki) - t159) - 0.18D2 * t32) t201 = t48 * (t169 * (t149 + t153 + t161 - 0.36D2 * Gi * Km * (-t1 #62 + Ki / 0.2D1)) + t186 + t199) t202 = t54 * t201 t205 = -t68 + 0.20D2 * Ki t207 = 0.9D1 * t61 t211 = 0.6D1 * t62 t224 = t48 * (t169 * (t149 + t1 * t205 + Gm * (-t207 - 0.20D2 * t1 #57) - t211) + t19 * (-t149 - t1 * t205 + Gm * (t207 + (0.20D2 * Gi # - 0.15D2 * Ki) * Km) + t211)) t225 = t54 * t224 t226 = 0.4D1 / 0.9D1 * t225 t228 = t79 * (0.8D1 / 0.9D1 * t202 - t226) t231 = 0.2D1 / 0.9D1 * t202 t232 = t231 + t226 t233 = Gi * t232 t234 = t202 - t225 t237 = t1 * (0.5D1 / 0.6D1 * t130 * t201 + t233 + 0.5D1 / 0.9D1 * #Ki * t234) t239 = t225 / 0.4D1 t241 = t61 * (t202 / 0.8D1 + t239) t243 = 0.5D1 / 0.9D1 * t225 t245 = Gi * (0.5D1 / 0.18D2 * t202 + t243) t250 = Gm * (t241 + Km * (t245 + 0.5D1 / 0.12D2 * Ki * t234)) t251 = t202 / 0.72D2 t252 = t225 / 0.36D2 t253 = -t251 - t252 t255 = 0.6D1 * t253 * t62 t259 = 0.2D1 / 0.3D1 * t225 t260 = 0.3D1 / 0.2D1 * Ki t261 = -Km + t260 t267 = t48 * (t169 * (-0.12D2 * t261 * t31 + t149 + t153 + t161) + # t186 + t199) t268 = t54 * t267 t269 = 0.2D1 / 0.9D1 * t268 t271 = t79 * (t259 - t269) t273 = 0.5D1 / 0.6D1 * t130 * t224 t274 = t226 + t269 t275 = Gi * t274 t276 = t225 - t268 t279 = t1 * (t273 + t275 + 0.5D1 / 0.18D2 * Ki * t276) t282 = t61 * (t239 + t268 / 0.8D1) t285 = Gi * (t243 + 0.5D1 / 0.18D2 * t268) t290 = Gm * (t282 + Km * (t285 + 0.5D1 / 0.24D2 * Ki * t276)) t291 = t268 / 0.72D2 t292 = -t252 - t291 t294 = 0.6D1 * t292 * t62 t297 = t125 * t116 * (t271 + t279 + t290 - t294) t301 = t1 * (0.48D2 * Km - t25 - 0.8D1 * Ki) t302 = 0.36D2 * t61 t305 = Gm * (t302 + 0.32D2 * t157 - t9) t310 = 0.36D2 * (-0.2D1 / 0.3D1 * Km + Ki) * Km * Gi t330 = t19 * (t1 * (-t150 + 0.160D3 * Gi + 0.88D2 * Ki) + Gm * (-0 #.72D2 * t61 + Km * (t25 + 0.114D3 * Ki) + 0.144D3 * t8) - 0.12D2 * # Gi * Km * (t12 - 0.6D1 * Ki)) t343 = f * (-0.128D3 * t79 + t1 * (-0.144D3 * Km - 0.128D3 * Gi - #0.80D2 * Ki) + Gm * (Km * (-t21 - 0.90D2 * Ki) - t9) - t33) t345 = t48 * (t169 * (t187 + t301 + t305 - t310) + t330 + t343) t346 = t54 * t345 t347 = 0.8D1 / 0.9D1 * t346 t349 = -t25 + 0.40D2 * Ki t350 = t1 * t349 t356 = 0.64D2 * t79 t361 = 0.48D2 * Ki t374 = t48 * (t169 * (t187 + t350 + Gm * (-t154 - 0.40D2 * t157) - # t63) + t19 * (-t356 + t1 * (t21 - 0.56D2 * Ki) + Gm * (t302 + Km #* (0.56D2 * Gi - t361)) + t63) + f * (-Gm * Km * t69 - t1 * t66)) t375 = t54 * t374 t376 = 0.4D1 / 0.9D1 * t375 t382 = 0.2D1 / 0.9D1 * t346 t383 = t382 + t376 + 0.16D2 t384 = Gi * t383 t385 = 0.5D1 / 0.9D1 * t346 t386 = 0.5D1 / 0.9D1 * t375 t392 = t375 / 0.4D1 t394 = t61 * (t346 / 0.8D1 + t392) t395 = 0.5D1 / 0.18D2 * t346 t397 = Gi * (t395 + t386 + 0.8D1) t398 = 0.5D1 / 0.12D2 * t346 t399 = 0.5D1 / 0.12D2 * t375 t406 = t346 / 0.72D2 t407 = t375 / 0.36D2 t408 = -t406 - t407 t410 = 0.6D1 * t408 * t62 t414 = 0.2D1 / 0.3D1 * t375 t416 = -0.12D2 * t14 * t31 t420 = t48 * (t169 * (t187 + t301 + t305 - t416) + t330 + t343) t421 = t54 * t420 t422 = 0.2D1 / 0.9D1 * t421 t424 = t79 * (t414 - t422) t426 = 0.5D1 / 0.6D1 * t130 * t374 t427 = t376 + t422 + 0.16D2 t428 = Gi * t427 t429 = 0.5D1 / 0.18D2 * t375 t430 = 0.5D1 / 0.18D2 * t421 t434 = t1 * (t426 + t428 + Ki * (-0.8D1 + t429 - t430)) t437 = t61 * (t392 + t421 / 0.8D1) t439 = Gi * (t386 + t430 + 0.8D1) t440 = 0.5D1 / 0.24D2 * t375 t447 = Gm * (t437 + Km * (t439 + Ki * (-0.9D1 + t440 - 0.5D1 / 0.2 #4D2 * t421))) t448 = t421 / 0.72D2 t449 = -t407 - t448 t451 = 0.6D1 * t449 * t62 t454 = t125 * t116 * (t424 + t434 + t447 - t451) t455 = 0.48D2 * t79 t456 = 0.72D2 * Km t458 = t1 * (-t456 + t4 + t5) t461 = 0.108D3 * t8 t463 = Gm * (-0.54D2 * t61 - 0.48D2 * t157 + t461) t487 = t19 * (-t149 + t1 * (0.84D2 * Km - 0.224D3 * Gi - 0.164D3 * # Ki) + Gm * (0.126D3 * t61 + Km * (-t21 - 0.177D3 * Ki) - 0.216D3 #* t8) - 0.12D2 * Gi * Km * (-t12 + 0.9D1 * Ki)) t488 = 0.224D3 * t79 t493 = t1 * (0.228D3 * Km + 0.176D3 * Gi + 0.152D3 * Ki) t499 = Gm * (t154 + (0.112D3 * Gi + 0.141D3 * Ki) * Km + t461) t500 = 0.9D1 / 0.2D1 * Ki t505 = f * (t488 + t493 + t499 - 0.12D2 * (-Km - t500) * t31) t508 = t54 * t48 * (t169 * (0.36D2 * Gi * Km * t261 - t455 + t458 #+ t463) + t487 + t505) t538 = 0.40D2 * Gi - 0.30D2 * Ki t545 = t48 * (t169 * (-t455 + t1 * (t4 - 0.60D2 * Ki) + Gm * (0.27 #D2 * t61 + 0.60D2 * t157) + 0.18D2 * t62) + t19 * (0.80D2 * t79 + #t1 * (-0.80D2 * Gi + 0.100D3 * Ki) + Gm * (-0.45D2 * t61 + (-0.100 #D3 * Gi + 0.75D2 * Ki) * Km) - 0.30D2 * t62) + f * (-t187 - t1 * t #349 + Gm * (Km * t538 + t154) + t63)) t546 = t54 * t545 t547 = 0.4D1 / 0.9D1 * t546 t549 = t79 * (0.8D1 / 0.9D1 * t508 - t547 - 0.64D2) t553 = 0.2D1 / 0.9D1 * t508 t554 = t553 - 0.16D2 + t547 t555 = Gi * t554 t557 = 0.5D1 / 0.9D1 * t546 t561 = t1 * (Km * (-0.60D2 + 0.5D1 / 0.6D1 * t508) + t555 + Ki * ( #-0.40D2 + 0.5D1 / 0.9D1 * t508 - t557)) t563 = t546 / 0.4D1 t565 = t61 * (t508 / 0.8D1 - 0.9D1 + t563) t566 = 0.5D1 / 0.18D2 * t508 t568 = Gi * (t566 - 0.20D2 + t557) t570 = 0.5D1 / 0.12D2 * t546 t576 = Gm * (t565 + Km * (t568 + Ki * (-0.30D2 + 0.5D1 / 0.12D2 * #t508 - t570))) t577 = t508 / 0.72D2 t578 = t546 / 0.36D2 t579 = 0.1D1 - t577 - t578 t581 = 0.6D1 * t579 * t62 t585 = 0.2D1 / 0.3D1 * t546 t591 = t169 * (-t455 + t458 + t463 - 0.12D2 * (0.3D1 * Km - t500) #* t31) t594 = t54 * t48 * (t591 + t487 + t505) t595 = 0.2D1 / 0.9D1 * t594 t597 = t79 * (t585 - t595 + 0.16D2) t599 = 0.5D1 / 0.6D1 * t130 * t545 t600 = t547 + t595 - 0.16D2 t601 = Gi * t600 t602 = 0.5D1 / 0.18D2 * t546 t603 = 0.5D1 / 0.18D2 * t594 t607 = t1 * (t599 + t601 + Ki * (0.20D2 + t602 - t603)) t610 = t61 * (t563 + t594 / 0.8D1 - 0.9D1) t612 = Gi * (t557 + t603 - 0.20D2) t613 = 0.5D1 / 0.24D2 * t546 t620 = Gm * (t610 + Km * (t612 + Ki * (0.15D2 + t613 - 0.5D1 / 0.2 #4D2 * t594))) t621 = t594 / 0.72D2 t622 = 0.1D1 - t578 - t621 t624 = 0.6D1 * t622 * t62 t627 = t125 * t116 * (t597 + t607 + t620 - t624) t636 = t54 * t48 * (t591 + t487 + f * (t488 + t493 + t499 - 0.36D2 # * Gi * Km * (-t162 - t260))) t637 = 0.2D1 / 0.9D1 * t636 t639 = t79 * (t585 - t637 + 0.16D2) t640 = t547 + t637 - 0.16D2 t641 = Gi * t640 t642 = 0.5D1 / 0.18D2 * t636 t646 = t1 * (t599 + t641 + Ki * (0.20D2 + t602 - t642)) t649 = t61 * (t563 + t636 / 0.8D1 - 0.9D1) t651 = Gi * (t557 + t642 - 0.20D2) t658 = Gm * (t649 + Km * (t651 + Ki * (0.15D2 + t613 - 0.5D1 / 0.2 #4D2 * t636))) t659 = t636 / 0.72D2 t660 = 0.1D1 - t578 - t659 t662 = 0.6D1 * t660 * t62 t665 = t125 * t116 * (t639 + t646 + t658 - t662) t668 = -Gi - Ki t683 = t54 * t48 * (t19 * (t1 * (t456 - t361) + Gm * (0.24D2 * Km #* t668 + t302) - t63) + f * (t356 + 0.48D2 * t1 * (Km + Ki) + Gm * # (-0.24D2 * Km * t668 + t154) + t63)) t684 = 0.8D1 / 0.9D1 * t683 t694 = t48 * (t19 * (-Gm * Km * t58 + t3 - t63) + f * (-t187 - t3 #+ Gm * (t154 + t59) + t63)) t695 = t54 * t694 t696 = 0.4D1 / 0.9D1 * t695 t699 = 0.5D1 / 0.6D1 * t683 t702 = 0.2D1 / 0.9D1 * t683 t703 = t702 + t696 t704 = Gi * t703 t705 = 0.5D1 / 0.9D1 * t683 t706 = 0.5D1 / 0.9D1 * t695 t714 = t61 * (t683 / 0.8D1 - 0.9D1 + t695 / 0.4D1) t715 = 0.5D1 / 0.18D2 * t683 t717 = Gi * (t715 - 0.12D2 + t706) t718 = 0.5D1 / 0.12D2 * t683 t719 = 0.5D1 / 0.12D2 * t695 t726 = t683 / 0.72D2 t727 = t695 / 0.36D2 t728 = 0.1D1 - t726 - t727 t730 = 0.6D1 * t728 * t62 t733 = t125 * t116 * (t79 * (t684 - t696 - 0.32D2) + t1 * (Km * (- #0.24D2 + t699) + t704 + Ki * (-0.24D2 + t705 - t706)) + Gm * (t714 # + Km * (t717 + Ki * (-0.12D2 + t718 - t719))) - t730) t734 = 0.2D1 / 0.3D1 * t695 t736 = t79 * (t734 - t702 + 0.16D2) t743 = t1 * (0.5D1 / 0.6D1 * t130 * t694 + t704 + Ki * (0.12D2 + 0 #.5D1 / 0.18D2 * t695 - t715)) t751 = Gm * (t714 + Km * (t717 + Ki * (0.6D1 + 0.5D1 / 0.24D2 * t6 #95 - 0.5D1 / 0.24D2 * t683))) t754 = t125 * t116 * (t736 + t743 + t751 - t730) t758 = Gm * (t302 + Km * (t25 - t189) - t9) t765 = 0.96D2 * (-Gm * t44 + t42 + t46) * (t155 + Gm) t767 = t48 * (f * (t187 + t301 + t758 - t416) + t765) t768 = t54 * t767 t774 = f * (t187 + t350 + Gm * (-Km * t538 - t154) - t63) t775 = t54 * t48 t776 = t775 * t774 t777 = 0.4D1 / 0.9D1 * t776 t779 = t79 * (0.8D1 / 0.9D1 * t768 - t777) t782 = 0.2D1 / 0.9D1 * t768 t783 = t782 + t777 t784 = Gi * t783 t785 = t768 - t776 t788 = t1 * (0.5D1 / 0.6D1 * t130 * t767 + t784 + 0.5D1 / 0.9D1 * #Ki * t785) t790 = t776 / 0.4D1 t792 = t61 * (t768 / 0.8D1 + t790) t794 = 0.5D1 / 0.9D1 * t776 t796 = Gi * (0.5D1 / 0.18D2 * t768 + t794) t801 = Gm * (t792 + Km * (t796 + 0.5D1 / 0.12D2 * Ki * t785)) t802 = t768 / 0.72D2 t803 = t776 / 0.36D2 t804 = -t802 - t803 t806 = 0.6D1 * t804 * t62 t809 = t125 * t116 * (t779 + t788 + t801 - t806) t810 = 0.2D1 / 0.3D1 * t776 t812 = t79 * (t810 - t782) t815 = 0.5D1 / 0.6D1 * Km * t775 * t774 t818 = t1 * (t815 + t784 - 0.5D1 / 0.18D2 * Ki * t785) t823 = Gm * (t792 + Km * (t796 - 0.5D1 / 0.24D2 * Ki * t785)) t826 = t125 * t116 * (t812 + t818 + t823 - t806) t830 = t48 * (f * (t187 + t301 + t758 - t310) + t765) t831 = t54 * t830 t832 = 0.2D1 / 0.9D1 * t831 t834 = t79 * (t810 - t832) t835 = t777 + t832 t836 = Gi * t835 t837 = t776 - t831 t840 = t1 * (t815 + t836 + 0.5D1 / 0.18D2 * Ki * t837) t843 = t61 * (t790 + t831 / 0.8D1) t846 = Gi * (t794 + 0.5D1 / 0.18D2 * t831) t851 = Gm * (t843 + Km * (t846 + 0.5D1 / 0.24D2 * Ki * t837)) t852 = t831 / 0.72D2 t853 = -t803 - t852 t855 = 0.6D1 * t853 * t62 t858 = t125 * t116 * (t834 + t840 + t851 - t855) t859 = 0.24D2 * t1 t860 = -t171 + t57 t862 = 0.12D2 * t31 t865 = 0.16D2 * t1 t867 = -0.18D2 * Km - t57 t872 = t48 * (t19 * (Gm * t860 - t859 + t862) + f * (Gm * t867 - t #862 - t865)) t873 = t79 * t54 t875 = 0.4D1 / 0.9D1 * t873 * t872 t876 = t130 * t872 t878 = Gi * t54 t879 = t878 * t872 t880 = 0.2D1 / 0.3D1 * t879 t882 = t1 * (0.5D1 / 0.6D1 * t876 + t880 + 0.8D1) t883 = t61 * t54 t884 = t883 * t872 t889 = 0.12D2 * Gi t891 = Gm * (0.3D1 / 0.8D1 * t884 + Km * (0.5D1 / 0.6D1 * t879 + 0 #.9D1) + t889) t892 = t876 / 0.24D2 t898 = t125 * t116 * (t875 + t882 + t891 - 0.6D1 * (-t892 - 0.1D1) # * t31) t899 = t1 * Mgi t900 = 0.24D2 * t899 t901 = Mgi * Gi t902 = 0.24D2 * t901 t903 = Km * Mgi t905 = t902 - 0.12D2 * t903 t908 = 0.12D2 * Gi * t903 t911 = 0.16D2 * t899 t913 = -t902 - 0.18D2 * t903 t918 = t48 * (t19 * (Gm * t905 - t900 + t908) + f * (Gm * t913 - t #908 - t911)) t920 = 0.4D1 / 0.9D1 * t873 * t918 t921 = t130 * t918 t923 = t878 * t918 t924 = 0.2D1 / 0.3D1 * t923 t925 = 0.8D1 * Mgi t927 = t1 * (0.5D1 / 0.6D1 * t921 + t924 + t925) t928 = t883 * t918 t931 = 0.9D1 * Mgi t934 = 0.12D2 * t901 t936 = Gm * (0.3D1 / 0.8D1 * t928 + Km * (0.5D1 / 0.6D1 * t923 + t #931) + t934) t937 = t921 / 0.24D2 t943 = t125 * t116 * (t920 + t927 + t936 - 0.6D1 * (-t937 - Mgi) * # t31) t951 = t48 * (t19 * (-Gm * t905 + t900 - t908) + f * (-Gm * t913 + # t908 + t911)) t953 = 0.4D1 / 0.9D1 * t873 * t951 t954 = t130 * t951 t956 = t878 * t951 t957 = 0.2D1 / 0.3D1 * t956 t959 = t1 * (0.5D1 / 0.6D1 * t954 + t957 - t925) t960 = t883 * t951 t966 = Gm * (0.3D1 / 0.8D1 * t960 + Km * (0.5D1 / 0.6D1 * t956 - t #931) - t934) t967 = t954 / 0.24D2 t973 = t125 * t116 * (t953 + t959 + t966 - 0.6D1 * (-t967 + Mgi) * # t31) t981 = t48 * (t19 * (-Gm * t860 + t859 - t862) + f * (-Gm * t867 + # t862 + t865)) t983 = 0.4D1 / 0.9D1 * t873 * t981 t984 = t130 * t981 t986 = t878 * t981 t987 = 0.2D1 / 0.3D1 * t986 t989 = t1 * (0.5D1 / 0.6D1 * t984 + t987 - 0.8D1) t990 = t883 * t981 t996 = Gm * (0.3D1 / 0.8D1 * t990 + Km * (0.5D1 / 0.6D1 * t986 - 0 #.9D1) - t889) t997 = t984 / 0.24D2 t1003 = t125 * t116 * (t983 + t989 + t996 - 0.6D1 * (-t997 + 0.1D1 #) * t31) t1004 = t1 * Mgm t1005 = 0.24D2 * t1004 t1006 = Mgm * Gi t1007 = 0.24D2 * t1006 t1008 = Km * Mgm t1010 = -t1007 + 0.12D2 * t1008 t1013 = 0.12D2 * Gi * t1008 t1016 = 0.16D2 * t1004 t1018 = t1007 + 0.18D2 * t1008 t1023 = t48 * (t19 * (Gm * t1010 + t1005 - t1013) + f * (Gm * t101 #8 + t1013 + t1016)) t1025 = 0.4D1 / 0.9D1 * t873 * t1023 t1026 = t130 * t1023 t1028 = t878 * t1023 t1029 = 0.2D1 / 0.3D1 * t1028 t1030 = 0.8D1 * Mgm t1032 = t1 * (0.5D1 / 0.6D1 * t1026 + t1029 - t1030) t1033 = t883 * t1023 t1036 = 0.9D1 * Mgm t1039 = 0.12D2 * t1006 t1041 = Gm * (0.3D1 / 0.8D1 * t1033 + Km * (0.5D1 / 0.6D1 * t1028 #- t1036) - t1039) t1042 = t1026 / 0.24D2 t1048 = t125 * t116 * (t1025 + t1032 + t1041 - 0.6D1 * (-t1042 + M #gm) * t31) t1056 = t48 * (t19 * (-Gm * t1010 - t1005 + t1013) + f * (-Gm * t1 #018 - t1013 - t1016)) t1058 = 0.4D1 / 0.9D1 * t873 * t1056 t1059 = t130 * t1056 t1061 = t878 * t1056 t1062 = 0.2D1 / 0.3D1 * t1061 t1064 = t1 * (0.5D1 / 0.6D1 * t1059 + t1062 + t1030) t1065 = t883 * t1056 t1071 = Gm * (0.3D1 / 0.8D1 * t1065 + Km * (0.5D1 / 0.6D1 * t1061 #+ t1036) + t1039) t1072 = t1059 / 0.24D2 t1078 = t125 * t116 * (t1058 + t1064 + t1071 - 0.6D1 * (-t1072 - M #gm) * t31) t1080 = t79 * (-t231 + t259) t1091 = t125 * t116 * (t1080 + t1 * (t273 - 0.5D1 / 0.18D2 * Ki * #t234 + t233) + Gm * (t241 + Km * (-0.5D1 / 0.24D2 * Ki * t234 + t2 #45)) - t255) t1107 = t125 * t116 * (t79 * (-t226 + 0.8D1 / 0.9D1 * t268) + t1 * # (0.5D1 / 0.6D1 * t130 * t267 - 0.5D1 / 0.9D1 * Ki * t276 + t275) #+ Gm * (t282 + Km * (-0.5D1 / 0.12D2 * Ki * t276 + t285)) - t294) t1109 = t79 * (-t382 + t414) t1123 = t125 * t116 * (t1109 + t1 * (t426 + Ki * (-t395 + t429 - 0 #.8D1) + t384) + Gm * (t394 + Km * (Ki * (-0.5D1 / 0.24D2 * t346 + #t440 - 0.9D1) + t397)) - t410) t1124 = 0.8D1 / 0.9D1 * t421 t1130 = 0.5D1 / 0.9D1 * t421 t1135 = 0.5D1 / 0.12D2 * t421 t1144 = t125 * t116 * (t79 * (-t376 + t1124 + 0.32D2) + t1 * (Km * # (0.5D1 / 0.6D1 * t421 + 0.36D2) + Ki * (-t386 + t1130 + 0.16D2) + # t428) + Gm * (t437 + Km * (Ki * (-t399 + t1135 + 0.18D2) + t439)) # - t451) t1146 = t79 * (-t553 + 0.16D2 + t585) t1160 = t125 * t116 * (t1146 + t1 * (t599 + Ki * (-t566 + 0.20D2 + # t602) + t555) + Gm * (t565 + Km * (Ki * (-0.5D1 / 0.24D2 * t508 + # 0.15D2 + t613) + t568)) - t581) t1181 = t125 * t116 * (t79 * (-t547 + 0.8D1 / 0.9D1 * t594 - 0.64D #2) + t1 * (Km * (0.5D1 / 0.6D1 * t594 - 0.60D2) + Ki * (-t557 + 0. #5D1 / 0.9D1 * t594 - 0.40D2) + t601) + Gm * (t610 + Km * (Ki * (-t #570 + 0.5D1 / 0.12D2 * t594 - 0.30D2) + t612)) - t624) t1202 = t125 * t116 * (t79 * (-t547 + 0.8D1 / 0.9D1 * t636 - 0.64D #2) + t1 * (Km * (-0.60D2 + 0.5D1 / 0.6D1 * t636) + t641 + Ki * (-0 #.40D2 - t557 + 0.5D1 / 0.9D1 * t636)) + Gm * (t649 + Km * (t651 + #Ki * (-0.30D2 - t570 + 0.5D1 / 0.12D2 * t636))) - t662) t1218 = t125 * t116 * (t79 * (-t777 + 0.8D1 / 0.9D1 * t831) + t1 * # (0.5D1 / 0.6D1 * t130 * t830 + t836 - 0.5D1 / 0.9D1 * Ki * t837) #+ Gm * (t843 + Km * (t846 - 0.5D1 / 0.12D2 * Ki * t837)) - t855) t1230 = t1 * (Gi * (-0.48D2 + 0.4D1 / 0.3D1 * t55 - t127) + Km * ( #t55 / 0.2D1 - t127) + Ki * (-0.16D2 + t55 / 0.18D2 + 0.17D2 / 0.18 #D2 * t76)) t1234 = t61 * (-t55 / 0.4D1 - t76 / 0.2D1) t1249 = Gm * (t1234 + Km * (-Gi * t85 + Ki * (-0.18D2 + 0.19D2 / 0 #.24D2 * t55 + t76 / 0.3D1)) - 0.36D2 * Gi * Ki * (0.5D1 / 0.3D1 - #t55 / 0.27D2 - t76 / 0.216D3)) t1253 = Km * (t55 / 0.108D3 + t76 / 0.54D2) t1262 = t55 / 0.3D1 t1264 = t79 * (t76 / 0.9D1 + t1262) t1275 = t1 * (Gi * (0.24D2 + t76 - t1262) + Km * (t76 / 0.6D1 - t1 #262) + Ki * (-0.16D2 + 0.19D2 / 0.36D2 * t76 + 0.17D2 / 0.36D2 * t #55)) t1291 = Gm * (t1234 + Km * (Gi * (-0.8D1 - t77 - t84) + Ki * (-0.1 #8D2 + 0.23D2 / 0.24D2 * t76 + t55 / 0.6D1)) - 0.36D2 * Gi * Ki * ( #0.1D1 / 0.6D1 - 0.17D2 / 0.432D3 * t76 - t55 / 0.432D3)) t1299 = t125 * t116 * (t1264 + t1275 + t1291 - 0.18D2 * (t1253 + K #i * (0.1D1 - t109 - t108)) * t31) t1307 = 0.17D2 / 0.18D2 * t225 t1313 = t225 / 0.2D1 t1315 = t61 * (-t202 / 0.4D1 - t1313) t1316 = -Gi * t232 t1318 = t225 / 0.3D1 t1324 = t225 / 0.216D3 t1332 = t225 / 0.54D2 t1342 = t225 / 0.9D1 t1343 = t268 / 0.3D1 t1345 = t79 * (t1342 + t1343) t1348 = t225 / 0.6D1 t1351 = 0.19D2 / 0.36D2 * t225 t1356 = t1 * (Gi * (t225 - t1343) + Km * (t1348 - t1343) + Ki * (t #1351 + 0.17D2 / 0.36D2 * t268)) t1359 = t61 * (-t1313 - t268 / 0.4D1) t1360 = -Gi * t274 t1361 = 0.23D2 / 0.24D2 * t225 t1367 = 0.17D2 / 0.432D3 * t225 t1374 = Gm * (t1359 + Km * (t1360 + Ki * (t1361 + t268 / 0.6D1)) - # 0.36D2 * Gi * Ki * (-t1367 - t268 / 0.432D3)) t1384 = t125 * t116 * (t1345 + t1356 + t1374 - 0.18D2 * (Km * (t13 #32 + t268 / 0.108D3) + Ki * t292) * t31) t1392 = 0.17D2 / 0.18D2 * t375 t1398 = t375 / 0.2D1 t1400 = t61 * (-t346 / 0.4D1 - t1398) t1403 = t375 / 0.3D1 t1409 = t375 / 0.216D3 t1417 = t375 / 0.54D2 t1427 = t375 / 0.9D1 t1428 = t421 / 0.3D1 t1430 = t79 * (t1427 + t1428) t1433 = t375 / 0.6D1 t1436 = 0.19D2 / 0.36D2 * t375 t1441 = t1 * (Gi * (t375 - t1428 - 0.24D2) + Km * (t1433 - t1428) #+ Ki * (t1436 + 0.17D2 / 0.36D2 * t421 + 0.16D2)) t1444 = t61 * (-t1398 - t421 / 0.4D1) t1447 = 0.23D2 / 0.24D2 * t375 t1453 = 0.17D2 / 0.432D3 * t375 t1460 = Gm * (t1444 + Km * (Gi * (-t376 - t422 + 0.8D1) + Ki * (0. #18D2 + t1447 + t421 / 0.6D1)) - 0.36D2 * Gi * Ki * (-0.1D1 / 0.6D1 # - t1453 - t421 / 0.432D3)) t1471 = t125 * t116 * (t1430 + t1441 + t1460 - 0.18D2 * (Km * (t14 #17 + t421 / 0.108D3) + Ki * (-t407 - t448 - 0.1D1)) * t31) t1479 = 0.17D2 / 0.18D2 * t546 t1485 = t546 / 0.2D1 t1487 = t61 * (0.18D2 - t508 / 0.4D1 - t1485) t1488 = -Gi * t554 t1490 = t546 / 0.3D1 t1496 = t546 / 0.216D3 t1504 = t546 / 0.54D2 t1514 = t546 / 0.9D1 t1515 = t594 / 0.3D1 t1517 = t79 * (t1514 + t1515 - 0.24D2) t1520 = t546 / 0.6D1 t1523 = 0.19D2 / 0.36D2 * t546 t1528 = t1 * (Gi * (t546 - t1515 + 0.24D2) + Km * (t1520 - t1515 + # 0.24D2) + Ki * (t1523 + 0.17D2 / 0.36D2 * t594 - 0.34D2)) t1531 = t61 * (0.18D2 - t1485 - t594 / 0.4D1) t1532 = -Gi * t600 t1533 = 0.23D2 / 0.24D2 * t546 t1539 = 0.17D2 / 0.432D3 * t546 t1546 = Gm * (t1531 + Km * (t1532 + Ki * (-0.12D2 + t1533 + t594 / # 0.6D1)) - 0.36D2 * Gi * Ki * (0.1D1 / 0.6D1 - t1539 - t594 / 0.43 #2D3)) t1557 = t636 / 0.3D1 t1559 = t79 * (t1514 + t1557 - 0.24D2) t1568 = t1 * (Gi * (t546 - t1557 + 0.24D2) + Km * (t1520 - t1557 + # 0.24D2) + Ki * (t1523 + 0.17D2 / 0.36D2 * t636 - 0.34D2)) t1571 = t61 * (0.18D2 - t1485 - t636 / 0.4D1) t1572 = -Gi * t640 t1584 = Gm * (t1571 + Km * (t1572 + Ki * (-0.12D2 + t1533 + t636 / # 0.6D1)) - 0.36D2 * Gi * Ki * (0.1D1 / 0.6D1 - t1539 - t636 / 0.43 #2D3)) t1606 = t1 * (Gi * (0.4D1 / 0.3D1 * t683 - 0.48D2 - t734) + Km * ( #t683 / 0.2D1 - 0.36D2 - t734) + Ki * (t683 / 0.18D2 + 0.12D2 + 0.1 #7D2 / 0.18D2 * t695)) t1610 = t61 * (0.18D2 - t683 / 0.4D1 - t695 / 0.2D1) t1625 = Gm * (t1610 + Km * (-Gi * t703 + Ki * (-0.39D2 + 0.19D2 / #0.24D2 * t683 + t695 / 0.3D1)) - 0.36D2 * Gi * Ki * (0.1D1 - t683 #/ 0.27D2 - t695 / 0.216D3)) t1629 = Km * (t683 / 0.108D3 - 0.2D1 / 0.3D1 + t695 / 0.54D2) t1638 = t683 / 0.3D1 t1640 = t79 * (t695 / 0.9D1 + t1638 - 0.24D2) t1651 = t1 * (Gi * (t695 - t1638) + Km * (t695 / 0.6D1 - t1638 + 0 #.24D2) + Ki * (0.19D2 / 0.36D2 * t695 + 0.17D2 / 0.36D2 * t683 - 0 #.18D2)) t1667 = Gm * (t1610 + Km * (Gi * (-t696 - t702 + 0.24D2) + Ki * (0 #.6D1 + 0.23D2 / 0.24D2 * t695 + t683 / 0.6D1)) - 0.36D2 * Gi * Ki #* (-0.17D2 / 0.432D3 * t695 - t683 / 0.432D3)) t1675 = t125 * t116 * (t1640 + t1651 + t1667 - 0.18D2 * (t1629 + K #i * (-t727 - t726)) * t31) t1683 = 0.17D2 / 0.18D2 * t776 t1687 = t1 * (Gi * (0.4D1 / 0.3D1 * t768 - t810) + Km * (t768 / 0. #2D1 - t810) + Ki * (t768 / 0.18D2 + t1683)) t1689 = t776 / 0.2D1 t1691 = t61 * (-t768 / 0.4D1 - t1689) t1692 = -Gi * t783 t1694 = t776 / 0.3D1 t1700 = t776 / 0.216D3 t1706 = Gm * (t1691 + Km * (t1692 + Ki * (0.19D2 / 0.24D2 * t768 + # t1694)) - 0.36D2 * Gi * Ki * (-t768 / 0.27D2 - t1700)) t1708 = t776 / 0.54D2 t1714 = 0.18D2 * (Km * (t768 / 0.108D3 + t1708) + Ki * t804) * t31 t1718 = t776 / 0.9D1 t1719 = t768 / 0.3D1 t1721 = t79 * (t1718 + t1719) t1724 = t776 / 0.6D1 t1727 = 0.19D2 / 0.36D2 * t776 t1732 = t1 * (Gi * (t776 - t1719) + Km * (t1724 - t1719) + Ki * (t #1727 + 0.17D2 / 0.36D2 * t768)) t1733 = 0.23D2 / 0.24D2 * t776 t1739 = 0.17D2 / 0.432D3 * t776 t1746 = Gm * (t1691 + Km * (t1692 + Ki * (t1733 + t768 / 0.6D1)) - # 0.36D2 * Gi * Ki * (-t1739 - t768 / 0.432D3)) t1750 = t831 / 0.3D1 t1752 = t79 * (t1718 + t1750) t1761 = t1 * (Gi * (t776 - t1750) + Km * (t1724 - t1750) + Ki * (t #1727 + 0.17D2 / 0.36D2 * t831)) t1764 = t61 * (-t1689 - t831 / 0.4D1) t1765 = -Gi * t835 t1777 = Gm * (t1764 + Km * (t1765 + Ki * (t1733 + t831 / 0.6D1)) - # 0.36D2 * Gi * Ki * (-t1739 - t831 / 0.432D3)) t1789 = Ki * t54 t1790 = t1789 * t872 t1792 = t1 * (t880 - t876 / 0.6D1 + t1790 - 0.16D2) t1797 = t1790 / 0.24D2 t1802 = Gm * (-0.3D1 / 0.4D1 * t884 + Km * (-t880 + 0.9D1 / 0.8D1 #* t1790 - 0.18D2) - 0.36D2 * Gi * (-t1797 + 0.2D1 / 0.3D1)) t1811 = t1789 * t918 t1812 = 0.16D2 * Mgi t1814 = t1 * (t924 - t921 / 0.6D1 + t1811 - t1812) t1817 = 0.18D2 * Mgi t1820 = t1811 / 0.24D2 t1821 = 0.2D1 / 0.3D1 * Mgi t1826 = Gm * (-0.3D1 / 0.4D1 * t928 + Km * (-t924 + 0.9D1 / 0.8D1 #* t1811 - t1817) - 0.36D2 * Gi * (-t1820 + t1821)) t1835 = t1789 * t951 t1837 = t1 * (t957 - t954 / 0.6D1 + t1835 + t1812) t1842 = t1835 / 0.24D2 t1847 = Gm * (-0.3D1 / 0.4D1 * t960 + Km * (-t957 + 0.9D1 / 0.8D1 #* t1835 + t1817) - 0.36D2 * Gi * (-t1842 - t1821)) t1856 = t1789 * t981 t1858 = t1 * (t987 - t984 / 0.6D1 + t1856 + 0.16D2) t1863 = t1856 / 0.24D2 t1868 = Gm * (-0.3D1 / 0.4D1 * t990 + Km * (-t987 + 0.9D1 / 0.8D1 #* t1856 + 0.18D2) - 0.36D2 * Gi * (-t1863 - 0.2D1 / 0.3D1)) t1877 = t1789 * t1023 t1878 = 0.16D2 * Mgm t1880 = t1 * (t1029 - t1026 / 0.6D1 + t1877 + t1878) t1883 = 0.18D2 * Mgm t1886 = t1877 / 0.24D2 t1887 = 0.2D1 / 0.3D1 * Mgm t1892 = Gm * (-0.3D1 / 0.4D1 * t1033 + Km * (-t1029 + 0.9D1 / 0.8D #1 * t1877 + t1883) - 0.36D2 * Gi * (-t1886 - t1887)) t1901 = t1789 * t1056 t1903 = t1 * (t1062 - t1059 / 0.6D1 + t1901 - t1878) t1908 = t1901 / 0.24D2 t1913 = Gm * (-0.3D1 / 0.4D1 * t1065 + Km * (-t1062 + 0.9D1 / 0.8D #1 * t1901 - t1883) - 0.36D2 * Gi * (-t1908 + t1887)) t1921 = Km * t110 t1922 = t55 / 0.48D2 t1923 = t76 / 0.24D2 t1932 = t125 * t116 * (t1264 + t1275 + t1291 + 0.12D2 * Gi * (t192 #1 + Ki * (-0.3D1 / 0.2D1 + t1922 + t1923)) * Km) t1941 = t125 * t116 * (t129 + t1230 + t1249 + 0.12D2 * Gi * (t1921 # + Ki * (t1923 + t1922)) * Km) t1942 = t202 / 0.3D1 t1968 = t225 / 0.24D2 t1977 = t125 * t116 * (t79 * (t1942 + t1342) + t1 * (Gi * (-t1942 #+ t225) + Km * (-t1942 + t1348) + Ki * (0.17D2 / 0.36D2 * t202 + t #1351)) + Gm * (t1315 + Km * (t1316 + Ki * (t202 / 0.6D1 + t1361)) #- 0.36D2 * Ki * (-t202 / 0.432D3 - t1367) * Gi) + 0.12D2 * Gi * (K #m * t253 + Ki * (t202 / 0.48D2 + t1968)) * Km) t2008 = 0.12D2 * Gi * (Km * t292 + Ki * (t1968 + t268 / 0.48D2)) * # Km t2011 = t125 * t116 * (t271 + t1 * (Gi * (-t259 + 0.4D1 / 0.3D1 * #t268) + Km * (-t259 + t268 / 0.2D1) + Ki * (t1307 + t268 / 0.18D2) #) + Gm * (t1359 + Km * (t1360 + Ki * (t1318 + 0.19D2 / 0.24D2 * t2 #68)) - 0.36D2 * Ki * (-t1324 - t268 / 0.27D2) * Gi) + t2008) t2014 = t125 * t116 * (t1345 + t1356 + t1374 + t2008) t2015 = t346 / 0.3D1 t2043 = t375 / 0.24D2 t2052 = t125 * t116 * (t79 * (t2015 + t1427) + t1 * (Gi * (-t2015 #+ t375 - 0.24D2) + Km * (-t2015 + t1433) + Ki * (0.17D2 / 0.36D2 * # t346 + t1436 + 0.16D2)) + Gm * (t1400 + Km * (Gi * (-t382 - t376 #+ 0.8D1) + Ki * (0.18D2 + t346 / 0.6D1 + t1447)) - 0.36D2 * Ki * ( #-0.1D1 / 0.6D1 - t346 / 0.432D3 - t1453) * Gi) + 0.12D2 * Gi * (Km # * t408 + Ki * (0.3D1 / 0.2D1 + t346 / 0.48D2 + t2043)) * Km) t2077 = Km * t449 t2078 = t421 / 0.48D2 t2087 = t125 * t116 * (t424 + t1 * (Gi * (-t414 + 0.4D1 / 0.3D1 * #t421 + 0.48D2) + Km * (-t414 + t421 / 0.2D1) + Ki * (t1392 + t421 #/ 0.18D2 + 0.16D2)) + Gm * (t1444 + Km * (-Gi * t427 + Ki * (0.18D #2 + t1403 + 0.19D2 / 0.24D2 * t421)) - 0.36D2 * Ki * (-0.5D1 / 0.3 #D1 - t1409 - t421 / 0.27D2) * Gi) + 0.12D2 * Gi * (t2077 + Ki * (t #2043 + t2078)) * Km) t2096 = t125 * t116 * (t1430 + t1441 + t1460 + 0.12D2 * Gi * (t207 #7 + Ki * (0.3D1 / 0.2D1 + t2043 + t2078)) * Km) t2097 = t508 / 0.3D1 t2123 = t546 / 0.24D2 t2132 = t125 * t116 * (t79 * (t2097 - 0.24D2 + t1514) + t1 * (Gi * # (-t2097 + 0.24D2 + t546) + Km * (-t2097 + 0.24D2 + t1520) + Ki * #(0.17D2 / 0.36D2 * t508 - 0.34D2 + t1523)) + Gm * (t1487 + Km * (t #1488 + Ki * (-0.12D2 + t508 / 0.6D1 + t1533)) - 0.36D2 * Ki * (0.1 #D1 / 0.6D1 - t508 / 0.432D3 - t1539) * Gi) + 0.12D2 * Gi * (Km * t #579 + Ki * (-0.3D1 / 0.2D1 + t508 / 0.48D2 + t2123)) * Km) t2163 = 0.12D2 * Gi * (Km * t622 + Ki * (-0.3D1 / 0.2D1 + t2123 + #t594 / 0.48D2)) * Km t2174 = 0.12D2 * Gi * (Km * t660 + Ki * (-0.3D1 / 0.2D1 + t2123 + #t636 / 0.48D2)) * Km t2178 = Km * t728 t2179 = t683 / 0.48D2 t2180 = t695 / 0.24D2 t2189 = t125 * t116 * (t1640 + t1651 + t1667 + 0.12D2 * Gi * (t217 #8 + Ki * (t2179 + t2180)) * Km) t2198 = t125 * t116 * (t736 + t1606 + t1625 + 0.12D2 * Gi * (t2178 # + Ki * (-0.3D1 / 0.2D1 + t2180 + t2179)) * Km) t2201 = t776 / 0.24D2 t2207 = 0.12D2 * Gi * (Km * t804 + Ki * (t768 / 0.48D2 + t2201)) * # Km t2210 = t125 * t116 * (t1721 + t1732 + t1746 + t2207) t2221 = 0.12D2 * Gi * (Km * t853 + Ki * (t2201 + t831 / 0.48D2)) * # Km t2232 = t125 * t116 * (t875 + t1792 + t1802 + 0.12D2 * Gi * (-t892 # + t1790 / 0.16D2 - 0.1D1) * Km) t2240 = t125 * t116 * (t920 + t1814 + t1826 + 0.12D2 * Gi * (-t937 # + t1811 / 0.16D2 - Mgi) * Km) t2248 = t125 * t116 * (t953 + t1837 + t1847 + 0.12D2 * Gi * (-t967 # + t1835 / 0.16D2 + Mgi) * Km) t2256 = t125 * t116 * (t983 + t1858 + t1868 + 0.12D2 * Gi * (-t997 # + t1856 / 0.16D2 + 0.1D1) * Km) t2264 = t125 * t116 * (t1025 + t1880 + t1892 + 0.12D2 * Gi * (-t10 #42 + t1877 / 0.16D2 + Mgm) * Km) t2272 = t125 * t116 * (t1058 + t1903 + t1913 + 0.12D2 * Gi * (-t10 #72 + t1901 / 0.16D2 - Mgm) * Km) t2329 = t79 * (t56 - t77) t2337 = t1 * (0.5D1 / 0.6D1 * t130 * t49 + Gi * (0.32D2 + t84 + t7 #7) + Ki * (0.8D1 + t87 - t88)) t2345 = Gm * (t96 + Km * (Gi * (0.16D2 + t97 + t88) + Ki * (0.9D1 #+ t100 - t101)) + t159) t2348 = -t55 / 0.216D3 - t76 / 0.108D3 t2362 = t125 * t116 * (-0.18D2 * Gi * t61 * t2348 + t129 + t137 + #t145) t2364 = t225 / 0.108D3 t2379 = t125 * t116 * (t271 + t279 + t290 - 0.18D2 * Gi * t61 * (- #t2364 - t268 / 0.216D3)) t2399 = t375 / 0.108D3 t2416 = t125 * t116 * (t424 + t434 + t447 - 0.18D2 * Gi * t61 * (- #t2399 - t421 / 0.216D3)) t2418 = t546 / 0.108D3 t2443 = t79 * (t684 - 0.64D2 - t696) t2451 = t1 * (Km * (-0.60D2 + t699) + Gi * (t702 - 0.48D2 + t696) #+ Ki * (t705 - 0.48D2 - t706)) t2459 = Gm * (t714 + Km * (Gi * (t715 - 0.36D2 + t706) + Ki * (t71 #8 - 0.39D2 - t719)) - t159) t2462 = 0.1D1 / 0.3D1 - t683 / 0.216D3 - t695 / 0.108D3 t2476 = t125 * t116 * (-0.18D2 * Gi * t61 * t2462 + t736 + t743 + #t751) t2478 = t776 / 0.108D3 t2482 = 0.18D2 * Gi * t61 * (-t768 / 0.216D3 - t2478) t2506 = Mgi / 0.3D1 t2531 = Mgm / 0.3D1 t2553 = t125 * t116 * (t2329 + t2337 + t2345 - 0.6D1 * Gi * (t1921 # - t13) * Km) t2578 = t125 * t116 * (t79 * (-t376 + t1124) + t1 * (0.5D1 / 0.6D1 # * t130 * t420 + Gi * (t376 + t422 - 0.32D2) + Ki * (-t386 + t1130 # - 0.8D1)) + Gm * (t437 + Km * (Gi * (t386 + t430 - 0.16D2) + Ki * # (-t399 + t1135 - 0.9D1)) - t159) - 0.6D1 * Gi * (t2077 + t13) * K #m) t2585 = t125 * t116 * (t2443 + t2451 + t2459 - 0.6D1 * Gi * (t2178 # + t13) * Km) t2586 = t76 / 0.72D2 t2587 = t225 / 0.72D2 t2588 = t375 / 0.72D2 t2590 = t546 / 0.72D2 t2591 = t726 - 0.1D1 t2592 = t695 / 0.72D2 t2593 = t776 / 0.72D2 t2595 = t54 * t872 / 0.72D2 t2597 = t54 * t918 / 0.72D2 t2599 = t54 * t951 / 0.72D2 t2601 = t54 * t981 / 0.72D2 t2603 = t54 * t1023 / 0.72D2 t2605 = t54 * t1056 / 0.72D2 JEA(1,1) = t126 JEA(1,2) = t148 JEA(1,3) = t148 JEA(1,4) = t125 * t116 * (t228 + t237 + t250 - t255) JEA(1,5) = t297 JEA(1,6) = t297 JEA(1,7) = t125 * t116 * (t79 * (t347 - t376 + 0.32D2) + t1 * (Km #* (0.36D2 + 0.5D1 / 0.6D1 * t346) + t384 + Ki * (0.16D2 + t385 - t #386)) + Gm * (t394 + Km * (t397 + Ki * (0.18D2 + t398 - t399))) - #t410) JEA(1,8) = t454 JEA(1,9) = t454 JEA(1,10) = t125 * t116 * (t549 + t561 + t576 - t581) JEA(1,11) = t627 JEA(1,12) = t665 JEA(1,13) = t733 JEA(1,14) = t754 JEA(1,15) = t754 JEA(1,16) = t809 JEA(1,17) = t826 JEA(1,18) = t858 JEA(1,19) = t898 JEA(1,20) = t943 JEA(1,21) = t973 JEA(1,22) = t1003 JEA(1,23) = t1048 JEA(1,24) = t1078 JEA(2,1) = t148 JEA(2,2) = t126 JEA(2,3) = t148 JEA(2,4) = t1091 JEA(2,5) = t1107 JEA(2,6) = t297 JEA(2,7) = t1123 JEA(2,8) = t1144 JEA(2,9) = t454 JEA(2,10) = t1160 JEA(2,11) = t1181 JEA(2,12) = t665 JEA(2,13) = t754 JEA(2,14) = t733 JEA(2,15) = t754 JEA(2,16) = t826 JEA(2,17) = t809 JEA(2,18) = t858 JEA(2,19) = t898 JEA(2,20) = t943 JEA(2,21) = t973 JEA(2,22) = t1003 JEA(2,23) = t1048 JEA(2,24) = t1078 JEA(3,1) = t148 JEA(3,2) = t148 JEA(3,3) = t126 JEA(3,4) = t1091 JEA(3,5) = t297 JEA(3,6) = t1107 JEA(3,7) = t1123 JEA(3,8) = t454 JEA(3,9) = t1144 JEA(3,10) = t1160 JEA(3,11) = t627 JEA(3,12) = t1202 JEA(3,13) = t754 JEA(3,14) = t754 JEA(3,15) = t733 JEA(3,16) = t826 JEA(3,17) = t826 JEA(3,18) = t1218 JEA(3,19) = t898 JEA(3,20) = t943 JEA(3,21) = t973 JEA(3,22) = t1003 JEA(3,23) = t1048 JEA(3,24) = t1078 JEA(4,1) = t125 * t116 * (t129 + t1230 + t1249 - 0.18D2 * (Ki * t1 #10 + t1253) * t31) JEA(4,2) = t1299 JEA(4,3) = t1299 JEA(4,4) = t125 * t116 * (t1080 + t1 * (Gi * (0.4D1 / 0.3D1 * t202 # - t259) + Km * (t202 / 0.2D1 - t259) + Ki * (t202 / 0.18D2 + t130 #7)) + Gm * (t1315 + Km * (t1316 + Ki * (0.19D2 / 0.24D2 * t202 + t #1318)) - 0.36D2 * Gi * Ki * (-t202 / 0.27D2 - t1324)) - 0.18D2 * ( #Km * (t202 / 0.108D3 + t1332) + Ki * t253) * t31) JEA(4,5) = t1384 JEA(4,6) = t1384 JEA(4,7) = t125 * t116 * (t1109 + t1 * (Gi * (0.4D1 / 0.3D1 * t346 # - t414 + 0.48D2) + Km * (t346 / 0.2D1 - t414) + Ki * (t346 / 0.18 #D2 + t1392 + 0.16D2)) + Gm * (t1400 + Km * (-Gi * t383 + Ki * (0.1 #8D2 + 0.19D2 / 0.24D2 * t346 + t1403)) - 0.36D2 * Gi * Ki * (-0.5D #1 / 0.3D1 - t346 / 0.27D2 - t1409)) - 0.18D2 * (Km * (t346 / 0.108 #D3 + t1417) + Ki * t408) * t31) JEA(4,8) = t1471 JEA(4,9) = t1471 JEA(4,10) = t125 * t116 * (t1146 + t1 * (Gi * (0.4D1 / 0.3D1 * t50 #8 - 0.96D2 - t585) + Km * (t508 / 0.2D1 - 0.36D2 - t585) + Ki * (t #508 / 0.18D2 - 0.4D1 + t1479)) + Gm * (t1487 + Km * (t1488 + Ki * #(-0.57D2 + 0.19D2 / 0.24D2 * t508 + t1490)) - 0.36D2 * Gi * Ki * ( #0.8D1 / 0.3D1 - t508 / 0.27D2 - t1496)) - 0.18D2 * (Km * (t508 / 0 #.108D3 - 0.2D1 / 0.3D1 + t1504) + Ki * t579) * t31) JEA(4,11) = t125 * t116 * (t1517 + t1528 + t1546 - 0.18D2 * (Km * #(t1504 + t594 / 0.108D3 - 0.2D1 / 0.3D1) + Ki * t622) * t31) JEA(4,12) = t125 * t116 * (t1559 + t1568 + t1584 - 0.18D2 * (Km * #(t1504 + t636 / 0.108D3 - 0.2D1 / 0.3D1) + Ki * t660) * t31) JEA(4,13) = t125 * t116 * (t736 + t1606 + t1625 - 0.18D2 * (Ki * t #728 + t1629) * t31) JEA(4,14) = t1675 JEA(4,15) = t1675 JEA(4,16) = t125 * t116 * (t812 + t1687 + t1706 - t1714) JEA(4,17) = t125 * t116 * (t1721 + t1732 + t1746 - t1714) JEA(4,18) = t125 * t116 * (t1752 + t1761 + t1777 - 0.18D2 * (Km * #(t1708 + t831 / 0.108D3) + Ki * t853) * t31) JEA(4,19) = t125 * t116 * (t875 + t1792 + t1802 - 0.18D2 * (t876 / # 0.36D2 - t1797 + 0.2D1 / 0.3D1) * t31) JEA(4,20) = t125 * t116 * (t920 + t1814 + t1826 - 0.18D2 * (t921 / # 0.36D2 - t1820 + t1821) * t31) JEA(4,21) = t125 * t116 * (t953 + t1837 + t1847 - 0.18D2 * (t954 / # 0.36D2 - t1842 - t1821) * t31) JEA(4,22) = t125 * t116 * (t983 + t1858 + t1868 - 0.18D2 * (t984 / # 0.36D2 - t1863 - 0.2D1 / 0.3D1) * t31) JEA(4,23) = t125 * t116 * (t1025 + t1880 + t1892 - 0.18D2 * (t1026 # / 0.36D2 - t1886 - t1887) * t31) JEA(4,24) = t125 * t116 * (t1058 + t1903 + t1913 - 0.18D2 * (t1059 # / 0.36D2 - t1908 + t1887) * t31) JEA(5,1) = t1932 JEA(5,2) = t1941 JEA(5,3) = t1932 JEA(5,4) = t1977 JEA(5,5) = t2011 JEA(5,6) = t2014 JEA(5,7) = t2052 JEA(5,8) = t2087 JEA(5,9) = t2096 JEA(5,10) = t2132 JEA(5,11) = t125 * t116 * (t597 + t1 * (Gi * (-t585 + 0.4D1 / 0.3D #1 * t594 - 0.96D2) + Km * (-t585 + t594 / 0.2D1 - 0.36D2) + Ki * ( #t1479 + t594 / 0.18D2 - 0.4D1)) + Gm * (t1531 + Km * (t1532 + Ki * # (-0.57D2 + t1490 + 0.19D2 / 0.24D2 * t594)) - 0.36D2 * Ki * (0.8D #1 / 0.3D1 - t1496 - t594 / 0.27D2) * Gi) + t2163) JEA(5,12) = t125 * t116 * (t1559 + t1568 + t1584 + t2174) JEA(5,13) = t2189 JEA(5,14) = t2198 JEA(5,15) = t2189 JEA(5,16) = t2210 JEA(5,17) = t125 * t116 * (t812 + t1687 + t1706 + t2207) JEA(5,18) = t125 * t116 * (t1752 + t1761 + t1777 + t2221) JEA(5,19) = t2232 JEA(5,20) = t2240 JEA(5,21) = t2248 JEA(5,22) = t2256 JEA(5,23) = t2264 JEA(5,24) = t2272 JEA(6,1) = t1932 JEA(6,2) = t1932 JEA(6,3) = t1941 JEA(6,4) = t1977 JEA(6,5) = t2014 JEA(6,6) = t2011 JEA(6,7) = t2052 JEA(6,8) = t2096 JEA(6,9) = t2087 JEA(6,10) = t2132 JEA(6,11) = t125 * t116 * (t1517 + t1528 + t1546 + t2163) JEA(6,12) = t125 * t116 * (t639 + t1 * (Gi * (-t585 + 0.4D1 / 0.3D #1 * t636 - 0.96D2) + Km * (-t585 + t636 / 0.2D1 - 0.36D2) + Ki * ( #t1479 + t636 / 0.18D2 - 0.4D1)) + Gm * (t1571 + Km * (t1572 + Ki * # (-0.57D2 + t1490 + 0.19D2 / 0.24D2 * t636)) - 0.36D2 * Ki * (0.8D #1 / 0.3D1 - t1496 - t636 / 0.27D2) * Gi) + t2174) JEA(6,13) = t2189 JEA(6,14) = t2189 JEA(6,15) = t2198 JEA(6,16) = t2210 JEA(6,17) = t2210 JEA(6,18) = t125 * t116 * (t834 + t1 * (Gi * (-t810 + 0.4D1 / 0.3D #1 * t831) + Km * (-t810 + t831 / 0.2D1) + Ki * (t1683 + t831 / 0.1 #8D2)) + Gm * (t1764 + Km * (t1765 + Ki * (t1694 + 0.19D2 / 0.24D2 #* t831)) - 0.36D2 * Ki * (-t1700 - t831 / 0.27D2) * Gi) + t2221) JEA(6,19) = t2232 JEA(6,20) = t2240 JEA(6,21) = t2248 JEA(6,22) = t2256 JEA(6,23) = t2264 JEA(6,24) = t2272 JEA(7,1) = t125 * t116 * (t2329 + t2337 + t2345 - 0.18D2 * Gi * Km # * (Km * t2348 - Ki)) JEA(7,2) = t2362 JEA(7,3) = t2362 JEA(7,4) = t125 * t116 * (t228 + t237 + t250 - 0.18D2 * Gi * t61 * # (-t202 / 0.216D3 - t2364)) JEA(7,5) = t2379 JEA(7,6) = t2379 JEA(7,7) = t125 * t116 * (t79 * (t347 - t376) + t1 * (0.5D1 / 0.6D #1 * t130 * t345 + Gi * (t382 + t376 - 0.32D2) + Ki * (t385 - t386 #- 0.8D1)) + Gm * (t394 + Km * (Gi * (t395 + t386 - 0.16D2) + Ki * #(t398 - t399 - 0.9D1)) - t159) - 0.18D2 * Gi * Km * (Km * (-t346 / # 0.216D3 - t2399) + Ki)) JEA(7,8) = t2416 JEA(7,9) = t2416 JEA(7,10) = t125 * t116 * (t549 + t561 + t576 - 0.18D2 * Gi * t61 #* (0.1D1 / 0.3D1 - t508 / 0.216D3 - t2418)) JEA(7,11) = t125 * t116 * (t597 + t607 + t620 - 0.18D2 * Gi * t61 #* (0.1D1 / 0.3D1 - t2418 - t594 / 0.216D3)) JEA(7,12) = t125 * t116 * (t639 + t646 + t658 - 0.18D2 * Gi * t61 #* (0.1D1 / 0.3D1 - t2418 - t636 / 0.216D3)) JEA(7,13) = t125 * t116 * (t2443 + t2451 + t2459 - 0.18D2 * Gi * K #m * (Km * t2462 + Ki)) JEA(7,14) = t2476 JEA(7,15) = t2476 JEA(7,16) = t125 * t116 * (t779 + t788 + t801 - t2482) JEA(7,17) = t125 * t116 * (t812 + t818 + t823 - t2482) JEA(7,18) = t125 * t116 * (t834 + t840 + t851 - 0.18D2 * Gi * t61 #* (-t2478 - t831 / 0.216D3)) JEA(7,19) = t125 * t116 * (t875 + t882 + t891 - 0.18D2 * Gi * Km * # (-t876 / 0.72D2 - 0.1D1 / 0.3D1)) JEA(7,20) = t125 * t116 * (t920 + t927 + t936 - 0.18D2 * Gi * Km * # (-t921 / 0.72D2 - t2506)) JEA(7,21) = t125 * t116 * (t953 + t959 + t966 - 0.18D2 * Gi * Km * # (-t954 / 0.72D2 + t2506)) JEA(7,22) = t125 * t116 * (t983 + t989 + t996 - 0.18D2 * Gi * Km * # (-t984 / 0.72D2 + 0.1D1 / 0.3D1)) JEA(7,23) = t125 * t116 * (t1025 + t1032 + t1041 - 0.18D2 * Gi * K #m * (-t1026 / 0.72D2 + t2531)) JEA(7,24) = t125 * t116 * (t1058 + t1064 + t1071 - 0.18D2 * Gi * K #m * (-t1059 / 0.72D2 - t2531)) JEA(8,1) = t148 JEA(8,2) = t2553 JEA(8,3) = t148 JEA(8,4) = t1091 JEA(8,5) = t1107 JEA(8,6) = t297 JEA(8,7) = t1123 JEA(8,8) = t2578 JEA(8,9) = t454 JEA(8,10) = t1160 JEA(8,11) = t1181 JEA(8,12) = t665 JEA(8,13) = t754 JEA(8,14) = t2585 JEA(8,15) = t754 JEA(8,16) = t826 JEA(8,17) = t809 JEA(8,18) = t858 JEA(8,19) = t898 JEA(8,20) = t943 JEA(8,21) = t973 JEA(8,22) = t1003 JEA(8,23) = t1048 JEA(8,24) = t1078 JEA(9,1) = t148 JEA(9,2) = t148 JEA(9,3) = t2553 JEA(9,4) = t1091 JEA(9,5) = t297 JEA(9,6) = t1107 JEA(9,7) = t1123 JEA(9,8) = t454 JEA(9,9) = t2578 JEA(9,10) = t1160 JEA(9,11) = t627 JEA(9,12) = t1202 JEA(9,13) = t754 JEA(9,14) = t754 JEA(9,15) = t2585 JEA(9,16) = t826 JEA(9,17) = t826 JEA(9,18) = t1218 JEA(9,19) = t898 JEA(9,20) = t943 JEA(9,21) = t973 JEA(9,22) = t1003 JEA(9,23) = t1048 JEA(9,24) = t1078 JEA(10,1) = t108 JEA(10,2) = t2586 JEA(10,3) = t2586 JEA(10,4) = t251 JEA(10,5) = t2587 JEA(10,6) = t2587 JEA(10,7) = t406 JEA(10,8) = t2588 JEA(10,9) = t2588 JEA(10,10) = t577 - 0.1D1 JEA(10,11) = t2590 JEA(10,12) = t2590 JEA(10,13) = t2591 JEA(10,14) = t2592 JEA(10,15) = t2592 JEA(10,16) = t802 JEA(10,17) = t2593 JEA(10,18) = t2593 JEA(10,19) = t2595 JEA(10,20) = t2597 JEA(10,21) = t2599 JEA(10,22) = t2601 JEA(10,23) = t2603 JEA(10,24) = t2605 JEA(11,1) = t2586 JEA(11,2) = t108 JEA(11,3) = t2586 JEA(11,4) = t2587 JEA(11,5) = t291 JEA(11,6) = t2587 JEA(11,7) = t2588 JEA(11,8) = t448 JEA(11,9) = t2588 JEA(11,10) = t2590 JEA(11,11) = t621 - 0.1D1 JEA(11,12) = t2590 JEA(11,13) = t2592 JEA(11,14) = t2591 JEA(11,15) = t2592 JEA(11,16) = t2593 JEA(11,17) = t802 JEA(11,18) = t2593 JEA(11,19) = t2595 JEA(11,20) = t2597 JEA(11,21) = t2599 JEA(11,22) = t2601 JEA(11,23) = t2603 JEA(11,24) = t2605 JEA(12,1) = t2586 JEA(12,2) = t2586 JEA(12,3) = t108 JEA(12,4) = t2587 JEA(12,5) = t2587 JEA(12,6) = t291 JEA(12,7) = t2588 JEA(12,8) = t2588 JEA(12,9) = t448 JEA(12,10) = t2590 JEA(12,11) = t2590 JEA(12,12) = t659 - 0.1D1 JEA(12,13) = t2592 JEA(12,14) = t2592 JEA(12,15) = t2591 JEA(12,16) = t2593 JEA(12,17) = t2593 JEA(12,18) = t852 JEA(12,19) = t2595 JEA(12,20) = t2597 JEA(12,21) = t2599 JEA(12,22) = t2601 JEA(12,23) = t2603 JEA(12,24) = t2605 if(affiche) then print*,'Dans JEA1inc3D matrice JEA' do i=1,12 write(*,'(24(E10.3))') (JEA(i,j),j=1,24) end do end if else print*,'Pb de dimansion dans jelc3d' err1=1 return end if return end