C
C There are a total of 42 entries in the algebraic variable array.
C There are a total of 10 entries in each of the rate and state variable arrays.
C There are a total of 47 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is Gi in component Gi (millimolar).
C ALGBRC(5) is Jglut in component Jglut (millimolar_per_millisecond).
C ALGBRC(7) is Jgk in component Jgk (millimolar_per_millisecond).
C STATES(2) is Ge in component Ge (millimolar).
C ALGBRC(2) is G_infinity in component Ge (millimolar).
C CONSTS(1) is tau_G in component Ge (millisecond).
C CONSTS(2) is Gmax in component Ge (millimolar).
C CONSTS(3) is Gmin in component Ge (millimolar).
C CONSTS(4) is SG in component Ge (dimensionless).
C CONSTS(5) is I_ in component Ge (dimensionless).
C STATES(3) is I in component I (dimensionless).
C STATES(4) is G6P in component G6P (millimolar).
C CONSTS(6) is kappa in component model_parameters (dimensionless).
C ALGBRC(42) is JPFK in component JPFK (micromolar_per_millisecond).
C STATES(5) is FBP in component FBP (micromolar).
C ALGBRC(8) is JGPDH in component JGPDH (micromolar_per_millisecond).
C ALGBRC(1) is F6P in component F6P (millimolar).
C CONSTS(7) is Kglut in component Jglut (millimolar).
C CONSTS(8) is Vglut in component Jglut (millimolar_per_millisecond).
C CONSTS(9) is Kgk in component Jgk (millimolar).
C CONSTS(10) is Vgk in component Jgk (millimolar_per_millisecond).
C CONSTS(11) is ngk in component Jgk (dimensionless).
C CONSTS(12) is lambda in component JPFK (dimensionless).
C CONSTS(13) is Vmax in component JPFK (micromolar_per_millisecond).
C CONSTS(47) is w0000 in component w (dimensionless).
C ALGBRC(26) is w1000 in component w (dimensionless).
C ALGBRC(9) is w0100 in component w (dimensionless).
C ALGBRC(10) is w0010 in component w (dimensionless).
C ALGBRC(28) is w0001 in component w (dimensionless).
C ALGBRC(30) is w1100 in component w (dimensionless).
C ALGBRC(32) is w1010 in component w (dimensionless).
C ALGBRC(34) is w1001 in component w (dimensionless).
C ALGBRC(35) is w0101 in component w (dimensionless).
C ALGBRC(36) is w0011 in component w (dimensionless).
C ALGBRC(11) is w0110 in component w (dimensionless).
C ALGBRC(37) is w1110 in component w (dimensionless).
C ALGBRC(38) is w0111 in component w (dimensionless).
C ALGBRC(39) is w1011 in component w (dimensionless).
C ALGBRC(40) is w1101 in component w (dimensionless).
C ALGBRC(41) is w1111 in component w (dimensionless).
C CONSTS(14) is f13 in component w (dimensionless).
C CONSTS(15) is f23 in component w (dimensionless).
C CONSTS(16) is f41 in component w (dimensionless).
C CONSTS(17) is f42 in component w (dimensionless).
C CONSTS(18) is f43 in component w (dimensionless).
C CONSTS(19) is K1 in component w (micromolar).
C CONSTS(20) is K2 in component w (micromolar).
C CONSTS(21) is K3 in component w (micromolar).
C CONSTS(22) is K4 in component w (micromolar).
C ALGBRC(24) is AMP in component AMP (micromolar).
C ALGBRC(25) is ATP in component ATP (micromolar).
C CONSTS(23) is Atot in component ATP (micromolar).
C STATES(6) is ADP in component ADP (micromolar).
C CONSTS(24) is tau_a in component ADP (millisecond).
C CONSTS(25) is r in component ADP (dimensionless).
C CONSTS(26) is r1 in component ADP (micromolar).
C ALGBRC(12) is gamma in component gamma (dimensionless).
C STATES(7) is Ca in component Ca (micromolar).
C CONSTS(27) is v_gamma in component gamma (dimensionless).
C CONSTS(28) is k_gamma in component gamma (micromolar_per_millisecond).
C STATES(8) is v in component membrane (millivolt).
C CONSTS(29) is cm in component membrane (femtofarad).
C ALGBRC(15) is I_Ca in component I_Ca (picoampere).
C ALGBRC(13) is I_K in component I_K (picoampere).
C ALGBRC(18) is I_K_Ca in component I_K_Ca (picoampere).
C ALGBRC(33) is I_K_ATP in component I_K_ATP (picoampere).
C CONSTS(30) is gK_ in component I_K (picosiemens).
C CONSTS(31) is vK in component model_parameters (millivolt).
C STATES(9) is n in component n (dimensionless).
C ALGBRC(6) is n_infinity in component n (dimensionless).
C ALGBRC(3) is tau_n in component n (millisecond).
C CONSTS(32) is gCa_ in component I_Ca (picosiemens).
C CONSTS(33) is vCa in component model_parameters (millivolt).
C ALGBRC(14) is m_infinity in component m (dimensionless).
C CONSTS(34) is gkCa_ in component I_K_Ca (picosiemens).
C ALGBRC(16) is gkCa in component I_K_Ca (picosiemens).
C CONSTS(35) is KD in component I_K_Ca (micromolar).
C CONSTS(36) is gkATP_ in component I_K_ATP (picosiemens).
C ALGBRC(31) is gkATP in component I_K_ATP (picosiemens).
C ALGBRC(29) is o_infinity in component I_K_ATP (dimensionless).
C ALGBRC(20) is MgADP in component I_K_ATP (micromolar).
C ALGBRC(22) is ADP3_ in component I_K_ATP (micromolar).
C ALGBRC(27) is ATP4_ in component I_K_ATP (micromolar).
C CONSTS(37) is fcyt in component Ca (dimensionless).
C ALGBRC(17) is Jmem in component Jmem (micromolar_per_millisecond).
C ALGBRC(23) is Jer in component Jer (micromolar_per_millisecond).
C STATES(10) is Caer in component Caer (micromolar).
C CONSTS(38) is fer in component Caer (dimensionless).
C CONSTS(39) is Vcyt_Ver in component Caer (dimensionless).
C CONSTS(40) is kPMCA in component Jmem (first_order_rate_constant).
C CONSTS(41) is alpha in component Jmem (micromolar_per_millisecond).
C ALGBRC(21) is Jleak in component Jleak (micromolar_per_millisecond).
C ALGBRC(19) is JSERCA in component JSERCA (micromolar_per_millisecond).
C CONSTS(42) is kSERCA in component JSERCA (first_order_rate_constant).
C CONSTS(43) is pleak in component Jleak (first_order_rate_constant).
C ALGBRC(4) is I_infinity in component I (dimensionless).
C CONSTS(44) is tau_I in component I (millisecond).
C CONSTS(45) is I_slope in component I (per_micromolar).
C CONSTS(46) is Canull in component I (micromolar).
C RATES(1) is d/dt Gi in component Gi (millimolar).
C RATES(2) is d/dt Ge in component Ge (millimolar).
C RATES(4) is d/dt G6P in component G6P (millimolar).
C RATES(5) is d/dt FBP in component FBP (micromolar).
C RATES(6) is d/dt ADP in component ADP (micromolar).
C RATES(8) is d/dt v in component membrane (millivolt).
C RATES(9) is d/dt n in component n (dimensionless).
C RATES(7) is d/dt Ca in component Ca (micromolar).
C RATES(10) is d/dt Caer in component Caer (micromolar).
C RATES(3) is d/dt I in component I (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0
      STATES(2) = 7
      CONSTS(1) = 450000
      CONSTS(2) = 15
      CONSTS(3) = 1
      CONSTS(4) = 1
      CONSTS(5) = 5
      STATES(3) = 0
      STATES(4) = 200
      CONSTS(6) = 0.005
      STATES(5) = 40
      CONSTS(7) = 7
      CONSTS(8) = 8
      CONSTS(9) = 7
      CONSTS(10) = 0.8
      CONSTS(11) = 4
      CONSTS(12) = 0.06
      CONSTS(13) = 2
      CONSTS(14) = 0.02
      CONSTS(15) = 0.2
      CONSTS(16) = 20
      CONSTS(17) = 20
      CONSTS(18) = 20
      CONSTS(19) = 30
      CONSTS(20) = 1
      CONSTS(21) = 50000
      CONSTS(22) = 1000
      CONSTS(23) = 3000
      STATES(6) = 780
      CONSTS(24) = 300000
      CONSTS(25) = 0.5
      CONSTS(26) = 0.35
      STATES(7) = 0.1
      CONSTS(27) = 2.2
      CONSTS(28) = 0.1
      STATES(8) = -60
      CONSTS(29) = 5300
      CONSTS(30) = 2700
      CONSTS(31) = -75
      STATES(9) = 0
      CONSTS(32) = 1000
      CONSTS(33) = 25
      CONSTS(34) = 400
      CONSTS(35) = 0.5
      CONSTS(36) = 40000
      CONSTS(37) = 0.01
      STATES(10) = 185
      CONSTS(38) = 0.01
      CONSTS(39) = 31
      CONSTS(40) = 0.18
      CONSTS(41) = 4.5e-6
      CONSTS(42) = 0.4
      CONSTS(43) = 0.0002
      CONSTS(44) = 10000
      CONSTS(45) = 210
      CONSTS(46) = 0.055
      CONSTS(47) = 1.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = CONSTS(3)+(CONSTS(2) - CONSTS(3))/(1.00000+EXP((STATES(3) - CONSTS(5))/CONSTS(4)))
      RATES(2) = (ALGBRC(2) - STATES(2))/CONSTS(1)
      ALGBRC(4) = TERNRY(STATES(7).GE.CONSTS(46),  CONSTS(45)*(STATES(7) - CONSTS(46)), TERNRY(STATES(7).LT.CONSTS(46), 0.00000, 0.0/0.0)
      RATES(3) = (ALGBRC(4) - STATES(3))/CONSTS(44)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(8)+16.0000)/5.60000))
      ALGBRC(3) = 1.00000/( 0.0350000*cosh((STATES(8)+16.0000)/22.4000))
      RATES(9) = (ALGBRC(6) - STATES(9))/ALGBRC(3)
      ALGBRC(5) = ( CONSTS(8)*(STATES(2) - STATES(1))*CONSTS(7))/( (CONSTS(7)+STATES(2))*(CONSTS(7)+STATES(1)))
      ALGBRC(7) = ( CONSTS(10)*STATES(1) ** CONSTS(11))/(CONSTS(9) ** CONSTS(11)+STATES(1) ** CONSTS(11))
      RATES(1) = ALGBRC(5) - ALGBRC(7)
      ALGBRC(14) = 1.00000/(1.00000+EXP(- (STATES(8)+20.0000)/12.0000))
      ALGBRC(15) =  CONSTS(32)*ALGBRC(14)*(STATES(8) - CONSTS(33))
      ALGBRC(17) = - ( CONSTS(41)*ALGBRC(15)+ CONSTS(40)*STATES(7))
      ALGBRC(21) =  CONSTS(43)*(STATES(10) - STATES(7))
      ALGBRC(19) =  CONSTS(42)*STATES(7)
      ALGBRC(23) = ALGBRC(21) - ALGBRC(19)
      RATES(7) =  CONSTS(37)*(ALGBRC(17)+ALGBRC(23))
      RATES(10) =  - CONSTS(38)*CONSTS(39)*ALGBRC(23)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      ALGBRC(8) =  0.200000* (STATES(5)/1.00000) ** (1.0 / 2)
      ALGBRC(12) = ( CONSTS(27)*ALGBRC(8))/(CONSTS(28)+ALGBRC(8))
      RATES(6) =  (1.00000/CONSTS(24))*(ALGBRC(25) -  STATES(6)*EXP( (CONSTS(25)+ALGBRC(12))*(1.00000 - STATES(7)/CONSTS(26))))
      ALGBRC(13) =  CONSTS(30)*STATES(9)*(STATES(8) - CONSTS(31))
      ALGBRC(16) = ( CONSTS(34)*STATES(7) ** 2.00000)/(CONSTS(35) ** 2.00000+STATES(7) ** 2.00000)
      ALGBRC(18) =  ALGBRC(16)*(STATES(8) - CONSTS(31))
      ALGBRC(20) =  0.165000*STATES(6)
      ALGBRC(22) =  0.135000*STATES(6)
      ALGBRC(27) =  0.00500000*ALGBRC(25)
      ALGBRC(29) = ( 0.0800000*(1.00000+( 2.00000*ALGBRC(20))/17.0000)+ 0.890000*ALGBRC(20)/17.0000 ** 2.00000)/( 1.00000+ALGBRC(20)/17.0000 ** 2.00000*(1.00000+ALGBRC(22)/26.0000+ALGBRC(27)/1.00000))
      ALGBRC(31) =  CONSTS(36)*ALGBRC(29)
      ALGBRC(33) =  ALGBRC(31)*(STATES(8) - CONSTS(31))
      RATES(8) = - (ALGBRC(13)+ALGBRC(15)+ALGBRC(18)+ALGBRC(33))/CONSTS(29)
      ALGBRC(26) = ALGBRC(24)/CONSTS(19) ** 1.00000
      ALGBRC(9) = STATES(5)/CONSTS(20) ** 1.00000
      ALGBRC(1) =  0.300000*STATES(4)
      ALGBRC(10) = ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(28) = ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(30) =  ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000
      ALGBRC(32) =  (1.00000/CONSTS(14) **  1.00000*1.00000)*ALGBRC(24)/CONSTS(19) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(34) =  (1.00000/CONSTS(16) **  1.00000*1.00000)*ALGBRC(24)/CONSTS(19) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(35) =  (1.00000/CONSTS(17) **  1.00000*1.00000)*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(36) =  (1.00000/CONSTS(18) **  1.00000*1.00000)*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(11) =  (1.00000/CONSTS(15) **  1.00000*1.00000)*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(37) =  (1.00000/( CONSTS(14) **  1.00000*1.00000*CONSTS(15) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(38) =  (1.00000/( CONSTS(15) **  1.00000*1.00000*CONSTS(17) **  1.00000*1.00000*CONSTS(18) **  1.00000*1.00000))*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(39) =  (1.00000/( CONSTS(14) **  1.00000*1.00000*CONSTS(16) **  1.00000*1.00000*CONSTS(18) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(40) =  (1.00000/( CONSTS(16) **  1.00000*1.00000*CONSTS(17) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(41) =  (1.00000/( CONSTS(14) **  1.00000*1.00000*CONSTS(15) **  1.00000*1.00000*CONSTS(16) **  1.00000*1.00000*CONSTS(17) **  1.00000*1.00000*CONSTS(18) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(42) = ( CONSTS(13)*( (1.00000 - CONSTS(12))*ALGBRC(37)+ CONSTS(12)*(ALGBRC(10)+ALGBRC(32)+ALGBRC(36)+ALGBRC(37)+ALGBRC(38)+ALGBRC(39)+ALGBRC(11)+ALGBRC(41))))/(CONSTS(47)+ALGBRC(26)+ALGBRC(9)+ALGBRC(10)+ALGBRC(28)+ALGBRC(30)+ALGBRC(32)+ALGBRC(34)+ALGBRC(35)+ALGBRC(36)+ALGBRC(11)+ALGBRC(37)+ALGBRC(38)+ALGBRC(39)+ALGBRC(40)+ALGBRC(41))
      RATES(4) =  CONSTS(6)*(ALGBRC(7) - ALGBRC(42))
      RATES(5) =  CONSTS(6)*(ALGBRC(42) -  0.500000*ALGBRC(8))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = CONSTS(3)+(CONSTS(2) - CONSTS(3))/(1.00000+EXP((STATES(3) - CONSTS(5))/CONSTS(4)))
      ALGBRC(4) = TERNRY(STATES(7).GE.CONSTS(46),  CONSTS(45)*(STATES(7) - CONSTS(46)), TERNRY(STATES(7).LT.CONSTS(46), 0.00000, 0.0/0.0)
      ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(8)+16.0000)/5.60000))
      ALGBRC(3) = 1.00000/( 0.0350000*cosh((STATES(8)+16.0000)/22.4000))
      ALGBRC(5) = ( CONSTS(8)*(STATES(2) - STATES(1))*CONSTS(7))/( (CONSTS(7)+STATES(2))*(CONSTS(7)+STATES(1)))
      ALGBRC(7) = ( CONSTS(10)*STATES(1) ** CONSTS(11))/(CONSTS(9) ** CONSTS(11)+STATES(1) ** CONSTS(11))
      ALGBRC(14) = 1.00000/(1.00000+EXP(- (STATES(8)+20.0000)/12.0000))
      ALGBRC(15) =  CONSTS(32)*ALGBRC(14)*(STATES(8) - CONSTS(33))
      ALGBRC(17) = - ( CONSTS(41)*ALGBRC(15)+ CONSTS(40)*STATES(7))
      ALGBRC(21) =  CONSTS(43)*(STATES(10) - STATES(7))
      ALGBRC(19) =  CONSTS(42)*STATES(7)
      ALGBRC(23) = ALGBRC(21) - ALGBRC(19)
      ALGBRC(8) =  0.200000* (STATES(5)/1.00000) ** (1.0 / 2)
      ALGBRC(12) = ( CONSTS(27)*ALGBRC(8))/(CONSTS(28)+ALGBRC(8))
      ALGBRC(13) =  CONSTS(30)*STATES(9)*(STATES(8) - CONSTS(31))
      ALGBRC(16) = ( CONSTS(34)*STATES(7) ** 2.00000)/(CONSTS(35) ** 2.00000+STATES(7) ** 2.00000)
      ALGBRC(18) =  ALGBRC(16)*(STATES(8) - CONSTS(31))
      ALGBRC(20) =  0.165000*STATES(6)
      ALGBRC(22) =  0.135000*STATES(6)
      ALGBRC(27) =  0.00500000*ALGBRC(25)
      ALGBRC(29) = ( 0.0800000*(1.00000+( 2.00000*ALGBRC(20))/17.0000)+ 0.890000*ALGBRC(20)/17.0000 ** 2.00000)/( 1.00000+ALGBRC(20)/17.0000 ** 2.00000*(1.00000+ALGBRC(22)/26.0000+ALGBRC(27)/1.00000))
      ALGBRC(31) =  CONSTS(36)*ALGBRC(29)
      ALGBRC(33) =  ALGBRC(31)*(STATES(8) - CONSTS(31))
      ALGBRC(26) = ALGBRC(24)/CONSTS(19) ** 1.00000
      ALGBRC(9) = STATES(5)/CONSTS(20) ** 1.00000
      ALGBRC(1) =  0.300000*STATES(4)
      ALGBRC(10) = ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(28) = ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(30) =  ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000
      ALGBRC(32) =  (1.00000/CONSTS(14) **  1.00000*1.00000)*ALGBRC(24)/CONSTS(19) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(34) =  (1.00000/CONSTS(16) **  1.00000*1.00000)*ALGBRC(24)/CONSTS(19) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(35) =  (1.00000/CONSTS(17) **  1.00000*1.00000)*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(36) =  (1.00000/CONSTS(18) **  1.00000*1.00000)*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(11) =  (1.00000/CONSTS(15) **  1.00000*1.00000)*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(37) =  (1.00000/( CONSTS(14) **  1.00000*1.00000*CONSTS(15) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000
      ALGBRC(38) =  (1.00000/( CONSTS(15) **  1.00000*1.00000*CONSTS(17) **  1.00000*1.00000*CONSTS(18) **  1.00000*1.00000))*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(39) =  (1.00000/( CONSTS(14) **  1.00000*1.00000*CONSTS(16) **  1.00000*1.00000*CONSTS(18) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(40) =  (1.00000/( CONSTS(16) **  1.00000*1.00000*CONSTS(17) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(41) =  (1.00000/( CONSTS(14) **  1.00000*1.00000*CONSTS(15) **  1.00000*1.00000*CONSTS(16) **  1.00000*1.00000*CONSTS(17) **  1.00000*1.00000*CONSTS(18) **  1.00000*1.00000))*ALGBRC(24)/CONSTS(19) ** 1.00000*STATES(5)/CONSTS(20) ** 1.00000*ALGBRC(1) ** 2.00000/CONSTS(21) ** 1.00000*ALGBRC(25) ** 2.00000/CONSTS(22) ** 1.00000
      ALGBRC(42) = ( CONSTS(13)*( (1.00000 - CONSTS(12))*ALGBRC(37)+ CONSTS(12)*(ALGBRC(10)+ALGBRC(32)+ALGBRC(36)+ALGBRC(37)+ALGBRC(38)+ALGBRC(39)+ALGBRC(11)+ALGBRC(41))))/(CONSTS(47)+ALGBRC(26)+ALGBRC(9)+ALGBRC(10)+ALGBRC(28)+ALGBRC(30)+ALGBRC(32)+ALGBRC(34)+ALGBRC(35)+ALGBRC(36)+ALGBRC(11)+ALGBRC(37)+ALGBRC(38)+ALGBRC(39)+ALGBRC(40)+ALGBRC(41))
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(25) - (CONSTS(23) - (STATES(6)+ALGBRC(24))))
      minfunc_0[2] = abs(ALGBRC(24) - STATES(6) ** 2.00000/ALGBRC(25))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END