Sources Esope | Cast3M

jelc3d
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