Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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