C
C There are a total of 45 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 64 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C ALGBRC(4) is SOVFThick in component sarcomere_geometry (dimensionless).
C ALGBRC(5) is SOVFThin in component sarcomere_geometry (dimensionless).
C ALGBRC(1) is sovr_ze in component sarcomere_geometry (micrometre).
C ALGBRC(2) is sovr_cle in component sarcomere_geometry (micrometre).
C ALGBRC(3) is len_sovr in component sarcomere_geometry (micrometre).
C CONSTS(1) is SLmax in component normalised_active_and_passive_force (micrometre).
C CONSTS(2) is SLmin in component normalised_active_and_passive_force (micrometre).
C CONSTS(3) is len_thin in component model_parameters (micrometre).
C CONSTS(4) is len_thick in component model_parameters (micrometre).
C CONSTS(5) is len_hbare in component model_parameters (micrometre).
C STATES(1) is SL in component normalised_active_and_passive_force (micrometre).
C STATES(2) is TRPNCaL in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C STATES(3) is TRPNCaH in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C ALGBRC(43) is dTRPNCaL in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C ALGBRC(44) is dTRPNCaH in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C ALGBRC(11) is kn_pT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C ALGBRC(15) is kp_nT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(56) is konT in component Ca_binding_to_troponin_to_thin_filament_regulation (second_order_rate_constant).
C CONSTS(57) is koffLT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(58) is koffHT in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(6) is Qkon in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(7) is Qkoff in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(8) is Qkn_p in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(9) is Qkp_n in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(10) is kon in component Ca_binding_to_troponin_to_thin_filament_regulation (second_order_rate_constant).
C CONSTS(11) is koffL in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(12) is koffH in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(13) is perm50 in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(14) is nperm in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(15) is kn_p in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(16) is kp_n in component Ca_binding_to_troponin_to_thin_filament_regulation (first_order_rate_constant).
C CONSTS(17) is koffmod in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C ALGBRC(7) is Tropreg in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C ALGBRC(9) is permtot in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C ALGBRC(13) is inprmt in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C CONSTS(18) is TmpC in component model_parameters (celsius).
C ALGBRC(41) is Cai in component equation_for_simulated_calcium_transient (micromolar).
C CONSTS(59) is fappT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C ALGBRC(18) is gappT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C ALGBRC(21) is hfT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C ALGBRC(22) is hbT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C ALGBRC(24) is gxbT in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C CONSTS(19) is fapp in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C CONSTS(20) is gapp in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C CONSTS(21) is hf in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C CONSTS(22) is hb in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C CONSTS(23) is gxb in component thin_filament_regulation_and_crossbridge_cycling_rates (first_order_rate_constant).
C CONSTS(24) is gslmod in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C ALGBRC(19) is hfmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C ALGBRC(20) is hbmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(25) is hfmdc in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(26) is hbmdc in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(27) is sigmap in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(28) is sigman in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(29) is xbmodsp in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(30) is Qfapp in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(31) is Qgapp in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(32) is Qhf in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(33) is Qhb in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(34) is Qgxb in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C ALGBRC(23) is gxbmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C ALGBRC(17) is gapslmd in component thin_filament_regulation_and_crossbridge_cycling_rates (dimensionless).
C CONSTS(35) is x_0 in component model_parameters (micrometre).
C STATES(4) is xXBpostr in component mean_strain_of_strongly_bound_states (micrometre).
C STATES(5) is xXBprer in component mean_strain_of_strongly_bound_states (micrometre).
C STATES(6) is XBpostr in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C STATES(7) is XBprer in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C ALGBRC(25) is dXBpostr in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
C ALGBRC(27) is dXBprer in component regulation_and_crossbridge_cycling_state_equations (first_order_rate_constant).
C STATES(8) is N_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C STATES(9) is P_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C ALGBRC(26) is P in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C STATES(10) is N in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C ALGBRC(32) is dxXBpostr in component mean_strain_of_strongly_bound_states (micrometre_per_millisecond).
C ALGBRC(31) is dxXBprer in component mean_strain_of_strongly_bound_states (micrometre_per_millisecond).
C CONSTS(36) is xPsi in component mean_strain_of_strongly_bound_states (dimensionless).
C ALGBRC(28) is dutyprer in component mean_strain_of_strongly_bound_states (dimensionless).
C ALGBRC(29) is dutypostr in component mean_strain_of_strongly_bound_states (dimensionless).
C ALGBRC(30) is dSL in component normalised_active_and_passive_force (micrometre_per_millisecond).
C CONSTS(62) is SSXBpostr in component normalised_active_and_passive_force (dimensionless).
C CONSTS(60) is SSXBprer in component normalised_active_and_passive_force (dimensionless).
C CONSTS(37) is kxb in component normalised_active_and_passive_force (millinewton_per_millimetre2).
C CONSTS(63) is Fnordv in component normalised_active_and_passive_force (millinewton_micrometre_per_millimetre2).
C ALGBRC(6) is force in component normalised_active_and_passive_force (millinewton_micrometre_per_millimetre2).
C ALGBRC(8) is active in component normalised_active_and_passive_force (unit_normalised_force).
C ALGBRC(14) is ppforce in component normalised_active_and_passive_force (unit_normalised_force).
C ALGBRC(10) is ppforce_t in component normalised_active_and_passive_force (unit_normalised_force).
C ALGBRC(12) is ppforce_c in component normalised_active_and_passive_force (unit_normalised_force).
C CONSTS(64) is preload in component normalised_active_and_passive_force (unit_normalised_force).
C ALGBRC(16) is afterload in component normalised_active_and_passive_force (unit_normalised_force).
C STATES(11) is intf in component normalised_active_and_passive_force (unit_normalised_force_millisecond).
C CONSTS(38) is SL_c in component normalised_active_and_passive_force (micrometre).
C CONSTS(39) is SLrest in component normalised_active_and_passive_force (micrometre).
C CONSTS(40) is SLset in component normalised_active_and_passive_force (micrometre).
C CONSTS(41) is PCon_t in component normalised_active_and_passive_force (unit_normalised_force).
C CONSTS(42) is PExp_t in component normalised_active_and_passive_force (per_micrometre).
C CONSTS(43) is PCon_c in component normalised_active_and_passive_force (unit_normalised_force).
C CONSTS(44) is PExp_c in component normalised_active_and_passive_force (per_micrometre).
C CONSTS(45) is massf in component normalised_active_and_passive_force (unit_normalised_force_millisecond2_per_micrometre).
C CONSTS(46) is visc in component normalised_active_and_passive_force (unit_normalised_force_millisecond_per_micrometre).
C CONSTS(47) is KSE in component normalised_active_and_passive_force (unit_normalised_force_per_micrometre).
C CONSTS(48) is SEon in component normalised_active_and_passive_force (dimensionless).
C ALGBRC(33) is FrSBXB in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (dimensionless).
C ALGBRC(34) is dFrSBXB in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (first_order_rate_constant).
C ALGBRC(36) is dsovr_ze in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micrometre_per_millisecond).
C ALGBRC(37) is dsovr_cle in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micrometre_per_millisecond).
C ALGBRC(38) is dlen_sovr in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micrometre_per_millisecond).
C ALGBRC(40) is dSOVFThick in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (first_order_rate_constant).
C ALGBRC(39) is dSOVFThin in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (first_order_rate_constant).
C CONSTS(49) is kxb in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (millinewton_per_millimetre2).
C ALGBRC(42) is dforce in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (millinewton_micrometre_per_millimetre2_per_millisecond).
C CONSTS(50) is Trop_conc in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micromolar).
C ALGBRC(35) is TropTot in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micromolar).
C ALGBRC(45) is dTropTot in component calculation_of_micromolar_per_millisecondes_of_Ca_for_apparent_Ca_binding (micromolar_per_millisecond).
C CONSTS(61) is beta in component equation_for_simulated_calcium_transient (dimensionless).
C CONSTS(51) is tau1 in component equation_for_simulated_calcium_transient (millisecond).
C CONSTS(52) is tau2 in component equation_for_simulated_calcium_transient (millisecond).
C CONSTS(53) is start_time in component equation_for_simulated_calcium_transient (millisecond).
C CONSTS(54) is Ca_amplitude in component equation_for_simulated_calcium_transient (micromolar).
C CONSTS(55) is Ca_diastolic in component equation_for_simulated_calcium_transient (micromolar).
C RATES(2) is d/dt TRPNCaL in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C RATES(3) is d/dt TRPNCaH in component Ca_binding_to_troponin_to_thin_filament_regulation (dimensionless).
C RATES(8) is d/dt N_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C RATES(9) is d/dt P_NoXB in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C RATES(10) is d/dt N in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C RATES(7) is d/dt XBprer in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C RATES(6) is d/dt XBpostr in component regulation_and_crossbridge_cycling_state_equations (dimensionless).
C RATES(5) is d/dt xXBprer in component mean_strain_of_strongly_bound_states (micrometre).
C RATES(4) is d/dt xXBpostr in component mean_strain_of_strongly_bound_states (micrometre).
C RATES(1) is d/dt SL in component normalised_active_and_passive_force (micrometre).
C RATES(11) is d/dt intf in component normalised_active_and_passive_force (unit_normalised_force_millisecond).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 2.4
      CONSTS(2) = 1.4
      CONSTS(3) = 1.2
      CONSTS(4) = 1.65
      CONSTS(5) = 0.1
      STATES(1) = 1.89999811516093
      STATES(2) = 0.0147730085063734
      STATES(3) = 0.13066096561522
      CONSTS(6) = 1.5
      CONSTS(7) = 1.3
      CONSTS(8) = 1.6
      CONSTS(9) = 1.6
      CONSTS(10) = 0.05
      CONSTS(11) = 0.25
      CONSTS(12) = 0.025
      CONSTS(13) = 0.5
      CONSTS(14) = 15
      CONSTS(15) = 0.5
      CONSTS(16) = 0.05
      CONSTS(17) = 1
      CONSTS(18) = 24
      CONSTS(19) = 0.5
      CONSTS(20) = 0.07
      CONSTS(21) = 2
      CONSTS(22) = 0.4
      CONSTS(23) = 0.07
      CONSTS(24) = 6
      CONSTS(25) = 5
      CONSTS(26) = 0
      CONSTS(27) = 8
      CONSTS(28) = 1
      CONSTS(29) = 1
      CONSTS(30) = 6.25
      CONSTS(31) = 2.5
      CONSTS(32) = 6.25
      CONSTS(33) = 6.25
      CONSTS(34) = 6.25
      CONSTS(35) = 0.007
      STATES(4) = 0.00700005394873882
      STATES(5) = 3.41212828972468e-8
      STATES(6) = 1.81017564383744e-6
      STATES(7) = 3.0494964880038e-7
      STATES(8) = 0.999999959256274
      STATES(9) = 4.07437173988636e-8
      STATES(10) = 0.999997834540066
      CONSTS(36) = 2
      CONSTS(37) = 120
      STATES(11) = -4.5113452510363e-6
      CONSTS(38) = 2.25
      CONSTS(39) = 1.85
      CONSTS(40) = 1.9
      CONSTS(41) = 0.002
      CONSTS(42) = 10
      CONSTS(43) = 0.02
      CONSTS(44) = 70
      CONSTS(45) = 50
      CONSTS(46) = 3
      CONSTS(47) = 1
      CONSTS(48) = 1
      CONSTS(49) = 120
      CONSTS(50) = 70
      CONSTS(51) = 20
      CONSTS(52) = 110
      CONSTS(53) = 5
      CONSTS(54) = 1.45
      CONSTS(55) = 0.09
      CONSTS(56) =  CONSTS(10)*CONSTS(6) ** (CONSTS(18) - 37.0000)/10.0000
      CONSTS(57) =  CONSTS(11)*CONSTS(17)*CONSTS(7) ** (CONSTS(18) - 37.0000)/10.0000
      CONSTS(58) =  CONSTS(12)*CONSTS(17)*CONSTS(7) ** (CONSTS(18) - 37.0000)/10.0000
      CONSTS(59) =  CONSTS(19)*CONSTS(29)*CONSTS(30) ** (CONSTS(18) - 37.0000)/10.0000
      CONSTS(60) = ( CONSTS(22)*CONSTS(19)+ CONSTS(23)*CONSTS(19))/( CONSTS(19)*CONSTS(21)+ CONSTS(23)*CONSTS(21)+ CONSTS(23)*CONSTS(20)+ CONSTS(22)*CONSTS(19)+ CONSTS(22)*CONSTS(20)+ CONSTS(23)*CONSTS(19))
      CONSTS(61) = CONSTS(51)/CONSTS(52) ** - 1.00000/(CONSTS(51)/CONSTS(52) - 1.00000) - CONSTS(51)/CONSTS(52) ** - 1.00000/(1.00000 - CONSTS(52)/CONSTS(51))
      CONSTS(62) = ( CONSTS(19)*CONSTS(21))/( CONSTS(19)*CONSTS(21)+ CONSTS(23)*CONSTS(21)+ CONSTS(23)*CONSTS(20)+ CONSTS(22)*CONSTS(19)+ CONSTS(22)*CONSTS(20)+ CONSTS(23)*CONSTS(19))
      CONSTS(63) =  CONSTS(37)*CONSTS(35)*CONSTS(62)
      CONSTS(64) =  (ABS(CONSTS(40) - CONSTS(39))/(CONSTS(40) - CONSTS(39)))*CONSTS(41)*(EXP( CONSTS(42)*ABS(CONSTS(40) - CONSTS(39))) - 1.00000)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(CONSTS(4)/2.00000.LT.STATES(1)/2.00000, CONSTS(4)/2.00000, STATES(1)/2.00000)
      ALGBRC(2) = TERNRY(STATES(1)/2.00000 - (STATES(1) - CONSTS(3)).GT.CONSTS(5)/2.00000, STATES(1)/2.00000 - (STATES(1) - CONSTS(3)), CONSTS(5)/2.00000)
      ALGBRC(3) = ALGBRC(1) - ALGBRC(2)
      ALGBRC(5) = ALGBRC(3)/CONSTS(3)
      ALGBRC(7) =  (1.00000 - ALGBRC(5))*STATES(2)+ ALGBRC(5)*STATES(3)
      ALGBRC(9) =  ABS(1.00000/(1.00000+CONSTS(13)/ALGBRC(7) ** CONSTS(14))) ** (1.0 / 2)
      ALGBRC(11) =  CONSTS(15)*ALGBRC(9)*CONSTS(8) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(13) = TERNRY(1.00000/ALGBRC(9).LT.100.000, 1.00000/ALGBRC(9), 100.000)
      ALGBRC(15) =  CONSTS(16)*ALGBRC(13)*CONSTS(9) ** (CONSTS(18) - 37.0000)/10.0000
      RATES(8) =  ALGBRC(15)*STATES(9) -  ALGBRC(11)*STATES(8)
      RATES(9) =  ALGBRC(11)*STATES(8) -  ALGBRC(15)*STATES(9)
      ALGBRC(4) = ( ALGBRC(3)*2.00000)/(CONSTS(4) - CONSTS(5))
      ALGBRC(6) =  CONSTS(37)*ALGBRC(4)*( STATES(4)*STATES(6)+ STATES(5)*STATES(7))
      ALGBRC(8) = ( 1.00000*ALGBRC(6))/CONSTS(63)
      ALGBRC(10) =  ((STATES(1) - CONSTS(39))/ABS(STATES(1) - CONSTS(39)))*CONSTS(41)*(EXP( CONSTS(42)*ABS(STATES(1) - CONSTS(39))) - 1.00000)
      ALGBRC(12) = TERNRY(STATES(1).GT.CONSTS(38),  CONSTS(43)*(EXP( CONSTS(44)*ABS(STATES(1) - CONSTS(38))) - 1.00000), 0.00000)
      ALGBRC(14) = ALGBRC(10)+ALGBRC(12)
      ALGBRC(16) = TERNRY(CONSTS(48).EQ.1.00000,  CONSTS(47)*(CONSTS(40) - STATES(1)), 0.00000)
      RATES(11) = (CONSTS(64)+ALGBRC(16)) - (ALGBRC(14)+ALGBRC(8))
      ALGBRC(19) = EXP( (- STATES(5)/ABS(STATES(5)))*CONSTS(25)*STATES(5)/CONSTS(35) ** 2.00000)
      ALGBRC(21) =  CONSTS(21)*ALGBRC(19)*CONSTS(29)*CONSTS(32) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(20) = EXP( ((STATES(4) - CONSTS(35))/ABS(STATES(4) - CONSTS(35)))*CONSTS(26)*(STATES(4) - CONSTS(35))/CONSTS(35) ** 2.00000)
      ALGBRC(22) =  CONSTS(22)*ALGBRC(20)*CONSTS(29)*CONSTS(33) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(23) = TERNRY(STATES(4).LT.CONSTS(35), EXP( CONSTS(27)*(CONSTS(35) - STATES(4))/CONSTS(35) ** 2.00000), EXP( CONSTS(28)*(STATES(4) - CONSTS(35))/CONSTS(35) ** 2.00000))
      ALGBRC(24) =  CONSTS(23)*ALGBRC(23)*CONSTS(29)*CONSTS(34) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(25) =  ALGBRC(21)*STATES(7) - ( ALGBRC(22)*STATES(6)+ ALGBRC(24)*STATES(6))
      RATES(6) = ALGBRC(25)
      ALGBRC(26) = ((1.00000 - STATES(10)) - STATES(7)) - STATES(6)
      RATES(10) =  ALGBRC(15)*ALGBRC(26) -  ALGBRC(11)*STATES(10)
      ALGBRC(17) = 1.00000+ (1.00000 - ALGBRC(4))*CONSTS(24)
      ALGBRC(18) =  CONSTS(20)*ALGBRC(17)*CONSTS(29)*CONSTS(31) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(27) = ( CONSTS(59)*ALGBRC(26)+ ALGBRC(22)*STATES(6)) - ( ALGBRC(18)*STATES(7)+ ALGBRC(21)*STATES(7))
      RATES(7) = ALGBRC(27)
      ALGBRC(30) = TERNRY(STATES(1).LE.CONSTS(1).AND.STATES(1).GT.CONSTS(2), (STATES(11)+ (CONSTS(40) - STATES(1))*CONSTS(46))/CONSTS(45), 0.00000)
      RATES(1) = ALGBRC(30)
      ALGBRC(28) = ( ALGBRC(22)*CONSTS(59)+ ALGBRC(24)*CONSTS(59))/( CONSTS(59)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(18)+ ALGBRC(22)*CONSTS(59)+ ALGBRC(22)*ALGBRC(18)+ ALGBRC(24)*CONSTS(59))
      ALGBRC(31) = ALGBRC(30)/2.00000+ (CONSTS(36)/ALGBRC(28))*( CONSTS(59)*- STATES(5)+ ALGBRC(22)*(STATES(4) - (CONSTS(35)+STATES(5))))
      RATES(5) = ALGBRC(31)
      ALGBRC(29) = ( CONSTS(59)*ALGBRC(21))/( CONSTS(59)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(18)+ ALGBRC(22)*CONSTS(59)+ ALGBRC(22)*ALGBRC(18)+ ALGBRC(24)*CONSTS(59))
      ALGBRC(32) = ALGBRC(30)/2.00000+ (CONSTS(36)/ALGBRC(29))*ALGBRC(21)*((STATES(5)+CONSTS(35)) - STATES(4))
      RATES(4) = ALGBRC(32)
      ALGBRC(41) = TERNRY(VOI.GT.CONSTS(53),  ((CONSTS(54) - CONSTS(55))/CONSTS(61))*(EXP(- (VOI - CONSTS(53))/CONSTS(51)) - EXP(- (VOI - CONSTS(53))/CONSTS(52)))+CONSTS(55), CONSTS(55))
      ALGBRC(43) =  CONSTS(56)*ALGBRC(41)*(1.00000 - STATES(2)) -  CONSTS(57)*STATES(2)
      RATES(2) = ALGBRC(43)
      ALGBRC(44) =  CONSTS(56)*ALGBRC(41)*(1.00000 - STATES(3)) -  CONSTS(58)*STATES(3)
      RATES(3) = ALGBRC(44)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(CONSTS(4)/2.00000.LT.STATES(1)/2.00000, CONSTS(4)/2.00000, STATES(1)/2.00000)
      ALGBRC(2) = TERNRY(STATES(1)/2.00000 - (STATES(1) - CONSTS(3)).GT.CONSTS(5)/2.00000, STATES(1)/2.00000 - (STATES(1) - CONSTS(3)), CONSTS(5)/2.00000)
      ALGBRC(3) = ALGBRC(1) - ALGBRC(2)
      ALGBRC(5) = ALGBRC(3)/CONSTS(3)
      ALGBRC(7) =  (1.00000 - ALGBRC(5))*STATES(2)+ ALGBRC(5)*STATES(3)
      ALGBRC(9) =  ABS(1.00000/(1.00000+CONSTS(13)/ALGBRC(7) ** CONSTS(14))) ** (1.0 / 2)
      ALGBRC(11) =  CONSTS(15)*ALGBRC(9)*CONSTS(8) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(13) = TERNRY(1.00000/ALGBRC(9).LT.100.000, 1.00000/ALGBRC(9), 100.000)
      ALGBRC(15) =  CONSTS(16)*ALGBRC(13)*CONSTS(9) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(4) = ( ALGBRC(3)*2.00000)/(CONSTS(4) - CONSTS(5))
      ALGBRC(6) =  CONSTS(37)*ALGBRC(4)*( STATES(4)*STATES(6)+ STATES(5)*STATES(7))
      ALGBRC(8) = ( 1.00000*ALGBRC(6))/CONSTS(63)
      ALGBRC(10) =  ((STATES(1) - CONSTS(39))/ABS(STATES(1) - CONSTS(39)))*CONSTS(41)*(EXP( CONSTS(42)*ABS(STATES(1) - CONSTS(39))) - 1.00000)
      ALGBRC(12) = TERNRY(STATES(1).GT.CONSTS(38),  CONSTS(43)*(EXP( CONSTS(44)*ABS(STATES(1) - CONSTS(38))) - 1.00000), 0.00000)
      ALGBRC(14) = ALGBRC(10)+ALGBRC(12)
      ALGBRC(16) = TERNRY(CONSTS(48).EQ.1.00000,  CONSTS(47)*(CONSTS(40) - STATES(1)), 0.00000)
      ALGBRC(19) = EXP( (- STATES(5)/ABS(STATES(5)))*CONSTS(25)*STATES(5)/CONSTS(35) ** 2.00000)
      ALGBRC(21) =  CONSTS(21)*ALGBRC(19)*CONSTS(29)*CONSTS(32) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(20) = EXP( ((STATES(4) - CONSTS(35))/ABS(STATES(4) - CONSTS(35)))*CONSTS(26)*(STATES(4) - CONSTS(35))/CONSTS(35) ** 2.00000)
      ALGBRC(22) =  CONSTS(22)*ALGBRC(20)*CONSTS(29)*CONSTS(33) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(23) = TERNRY(STATES(4).LT.CONSTS(35), EXP( CONSTS(27)*(CONSTS(35) - STATES(4))/CONSTS(35) ** 2.00000), EXP( CONSTS(28)*(STATES(4) - CONSTS(35))/CONSTS(35) ** 2.00000))
      ALGBRC(24) =  CONSTS(23)*ALGBRC(23)*CONSTS(29)*CONSTS(34) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(25) =  ALGBRC(21)*STATES(7) - ( ALGBRC(22)*STATES(6)+ ALGBRC(24)*STATES(6))
      ALGBRC(26) = ((1.00000 - STATES(10)) - STATES(7)) - STATES(6)
      ALGBRC(17) = 1.00000+ (1.00000 - ALGBRC(4))*CONSTS(24)
      ALGBRC(18) =  CONSTS(20)*ALGBRC(17)*CONSTS(29)*CONSTS(31) ** (CONSTS(18) - 37.0000)/10.0000
      ALGBRC(27) = ( CONSTS(59)*ALGBRC(26)+ ALGBRC(22)*STATES(6)) - ( ALGBRC(18)*STATES(7)+ ALGBRC(21)*STATES(7))
      ALGBRC(30) = TERNRY(STATES(1).LE.CONSTS(1).AND.STATES(1).GT.CONSTS(2), (STATES(11)+ (CONSTS(40) - STATES(1))*CONSTS(46))/CONSTS(45), 0.00000)
      ALGBRC(28) = ( ALGBRC(22)*CONSTS(59)+ ALGBRC(24)*CONSTS(59))/( CONSTS(59)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(18)+ ALGBRC(22)*CONSTS(59)+ ALGBRC(22)*ALGBRC(18)+ ALGBRC(24)*CONSTS(59))
      ALGBRC(31) = ALGBRC(30)/2.00000+ (CONSTS(36)/ALGBRC(28))*( CONSTS(59)*- STATES(5)+ ALGBRC(22)*(STATES(4) - (CONSTS(35)+STATES(5))))
      ALGBRC(29) = ( CONSTS(59)*ALGBRC(21))/( CONSTS(59)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(21)+ ALGBRC(24)*ALGBRC(18)+ ALGBRC(22)*CONSTS(59)+ ALGBRC(22)*ALGBRC(18)+ ALGBRC(24)*CONSTS(59))
      ALGBRC(32) = ALGBRC(30)/2.00000+ (CONSTS(36)/ALGBRC(29))*ALGBRC(21)*((STATES(5)+CONSTS(35)) - STATES(4))
      ALGBRC(41) = TERNRY(VOI.GT.CONSTS(53),  ((CONSTS(54) - CONSTS(55))/CONSTS(61))*(EXP(- (VOI - CONSTS(53))/CONSTS(51)) - EXP(- (VOI - CONSTS(53))/CONSTS(52)))+CONSTS(55), CONSTS(55))
      ALGBRC(43) =  CONSTS(56)*ALGBRC(41)*(1.00000 - STATES(2)) -  CONSTS(57)*STATES(2)
      ALGBRC(44) =  CONSTS(56)*ALGBRC(41)*(1.00000 - STATES(3)) -  CONSTS(58)*STATES(3)
      ALGBRC(33) = (STATES(6)+STATES(7))/(CONSTS(62)+CONSTS(60))
      ALGBRC(34) = (ALGBRC(25)+ALGBRC(27))/(CONSTS(62)+CONSTS(60))
      ALGBRC(35) =  CONSTS(50)*( (1.00000 - ALGBRC(5))*STATES(2)+ ALGBRC(5)*( ALGBRC(33)*STATES(3)+ (1.00000 - ALGBRC(33))*STATES(2)))
      ALGBRC(36) = TERNRY(STATES(1).LT.CONSTS(4),  - 0.500000*ALGBRC(30), 0.00000)
      ALGBRC(37) = TERNRY( 2.00000*CONSTS(3) - STATES(1).GT.CONSTS(5),  - 0.500000*ALGBRC(30), 0.00000)
      ALGBRC(38) = ALGBRC(36) - ALGBRC(37)
      ALGBRC(39) = ALGBRC(38)/CONSTS(3)
      ALGBRC(40) = ( 2.00000*ALGBRC(38))/(CONSTS(4) - CONSTS(5))
      ALGBRC(42) =  CONSTS(49)*ALGBRC(40)*( STATES(4)*STATES(6)+ STATES(5)*STATES(7))+ CONSTS(49)*ALGBRC(4)*( ALGBRC(32)*STATES(6)+ STATES(4)*ALGBRC(25)+ ALGBRC(31)*STATES(7)+ STATES(5)*ALGBRC(27))
      ALGBRC(45) =  CONSTS(50)*( - ALGBRC(39)*STATES(2)+ (1.00000 - ALGBRC(5))*ALGBRC(43)+ ALGBRC(39)*( ALGBRC(33)*STATES(3)+ (1.00000 - ALGBRC(33))*STATES(2))+ ALGBRC(5)*(( ALGBRC(34)*STATES(3)+ ALGBRC(33)*ALGBRC(44)+ (1.00000 - ALGBRC(33))*ALGBRC(43)) -  ALGBRC(34)*STATES(2)))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END