C
C There are a total of 56 entries in the algebraic variable array.
C There are a total of 6 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 (ms).
C CONSTS(1) is V_myo in component cell_geometry (um3).
C CONSTS(2) is V_SR in component cell_geometry (um3).
C CONSTS(3) is A_cap in component cell_geometry (um2).
C ALGBRC(1) is V in component membrane (mV).
C CONSTS(4) is R in component membrane (mJ_per_mole_K).
C CONSTS(5) is T in component membrane (kelvin).
C CONSTS(6) is F in component membrane (C_per_mole).
C ALGBRC(4) is FVRT in component membrane (dimensionless).
C ALGBRC(5) is FVRT_Ca in component membrane (dimensionless).
C CONSTS(7) is g_D in component CaRU (um3_per_ms).
C CONSTS(8) is J_R in component CaRU (um3_per_ms).
C CONSTS(9) is J_L in component CaRU (um3_per_ms).
C CONSTS(10) is N in component CaRU (dimensionless).
C STATES(1) is Ca_i in component intracellular_ion_concentrations (mM).
C CONSTS(11) is Ca_o in component extracellular_ion_concentrations (mM).
C STATES(2) is Ca_SR in component intracellular_ion_concentrations (mM).
C ALGBRC(35) is I_RyR in component RyR_current (mM_per_ms).
C ALGBRC(44) is I_LCC in component LCC_current (mM_per_ms).
C ALGBRC(16) is C_oc in component DS_Calcium_Concentrations (mM).
C ALGBRC(14) is C_co in component DS_Calcium_Concentrations (mM).
C CONSTS(12) is V_L in component CaRU_Transitions (mV).
C CONSTS(13) is del_VL in component CaRU_Transitions (mV).
C CONSTS(14) is phi_L in component CaRU_Transitions (dimensionless).
C CONSTS(15) is t_L in component CaRU_Transitions (ms).
C CONSTS(16) is tau_L in component CaRU_Transitions (ms).
C CONSTS(44) is t_R in component CaRU_Transitions (ms).
C CONSTS(17) is tau_R in component CaRU_Transitions (ms).
C CONSTS(18) is phi_R in component CaRU_Transitions (dimensionless).
C CONSTS(19) is theta_R in component CaRU_Transitions (dimensionless).
C CONSTS(20) is K_RyR in component CaRU_Transitions (mM).
C CONSTS(21) is K_L in component CaRU_Transitions (mM).
C CONSTS(22) is a in component CaRU_Transitions (dimensionless).
C CONSTS(23) is b in component CaRU_Transitions (dimensionless).
C CONSTS(24) is c in component CaRU_Transitions (dimensionless).
C CONSTS(25) is d in component CaRU_Transitions (dimensionless).
C ALGBRC(6) is expVL in component CaRU_Transitions (dimensionless).
C ALGBRC(8) is alpha_p in component CaRU_Transitions (per_ms).
C CONSTS(45) is alpha_m in component CaRU_Transitions (per_ms).
C ALGBRC(17) is beta_poc in component CaRU_Transitions (per_ms).
C ALGBRC(9) is beta_pcc in component CaRU_Transitions (per_ms).
C CONSTS(46) is beta_m in component CaRU_Transitions (per_ms).
C ALGBRC(15) is epsilon_pco in component CaRU_Transitions (per_ms).
C ALGBRC(10) is epsilon_pcc in component CaRU_Transitions (per_ms).
C ALGBRC(11) is epsilon_m in component CaRU_Transitions (per_ms).
C ALGBRC(18) is mu_poc in component CaRU_Transitions (per_ms).
C ALGBRC(12) is mu_pcc in component CaRU_Transitions (per_ms).
C ALGBRC(19) is mu_moc in component CaRU_Transitions (per_ms).
C ALGBRC(13) is mu_mcc in component CaRU_Transitions (per_ms).
C ALGBRC(2) is C_cc in component DS_Calcium_Concentrations (mM).
C ALGBRC(7) is C_oo in component DS_Calcium_Concentrations (mM).
C ALGBRC(23) is J_Loo in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(22) is J_Loc in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(20) is J_Rco in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(21) is J_Roo in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(24) is denom in component CaRU_states (per_ms3).
C ALGBRC(25) is y_oc in component CaRU_states (dimensionless).
C ALGBRC(26) is y_co in component CaRU_states (dimensionless).
C ALGBRC(27) is y_oo in component CaRU_states (dimensionless).
C ALGBRC(28) is y_cc in component CaRU_states (dimensionless).
C ALGBRC(30) is y_ci in component CaRU_states (dimensionless).
C ALGBRC(33) is y_oi in component CaRU_states (dimensionless).
C ALGBRC(36) is y_ic in component CaRU_states (dimensionless).
C ALGBRC(39) is y_io in component CaRU_states (dimensionless).
C ALGBRC(42) is y_ii in component CaRU_states (dimensionless).
C ALGBRC(31) is r_1 in component CaRU_reduced_states (per_ms).
C ALGBRC(34) is r_2 in component CaRU_reduced_states (per_ms).
C ALGBRC(37) is r_3 in component CaRU_reduced_states (per_ms).
C ALGBRC(40) is r_4 in component CaRU_reduced_states (per_ms).
C ALGBRC(43) is r_5 in component CaRU_reduced_states (per_ms).
C ALGBRC(45) is r_6 in component CaRU_reduced_states (per_ms).
C ALGBRC(47) is r_7 in component CaRU_reduced_states (per_ms).
C ALGBRC(49) is r_8 in component CaRU_reduced_states (per_ms).
C STATES(3) is z_1 in component CaRU_reduced_states (dimensionless).
C STATES(4) is z_2 in component CaRU_reduced_states (dimensionless).
C STATES(5) is z_3 in component CaRU_reduced_states (dimensionless).
C ALGBRC(51) is z_4 in component CaRU_reduced_states (dimensionless).
C ALGBRC(29) is J_R1 in component RyR_current (um3_mM_per_ms).
C ALGBRC(32) is J_R3 in component RyR_current (um3_mM_per_ms).
C ALGBRC(38) is J_L1 in component LCC_current (um3_mM_per_ms).
C ALGBRC(41) is J_L2 in component LCC_current (um3_mM_per_ms).
C CONSTS(26) is K_mNa in component Na_Ca_Exchanger (mM).
C CONSTS(27) is K_mCa in component Na_Ca_Exchanger (mM).
C CONSTS(28) is eta in component Na_Ca_Exchanger (dimensionless).
C CONSTS(29) is k_sat in component Na_Ca_Exchanger (dimensionless).
C CONSTS(30) is g_NCX in component Na_Ca_Exchanger (mM_per_ms).
C CONSTS(31) is Na_i in component intracellular_ion_concentrations (mM).
C CONSTS(32) is Na_o in component extracellular_ion_concentrations (mM).
C ALGBRC(46) is I_NaCa in component Na_Ca_Exchanger (mM_per_ms).
C CONSTS(33) is g_SERCA in component SERCA (mM_per_ms).
C CONSTS(34) is K_SERCA in component SERCA (mM).
C ALGBRC(48) is I_SERCA in component SERCA (mM_per_ms).
C CONSTS(35) is g_pCa in component Sarcolemmal_Ca_pump (mM_per_ms).
C CONSTS(36) is K_mpCa in component Sarcolemmal_Ca_pump (mM).
C ALGBRC(50) is I_pCa in component Sarcolemmal_Ca_pump (mM_per_ms).
C ALGBRC(52) is E_Ca in component Background_Ca_current (mV).
C CONSTS(37) is g_CaB in component Background_Ca_current (mM_per_mV_ms).
C ALGBRC(53) is I_CaB in component Background_Ca_current (mM_per_ms).
C CONSTS(38) is g_SRl in component SR_Ca_leak_current (per_ms).
C ALGBRC(54) is I_SR in component SR_Ca_leak_current (mM_per_ms).
C CONSTS(39) is k_m_TRPN in component troponin_Ca_buffer (per_ms).
C CONSTS(40) is k_p_TRPN in component troponin_Ca_buffer (per_mM_ms).
C CONSTS(41) is B_TRPN in component troponin_Ca_buffer (mM).
C STATES(6) is TRPN in component intracellular_ion_concentrations (mM).
C ALGBRC(55) is I_TRPN in component troponin_Ca_buffer (mM_per_ms).
C CONSTS(42) is k_CMDN in component calmodulin_Ca_buffer (mM).
C CONSTS(43) is B_CMDN in component calmodulin_Ca_buffer (mM).
C ALGBRC(56) is beta_CMDN in component calmodulin_Ca_buffer (dimensionless).
C ALGBRC(3) is CaSR_plot in component intracellular_ion_concentrations (mM).
C RATES(3) is d/dt z_1 in component CaRU_reduced_states (dimensionless).
C RATES(4) is d/dt z_2 in component CaRU_reduced_states (dimensionless).
C RATES(5) is d/dt z_3 in component CaRU_reduced_states (dimensionless).
C RATES(6) is d/dt TRPN in component intracellular_ion_concentrations (mM).
C RATES(1) is d/dt Ca_i in component intracellular_ion_concentrations (mM).
C RATES(2) is d/dt Ca_SR in component intracellular_ion_concentrations (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 25.84e3
      CONSTS(2) = 2.098e3
      CONSTS(3) = 1.534e4
      CONSTS(4) = 8314.5
      CONSTS(5) = 295
      CONSTS(6) = 96487
      CONSTS(7) = 0.065
      CONSTS(8) = 0.02
      CONSTS(9) = 9.13e-4
      CONSTS(10) = 50000
      STATES(1) = 0.0001
      CONSTS(11) = 1
      STATES(2) = 700e-3
      CONSTS(12) = -2
      CONSTS(13) = 7
      CONSTS(14) = 2.35
      CONSTS(15) = 1
      CONSTS(16) = 650
      CONSTS(17) = 2.43
      CONSTS(18) = 0.05
      CONSTS(19) = 0.012
      CONSTS(20) = 41e-3
      CONSTS(21) = 0.22e-3
      CONSTS(22) = 0.0625
      CONSTS(23) = 14
      CONSTS(24) = 0.01
      CONSTS(25) = 100
      STATES(3) = 0.98859
      STATES(4) = 0.0087302
      STATES(5) = 0.0026566
      CONSTS(26) = 87.5
      CONSTS(27) = 1.38
      CONSTS(28) = 0.35
      CONSTS(29) = 0.1
      CONSTS(30) = 38.5e-3
      CONSTS(31) = 10
      CONSTS(32) = 140
      CONSTS(33) = 0.45e-3
      CONSTS(34) = 0.5e-3
      CONSTS(35) = 0.0035e-3
      CONSTS(36) = 0.5e-3
      CONSTS(37) = 2.6875e-8
      CONSTS(38) = 1.8951e-5
      CONSTS(39) = 0.04
      CONSTS(40) = 0.04e3
      CONSTS(41) = 70e-3
      STATES(6) = 0.0636364
      CONSTS(42) = 2.382e-3
      CONSTS(43) = 50e-3
      CONSTS(44) =  1.17000*CONSTS(15)
      CONSTS(45) = CONSTS(14)/CONSTS(15)
      CONSTS(46) = CONSTS(18)/CONSTS(44)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(VOI.GE.0.00000.AND.VOI.LE.200.000, 0.00000, - 80.0000)
      ALGBRC(4) = ( CONSTS(6)*ALGBRC(1))/( CONSTS(4)*CONSTS(5))
      ALGBRC(5) =  2.00000*ALGBRC(4)
      ALGBRC(16) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-09, (STATES(1)+( (CONSTS(9)/CONSTS(7))*CONSTS(11)*ALGBRC(5)*EXP(- ALGBRC(5)))/(1.00000 - EXP(- ALGBRC(5))))/(1.00000+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), (STATES(1)+ (CONSTS(9)/CONSTS(7))*CONSTS(11))/(1.00000+CONSTS(9)/CONSTS(7)))
      ALGBRC(18) = (ALGBRC(16) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000)/( CONSTS(17)*(ALGBRC(16) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(12) = (STATES(1) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000)/( CONSTS(17)*(STATES(1) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(6) = EXP((ALGBRC(1) - CONSTS(12))/CONSTS(13))
      ALGBRC(8) = ALGBRC(6)/( CONSTS(15)*(ALGBRC(6)+1.00000))
      ALGBRC(9) = STATES(1) ** 2.00000/( CONSTS(44)*(STATES(1) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(17) = ALGBRC(16) ** 2.00000/( CONSTS(44)*(ALGBRC(16) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(24) =  (ALGBRC(8)+CONSTS(45))*( (CONSTS(45)+CONSTS(46)+ALGBRC(17))*(CONSTS(46)+ALGBRC(9))+ ALGBRC(8)*(CONSTS(46)+ALGBRC(17)))
      ALGBRC(25) = ( ALGBRC(8)*CONSTS(46)*(ALGBRC(8)+CONSTS(45)+CONSTS(46)+ALGBRC(9)))/ALGBRC(24)
      ALGBRC(28) = ( CONSTS(45)*CONSTS(46)*(CONSTS(45)+ALGBRC(8)+CONSTS(46)+ALGBRC(17)))/ALGBRC(24)
      ALGBRC(31) =  ALGBRC(25)*ALGBRC(18)+ ALGBRC(28)*ALGBRC(12)
      ALGBRC(19) = ( CONSTS(19)*CONSTS(25)*(ALGBRC(16) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))/( CONSTS(17)*( CONSTS(25)*ALGBRC(16) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))
      ALGBRC(13) = ( CONSTS(19)*CONSTS(25)*(STATES(1) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))/( CONSTS(17)*( CONSTS(25)*STATES(1) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))
      ALGBRC(34) = ( ALGBRC(8)*ALGBRC(19)+ CONSTS(45)*ALGBRC(13))/(ALGBRC(8)+CONSTS(45))
      ALGBRC(14) = (STATES(1)+ (CONSTS(8)/CONSTS(7))*STATES(2))/(1.00000+CONSTS(8)/CONSTS(7))
      ALGBRC(15) = ( ALGBRC(14)*(ALGBRC(6)+CONSTS(22)))/( CONSTS(16)*CONSTS(21)*(ALGBRC(6)+1.00000))
      ALGBRC(10) = ( STATES(1)*(ALGBRC(6)+CONSTS(22)))/( CONSTS(16)*CONSTS(21)*(ALGBRC(6)+1.00000))
      ALGBRC(26) = ( CONSTS(45)*( ALGBRC(9)*(CONSTS(45)+CONSTS(46)+ALGBRC(17))+ ALGBRC(17)*ALGBRC(8)))/ALGBRC(24)
      ALGBRC(43) =  ALGBRC(26)*ALGBRC(15)+ ALGBRC(28)*ALGBRC(10)
      ALGBRC(11) = ( CONSTS(23)*(ALGBRC(6)+CONSTS(22)))/( CONSTS(16)*( CONSTS(23)*ALGBRC(6)+CONSTS(22)))
      ALGBRC(45) = ALGBRC(11)
      RATES(3) =  - (ALGBRC(31)+ALGBRC(43))*STATES(3)+ ALGBRC(34)*STATES(4)+ ALGBRC(45)*STATES(5)
      ALGBRC(47) = ( CONSTS(45)*ALGBRC(10))/(ALGBRC(8)+CONSTS(45))
      ALGBRC(49) = ALGBRC(11)
      ALGBRC(51) = ((1.00000 - STATES(3)) - STATES(4)) - STATES(5)
      RATES(4) = ( ALGBRC(31)*STATES(3) -  (ALGBRC(34)+ALGBRC(47))*STATES(4))+ ALGBRC(49)*ALGBRC(51)
      ALGBRC(37) = ( CONSTS(46)*ALGBRC(12))/(CONSTS(46)+ALGBRC(9))
      ALGBRC(40) = ALGBRC(13)
      RATES(5) = ( ALGBRC(43)*STATES(3) -  (ALGBRC(45)+ALGBRC(37))*STATES(5))+ ALGBRC(40)*ALGBRC(51)
      ALGBRC(20) = ( CONSTS(8)*(STATES(2) - STATES(1)))/(1.00000+CONSTS(8)/CONSTS(7))
      ALGBRC(21) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-05, ( CONSTS(8)*((STATES(2) - STATES(1))+ (( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5))))*(STATES(2) -  CONSTS(11)*EXP(- ALGBRC(5)))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), ( CONSTS(8)*((STATES(2) - STATES(1))+ (( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(STATES(2) -  CONSTS(11)*EXP(- 1.00000e-05))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
      ALGBRC(27) = ( ALGBRC(8)*( ALGBRC(17)*(ALGBRC(8)+CONSTS(46)+ALGBRC(9))+ ALGBRC(9)*CONSTS(45)))/ALGBRC(24)
      ALGBRC(29) =  ALGBRC(27)*ALGBRC(21)+ ALGBRC(20)*ALGBRC(26)
      ALGBRC(32) = ( ALGBRC(20)*ALGBRC(9))/(CONSTS(46)+ALGBRC(9))
      ALGBRC(35) = ( ( STATES(3)*ALGBRC(29)+ STATES(5)*ALGBRC(32))*CONSTS(10))/CONSTS(1)
      ALGBRC(48) = ( CONSTS(33)*STATES(1) ** 2.00000)/(CONSTS(34) ** 2.00000+STATES(1) ** 2.00000)
      ALGBRC(54) =  CONSTS(38)*(STATES(2) - STATES(1))
      RATES(2) =  (CONSTS(1)/CONSTS(2))*((- ALGBRC(35)+ALGBRC(48)) - ALGBRC(54))
      ALGBRC(55) =  CONSTS(39)*(CONSTS(41) - STATES(6)) -  CONSTS(40)*STATES(6)*STATES(1)
      RATES(6) = ALGBRC(55)
      ALGBRC(23) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-05, ( (( CONSTS(9)*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5))))*(( CONSTS(11)*EXP(- ALGBRC(5)) - STATES(1))+ (CONSTS(8)/CONSTS(7))*( CONSTS(11)*EXP(- ALGBRC(5)) - STATES(2))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(ALGBRC(5)))), ( (( CONSTS(9)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(( CONSTS(11)*EXP(- 1.00000e-05) - STATES(1))+ (CONSTS(8)/CONSTS(7))*( CONSTS(11)*EXP(- 1.00000e-05) - STATES(2))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
      ALGBRC(22) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-05, ( (( CONSTS(9)*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5))))*( CONSTS(11)*EXP(- ALGBRC(5)) - STATES(1)))/(1.00000+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), ( (( CONSTS(9)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*( CONSTS(11)*EXP(- 1.00000e-05) - STATES(1)))/(1.00000+( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
      ALGBRC(38) =  ALGBRC(23)*ALGBRC(27)+ ALGBRC(22)*ALGBRC(25)
      ALGBRC(41) = ( ALGBRC(22)*ALGBRC(8))/(ALGBRC(8)+CONSTS(45))
      ALGBRC(44) = ( ( STATES(3)*ALGBRC(38)+ STATES(4)*ALGBRC(41))*CONSTS(10))/CONSTS(1)
      ALGBRC(46) = ( CONSTS(30)*( EXP( CONSTS(28)*ALGBRC(4))*CONSTS(31) ** 3.00000*CONSTS(11) -  EXP( (CONSTS(28) - 1.00000)*ALGBRC(4))*CONSTS(32) ** 3.00000*STATES(1)))/( (CONSTS(32) ** 3.00000+CONSTS(26) ** 3.00000)*(CONSTS(11)+CONSTS(27))*(1.00000+ CONSTS(29)*EXP( (CONSTS(28) - 1.00000)*ALGBRC(4))))
      ALGBRC(50) = ( CONSTS(35)*STATES(1))/(CONSTS(36)+STATES(1))
      ALGBRC(52) =  (( CONSTS(4)*CONSTS(5))/( 2.00000*CONSTS(6)))*log(CONSTS(11)/STATES(1))
      ALGBRC(53) =  CONSTS(37)*(ALGBRC(52) - ALGBRC(1))
      ALGBRC(56) = 1.00000+( CONSTS(42)*CONSTS(43))/CONSTS(42)+STATES(1) ** 2.00000 ** - 1.00000
      RATES(1) =  ALGBRC(56)*(((((ALGBRC(44)+ALGBRC(35)) - ALGBRC(48))+ALGBRC(54)+ALGBRC(46)) - ALGBRC(50))+ALGBRC(53)+ALGBRC(55))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(VOI.GE.0.00000.AND.VOI.LE.200.000, 0.00000, - 80.0000)
      ALGBRC(4) = ( CONSTS(6)*ALGBRC(1))/( CONSTS(4)*CONSTS(5))
      ALGBRC(5) =  2.00000*ALGBRC(4)
      ALGBRC(16) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-09, (STATES(1)+( (CONSTS(9)/CONSTS(7))*CONSTS(11)*ALGBRC(5)*EXP(- ALGBRC(5)))/(1.00000 - EXP(- ALGBRC(5))))/(1.00000+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), (STATES(1)+ (CONSTS(9)/CONSTS(7))*CONSTS(11))/(1.00000+CONSTS(9)/CONSTS(7)))
      ALGBRC(18) = (ALGBRC(16) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000)/( CONSTS(17)*(ALGBRC(16) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(12) = (STATES(1) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000)/( CONSTS(17)*(STATES(1) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(6) = EXP((ALGBRC(1) - CONSTS(12))/CONSTS(13))
      ALGBRC(8) = ALGBRC(6)/( CONSTS(15)*(ALGBRC(6)+1.00000))
      ALGBRC(9) = STATES(1) ** 2.00000/( CONSTS(44)*(STATES(1) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(17) = ALGBRC(16) ** 2.00000/( CONSTS(44)*(ALGBRC(16) ** 2.00000+CONSTS(20) ** 2.00000))
      ALGBRC(24) =  (ALGBRC(8)+CONSTS(45))*( (CONSTS(45)+CONSTS(46)+ALGBRC(17))*(CONSTS(46)+ALGBRC(9))+ ALGBRC(8)*(CONSTS(46)+ALGBRC(17)))
      ALGBRC(25) = ( ALGBRC(8)*CONSTS(46)*(ALGBRC(8)+CONSTS(45)+CONSTS(46)+ALGBRC(9)))/ALGBRC(24)
      ALGBRC(28) = ( CONSTS(45)*CONSTS(46)*(CONSTS(45)+ALGBRC(8)+CONSTS(46)+ALGBRC(17)))/ALGBRC(24)
      ALGBRC(31) =  ALGBRC(25)*ALGBRC(18)+ ALGBRC(28)*ALGBRC(12)
      ALGBRC(19) = ( CONSTS(19)*CONSTS(25)*(ALGBRC(16) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))/( CONSTS(17)*( CONSTS(25)*ALGBRC(16) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))
      ALGBRC(13) = ( CONSTS(19)*CONSTS(25)*(STATES(1) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))/( CONSTS(17)*( CONSTS(25)*STATES(1) ** 2.00000+ CONSTS(24)*CONSTS(20) ** 2.00000))
      ALGBRC(34) = ( ALGBRC(8)*ALGBRC(19)+ CONSTS(45)*ALGBRC(13))/(ALGBRC(8)+CONSTS(45))
      ALGBRC(14) = (STATES(1)+ (CONSTS(8)/CONSTS(7))*STATES(2))/(1.00000+CONSTS(8)/CONSTS(7))
      ALGBRC(15) = ( ALGBRC(14)*(ALGBRC(6)+CONSTS(22)))/( CONSTS(16)*CONSTS(21)*(ALGBRC(6)+1.00000))
      ALGBRC(10) = ( STATES(1)*(ALGBRC(6)+CONSTS(22)))/( CONSTS(16)*CONSTS(21)*(ALGBRC(6)+1.00000))
      ALGBRC(26) = ( CONSTS(45)*( ALGBRC(9)*(CONSTS(45)+CONSTS(46)+ALGBRC(17))+ ALGBRC(17)*ALGBRC(8)))/ALGBRC(24)
      ALGBRC(43) =  ALGBRC(26)*ALGBRC(15)+ ALGBRC(28)*ALGBRC(10)
      ALGBRC(11) = ( CONSTS(23)*(ALGBRC(6)+CONSTS(22)))/( CONSTS(16)*( CONSTS(23)*ALGBRC(6)+CONSTS(22)))
      ALGBRC(45) = ALGBRC(11)
      ALGBRC(47) = ( CONSTS(45)*ALGBRC(10))/(ALGBRC(8)+CONSTS(45))
      ALGBRC(49) = ALGBRC(11)
      ALGBRC(51) = ((1.00000 - STATES(3)) - STATES(4)) - STATES(5)
      ALGBRC(37) = ( CONSTS(46)*ALGBRC(12))/(CONSTS(46)+ALGBRC(9))
      ALGBRC(40) = ALGBRC(13)
      ALGBRC(20) = ( CONSTS(8)*(STATES(2) - STATES(1)))/(1.00000+CONSTS(8)/CONSTS(7))
      ALGBRC(21) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-05, ( CONSTS(8)*((STATES(2) - STATES(1))+ (( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5))))*(STATES(2) -  CONSTS(11)*EXP(- ALGBRC(5)))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), ( CONSTS(8)*((STATES(2) - STATES(1))+ (( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(STATES(2) -  CONSTS(11)*EXP(- 1.00000e-05))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
      ALGBRC(27) = ( ALGBRC(8)*( ALGBRC(17)*(ALGBRC(8)+CONSTS(46)+ALGBRC(9))+ ALGBRC(9)*CONSTS(45)))/ALGBRC(24)
      ALGBRC(29) =  ALGBRC(27)*ALGBRC(21)+ ALGBRC(20)*ALGBRC(26)
      ALGBRC(32) = ( ALGBRC(20)*ALGBRC(9))/(CONSTS(46)+ALGBRC(9))
      ALGBRC(35) = ( ( STATES(3)*ALGBRC(29)+ STATES(5)*ALGBRC(32))*CONSTS(10))/CONSTS(1)
      ALGBRC(48) = ( CONSTS(33)*STATES(1) ** 2.00000)/(CONSTS(34) ** 2.00000+STATES(1) ** 2.00000)
      ALGBRC(54) =  CONSTS(38)*(STATES(2) - STATES(1))
      ALGBRC(55) =  CONSTS(39)*(CONSTS(41) - STATES(6)) -  CONSTS(40)*STATES(6)*STATES(1)
      ALGBRC(23) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-05, ( (( CONSTS(9)*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5))))*(( CONSTS(11)*EXP(- ALGBRC(5)) - STATES(1))+ (CONSTS(8)/CONSTS(7))*( CONSTS(11)*EXP(- ALGBRC(5)) - STATES(2))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(ALGBRC(5)))), ( (( CONSTS(9)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(( CONSTS(11)*EXP(- 1.00000e-05) - STATES(1))+ (CONSTS(8)/CONSTS(7))*( CONSTS(11)*EXP(- 1.00000e-05) - STATES(2))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
      ALGBRC(22) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-05, ( (( CONSTS(9)*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5))))*( CONSTS(11)*EXP(- ALGBRC(5)) - STATES(1)))/(1.00000+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), ( (( CONSTS(9)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*( CONSTS(11)*EXP(- 1.00000e-05) - STATES(1)))/(1.00000+( (CONSTS(9)/CONSTS(7))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
      ALGBRC(38) =  ALGBRC(23)*ALGBRC(27)+ ALGBRC(22)*ALGBRC(25)
      ALGBRC(41) = ( ALGBRC(22)*ALGBRC(8))/(ALGBRC(8)+CONSTS(45))
      ALGBRC(44) = ( ( STATES(3)*ALGBRC(38)+ STATES(4)*ALGBRC(41))*CONSTS(10))/CONSTS(1)
      ALGBRC(46) = ( CONSTS(30)*( EXP( CONSTS(28)*ALGBRC(4))*CONSTS(31) ** 3.00000*CONSTS(11) -  EXP( (CONSTS(28) - 1.00000)*ALGBRC(4))*CONSTS(32) ** 3.00000*STATES(1)))/( (CONSTS(32) ** 3.00000+CONSTS(26) ** 3.00000)*(CONSTS(11)+CONSTS(27))*(1.00000+ CONSTS(29)*EXP( (CONSTS(28) - 1.00000)*ALGBRC(4))))
      ALGBRC(50) = ( CONSTS(35)*STATES(1))/(CONSTS(36)+STATES(1))
      ALGBRC(52) =  (( CONSTS(4)*CONSTS(5))/( 2.00000*CONSTS(6)))*log(CONSTS(11)/STATES(1))
      ALGBRC(53) =  CONSTS(37)*(ALGBRC(52) - ALGBRC(1))
      ALGBRC(56) = 1.00000+( CONSTS(42)*CONSTS(43))/CONSTS(42)+STATES(1) ** 2.00000 ** - 1.00000
      ALGBRC(2) = STATES(1)
      ALGBRC(3) = ( STATES(2)*CONSTS(2))/CONSTS(1)
      ALGBRC(7) = TERNRY(ABS(ALGBRC(5)).GT.1.00000e-09, (STATES(1)+ (CONSTS(8)/CONSTS(7))*STATES(2)+( (CONSTS(9)/CONSTS(7))*CONSTS(11)*ALGBRC(5)*EXP(- ALGBRC(5)))/(1.00000 - EXP(- ALGBRC(5))))/(1.00000+CONSTS(8)/CONSTS(7)+( (CONSTS(9)/CONSTS(7))*ALGBRC(5))/(1.00000 - EXP(- ALGBRC(5)))), (STATES(1)+ (CONSTS(8)/CONSTS(7))*STATES(2)+ (CONSTS(9)/CONSTS(7))*CONSTS(11))/(1.00000+CONSTS(8)/CONSTS(7)+CONSTS(9)/CONSTS(7)))
      ALGBRC(30) = CONSTS(45)/(ALGBRC(8)+CONSTS(45))
      ALGBRC(33) = ALGBRC(8)/(ALGBRC(8)+CONSTS(45))
      ALGBRC(36) = CONSTS(46)/(ALGBRC(9)+CONSTS(46))
      ALGBRC(39) = ALGBRC(9)/(ALGBRC(9)+CONSTS(46))
      ALGBRC(42) = (((((((1.00000 - ALGBRC(25)) - ALGBRC(26)) - ALGBRC(27)) - ALGBRC(28)) - ALGBRC(30)) - ALGBRC(36)) - ALGBRC(33)) - ALGBRC(39)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END