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 79 entries in the algebraic variable array.
C There are a total of 11 entries in each of the rate and state variable arrays.
C There are a total of 75 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is Vm in component membrane (millivolt).
C CONSTS(1) is cm in component membrane (femtofarad).
C ALGBRC(4) is Ica in component Ica (femtoampere).
C ALGBRC(2) is Ik in component Ik (femtoampere).
C ALGBRC(5) is Ikca in component Ikca (femtoampere).
C ALGBRC(44) is Ikatp in component Ikatp (femtoampere).
C CONSTS(2) is gK in component Ik (picosiemens).
C CONSTS(3) is VK in component model_parameters (millivolt).
C STATES(2) is n in component n (dimensionless).
C ALGBRC(1) is n_infinity in component n (dimensionless).
C CONSTS(4) is tau_n in component n (millisecond).
C CONSTS(5) is vn in component n (millivolt).
C CONSTS(6) is sn in component n (millivolt).
C CONSTS(7) is gCa in component Ica (picosiemens).
C CONSTS(8) is VCa in component model_parameters (millivolt).
C ALGBRC(3) is m_infinity in component m (dimensionless).
C CONSTS(9) is v in component m (millivolt).
C CONSTS(10) is sm in component m (millivolt).
C CONSTS(11) is gkCa in component Ikca (picosiemens).
C CONSTS(12) is kd in component Ikca (micromolar).
C STATES(3) is c in component c (micromolar).
C CONSTS(13) is gkATP_ in component Ikatp (picosiemens).
C ALGBRC(42) is katpo in component Ikatp (dimensionless).
C ALGBRC(7) is topo in component Ikatp (dimensionless).
C ALGBRC(40) is bottomo in component Ikatp (dimensionless).
C ALGBRC(6) is mgadp in component Ikatp (micromolar).
C ALGBRC(8) is adp3m in component Ikatp (micromolar).
C ALGBRC(37) is atp4m in component Ikatp (micromolar).
C ALGBRC(36) is atp in component atp (micromolar).
C STATES(4) is adp in component adp (micromolar).
C ALGBRC(9) is JGPDH in component JGPDH (micromolar_millisecond).
C CONSTS(14) is kGPDH in component JGPDH (micromolar_millisecond).
C STATES(5) is FBP in component FBP (micromolar).
C STATES(6) is G6P in component G6P (micromolar).
C CONSTS(75) is JGK_ms in component JGK (micromolar_millisecond).
C ALGBRC(79) is JPFK_ms in component JPFK (micromolar_millisecond).
C ALGBRC(10) is F6P in component F6P (micromolar).
C ALGBRC(78) is JPFK in component JPFK (micromolar_second).
C CONSTS(15) is bottom1 in component JPFK (dimensionless).
C CONSTS(16) is topa1 in component JPFK (dimensionless).
C CONSTS(17) is k1 in component JPFK (micromolar).
C CONSTS(18) is k2 in component JPFK (micromolar).
C CONSTS(19) is k3 in component JPFK (micromolar).
C CONSTS(20) is k4 in component JPFK (micromolar).
C CONSTS(21) is VmaxPFK in component JPFK (micromolar_millisecond).
C ALGBRC(38) is weight2 in component JPFK (dimensionless).
C CONSTS(72) is topa2 in component JPFK (dimensionless).
C ALGBRC(41) is bottom2 in component JPFK (dimensionless).
C ALGBRC(14) is topa3 in component JPFK (dimensionless).
C ALGBRC(12) is weight3 in component JPFK (dimensionless).
C ALGBRC(43) is bottom3 in component JPFK (dimensionless).
C CONSTS(22) is f13 in component JPFK (dimensionless).
C CONSTS(23) is f43 in component JPFK (dimensionless).
C CONSTS(24) is f23 in component JPFK (dimensionless).
C CONSTS(25) is f42 in component JPFK (dimensionless).
C CONSTS(26) is f41 in component JPFK (dimensionless).
C ALGBRC(45) is weight4 in component JPFK (dimensionless).
C ALGBRC(46) is topa4 in component JPFK (dimensionless).
C ALGBRC(47) is bottom4 in component JPFK (dimensionless).
C ALGBRC(16) is weight5 in component JPFK (dimensionless).
C ALGBRC(48) is topa5 in component JPFK (dimensionless).
C ALGBRC(49) is bottom5 in component JPFK (dimensionless).
C ALGBRC(50) is weight6 in component JPFK (dimensionless).
C ALGBRC(51) is topa6 in component JPFK (dimensionless).
C ALGBRC(52) is bottom6 in component JPFK (dimensionless).
C ALGBRC(17) is weight7 in component JPFK (dimensionless).
C ALGBRC(53) is topa7 in component JPFK (dimensionless).
C ALGBRC(54) is bottom7 in component JPFK (dimensionless).
C ALGBRC(55) is weight8 in component JPFK (dimensionless).
C ALGBRC(56) is topa8 in component JPFK (dimensionless).
C ALGBRC(57) is bottom8 in component JPFK (dimensionless).
C CONSTS(74) is weight9 in component JPFK (dimensionless).
C ALGBRC(58) is topa9 in component JPFK (dimensionless).
C ALGBRC(59) is bottom9 in component JPFK (dimensionless).
C ALGBRC(60) is weight10 in component JPFK (dimensionless).
C ALGBRC(61) is topa10 in component JPFK (dimensionless).
C ALGBRC(62) is bottom10 in component JPFK (dimensionless).
C ALGBRC(18) is weight11 in component JPFK (dimensionless).
C ALGBRC(63) is topa11 in component JPFK (dimensionless).
C ALGBRC(64) is bottom11 in component JPFK (dimensionless).
C ALGBRC(65) is weight12 in component JPFK (dimensionless).
C ALGBRC(66) is topa12 in component JPFK (dimensionless).
C ALGBRC(67) is bottom12 in component JPFK (dimensionless).
C ALGBRC(19) is weight13 in component JPFK (dimensionless).
C ALGBRC(68) is topa13 in component JPFK (dimensionless).
C ALGBRC(69) is bottom13 in component JPFK (dimensionless).
C ALGBRC(70) is weight14 in component JPFK (dimensionless).
C ALGBRC(71) is topa14 in component JPFK (dimensionless).
C ALGBRC(72) is bottom14 in component JPFK (dimensionless).
C ALGBRC(20) is weight15 in component JPFK (dimensionless).
C ALGBRC(73) is topa15 in component JPFK (dimensionless).
C ALGBRC(74) is bottom15 in component JPFK (dimensionless).
C ALGBRC(75) is weight16 in component JPFK (dimensionless).
C ALGBRC(76) is topa16 in component JPFK (dimensionless).
C ALGBRC(77) is bottom16 in component JPFK (dimensionless).
C ALGBRC(21) is topb in component JPFK (dimensionless).
C CONSTS(27) is AMP in component model_parameters (micromolar).
C CONSTS(28) is lambda in component JPFK (dimensionless).
C CONSTS(29) is kappa in component model_parameters (dimensionless).
C ALGBRC(15) is JPDH in component JPDH (micromolar_millisecond).
C CONSTS(30) is p1 in component JPDH (dimensionless).
C CONSTS(31) is p2 in component JPDH (dimensionless).
C CONSTS(32) is p3 in component JPDH (micromolar).
C CONSTS(33) is JGPDHbas in component JPDH (micromolar_millisecond).
C ALGBRC(13) is NADm in component NADm (millimolar).
C STATES(7) is Cam in component Cam (micromolar).
C STATES(8) is NADHm in component NADHm (millimolar).
C CONSTS(34) is gamma in component model_parameters (dimensionless).
C ALGBRC(11) is JO in component JO (micromolar_millisecond).
C CONSTS(35) is p4 in component JO (micromolar_millisecond).
C CONSTS(36) is p5 in component JO (millimolar).
C CONSTS(37) is p6 in component JO (millivolt).
C CONSTS(38) is p7 in component JO (millivolt).
C STATES(9) is delta_psi in component delta_psi (millivolt).
C CONSTS(39) is NADmtot in component NADm (millimolar).
C CONSTS(40) is Cmito in component delta_psi (micromolar_millivolt).
C ALGBRC(22) is JHres in component JHres (micromolar_millisecond).
C ALGBRC(30) is JHatp in component JHatp (micromolar_millisecond).
C ALGBRC(34) is JANT in component JANT (micromolar_millisecond).
C ALGBRC(23) is JHleak in component JHleak (micromolar_millisecond).
C ALGBRC(25) is JNaCa in component JNaCa (micromolar_millisecond).
C ALGBRC(24) is Juni in component Juni (micromolar_millisecond).
C CONSTS(41) is p8 in component JHres (micromolar_millisecond).
C CONSTS(42) is p9 in component JHres (millimolar).
C CONSTS(43) is p10 in component JHres (millivolt).
C CONSTS(44) is p11 in component JHres (millivolt).
C ALGBRC(28) is JF1F0 in component JF1F0 (micromolar_millisecond).
C CONSTS(45) is p13 in component JF1F0 (millimolar).
C CONSTS(46) is p14 in component JF1F0 (millivolt).
C CONSTS(47) is p15 in component JF1F0 (millivolt).
C CONSTS(48) is p16 in component JF1F0 (micromolar_millisecond).
C ALGBRC(26) is ATPm in component ATPm (millimolar).
C CONSTS(49) is JGK in component JGK (micromolar_second).
C CONSTS(50) is p17 in component JHleak (micromolar_millisecond_millivolt).
C CONSTS(51) is p18 in component JHleak (micromolar_millisecond).
C CONSTS(52) is p19 in component JANT (micromolar_millisecond).
C CONSTS(53) is p20 in component JANT (dimensionless).
C CONSTS(54) is FRT in component JANT (per_millivolt).
C ALGBRC(32) is RATm in component RATm (dimensionless).
C STATES(10) is ADPm in component ADPm (millimolar).
C CONSTS(55) is p21 in component Juni (per_micromolar_millisecond_millivolt).
C CONSTS(56) is p22 in component Juni (second_order_rate_constant).
C CONSTS(57) is p23 in component JNaCa (micromolar_millisecond).
C CONSTS(58) is p24 in component JNaCa (per_millivolt).
C CONSTS(59) is fmito in component Cam (dimensionless).
C ALGBRC(27) is Jmito in component Jmito (micromolar_millisecond).
C CONSTS(60) is Amtot in component ATPm (millimolar).
C CONSTS(73) is delta in component model_parameters (dimensionless).
C ALGBRC(39) is Jhyd in component Jhyd (micromolar_millisecond).
C CONSTS(61) is khyd in component Jhyd (second_order_rate_constant).
C CONSTS(62) is khydbas in component Jhyd (first_order_rate_constant).
C CONSTS(63) is atot in component atp (micromolar).
C CONSTS(64) is fcyt in component c (dimensionless).
C ALGBRC(35) is Jer in component Jer (micromolar_millisecond).
C ALGBRC(29) is Jmem in component Jmem (micromolar_millisecond).
C CONSTS(65) is kPMCA in component Jmem (first_order_rate_constant).
C CONSTS(66) is alpha in component Jmem (micromolar_millisecond_femtoampere).
C CONSTS(67) is Cbas in component Jmem (micromolar).
C ALGBRC(31) is Jleak in component Jleak (micromolar_millisecond).
C CONSTS(68) is pleak in component Jleak (first_order_rate_constant).
C STATES(11) is Caer in component Caer (micromolar).
C ALGBRC(33) is JSERCA in component JSERCA (micromolar_millisecond).
C CONSTS(69) is kSERCA in component JSERCA (first_order_rate_constant).
C CONSTS(70) is fer in component Caer (dimensionless).
C CONSTS(71) is Vc_Ver in component Caer (dimensionless).
C RATES(1) is d/dt Vm in component membrane (millivolt).
C RATES(2) is d/dt n in component n (dimensionless).
C RATES(6) is d/dt G6P in component G6P (micromolar).
C RATES(5) is d/dt FBP in component FBP (micromolar).
C RATES(8) is d/dt NADHm in component NADHm (millimolar).
C RATES(9) is d/dt delta_psi in component delta_psi (millivolt).
C RATES(7) is d/dt Cam in component Cam (micromolar).
C RATES(10) is d/dt ADPm in component ADPm (millimolar).
C RATES(4) is d/dt adp in component adp (micromolar).
C RATES(3) is d/dt c in component c (micromolar).
C RATES(11) is d/dt Caer in component Caer (micromolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -60.0
CONSTS(1) = 5300.0
CONSTS(2) = 2700.0
CONSTS(3) = -75.0
STATES(2) = 0.0
CONSTS(4) = 20.0
CONSTS(5) = -16.0
CONSTS(6) = 5.0
CONSTS(7) = 1000.0
CONSTS(8) = 25.0
CONSTS(9) = -20.0
CONSTS(10) = 12.0
CONSTS(11) = 300.0
CONSTS(12) = 0.5
STATES(3) = 0.17
CONSTS(13) = 16000.0
STATES(4) = 1137.0
CONSTS(14) = 0.0005
STATES(5) = 2.16
STATES(6) = 301.0
CONSTS(15) = 1
CONSTS(16) = 0
CONSTS(17) = 30
CONSTS(18) = 1
CONSTS(19) = 50000
CONSTS(20) = 1000
CONSTS(21) = 5.0
CONSTS(22) = 0.02
CONSTS(23) = 20
CONSTS(24) = 0.2
CONSTS(25) = 20
CONSTS(26) = 20
CONSTS(27) = 500.0
CONSTS(28) = 0.06
CONSTS(29) = 0.001
CONSTS(30) = 400.0
CONSTS(31) = 1.0
CONSTS(32) = 0.01
CONSTS(33) = 0.0005
STATES(7) = 0.2
STATES(8) = 0.4
CONSTS(34) = 0.001
CONSTS(35) = 0.6
CONSTS(36) = 0.1
CONSTS(37) = 177.0
CONSTS(38) = 5.0
STATES(9) = 164.0
CONSTS(39) = 10.0
CONSTS(40) = 1.8
CONSTS(41) = 7.0
CONSTS(42) = 0.1
CONSTS(43) = 177.0
CONSTS(44) = 5.0
CONSTS(45) = 10.0
CONSTS(46) = 190.0
CONSTS(47) = 8.5
CONSTS(48) = 35.0
CONSTS(49) = 0.4
CONSTS(50) = 0.002
CONSTS(51) = -0.03
CONSTS(52) = 0.35
CONSTS(53) = 2.0
CONSTS(54) = 0.037410133
STATES(10) = 11.1
CONSTS(55) = 0.04
CONSTS(56) = 1.1
CONSTS(57) = 0.01
CONSTS(58) = 0.016
CONSTS(59) = 0.01
CONSTS(60) = 15.0
CONSTS(61) = 0.00005
CONSTS(62) = 0.00005
CONSTS(63) = 2500.0
CONSTS(64) = 0.01
CONSTS(65) = 0.1
CONSTS(66) = 4.5E-6
CONSTS(67) = 0.05
CONSTS(68) = 0.0002
STATES(11) = 345.0
CONSTS(69) = 0.4
CONSTS(70) = 0.01
CONSTS(71) = 31.0
CONSTS(72) = CONSTS(16)
CONSTS(73) = 3.90000/53.2000
CONSTS(74) = CONSTS(27)/CONSTS(17)
CONSTS(75) = CONSTS(29)*CONSTS(49)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = 1.00000/(1.00000+EXP((CONSTS(5) - STATES(1))/CONSTS(6)))
RATES(2) = (ALGBRC(1) - STATES(2))/CONSTS(4)
ALGBRC(9) = CONSTS(14)* ABS(STATES(5)/1.00000) ** (1.0 / 2)
ALGBRC(13) = CONSTS(39) - STATES(8)
ALGBRC(15) = (CONSTS(30)/(CONSTS(31)+STATES(8)/ALGBRC(13)))*(STATES(7)/(CONSTS(32)+STATES(7)))*(ALGBRC(9)+CONSTS(33))
ALGBRC(11) = ( CONSTS(35)*(STATES(8)/(CONSTS(36)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(37))/CONSTS(38)))
RATES(8) = CONSTS(34)*(ALGBRC(15) - ALGBRC(11))
ALGBRC(25) = (CONSTS(57)/STATES(3))*STATES(7)*EXP( CONSTS(58)*STATES(9))
ALGBRC(24) = ( CONSTS(55)*STATES(9) - CONSTS(56))*STATES(3) ** 2.00000
ALGBRC(27) = ALGBRC(25) - ALGBRC(24)
RATES(7) = - CONSTS(59)*ALGBRC(27)
ALGBRC(22) = ( CONSTS(41)*(STATES(8)/(CONSTS(42)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(43))/CONSTS(44)))
ALGBRC(26) = CONSTS(60) - STATES(10)
ALGBRC(28) = (( CONSTS(48)*CONSTS(45))/(CONSTS(45)+ALGBRC(26)))/(1.00000+EXP((CONSTS(46) - STATES(9))/CONSTS(47)))
ALGBRC(30) = 3.00000*ALGBRC(28)
ALGBRC(32) = ALGBRC(26)/STATES(10)
ALGBRC(34) = CONSTS(52)*((ALGBRC(32)/(ALGBRC(32)+CONSTS(53)))/EXP( -0.500000*CONSTS(54)*STATES(9)))
ALGBRC(23) = CONSTS(50)*STATES(9)+CONSTS(51)
RATES(9) = (ALGBRC(22) - (ALGBRC(30)+ALGBRC(34)+ALGBRC(23)+ALGBRC(25)+ 2.00000*ALGBRC(24)))/CONSTS(40)
RATES(10) = CONSTS(34)*(ALGBRC(34) - ALGBRC(28))
ALGBRC(31) = CONSTS(68)*(STATES(11) - STATES(3))
ALGBRC(33) = CONSTS(69)*STATES(3)
ALGBRC(35) = ALGBRC(31) - ALGBRC(33)
ALGBRC(3) = 1.00000/(1.00000+EXP((CONSTS(9) - STATES(1))/CONSTS(10)))
ALGBRC(4) = CONSTS(7)*ALGBRC(3)*(STATES(1) - CONSTS(8))
ALGBRC(29) = - ( CONSTS(66)*ALGBRC(4)+ CONSTS(65)*(STATES(3) - CONSTS(67)))
RATES(3) = CONSTS(64)*(ALGBRC(29)+ALGBRC(35)+ CONSTS(73)*ALGBRC(27))
RATES(11) = - CONSTS(70)*CONSTS(71)*ALGBRC(35)
ALGBRC(36) = CONSTS(63) - STATES(4)
ALGBRC(39) = ( CONSTS(61)*STATES(3)+CONSTS(62))*ALGBRC(36)
RATES(4) = - CONSTS(73)*ALGBRC(34)+ALGBRC(39)
ALGBRC(2) = CONSTS(2)*STATES(2)*(STATES(1) - CONSTS(3))
ALGBRC(5) = (CONSTS(11)/(1.00000+CONSTS(12)/STATES(3) ** 2.00000))*(STATES(1) - CONSTS(3))
ALGBRC(6) = 0.165000*STATES(4)
ALGBRC(7) = 0.0800000*(1.00000+( 2.00000*ALGBRC(6))/17.0000)+ 0.890000*ALGBRC(6)/17.0000 ** 2.00000
ALGBRC(8) = 0.135000*STATES(4)
ALGBRC(37) = 0.0500000*ALGBRC(36)
ALGBRC(40) = 1.00000+ALGBRC(6)/17.0000 ** 2.00000*(1.00000+ALGBRC(8)/26.0000+ALGBRC(37)/1.00000)
ALGBRC(42) = ALGBRC(7)/ALGBRC(40)
ALGBRC(44) = CONSTS(13)*ALGBRC(42)*(STATES(1) - CONSTS(3))
RATES(1) = - (ALGBRC(2)+ALGBRC(4)+ALGBRC(5)+ALGBRC(44))/CONSTS(1)
ALGBRC(10) = 0.300000*STATES(6)
ALGBRC(12) = ALGBRC(10) ** 2.00000/( CONSTS(19)*1.00000)
ALGBRC(14) = CONSTS(72)+ALGBRC(12)
ALGBRC(45) = ALGBRC(10)*ALGBRC(36) ** 2.00000/( CONSTS(23)*CONSTS(19)*CONSTS(20)*1.00000 ** 2.00000)
ALGBRC(46) = ALGBRC(14)+ALGBRC(45)
ALGBRC(48) = ALGBRC(46)
ALGBRC(51) = ALGBRC(48)
ALGBRC(17) = ( STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(24)*1.00000)
ALGBRC(53) = ALGBRC(51)+ALGBRC(17)
ALGBRC(55) = ( STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(25)*CONSTS(23)*1.00000 ** 2.00000)
ALGBRC(56) = ALGBRC(53)+ALGBRC(55)
ALGBRC(58) = ALGBRC(56)
ALGBRC(61) = ALGBRC(58)
ALGBRC(18) = ( CONSTS(27)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(22)*1.00000)
ALGBRC(63) = ALGBRC(61)+ALGBRC(18)
ALGBRC(65) = ( CONSTS(27)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(20)*CONSTS(22)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
ALGBRC(66) = ALGBRC(63)+ALGBRC(65)
ALGBRC(68) = ALGBRC(66)
ALGBRC(71) = ALGBRC(68)
ALGBRC(73) = ALGBRC(71)
ALGBRC(75) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(22)*CONSTS(25)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
ALGBRC(76) = ALGBRC(73)+ALGBRC(75)
ALGBRC(38) = ALGBRC(36) ** 2.00000/( CONSTS(20)*1.00000)
ALGBRC(41) = CONSTS(15)+ALGBRC(38)
ALGBRC(43) = ALGBRC(41)+ALGBRC(12)
ALGBRC(47) = ALGBRC(43)+ALGBRC(45)
ALGBRC(16) = STATES(5)/CONSTS(18)
ALGBRC(49) = ALGBRC(47)+ALGBRC(16)
ALGBRC(50) = ( STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(20)*CONSTS(25)*1.00000)
ALGBRC(52) = ALGBRC(49)+ALGBRC(50)
ALGBRC(54) = ALGBRC(52)+ALGBRC(17)
ALGBRC(57) = ALGBRC(54)+ALGBRC(55)
ALGBRC(59) = ALGBRC(57)+CONSTS(74)
ALGBRC(60) = ( CONSTS(27)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(20)*CONSTS(26)*1.00000)
ALGBRC(62) = ALGBRC(59)+ALGBRC(60)
ALGBRC(64) = ALGBRC(62)+ALGBRC(18)
ALGBRC(67) = ALGBRC(64)+ALGBRC(65)
ALGBRC(19) = ( CONSTS(27)*STATES(5))/( CONSTS(17)*CONSTS(18))
ALGBRC(69) = ALGBRC(67)+ALGBRC(19)
ALGBRC(70) = ( CONSTS(27)*STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(20)*CONSTS(25)*CONSTS(26)*1.00000)
ALGBRC(72) = ALGBRC(69)+ALGBRC(70)
ALGBRC(20) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(24)*CONSTS(22)*1.00000)
ALGBRC(74) = ALGBRC(72)+ALGBRC(20)
ALGBRC(77) = ALGBRC(74)+ALGBRC(75)
ALGBRC(21) = ALGBRC(20)
ALGBRC(78) = ( CONSTS(28)*CONSTS(21)*ALGBRC(76)+ CONSTS(21)*ALGBRC(21))/ALGBRC(77)
ALGBRC(79) = CONSTS(29)*ALGBRC(78)
RATES(6) = CONSTS(75) - ALGBRC(79)
RATES(5) = ALGBRC(79) - 0.500000*ALGBRC(9)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = 1.00000/(1.00000+EXP((CONSTS(5) - STATES(1))/CONSTS(6)))
ALGBRC(9) = CONSTS(14)* ABS(STATES(5)/1.00000) ** (1.0 / 2)
ALGBRC(13) = CONSTS(39) - STATES(8)
ALGBRC(15) = (CONSTS(30)/(CONSTS(31)+STATES(8)/ALGBRC(13)))*(STATES(7)/(CONSTS(32)+STATES(7)))*(ALGBRC(9)+CONSTS(33))
ALGBRC(11) = ( CONSTS(35)*(STATES(8)/(CONSTS(36)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(37))/CONSTS(38)))
ALGBRC(25) = (CONSTS(57)/STATES(3))*STATES(7)*EXP( CONSTS(58)*STATES(9))
ALGBRC(24) = ( CONSTS(55)*STATES(9) - CONSTS(56))*STATES(3) ** 2.00000
ALGBRC(27) = ALGBRC(25) - ALGBRC(24)
ALGBRC(22) = ( CONSTS(41)*(STATES(8)/(CONSTS(42)+STATES(8))))/(1.00000+EXP((STATES(9) - CONSTS(43))/CONSTS(44)))
ALGBRC(26) = CONSTS(60) - STATES(10)
ALGBRC(28) = (( CONSTS(48)*CONSTS(45))/(CONSTS(45)+ALGBRC(26)))/(1.00000+EXP((CONSTS(46) - STATES(9))/CONSTS(47)))
ALGBRC(30) = 3.00000*ALGBRC(28)
ALGBRC(32) = ALGBRC(26)/STATES(10)
ALGBRC(34) = CONSTS(52)*((ALGBRC(32)/(ALGBRC(32)+CONSTS(53)))/EXP( -0.500000*CONSTS(54)*STATES(9)))
ALGBRC(23) = CONSTS(50)*STATES(9)+CONSTS(51)
ALGBRC(31) = CONSTS(68)*(STATES(11) - STATES(3))
ALGBRC(33) = CONSTS(69)*STATES(3)
ALGBRC(35) = ALGBRC(31) - ALGBRC(33)
ALGBRC(3) = 1.00000/(1.00000+EXP((CONSTS(9) - STATES(1))/CONSTS(10)))
ALGBRC(4) = CONSTS(7)*ALGBRC(3)*(STATES(1) - CONSTS(8))
ALGBRC(29) = - ( CONSTS(66)*ALGBRC(4)+ CONSTS(65)*(STATES(3) - CONSTS(67)))
ALGBRC(36) = CONSTS(63) - STATES(4)
ALGBRC(39) = ( CONSTS(61)*STATES(3)+CONSTS(62))*ALGBRC(36)
ALGBRC(2) = CONSTS(2)*STATES(2)*(STATES(1) - CONSTS(3))
ALGBRC(5) = (CONSTS(11)/(1.00000+CONSTS(12)/STATES(3) ** 2.00000))*(STATES(1) - CONSTS(3))
ALGBRC(6) = 0.165000*STATES(4)
ALGBRC(7) = 0.0800000*(1.00000+( 2.00000*ALGBRC(6))/17.0000)+ 0.890000*ALGBRC(6)/17.0000 ** 2.00000
ALGBRC(8) = 0.135000*STATES(4)
ALGBRC(37) = 0.0500000*ALGBRC(36)
ALGBRC(40) = 1.00000+ALGBRC(6)/17.0000 ** 2.00000*(1.00000+ALGBRC(8)/26.0000+ALGBRC(37)/1.00000)
ALGBRC(42) = ALGBRC(7)/ALGBRC(40)
ALGBRC(44) = CONSTS(13)*ALGBRC(42)*(STATES(1) - CONSTS(3))
ALGBRC(10) = 0.300000*STATES(6)
ALGBRC(12) = ALGBRC(10) ** 2.00000/( CONSTS(19)*1.00000)
ALGBRC(14) = CONSTS(72)+ALGBRC(12)
ALGBRC(45) = ALGBRC(10)*ALGBRC(36) ** 2.00000/( CONSTS(23)*CONSTS(19)*CONSTS(20)*1.00000 ** 2.00000)
ALGBRC(46) = ALGBRC(14)+ALGBRC(45)
ALGBRC(48) = ALGBRC(46)
ALGBRC(51) = ALGBRC(48)
ALGBRC(17) = ( STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(24)*1.00000)
ALGBRC(53) = ALGBRC(51)+ALGBRC(17)
ALGBRC(55) = ( STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(25)*CONSTS(23)*1.00000 ** 2.00000)
ALGBRC(56) = ALGBRC(53)+ALGBRC(55)
ALGBRC(58) = ALGBRC(56)
ALGBRC(61) = ALGBRC(58)
ALGBRC(18) = ( CONSTS(27)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(22)*1.00000)
ALGBRC(63) = ALGBRC(61)+ALGBRC(18)
ALGBRC(65) = ( CONSTS(27)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(19)*CONSTS(20)*CONSTS(22)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
ALGBRC(66) = ALGBRC(63)+ALGBRC(65)
ALGBRC(68) = ALGBRC(66)
ALGBRC(71) = ALGBRC(68)
ALGBRC(73) = ALGBRC(71)
ALGBRC(75) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(20)*CONSTS(24)*CONSTS(22)*CONSTS(25)*CONSTS(26)*CONSTS(23)*1.00000 ** 2.00000)
ALGBRC(76) = ALGBRC(73)+ALGBRC(75)
ALGBRC(38) = ALGBRC(36) ** 2.00000/( CONSTS(20)*1.00000)
ALGBRC(41) = CONSTS(15)+ALGBRC(38)
ALGBRC(43) = ALGBRC(41)+ALGBRC(12)
ALGBRC(47) = ALGBRC(43)+ALGBRC(45)
ALGBRC(16) = STATES(5)/CONSTS(18)
ALGBRC(49) = ALGBRC(47)+ALGBRC(16)
ALGBRC(50) = ( STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(18)*CONSTS(20)*CONSTS(25)*1.00000)
ALGBRC(52) = ALGBRC(49)+ALGBRC(50)
ALGBRC(54) = ALGBRC(52)+ALGBRC(17)
ALGBRC(57) = ALGBRC(54)+ALGBRC(55)
ALGBRC(59) = ALGBRC(57)+CONSTS(74)
ALGBRC(60) = ( CONSTS(27)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(20)*CONSTS(26)*1.00000)
ALGBRC(62) = ALGBRC(59)+ALGBRC(60)
ALGBRC(64) = ALGBRC(62)+ALGBRC(18)
ALGBRC(67) = ALGBRC(64)+ALGBRC(65)
ALGBRC(19) = ( CONSTS(27)*STATES(5))/( CONSTS(17)*CONSTS(18))
ALGBRC(69) = ALGBRC(67)+ALGBRC(19)
ALGBRC(70) = ( CONSTS(27)*STATES(5)*ALGBRC(36) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(20)*CONSTS(25)*CONSTS(26)*1.00000)
ALGBRC(72) = ALGBRC(69)+ALGBRC(70)
ALGBRC(20) = ( CONSTS(27)*STATES(5)*ALGBRC(10) ** 2.00000)/( CONSTS(17)*CONSTS(18)*CONSTS(19)*CONSTS(24)*CONSTS(22)*1.00000)
ALGBRC(74) = ALGBRC(72)+ALGBRC(20)
ALGBRC(77) = ALGBRC(74)+ALGBRC(75)
ALGBRC(21) = ALGBRC(20)
ALGBRC(78) = ( CONSTS(28)*CONSTS(21)*ALGBRC(76)+ CONSTS(21)*ALGBRC(21))/ALGBRC(77)
ALGBRC(79) = CONSTS(29)*ALGBRC(78)
RETURN
END