C
C There are a total of 49 entries in the algebraic variable array.
C There are a total of 18 entries in each of the rate and state variable arrays.
C There are a total of 46 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is R in component constants (joule_per_kilomole_kelvin).
C CONSTS(2) is T in component constants (kelvin).
C CONSTS(3) is F in component constants (coulomb_per_mole).
C STATES(1) is E in component membrane (millivolt).
C CONSTS(4) is C in component membrane (nanoF).
C ALGBRC(48) is i_tot in component membrane (picoA).
C ALGBRC(19) is i_CaL in component L_type_calcium_current (picoA).
C ALGBRC(20) is i_CaT in component T_type_calcium_current (picoA).
C ALGBRC(21) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(24) is i_K in component delayed_rectifying_potassium_current (picoA).
C ALGBRC(27) is i_f in component hyperpolarising_activated_current (picoA).
C ALGBRC(28) is i_p in component sodium_potassium_pump (picoA).
C ALGBRC(43) is i_NaCa in component sodium_calcium_exchange_current (picoA).
C ALGBRC(44) is i_bNa in component background_sodium_current (picoA).
C ALGBRC(46) is i_bK in component background_potassium_current (picoA).
C ALGBRC(9) is E_Ca in component reversal_potentials (millivolt).
C ALGBRC(17) is E_Na in component reversal_potentials (millivolt).
C ALGBRC(18) is E_K in component reversal_potentials (millivolt).
C STATES(2) is Cai in component ion_concentrations (millimolar).
C STATES(3) is Cao in component ion_concentrations (millimolar).
C STATES(4) is Nai in component ion_concentrations (millimolar).
C STATES(5) is Nao in component ion_concentrations (millimolar).
C STATES(6) is Ki in component ion_concentrations (millimolar).
C STATES(7) is Ko in component ion_concentrations (millimolar).
C CONSTS(5) is g_CaL in component L_type_calcium_current (nanoS).
C STATES(8) is dL in component L_type_calcium_current_d_gate (dimensionless).
C STATES(9) is fL in component L_type_calcium_current_f_gate (dimensionless).
C STATES(10) is fL2 in component L_type_calcium_current_f2_gate (dimensionless).
C ALGBRC(1) is dL_infinity in component L_type_calcium_current_d_gate (dimensionless).
C CONSTS(6) is tau_dL in component L_type_calcium_current_d_gate (second).
C ALGBRC(2) is fL_infinity in component L_type_calcium_current_f_gate (dimensionless).
C ALGBRC(10) is tau_fL in component L_type_calcium_current_f_gate (second).
C CONSTS(7) is alpha_fL2 in component L_type_calcium_current_f2_gate (per_second).
C CONSTS(8) is beta_fL2 in component L_type_calcium_current_f2_gate (per_millimolar_second).
C CONSTS(9) is g_CaT in component T_type_calcium_current (nanoS).
C STATES(11) is dT in component T_type_calcium_current_d_gate (dimensionless).
C STATES(12) is fT in component T_type_calcium_current_f_gate (dimensionless).
C ALGBRC(3) is dT_infinity in component T_type_calcium_current_d_gate (dimensionless).
C ALGBRC(11) is tau_dT in component T_type_calcium_current_d_gate (second).
C ALGBRC(4) is fT_infinity in component T_type_calcium_current_f_gate (dimensionless).
C ALGBRC(12) is tau_fT in component T_type_calcium_current_f_gate (second).
C CONSTS(10) is g_Na in component fast_sodium_current (nanoS).
C STATES(13) is m in component fast_sodium_current_m_gate (dimensionless).
C STATES(14) is h in component fast_sodium_current_h_gate (dimensionless).
C ALGBRC(5) is alpha_m in component fast_sodium_current_m_gate (per_second).
C ALGBRC(13) is beta_m in component fast_sodium_current_m_gate (per_second).
C ALGBRC(6) is alpha_h in component fast_sodium_current_h_gate (per_second).
C ALGBRC(14) is beta_h in component fast_sodium_current_h_gate (per_second).
C ALGBRC(22) is i_KK in component delayed_rectifying_potassium_current (picoA).
C ALGBRC(23) is i_KNa in component delayed_rectifying_potassium_current (picoA).
C CONSTS(11) is Kk in component delayed_rectifying_potassium_current (picoA_per_millimolar).
C CONSTS(12) is P_KNa in component delayed_rectifying_potassium_current (dimensionless).
C STATES(15) is x in component delayed_rectifying_potassium_current_x_gate (dimensionless).
C ALGBRC(7) is x_infinity in component delayed_rectifying_potassium_current_x_gate (dimensionless).
C ALGBRC(15) is tau_x in component delayed_rectifying_potassium_current_x_gate (second).
C ALGBRC(25) is i_fNa in component hyperpolarising_activated_current (picoA).
C ALGBRC(26) is i_fK in component hyperpolarising_activated_current (picoA).
C CONSTS(13) is Kmf in component hyperpolarising_activated_current (millimolar).
C CONSTS(14) is g_fNa in component hyperpolarising_activated_current (nanoS).
C CONSTS(15) is g_fK in component hyperpolarising_activated_current (nanoS).
C STATES(16) is y in component hyperpolarising_activated_current_y_gate (dimensionless).
C ALGBRC(8) is alpha_y in component hyperpolarising_activated_current_y_gate (per_second).
C ALGBRC(16) is beta_y in component hyperpolarising_activated_current_y_gate (per_second).
C CONSTS(16) is KmNa in component sodium_potassium_pump (millimolar).
C CONSTS(17) is KmK in component sodium_potassium_pump (millimolar).
C CONSTS(18) is i_pmax in component sodium_potassium_pump (picoA).
C CONSTS(19) is kNaCa in component sodium_calcium_exchange_current (picoA).
C ALGBRC(39) is x1 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(40) is x2 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(41) is x3 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(42) is x4 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(35) is k41 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(33) is k34 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(31) is k23 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(32) is k21 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(30) is k32 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(38) is k43 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(37) is k12 in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(36) is k14 in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(20) is Qci in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(21) is Qn in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(22) is Qco in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(23) is K3ni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(24) is Kci in component sodium_calcium_exchange_current (millimolar).
C CONSTS(25) is K1ni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(26) is K2ni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(27) is Kcni in component sodium_calcium_exchange_current (millimolar).
C CONSTS(28) is K3no in component sodium_calcium_exchange_current (millimolar).
C CONSTS(29) is K1no in component sodium_calcium_exchange_current (millimolar).
C CONSTS(30) is K2no in component sodium_calcium_exchange_current (millimolar).
C CONSTS(31) is Kco in component sodium_calcium_exchange_current (millimolar).
C ALGBRC(29) is do in component sodium_calcium_exchange_current (dimensionless).
C ALGBRC(34) is di in component sodium_calcium_exchange_current (dimensionless).
C CONSTS(32) is g_Nab in component background_sodium_current (nanoS).
C CONSTS(33) is KbK in component background_potassium_current (picoA_per_millimolar).
C ALGBRC(45) is i_up in component sarcoplasmic_reticulum_kinetics (picoA).
C ALGBRC(47) is i_tr in component sarcoplasmic_reticulum_kinetics (picoA).
C ALGBRC(49) is i_rel in component sarcoplasmic_reticulum_kinetics (picoA).
C CONSTS(34) is V_i in component ion_concentrations (microlitre).
C CONSTS(44) is V_rel in component sarcoplasmic_reticulum_kinetics (microlitre).
C CONSTS(46) is V_up in component sarcoplasmic_reticulum_kinetics (microlitre).
C CONSTS(35) is i_up_max in component sarcoplasmic_reticulum_kinetics (picoA).
C CONSTS(36) is KmCaup in component sarcoplasmic_reticulum_kinetics (millimolar).
C CONSTS(37) is KmCarel in component sarcoplasmic_reticulum_kinetics (millimolar).
C CONSTS(38) is tau_rel in component sarcoplasmic_reticulum_kinetics (second).
C CONSTS(39) is tau_tr in component sarcoplasmic_reticulum_kinetics (second).
C STATES(17) is Caup in component ion_concentrations (millimolar).
C STATES(18) is Carel in component ion_concentrations (millimolar).
C CONSTS(45) is V_e in component ion_concentrations (microlitre).
C CONSTS(40) is tau_b in component ion_concentrations (second).
C CONSTS(41) is Nab in component ion_concentrations (millimolar).
C CONSTS(42) is Cab in component ion_concentrations (millimolar).
C CONSTS(43) is Kb in component ion_concentrations (millimolar).
C RATES(1) is d/dt E in component membrane (millivolt).
C RATES(8) is d/dt dL in component L_type_calcium_current_d_gate (dimensionless).
C RATES(9) is d/dt fL in component L_type_calcium_current_f_gate (dimensionless).
C RATES(10) is d/dt fL2 in component L_type_calcium_current_f2_gate (dimensionless).
C RATES(11) is d/dt dT in component T_type_calcium_current_d_gate (dimensionless).
C RATES(12) is d/dt fT in component T_type_calcium_current_f_gate (dimensionless).
C RATES(13) is d/dt m in component fast_sodium_current_m_gate (dimensionless).
C RATES(14) is d/dt h in component fast_sodium_current_h_gate (dimensionless).
C RATES(15) is d/dt x in component delayed_rectifying_potassium_current_x_gate (dimensionless).
C RATES(16) is d/dt y in component hyperpolarising_activated_current_y_gate (dimensionless).
C RATES(4) is d/dt Nai in component ion_concentrations (millimolar).
C RATES(5) is d/dt Nao in component ion_concentrations (millimolar).
C RATES(6) is d/dt Ki in component ion_concentrations (millimolar).
C RATES(7) is d/dt Ko in component ion_concentrations (millimolar).
C RATES(2) is d/dt Cai in component ion_concentrations (millimolar).
C RATES(3) is d/dt Cao in component ion_concentrations (millimolar).
C RATES(17) is d/dt Caup in component ion_concentrations (millimolar).
C RATES(18) is d/dt Carel in component ion_concentrations (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 8314.472
      CONSTS(2) = 310
      CONSTS(3) = 96485.3415
      STATES(1) = -64.9
      CONSTS(4) = 3.2e-5
      STATES(2) = 0.000034
      STATES(3) = 2.0004
      STATES(4) = 7.4994
      STATES(5) = 139.9929
      STATES(6) = 140.0073
      STATES(7) = 5.4243
      CONSTS(5) = 0.4
      STATES(8) = 0.0001
      STATES(9) = 0.1505
      STATES(10) = 0.219
      CONSTS(6) = 0.002
      CONSTS(7) = 3
      CONSTS(8) = 40000
      CONSTS(9) = 0.085
      STATES(11) = 0.001
      STATES(12) = 0.1328
      CONSTS(10) = 0.25
      STATES(13) = 0.0139
      STATES(14) = 0.0087
      CONSTS(11) = 0.00026
      CONSTS(12) = 0.035
      STATES(15) = 0.5682
      CONSTS(13) = 10.3
      CONSTS(14) = 0.0081
      CONSTS(15) = 0.0135
      STATES(16) = 0.0287
      CONSTS(16) = 40
      CONSTS(17) = 1
      CONSTS(18) = 0.226
      CONSTS(19) = 4
      CONSTS(20) = 0.1369
      CONSTS(21) = 0.4315
      CONSTS(22) = 0
      CONSTS(23) = 26.44
      CONSTS(24) = 0.0207
      CONSTS(25) = 395.3
      CONSTS(26) = 2.289
      CONSTS(27) = 26.44
      CONSTS(28) = 4.663
      CONSTS(29) = 1628
      CONSTS(30) = 561.4
      CONSTS(31) = 3.663
      CONSTS(32) = 0.00024
      CONSTS(33) = 0.00007
      CONSTS(34) = 2.5e-6
      CONSTS(35) = 0.0212
      CONSTS(36) = 0.0005
      CONSTS(37) = 0.001
      CONSTS(38) = 0.005
      CONSTS(39) = 0.4
      STATES(17) = 0.5832
      STATES(18) = 0.1101
      CONSTS(40) = 0.1
      CONSTS(41) = 140
      CONSTS(42) = 2
      CONSTS(43) = 5.4
      CONSTS(44) =  0.00600000*CONSTS(34)
      CONSTS(45) =  0.200000*CONSTS(34)
      CONSTS(46) =  0.0140000*CONSTS(34)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(10) =  CONSTS(7)*(1.00000 - STATES(10)) -  CONSTS(8)*STATES(2)*STATES(10)
      ALGBRC(1) = 1.00000/(1.00000+EXP((STATES(1)+6.60000)/- 6.60000))
      RATES(8) = (ALGBRC(1) - STATES(8))/CONSTS(6)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+25.0000)/6.00000))
      ALGBRC(10) = 0.0310000+1.00000/(1.00000+EXP((STATES(1)+37.6000)/8.10000))
      RATES(9) = (ALGBRC(2) - STATES(9))/ALGBRC(10)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+23.0000)/- 6.10000))
      ALGBRC(11) = 0.000600000+0.00540000/(1.00000+EXP( 0.0300000*(STATES(1)+100.000)))
      RATES(11) = (ALGBRC(3) - STATES(11))/ALGBRC(11)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+75.0000)/6.60000))
      ALGBRC(12) = 0.00100000+0.0400000/(1.00000+EXP( 0.0800000*(STATES(1)+65.0000)))
      RATES(12) = (ALGBRC(4) - STATES(12))/ALGBRC(12)
      ALGBRC(5) = ( 200.000*(STATES(1)+34.3000))/(1.00000 - EXP( - 0.0900000*(STATES(1)+34.3000)))
      ALGBRC(13) =  8000.00*EXP( - 0.150000*(STATES(1)+56.2000))
      RATES(13) =  ALGBRC(5)*(1.00000 - STATES(13)) -  ALGBRC(13)*STATES(13)
      ALGBRC(6) =  32.4000*EXP( - 0.140000*(STATES(1)+93.4000))
      ALGBRC(14) = 709.000/(1.00000+ 4.20000*EXP( - 0.0600000*(STATES(1)+45.4000)))
      RATES(14) =  ALGBRC(6)*(1.00000 - STATES(14)) -  ALGBRC(14)*STATES(14)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+25.1000)/- 7.40000))
      ALGBRC(15) = 1.00000/( 17.0000*EXP( 0.0398000*STATES(1))+ 0.211000*EXP( - 0.0510000*STATES(1)))
      RATES(15) = (ALGBRC(7) - STATES(15))/ALGBRC(15)
      ALGBRC(8) = ( 0.360000*(STATES(1)+137.800))/(EXP( 0.0660000*(STATES(1)+137.800)) - 1.00000)
      ALGBRC(16) = ( 0.100000*(STATES(1)+76.3000))/(1.00000 - EXP( - 0.210000*(STATES(1)+76.3000)))
      RATES(16) =  ALGBRC(8)*(1.00000 - STATES(16)) -  ALGBRC(16)*STATES(16)
      ALGBRC(9) =  (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(STATES(3)/STATES(2))
      ALGBRC(19) =  CONSTS(5)*STATES(8)*STATES(9)*STATES(10)*((STATES(1) - ALGBRC(9))+75.0000)
      ALGBRC(20) =  CONSTS(9)*STATES(11)*STATES(12)*((STATES(1) - ALGBRC(9))+75.0000)
      ALGBRC(35) = EXP(( - CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(33) = STATES(5)/(CONSTS(28)+STATES(5))
      ALGBRC(29) = 1.00000+STATES(3)/CONSTS(31)+ (STATES(3)/CONSTS(31))*EXP(( CONSTS(22)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+STATES(5)/CONSTS(29)+STATES(5) ** 2.00000/( CONSTS(29)*CONSTS(30))+STATES(5) ** 3.00000/( CONSTS(29)*CONSTS(30)*CONSTS(28))
      ALGBRC(31) = ( (STATES(5) ** 2.00000/( CONSTS(29)*CONSTS(30))+STATES(5) ** 3.00000/( CONSTS(29)*CONSTS(30)*CONSTS(28)))*EXP(( - CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(29)
      ALGBRC(32) = ( (STATES(3)/CONSTS(31))*EXP(( - CONSTS(22)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(29)
      ALGBRC(30) = EXP(( CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(38) = STATES(4)/(CONSTS(23)+STATES(4))
      ALGBRC(39) =  ALGBRC(35)*ALGBRC(33)*(ALGBRC(31)+ALGBRC(32))+ ALGBRC(32)*ALGBRC(30)*(ALGBRC(38)+ALGBRC(35))
      ALGBRC(34) = 1.00000+STATES(2)/CONSTS(24)+ (STATES(2)/CONSTS(24))*EXP(( - CONSTS(20)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+( STATES(2)*STATES(4))/( CONSTS(24)*CONSTS(27))+STATES(4)/CONSTS(25)+STATES(4) ** 2.00000/( CONSTS(25)*CONSTS(26))+STATES(4) ** 3.00000/( CONSTS(25)*CONSTS(26)*CONSTS(23))
      ALGBRC(37) = ( (STATES(2)/CONSTS(24))*EXP(( - CONSTS(20)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(34)
      ALGBRC(36) = ( (STATES(4) ** 2.00000/( CONSTS(25)*CONSTS(26))+STATES(4) ** 3.00000/( CONSTS(25)*CONSTS(26)*CONSTS(23)))*EXP(( CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(34)
      ALGBRC(40) =  ALGBRC(30)*ALGBRC(38)*(ALGBRC(36)+ALGBRC(37))+ ALGBRC(35)*ALGBRC(37)*(ALGBRC(33)+ALGBRC(30))
      ALGBRC(41) =  ALGBRC(36)*ALGBRC(38)*(ALGBRC(31)+ALGBRC(32))+ ALGBRC(37)*ALGBRC(31)*(ALGBRC(38)+ALGBRC(35))
      ALGBRC(42) =  ALGBRC(31)*ALGBRC(33)*(ALGBRC(36)+ALGBRC(37))+ ALGBRC(36)*ALGBRC(32)*(ALGBRC(33)+ALGBRC(30))
      ALGBRC(43) = ( CONSTS(19)*( ALGBRC(40)*ALGBRC(32) -  ALGBRC(39)*ALGBRC(37)))/(ALGBRC(39)+ALGBRC(40)+ALGBRC(41)+ALGBRC(42))
      RATES(3) = ( 1.00000*((ALGBRC(19)+ALGBRC(20)) -  2.00000*ALGBRC(43)))/( 2.00000*CONSTS(3)*1.00000*CONSTS(45))+(CONSTS(42) - STATES(3))/CONSTS(40)
      ALGBRC(17) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(5)/STATES(4))
      ALGBRC(21) =  CONSTS(10)*STATES(13) ** 3.00000*STATES(14)*(STATES(1) - ALGBRC(17))
      ALGBRC(28) =  (( (( CONSTS(18)*STATES(4))/(STATES(4)+CONSTS(16)))*STATES(7))/(STATES(7)+CONSTS(17)))*(1.00000 - (STATES(1) - 40.0000)/211.000 ** 2.00000)
      ALGBRC(44) =  CONSTS(32)*(STATES(1) - ALGBRC(17))
      ALGBRC(23) =  STATES(15)*CONSTS(11)*CONSTS(12)*STATES(7)/1.00000 ** 0.590000*(STATES(4) -  STATES(5)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))
      ALGBRC(25) =  (( STATES(16)*STATES(7) ** 1.83000)/(STATES(7) ** 1.83000+CONSTS(13) ** 1.83000))*CONSTS(14)*(STATES(1) - ALGBRC(17))
      RATES(4) = ( - 1.00000*(ALGBRC(44)+ALGBRC(25)+ALGBRC(21)+ 3.00000*ALGBRC(28)+ 3.00000*ALGBRC(43)+ALGBRC(23)))/( CONSTS(3)*1.00000*CONSTS(34))
      RATES(5) = ( 1.00000*(ALGBRC(44)+ALGBRC(25)+ALGBRC(21)+ 3.00000*ALGBRC(28)+ 3.00000*ALGBRC(43)+ALGBRC(23)))/( CONSTS(3)*1.00000*CONSTS(45))+(CONSTS(41) - STATES(5))/CONSTS(40)
      ALGBRC(46) =  CONSTS(33)*STATES(7)/1.00000 ** 0.410000*(STATES(6) -  STATES(7)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))
      ALGBRC(22) =  STATES(15)*CONSTS(11)*STATES(7)/1.00000 ** 0.590000*(STATES(6) -  STATES(7)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))
      ALGBRC(18) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(7)/STATES(6))
      ALGBRC(26) =  (( STATES(16)*STATES(7) ** 1.83000)/(STATES(7) ** 1.83000+CONSTS(13) ** 1.83000))*CONSTS(15)*(STATES(1) - ALGBRC(18))
      RATES(6) = ( - 1.00000*(((ALGBRC(22)+ALGBRC(26)) -  2.00000*ALGBRC(28))+ALGBRC(46)))/( CONSTS(3)*1.00000*CONSTS(34))
      RATES(7) = ( 1.00000*(((ALGBRC(22)+ALGBRC(26)) -  2.00000*ALGBRC(28))+ALGBRC(46)))/( CONSTS(3)*1.00000*CONSTS(45))+(CONSTS(43) - STATES(7))/CONSTS(40)
      ALGBRC(45) = ( CONSTS(35)*STATES(2) ** 2.00000)/(STATES(2) ** 2.00000+CONSTS(36) ** 2.00000)
      ALGBRC(47) =  (( 2.00000*1.00000*CONSTS(44)*CONSTS(3))/( 1.00000*CONSTS(39)))*STATES(17)
      RATES(17) = ( 1.00000*(ALGBRC(45) - ALGBRC(47)))/( 2.00000*1.00000*CONSTS(46)*CONSTS(3))
      ALGBRC(24) = ALGBRC(22)+ALGBRC(23)
      ALGBRC(27) = ALGBRC(26)+ALGBRC(25)
      ALGBRC(48) = ALGBRC(19)+ALGBRC(20)+ALGBRC(21)+ALGBRC(24)+ALGBRC(27)+ALGBRC(28)+ALGBRC(43)+ALGBRC(44)+ALGBRC(46)
      RATES(1) = - ALGBRC(48)/CONSTS(4)
      ALGBRC(49) = ( (( 2.00000*1.00000*CONSTS(44)*CONSTS(3))/( 1.00000*CONSTS(38)))*STATES(18)*STATES(2) ** 2.00000)/(STATES(2) ** 2.00000+CONSTS(37) ** 2.00000)
      RATES(2) = ( - 1.00000*(((ALGBRC(19)+ALGBRC(20)) -  2.00000*ALGBRC(43))+ALGBRC(45)+- ALGBRC(49)))/( 2.00000*CONSTS(3)*1.00000*CONSTS(34))
      RATES(18) = ( 1.00000*(ALGBRC(47) - ALGBRC(49)))/( 2.00000*1.00000*CONSTS(44)*CONSTS(3))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = 1.00000/(1.00000+EXP((STATES(1)+6.60000)/- 6.60000))
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+25.0000)/6.00000))
      ALGBRC(10) = 0.0310000+1.00000/(1.00000+EXP((STATES(1)+37.6000)/8.10000))
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+23.0000)/- 6.10000))
      ALGBRC(11) = 0.000600000+0.00540000/(1.00000+EXP( 0.0300000*(STATES(1)+100.000)))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+75.0000)/6.60000))
      ALGBRC(12) = 0.00100000+0.0400000/(1.00000+EXP( 0.0800000*(STATES(1)+65.0000)))
      ALGBRC(5) = ( 200.000*(STATES(1)+34.3000))/(1.00000 - EXP( - 0.0900000*(STATES(1)+34.3000)))
      ALGBRC(13) =  8000.00*EXP( - 0.150000*(STATES(1)+56.2000))
      ALGBRC(6) =  32.4000*EXP( - 0.140000*(STATES(1)+93.4000))
      ALGBRC(14) = 709.000/(1.00000+ 4.20000*EXP( - 0.0600000*(STATES(1)+45.4000)))
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+25.1000)/- 7.40000))
      ALGBRC(15) = 1.00000/( 17.0000*EXP( 0.0398000*STATES(1))+ 0.211000*EXP( - 0.0510000*STATES(1)))
      ALGBRC(8) = ( 0.360000*(STATES(1)+137.800))/(EXP( 0.0660000*(STATES(1)+137.800)) - 1.00000)
      ALGBRC(16) = ( 0.100000*(STATES(1)+76.3000))/(1.00000 - EXP( - 0.210000*(STATES(1)+76.3000)))
      ALGBRC(9) =  (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(STATES(3)/STATES(2))
      ALGBRC(19) =  CONSTS(5)*STATES(8)*STATES(9)*STATES(10)*((STATES(1) - ALGBRC(9))+75.0000)
      ALGBRC(20) =  CONSTS(9)*STATES(11)*STATES(12)*((STATES(1) - ALGBRC(9))+75.0000)
      ALGBRC(35) = EXP(( - CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(33) = STATES(5)/(CONSTS(28)+STATES(5))
      ALGBRC(29) = 1.00000+STATES(3)/CONSTS(31)+ (STATES(3)/CONSTS(31))*EXP(( CONSTS(22)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+STATES(5)/CONSTS(29)+STATES(5) ** 2.00000/( CONSTS(29)*CONSTS(30))+STATES(5) ** 3.00000/( CONSTS(29)*CONSTS(30)*CONSTS(28))
      ALGBRC(31) = ( (STATES(5) ** 2.00000/( CONSTS(29)*CONSTS(30))+STATES(5) ** 3.00000/( CONSTS(29)*CONSTS(30)*CONSTS(28)))*EXP(( - CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(29)
      ALGBRC(32) = ( (STATES(3)/CONSTS(31))*EXP(( - CONSTS(22)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(29)
      ALGBRC(30) = EXP(( CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2)))
      ALGBRC(38) = STATES(4)/(CONSTS(23)+STATES(4))
      ALGBRC(39) =  ALGBRC(35)*ALGBRC(33)*(ALGBRC(31)+ALGBRC(32))+ ALGBRC(32)*ALGBRC(30)*(ALGBRC(38)+ALGBRC(35))
      ALGBRC(34) = 1.00000+STATES(2)/CONSTS(24)+ (STATES(2)/CONSTS(24))*EXP(( - CONSTS(20)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+( STATES(2)*STATES(4))/( CONSTS(24)*CONSTS(27))+STATES(4)/CONSTS(25)+STATES(4) ** 2.00000/( CONSTS(25)*CONSTS(26))+STATES(4) ** 3.00000/( CONSTS(25)*CONSTS(26)*CONSTS(23))
      ALGBRC(37) = ( (STATES(2)/CONSTS(24))*EXP(( - CONSTS(20)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))/ALGBRC(34)
      ALGBRC(36) = ( (STATES(4) ** 2.00000/( CONSTS(25)*CONSTS(26))+STATES(4) ** 3.00000/( CONSTS(25)*CONSTS(26)*CONSTS(23)))*EXP(( CONSTS(21)*STATES(1)*CONSTS(3))/( 2.00000*CONSTS(1)*CONSTS(2))))/ALGBRC(34)
      ALGBRC(40) =  ALGBRC(30)*ALGBRC(38)*(ALGBRC(36)+ALGBRC(37))+ ALGBRC(35)*ALGBRC(37)*(ALGBRC(33)+ALGBRC(30))
      ALGBRC(41) =  ALGBRC(36)*ALGBRC(38)*(ALGBRC(31)+ALGBRC(32))+ ALGBRC(37)*ALGBRC(31)*(ALGBRC(38)+ALGBRC(35))
      ALGBRC(42) =  ALGBRC(31)*ALGBRC(33)*(ALGBRC(36)+ALGBRC(37))+ ALGBRC(36)*ALGBRC(32)*(ALGBRC(33)+ALGBRC(30))
      ALGBRC(43) = ( CONSTS(19)*( ALGBRC(40)*ALGBRC(32) -  ALGBRC(39)*ALGBRC(37)))/(ALGBRC(39)+ALGBRC(40)+ALGBRC(41)+ALGBRC(42))
      ALGBRC(17) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(5)/STATES(4))
      ALGBRC(21) =  CONSTS(10)*STATES(13) ** 3.00000*STATES(14)*(STATES(1) - ALGBRC(17))
      ALGBRC(28) =  (( (( CONSTS(18)*STATES(4))/(STATES(4)+CONSTS(16)))*STATES(7))/(STATES(7)+CONSTS(17)))*(1.00000 - (STATES(1) - 40.0000)/211.000 ** 2.00000)
      ALGBRC(44) =  CONSTS(32)*(STATES(1) - ALGBRC(17))
      ALGBRC(23) =  STATES(15)*CONSTS(11)*CONSTS(12)*STATES(7)/1.00000 ** 0.590000*(STATES(4) -  STATES(5)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))
      ALGBRC(25) =  (( STATES(16)*STATES(7) ** 1.83000)/(STATES(7) ** 1.83000+CONSTS(13) ** 1.83000))*CONSTS(14)*(STATES(1) - ALGBRC(17))
      ALGBRC(46) =  CONSTS(33)*STATES(7)/1.00000 ** 0.410000*(STATES(6) -  STATES(7)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))
      ALGBRC(22) =  STATES(15)*CONSTS(11)*STATES(7)/1.00000 ** 0.590000*(STATES(6) -  STATES(7)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))))
      ALGBRC(18) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(STATES(7)/STATES(6))
      ALGBRC(26) =  (( STATES(16)*STATES(7) ** 1.83000)/(STATES(7) ** 1.83000+CONSTS(13) ** 1.83000))*CONSTS(15)*(STATES(1) - ALGBRC(18))
      ALGBRC(45) = ( CONSTS(35)*STATES(2) ** 2.00000)/(STATES(2) ** 2.00000+CONSTS(36) ** 2.00000)
      ALGBRC(47) =  (( 2.00000*1.00000*CONSTS(44)*CONSTS(3))/( 1.00000*CONSTS(39)))*STATES(17)
      ALGBRC(24) = ALGBRC(22)+ALGBRC(23)
      ALGBRC(27) = ALGBRC(26)+ALGBRC(25)
      ALGBRC(48) = ALGBRC(19)+ALGBRC(20)+ALGBRC(21)+ALGBRC(24)+ALGBRC(27)+ALGBRC(28)+ALGBRC(43)+ALGBRC(44)+ALGBRC(46)
      ALGBRC(49) = ( (( 2.00000*1.00000*CONSTS(44)*CONSTS(3))/( 1.00000*CONSTS(38)))*STATES(18)*STATES(2) ** 2.00000)/(STATES(2) ** 2.00000+CONSTS(37) ** 2.00000)
      RETURN
      END