C
C There are a total of 233 entries in the algebraic variable array.
C There are a total of 44 entries in each of the rate and state variable arrays.
C There are a total of 186 entries in the constant variable array.
C
C
C CONSTS(1) is kappa_R_TRPNCa in component BG_parameters (fmol_per_sec).
C CONSTS(2) is kappa_R_BU in component BG_parameters (fmol_per_sec).
C CONSTS(3) is kappa_R_UW in component BG_parameters (fmol_per_sec).
C CONSTS(4) is kappa_R_WS in component BG_parameters (fmol_per_sec).
C CONSTS(5) is kappa_R_SU in component BG_parameters (fmol_per_sec).
C CONSTS(6) is kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec).
C CONSTS(7) is kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec).
C CONSTS(8) is kappa_LCC_K1 in component BG_parameters (fmol_per_sec).
C CONSTS(9) is kappa_LCC_K2 in component BG_parameters (fmol_per_sec).
C CONSTS(10) is kappa_d000 in component BG_parameters (fmol_per_sec).
C CONSTS(11) is kappa_d010 in component BG_parameters (fmol_per_sec).
C CONSTS(12) is kappa_d020 in component BG_parameters (fmol_per_sec).
C CONSTS(13) is kappa_d001 in component BG_parameters (fmol_per_sec).
C CONSTS(14) is kappa_d011 in component BG_parameters (fmol_per_sec).
C CONSTS(15) is kappa_d021 in component BG_parameters (fmol_per_sec).
C CONSTS(16) is kappa_f1_000 in component BG_parameters (fmol_per_sec).
C CONSTS(17) is kappa_f1_100 in component BG_parameters (fmol_per_sec).
C CONSTS(18) is kappa_f1_001 in component BG_parameters (fmol_per_sec).
C CONSTS(19) is kappa_f1_101 in component BG_parameters (fmol_per_sec).
C CONSTS(20) is kappa_f2_000 in component BG_parameters (fmol_per_sec).
C CONSTS(21) is kappa_f2_100 in component BG_parameters (fmol_per_sec).
C CONSTS(22) is kappa_f2_001 in component BG_parameters (fmol_per_sec).
C CONSTS(23) is kappa_f2_101 in component BG_parameters (fmol_per_sec).
C CONSTS(24) is kappa_f3_010 in component BG_parameters (fmol_per_sec).
C CONSTS(25) is kappa_f3_110 in component BG_parameters (fmol_per_sec).
C CONSTS(26) is kappa_f3_011 in component BG_parameters (fmol_per_sec).
C CONSTS(27) is kappa_f3_111 in component BG_parameters (fmol_per_sec).
C CONSTS(28) is kappa_fCa000 in component BG_parameters (fmol_per_sec).
C CONSTS(29) is kappa_fCa100 in component BG_parameters (fmol_per_sec).
C CONSTS(30) is kappa_fCa010 in component BG_parameters (fmol_per_sec).
C CONSTS(31) is kappa_fCa110 in component BG_parameters (fmol_per_sec).
C CONSTS(32) is kappa_fCa020 in component BG_parameters (fmol_per_sec).
C CONSTS(33) is kappa_fCa120 in component BG_parameters (fmol_per_sec).
C CONSTS(34) is kappa_PLBph1 in component BG_parameters (fmol_per_sec).
C CONSTS(35) is kappa_PLBph2 in component BG_parameters (fmol_per_sec).
C CONSTS(36) is kappa_PLBd1 in component BG_parameters (fmol_per_sec).
C CONSTS(37) is kappa_PLBd2 in component BG_parameters (fmol_per_sec).
C CONSTS(38) is kappa_Inh in component BG_parameters (fmol_per_sec).
C CONSTS(39) is kappa_RyR in component BG_parameters (fmol_per_sec).
C CONSTS(40) is kappa_OC in component BG_parameters (fmol_per_sec).
C CONSTS(41) is kappa_CCI in component BG_parameters (fmol_per_sec).
C CONSTS(42) is kappa_CII in component BG_parameters (fmol_per_sec).
C CONSTS(43) is kappa_IO in component BG_parameters (fmol_per_sec).
C CONSTS(44) is kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec).
C CONSTS(45) is kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec).
C CONSTS(46) is kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec).
C CONSTS(47) is kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec).
C CONSTS(48) is kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec).
C CONSTS(49) is kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec).
C CONSTS(50) is kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec).
C CONSTS(51) is kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec).
C CONSTS(52) is kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec).
C CONSTS(53) is K_TRPN in component BG_parameters (per_fmol).
C CONSTS(54) is K_Cai in component BG_parameters (per_fmol).
C CONSTS(55) is K_Ca_TRPN in component BG_parameters (per_fmol).
C CONSTS(56) is K_B in component BG_parameters (per_fmol).
C CONSTS(57) is K_U in component BG_parameters (per_fmol).
C CONSTS(58) is K_W in component BG_parameters (per_fmol).
C CONSTS(59) is K_S in component BG_parameters (per_fmol).
C CONSTS(60) is K_Ca_i in component BG_parameters (per_fmol).
C CONSTS(61) is K_Ca_o in component BG_parameters (per_fmol).
C CONSTS(62) is K_K_i in component BG_parameters (per_fmol).
C CONSTS(63) is K_K_o in component BG_parameters (per_fmol).
C CONSTS(64) is K_000_LCC in component BG_parameters (per_fmol).
C CONSTS(65) is K_010_LCC in component BG_parameters (per_fmol).
C CONSTS(66) is K_020_LCC in component BG_parameters (per_fmol).
C CONSTS(67) is K_100_LCC in component BG_parameters (per_fmol).
C CONSTS(68) is K_110_LCC in component BG_parameters (per_fmol).
C CONSTS(69) is K_120_LCC in component BG_parameters (per_fmol).
C CONSTS(70) is K_001_LCC in component BG_parameters (per_fmol).
C CONSTS(71) is K_011_LCC in component BG_parameters (per_fmol).
C CONSTS(72) is K_021_LCC in component BG_parameters (per_fmol).
C CONSTS(73) is K_101_LCC in component BG_parameters (per_fmol).
C CONSTS(74) is K_111_LCC in component BG_parameters (per_fmol).
C CONSTS(75) is K_121_LCC in component BG_parameters (per_fmol).
C CONSTS(76) is K_PLB in component BG_parameters (per_fmol).
C CONSTS(77) is K_PKACI in component BG_parameters (per_fmol).
C CONSTS(78) is K_PLB_PKACI in component BG_parameters (per_fmol).
C CONSTS(79) is K_PP1 in component BG_parameters (per_fmol).
C CONSTS(80) is K_PLBp_PP1 in component BG_parameters (per_fmol).
C CONSTS(81) is K_PLBp in component BG_parameters (per_fmol).
C CONSTS(82) is K_Ip in component BG_parameters (per_fmol).
C CONSTS(83) is K_Ip_PP1 in component BG_parameters (per_fmol).
C CONSTS(84) is K_Ca_SR in component BG_parameters (per_fmol).
C CONSTS(85) is K_C_RyR in component BG_parameters (per_fmol).
C CONSTS(86) is K_CI_RyR in component BG_parameters (per_fmol).
C CONSTS(87) is K_I_RyR in component BG_parameters (per_fmol).
C CONSTS(88) is K_O_RyR in component BG_parameters (per_fmol).
C CONSTS(89) is K_P1_SERCA in component BG_parameters (per_fmol).
C CONSTS(90) is K_P2_SERCA in component BG_parameters (per_fmol).
C CONSTS(91) is K_P2a_SERCA in component BG_parameters (per_fmol).
C CONSTS(92) is K_P4_SERCA in component BG_parameters (per_fmol).
C CONSTS(93) is K_P5_SERCA in component BG_parameters (per_fmol).
C CONSTS(94) is K_P6_SERCA in component BG_parameters (per_fmol).
C CONSTS(95) is K_P8_SERCA in component BG_parameters (per_fmol).
C CONSTS(96) is K_P9_SERCA in component BG_parameters (per_fmol).
C CONSTS(97) is K_P10_SERCA in component BG_parameters (per_fmol).
C CONSTS(98) is K_H in component BG_parameters (per_fmol).
C CONSTS(99) is K_MgATP in component BG_parameters (per_fmol).
C CONSTS(100) is K_MgADP in component BG_parameters (per_fmol).
C CONSTS(101) is K_Pi in component BG_parameters (per_fmol).
C CONSTS(102) is z_rCa in component BG_parameters (dimensionless).
C CONSTS(103) is z_fd in component BG_parameters (dimensionless).
C CONSTS(104) is z_ff1 in component BG_parameters (dimensionless).
C CONSTS(105) is z_ff2 in component BG_parameters (dimensionless).
C CONSTS(106) is z_ff3 in component BG_parameters (dimensionless).
C CONSTS(107) is z_rd in component BG_parameters (dimensionless).
C CONSTS(108) is z_rf1 in component BG_parameters (dimensionless).
C CONSTS(109) is z_rf2 in component BG_parameters (dimensionless).
C CONSTS(110) is z_rf3 in component BG_parameters (dimensionless).
C VOI is time in component environment (second).
C CONSTS(111) is vol_myo in component environment (pL).
C CONSTS(112) is freq in component environment (dimensionless).
C CONSTS(113) is C_m in component environment (fF).
C CONSTS(114) is F in component constants (C_per_mol).
C STATES(1) is q_membrane in component environment (fC).
C ALGBRC(1) is V_m in component environment (volt).
C STATES(2) is q_TRPN in component environment (fmol).
C STATES(3) is q_Ca_TRPN in component environment (fmol).
C STATES(4) is q_Cai in component environment (fmol).
C STATES(5) is q_B in component environment (fmol).
C STATES(6) is q_U in component environment (fmol).
C STATES(7) is q_W in component environment (fmol).
C STATES(8) is q_S in component environment (fmol).
C STATES(9) is q_Cao in component environment (fmol).
C STATES(10) is q_Ki in component environment (fmol).
C STATES(11) is q_Ko in component environment (fmol).
C STATES(12) is q_PLB in component environment (fmol).
C STATES(13) is q_PKACI in component environment (fmol).
C STATES(14) is q_PLB_PKACI in component environment (fmol).
C STATES(15) is q_PP1 in component environment (fmol).
C STATES(16) is q_PLBp_PP1 in component environment (fmol).
C STATES(17) is q_PLBp in component environment (fmol).
C STATES(18) is q_Ip in component environment (fmol).
C STATES(19) is q_Ip_PP1 in component environment (fmol).
C STATES(20) is q_Ca_SR in component environment (fmol).
C STATES(21) is q_H in component environment (fmol).
C STATES(22) is q_P in component environment (fmol).
C STATES(23) is q_MgADP in component environment (fmol).
C STATES(24) is q_MgATP in component environment (fmol).
C STATES(25) is q_000_LCC in component environment (fmol).
C STATES(26) is q_010_LCC in component environment (fmol).
C STATES(27) is q_020_LCC in component environment (fmol).
C STATES(28) is q_100_LCC in component environment (fmol).
C STATES(29) is q_110_LCC in component environment (fmol).
C STATES(30) is q_120_LCC in component environment (fmol).
C STATES(31) is q_001_LCC in component environment (fmol).
C STATES(32) is q_011_LCC in component environment (fmol).
C STATES(33) is q_021_LCC in component environment (fmol).
C STATES(34) is q_101_LCC in component environment (fmol).
C STATES(35) is q_111_LCC in component environment (fmol).
C STATES(36) is q_121_LCC in component environment (fmol).
C STATES(37) is q_C_RyR in component environment (fmol).
C STATES(38) is q_CI_RyR in component environment (fmol).
C STATES(39) is q_I_RyR in component environment (fmol).
C STATES(40) is q_O_RyR in component environment (fmol).
C CONSTS(115) is q_P1_SERCA in component environment (fmol).
C CONSTS(116) is q_P2_SERCA in component environment (fmol).
C CONSTS(117) is q_P2a_SERCA in component environment (fmol).
C CONSTS(118) is q_P4_SERCA in component environment (fmol).
C CONSTS(119) is q_P5_SERCA in component environment (fmol).
C CONSTS(120) is q_P6_SERCA in component environment (fmol).
C CONSTS(121) is q_P8_SERCA in component environment (fmol).
C CONSTS(122) is q_P9_SERCA in component environment (fmol).
C CONSTS(123) is q_P10_SERCA in component environment (fmol).
C ALGBRC(233) is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(44) is v_R_BU in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(47) is v_R_UW in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(50) is v_R_WS in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(53) is v_R_SU in component crossbridge_TRPN (fmol_per_sec).
C ALGBRC(230) is T_total in component crossbridge_TRPN (kPa).
C STATES(41) is SL in component crossbridge_TRPN (metre).
C ALGBRC(164) is v_Ca_i_LCC in component LCC (fmol_per_sec).
C ALGBRC(162) is v_Ca_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(168) is v_K_i_LCC in component LCC (fmol_per_sec).
C ALGBRC(169) is v_K_o_LCC in component LCC (fmol_per_sec).
C ALGBRC(170) is I_mem_LCC in component LCC (fA).
C ALGBRC(45) is v_PLBph1 in component PLB (fmol_per_sec).
C ALGBRC(48) is v_PLBph2 in component PLB (fmol_per_sec).
C ALGBRC(51) is v_PLBd1 in component PLB (fmol_per_sec).
C ALGBRC(54) is v_PLBd2 in component PLB (fmol_per_sec).
C ALGBRC(56) is v_Inh in component PLB (fmol_per_sec).
C ALGBRC(178) is v_RyR in component RyR (fmol_per_sec).
C ALGBRC(183) is v_RyRgate_Ca_i in component RyR (fmol_per_sec).
C ALGBRC(210) is v_Ca_i_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(220) is v_Ca_SR_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(227) is v_H_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(205) is v_MgADP_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(202) is v_MgATP_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(200) is v_Pi_SERCA in component SERCA (fmol_per_sec).
C ALGBRC(10) is I_pulse in component environment (fA).
C CONSTS(124) is pulse_start in component environment (second).
C CONSTS(125) is pulse_end in component environment (second).
C CONSTS(126) is pulseMag in component environment (fA).
C CONSTS(127) is pulseHolding in component environment (fA).
C ALGBRC(173) is sum_I in component environment (fA).
C CONSTS(128) is zCa in component ion_valences (dimensionless).
C CONSTS(129) is zK in component ion_valences (dimensionless).
C ALGBRC(3) is PLB_T in component environment (fmol).
C ALGBRC(4) is PKACI_T in component environment (fmol).
C ALGBRC(5) is Ip_T in component environment (fmol).
C ALGBRC(6) is K_i_T in component environment (fmol).
C ALGBRC(7) is K_o_T in component environment (fmol).
C CONSTS(130) is R in component constants (J_per_K_per_mol).
C CONSTS(131) is T in component constants (kelvin).
C CONSTS(132) is zNa in component ion_valences (dimensionless).
C CONSTS(133) is zCl in component ion_valences (dimensionless).
C CONSTS(134) is n_Tm in component crossbridge_TRPN (dimensionless).
C ALGBRC(25) is mu_TRPN in component crossbridge_TRPN (J_per_mol).
C ALGBRC(29) is mu_Cai in component crossbridge_TRPN (J_per_mol).
C ALGBRC(31) is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol).
C ALGBRC(32) is mu_B in component crossbridge_TRPN (J_per_mol).
C ALGBRC(35) is mu_U in component crossbridge_TRPN (J_per_mol).
C ALGBRC(38) is mu_W in component crossbridge_TRPN (J_per_mol).
C ALGBRC(41) is mu_S in component crossbridge_TRPN (J_per_mol).
C ALGBRC(231) is tension in component crossbridge_TRPN (N_per_mm2).
C ALGBRC(232) is mu_tension in component crossbridge_TRPN (J_per_mol).
C CONSTS(135) is kf_coeff in component crossbridge_TRPN (dimensionless).
C CONSTS(161) is n in component crossbridge_TRPN (J_per_mol).
C CONSTS(136) is hh in component crossbridge_TRPN (mm2_per_N).
C CONSTS(137) is SL_0 in component crossbridge_TRPN (metre).
C CONSTS(138) is q_MS in component crossbridge_TRPN (fmol).
C CONSTS(139) is r_s in component crossbridge_TRPN (dimensionless).
C CONSTS(140) is r_w in component crossbridge_TRPN (dimensionless).
C CONSTS(141) is A_eff in component crossbridge_TRPN (dimensionless).
C CONSTS(142) is phi in component crossbridge_TRPN (dimensionless).
C CONSTS(143) is k_uw in component crossbridge_TRPN (per_sec).
C CONSTS(144) is k_ws in component crossbridge_TRPN (per_sec).
C STATES(42) is G_w in component crossbridge_TRPN (metre).
C STATES(43) is G_s in component crossbridge_TRPN (metre).
C CONSTS(158) is c_w in component crossbridge_TRPN (per_sec).
C CONSTS(159) is c_s in component crossbridge_TRPN (per_sec).
C CONSTS(160) is A_w in component crossbridge_TRPN (dimensionless).
C CONSTS(164) is A_s in component crossbridge_TRPN (dimensionless).
C ALGBRC(23) is mu_1 in component crossbridge_TRPN (J_per_m).
C CONSTS(165) is mu_2 in component crossbridge_TRPN (J_per_m).
C ALGBRC(8) is mu_3 in component crossbridge_TRPN (J_per_m).
C ALGBRC(24) is mu_4 in component crossbridge_TRPN (J_per_m).
C CONSTS(168) is mu_5 in component crossbridge_TRPN (J_per_m).
C ALGBRC(9) is mu_6 in component crossbridge_TRPN (J_per_m).
C ALGBRC(27) is v_1 in component crossbridge_TRPN (m_per_s).
C CONSTS(145) is v_2 in component crossbridge_TRPN (m_per_s).
C ALGBRC(28) is v_4 in component crossbridge_TRPN (m_per_s).
C CONSTS(146) is v_to_mu in component crossbridge_TRPN (Js_per_m2).
C ALGBRC(43) is T_active in component crossbridge_TRPN (kPa).
C CONSTS(147) is T_ref in component crossbridge_TRPN (kPa).
C ALGBRC(40) is mu_T_a in component crossbridge_TRPN (kPa).
C ALGBRC(34) is mu_T_S in component crossbridge_TRPN (kPa).
C ALGBRC(37) is mu_T_W in component crossbridge_TRPN (kPa).
C STATES(44) is Cdd in component crossbridge_TRPN (metre).
C ALGBRC(213) is eta in component crossbridge_TRPN (per_sec).
C CONSTS(148) is eta_l in component crossbridge_TRPN (per_sec).
C CONSTS(149) is eta_s in component crossbridge_TRPN (per_sec).
C CONSTS(150) is k in component crossbridge_TRPN (dimensionless).
C CONSTS(151) is alpha in component crossbridge_TRPN (kPa).
C ALGBRC(219) is v_Cdd in component crossbridge_TRPN (m_per_s).
C ALGBRC(214) is mu_d in component crossbridge_TRPN (J_per_m).
C ALGBRC(215) is mu_k in component crossbridge_TRPN (J_per_m).
C ALGBRC(216) is mu_f in component crossbridge_TRPN (J_per_m).
C ALGBRC(223) is mu_T_passive in component crossbridge_TRPN (kPa).
C ALGBRC(226) is T_passive in component crossbridge_TRPN (kPa).
C ALGBRC(49) is mu_Ca_o in component LCC (J_per_mol).
C ALGBRC(52) is mu_Ca_i in component LCC (J_per_mol).
C ALGBRC(55) is mu_K_o in component LCC (J_per_mol).
C ALGBRC(57) is mu_K_i in component LCC (J_per_mol).
C ALGBRC(46) is V_mem in component LCC (volt).
C ALGBRC(66) is Af_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(98) is Ar_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(94) is Am_LCC_Ca1 in component LCC (J_per_mol).
C ALGBRC(159) is v_LCC_Ca1 in component LCC (fmol_per_sec).
C ALGBRC(67) is Af_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(99) is Ar_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(95) is Am_LCC_Ca2 in component LCC (J_per_mol).
C ALGBRC(161) is v_LCC_Ca2 in component LCC (fmol_per_sec).
C ALGBRC(68) is Af_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(100) is Ar_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(96) is Am_LCC_K1 in component LCC (J_per_mol).
C ALGBRC(163) is v_LCC_K1 in component LCC (fmol_per_sec).
C ALGBRC(69) is Af_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(101) is Ar_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(97) is Am_LCC_K2 in component LCC (J_per_mol).
C ALGBRC(165) is v_LCC_K2 in component LCC (fmol_per_sec).
C ALGBRC(88) is Af_fCa00 in component LCC (J_per_mol).
C ALGBRC(138) is Ar_fCa00 in component LCC (J_per_mol).
C ALGBRC(139) is v_fCa00 in component LCC (fmol_per_sec).
C ALGBRC(89) is Af_fCa01 in component LCC (J_per_mol).
C ALGBRC(140) is Ar_fCa01 in component LCC (J_per_mol).
C ALGBRC(143) is v_fCa01 in component LCC (fmol_per_sec).
C ALGBRC(90) is Af_fCa02 in component LCC (J_per_mol).
C ALGBRC(144) is Ar_fCa02 in component LCC (J_per_mol).
C ALGBRC(147) is v_fCa02 in component LCC (fmol_per_sec).
C ALGBRC(91) is Af_fCa10 in component LCC (J_per_mol).
C ALGBRC(148) is Ar_fCa10 in component LCC (J_per_mol).
C ALGBRC(151) is v_fCa10 in component LCC (fmol_per_sec).
C ALGBRC(92) is Af_fCa11 in component LCC (J_per_mol).
C ALGBRC(152) is Ar_fCa11 in component LCC (J_per_mol).
C ALGBRC(155) is v_fCa11 in component LCC (fmol_per_sec).
C ALGBRC(93) is Af_fCa12 in component LCC (J_per_mol).
C ALGBRC(156) is Ar_fCa12 in component LCC (J_per_mol).
C ALGBRC(158) is v_fCa12 in component LCC (fmol_per_sec).
C ALGBRC(58) is mu_S000_LCC in component LCC (J_per_mol).
C ALGBRC(141) is v_S000_LCC in component LCC (fmol_per_sec).
C ALGBRC(59) is mu_S010_LCC in component LCC (J_per_mol).
C ALGBRC(145) is v_S010_LCC in component LCC (fmol_per_sec).
C ALGBRC(60) is mu_S020_LCC in component LCC (J_per_mol).
C ALGBRC(149) is v_S020_LCC in component LCC (fmol_per_sec).
C ALGBRC(61) is mu_S100_LCC in component LCC (J_per_mol).
C ALGBRC(153) is v_S100_LCC in component LCC (fmol_per_sec).
C ALGBRC(62) is mu_S110_LCC in component LCC (J_per_mol).
C ALGBRC(157) is v_S110_LCC in component LCC (fmol_per_sec).
C ALGBRC(64) is mu_S120_LCC in component LCC (J_per_mol).
C ALGBRC(160) is v_S120_LCC in component LCC (fmol_per_sec).
C ALGBRC(11) is mu_S001_LCC in component LCC (J_per_mol).
C ALGBRC(142) is v_S001_LCC in component LCC (fmol_per_sec).
C ALGBRC(12) is mu_S011_LCC in component LCC (J_per_mol).
C ALGBRC(146) is v_S011_LCC in component LCC (fmol_per_sec).
C ALGBRC(13) is mu_S021_LCC in component LCC (J_per_mol).
C ALGBRC(150) is v_S021_LCC in component LCC (fmol_per_sec).
C ALGBRC(14) is mu_S101_LCC in component LCC (J_per_mol).
C ALGBRC(154) is v_S101_LCC in component LCC (fmol_per_sec).
C ALGBRC(63) is mu_S111_LCC in component LCC (J_per_mol).
C ALGBRC(167) is v_S111_LCC in component LCC (fmol_per_sec).
C ALGBRC(65) is mu_S121_LCC in component LCC (J_per_mol).
C ALGBRC(172) is v_S121_LCC in component LCC (fmol_per_sec).
C ALGBRC(70) is Af_d000 in component LCC (J_per_mol).
C ALGBRC(102) is Ar_d000 in component LCC (J_per_mol).
C ALGBRC(103) is v_d000 in component LCC (fmol_per_sec).
C ALGBRC(72) is Af_d010 in component LCC (J_per_mol).
C ALGBRC(106) is Ar_d010 in component LCC (J_per_mol).
C ALGBRC(107) is v_d010 in component LCC (fmol_per_sec).
C ALGBRC(74) is Af_d020 in component LCC (J_per_mol).
C ALGBRC(110) is Ar_d020 in component LCC (J_per_mol).
C ALGBRC(111) is v_d020 in component LCC (fmol_per_sec).
C ALGBRC(76) is Af_f1_000 in component LCC (J_per_mol).
C ALGBRC(114) is Ar_f1_000 in component LCC (J_per_mol).
C ALGBRC(115) is v_f1_000 in component LCC (fmol_per_sec).
C ALGBRC(80) is Af_f2_000 in component LCC (J_per_mol).
C ALGBRC(122) is Ar_f2_000 in component LCC (J_per_mol).
C ALGBRC(123) is v_f2_000 in component LCC (fmol_per_sec).
C ALGBRC(84) is Af_f3_010 in component LCC (J_per_mol).
C ALGBRC(130) is Ar_f3_010 in component LCC (J_per_mol).
C ALGBRC(131) is v_f3_010 in component LCC (fmol_per_sec).
C ALGBRC(78) is Af_f1_100 in component LCC (J_per_mol).
C ALGBRC(118) is Ar_f1_100 in component LCC (J_per_mol).
C ALGBRC(119) is v_f1_100 in component LCC (fmol_per_sec).
C ALGBRC(82) is Af_f2_100 in component LCC (J_per_mol).
C ALGBRC(126) is Ar_f2_100 in component LCC (J_per_mol).
C ALGBRC(127) is v_f2_100 in component LCC (fmol_per_sec).
C ALGBRC(86) is Af_f3_110 in component LCC (J_per_mol).
C ALGBRC(134) is Ar_f3_110 in component LCC (J_per_mol).
C ALGBRC(135) is v_f3_110 in component LCC (fmol_per_sec).
C ALGBRC(71) is Af_d001 in component LCC (J_per_mol).
C ALGBRC(104) is Ar_d001 in component LCC (J_per_mol).
C ALGBRC(105) is v_d001 in component LCC (fmol_per_sec).
C ALGBRC(73) is Af_d011 in component LCC (J_per_mol).
C ALGBRC(108) is Ar_d011 in component LCC (J_per_mol).
C ALGBRC(109) is v_d011 in component LCC (fmol_per_sec).
C ALGBRC(75) is Af_d021 in component LCC (J_per_mol).
C ALGBRC(112) is Ar_d021 in component LCC (J_per_mol).
C ALGBRC(113) is v_d021 in component LCC (fmol_per_sec).
C ALGBRC(77) is Af_f1_001 in component LCC (J_per_mol).
C ALGBRC(116) is Ar_f1_001 in component LCC (J_per_mol).
C ALGBRC(117) is v_f1_001 in component LCC (fmol_per_sec).
C ALGBRC(81) is Af_f2_001 in component LCC (J_per_mol).
C ALGBRC(124) is Ar_f2_001 in component LCC (J_per_mol).
C ALGBRC(125) is v_f2_001 in component LCC (fmol_per_sec).
C ALGBRC(85) is Af_f3_011 in component LCC (J_per_mol).
C ALGBRC(132) is Ar_f3_011 in component LCC (J_per_mol).
C ALGBRC(133) is v_f3_011 in component LCC (fmol_per_sec).
C ALGBRC(79) is Af_f1_101 in component LCC (J_per_mol).
C ALGBRC(120) is Ar_f1_101 in component LCC (J_per_mol).
C ALGBRC(121) is v_f1_101 in component LCC (fmol_per_sec).
C ALGBRC(83) is Af_f2_101 in component LCC (J_per_mol).
C ALGBRC(128) is Ar_f2_101 in component LCC (J_per_mol).
C ALGBRC(129) is v_f2_101 in component LCC (fmol_per_sec).
C ALGBRC(87) is Af_f3_111 in component LCC (J_per_mol).
C ALGBRC(136) is Ar_f3_111 in component LCC (J_per_mol).
C ALGBRC(137) is v_f3_111 in component LCC (fmol_per_sec).
C ALGBRC(2) is mu_PLB in component PLB (J_per_mol).
C ALGBRC(22) is mu_PKACI in component PLB (J_per_mol).
C ALGBRC(26) is mu_PLB_PKACI in component PLB (J_per_mol).
C ALGBRC(30) is mu_PP1 in component PLB (J_per_mol).
C ALGBRC(33) is mu_PLBp_PP1 in component PLB (J_per_mol).
C ALGBRC(36) is mu_PLBp in component PLB (J_per_mol).
C ALGBRC(39) is mu_Ip in component PLB (J_per_mol).
C ALGBRC(42) is mu_Ip_PP1 in component PLB (J_per_mol).
C CONSTS(152) is nCa_1 in component RyR (dimensionless).
C CONSTS(153) is nCa_2 in component RyR (dimensionless).
C ALGBRC(166) is mu_Ca_SR in component RyR (J_per_mol).
C ALGBRC(171) is mu_Ca_i in component RyR (J_per_mol).
C ALGBRC(15) is q_Ca_gate_complexes in component RyR (fmol).
C ALGBRC(174) is mu_C_RyR in component RyR (J_per_mol).
C ALGBRC(175) is mu_CI_RyR in component RyR (J_per_mol).
C ALGBRC(176) is mu_I_RyR in component RyR (J_per_mol).
C ALGBRC(177) is mu_O_RyR in component RyR (J_per_mol).
C ALGBRC(179) is v_OC in component RyR (fmol_per_sec).
C ALGBRC(180) is v_CCI in component RyR (fmol_per_sec).
C ALGBRC(181) is v_CII in component RyR (fmol_per_sec).
C ALGBRC(182) is v_IO in component RyR (fmol_per_sec).
C CONSTS(154) is n_Ca_i in component SERCA (dimensionless).
C CONSTS(155) is n_Ca_SR in component SERCA (dimensionless).
C CONSTS(156) is n_H in component SERCA (dimensionless).
C CONSTS(162) is q_SERCA_Ca_complexes in component SERCA (fmol).
C ALGBRC(16) is c_Ca_i in component SERCA (mM).
C ALGBRC(17) is c_Ca_SR in component SERCA (mM).
C ALGBRC(18) is c_H in component SERCA (mM).
C ALGBRC(19) is c_MgADP in component SERCA (mM).
C ALGBRC(20) is c_MgATP in component SERCA (mM).
C ALGBRC(21) is c_Pi in component SERCA (mM).
C CONSTS(157) is vol_i in component SERCA (pL).
C CONSTS(163) is vol_sr in component SERCA (pL).
C CONSTS(167) is vol_isr in component SERCA (pL).
C ALGBRC(184) is mu_Ca_i in component SERCA (J_per_mol).
C ALGBRC(185) is mu_Ca_SR in component SERCA (J_per_mol).
C ALGBRC(186) is mu_H in component SERCA (J_per_mol).
C ALGBRC(187) is mu_MgADP in component SERCA (J_per_mol).
C ALGBRC(188) is mu_MgATP in component SERCA (J_per_mol).
C CONSTS(166) is mu_P1 in component SERCA (J_per_mol).
C ALGBRC(204) is v_P1 in component SERCA (fmol_per_sec).
C ALGBRC(189) is mu_Pi in component SERCA (J_per_mol).
C CONSTS(169) is mu_P2 in component SERCA (J_per_mol).
C ALGBRC(211) is v_P2 in component SERCA (fmol_per_sec).
C CONSTS(170) is mu_P2a in component SERCA (J_per_mol).
C ALGBRC(208) is v_P2a in component SERCA (fmol_per_sec).
C CONSTS(171) is mu_P4 in component SERCA (J_per_mol).
C ALGBRC(217) is v_P4 in component SERCA (fmol_per_sec).
C CONSTS(172) is mu_P5 in component SERCA (J_per_mol).
C ALGBRC(218) is v_P5 in component SERCA (fmol_per_sec).
C CONSTS(174) is mu_P6 in component SERCA (J_per_mol).
C ALGBRC(222) is v_P6 in component SERCA (fmol_per_sec).
C CONSTS(175) is mu_P8 in component SERCA (J_per_mol).
C ALGBRC(225) is v_P8 in component SERCA (fmol_per_sec).
C CONSTS(176) is mu_P9 in component SERCA (J_per_mol).
C ALGBRC(228) is v_P9 in component SERCA (fmol_per_sec).
C CONSTS(173) is mu_P10 in component SERCA (J_per_mol).
C ALGBRC(229) is v_P10 in component SERCA (fmol_per_sec).
C ALGBRC(191) is Af_R1_2 in component SERCA (J_per_mol).
C CONSTS(178) is Ar_R1_2 in component SERCA (J_per_mol).
C ALGBRC(201) is v_SERCA_R1_2 in component SERCA (fmol_per_sec).
C CONSTS(179) is Af_R5_6 in component SERCA (J_per_mol).
C ALGBRC(192) is Ar_R5_6 in component SERCA (J_per_mol).
C ALGBRC(203) is v_SERCA_R5_6 in component SERCA (fmol_per_sec).
C ALGBRC(193) is Af_R2_2a in component SERCA (J_per_mol).
C CONSTS(180) is Ar_R2_2a in component SERCA (J_per_mol).
C ALGBRC(206) is v_SERCA_R2_2a in component SERCA (fmol_per_sec).
C ALGBRC(194) is Af_R2_4 in component SERCA (J_per_mol).
C CONSTS(181) is Ar_R2_4 in component SERCA (J_per_mol).
C ALGBRC(207) is v_SERCA_R2_4 in component SERCA (fmol_per_sec).
C CONSTS(182) is Af_R4_5 in component SERCA (J_per_mol).
C ALGBRC(195) is Ar_R4_5 in component SERCA (J_per_mol).
C ALGBRC(209) is v_SERCA_R4_5 in component SERCA (fmol_per_sec).
C CONSTS(183) is Af_R6_8 in component SERCA (J_per_mol).
C ALGBRC(196) is Ar_R6_8 in component SERCA (J_per_mol).
C ALGBRC(212) is v_SERCA_R6_8 in component SERCA (fmol_per_sec).
C ALGBRC(197) is Af_R8_9 in component SERCA (J_per_mol).
C CONSTS(184) is Ar_R8_9 in component SERCA (J_per_mol).
C ALGBRC(221) is v_SERCA_R8_9 in component SERCA (fmol_per_sec).
C CONSTS(185) is Af_R9_10 in component SERCA (J_per_mol).
C ALGBRC(198) is Ar_R9_10 in component SERCA (J_per_mol).
C ALGBRC(224) is v_SERCA_R9_10 in component SERCA (fmol_per_sec).
C CONSTS(177) is Af_R10_1 in component SERCA (J_per_mol).
C ALGBRC(190) is Ar_R10_1 in component SERCA (J_per_mol).
C ALGBRC(199) is v_SERCA_R10_1 in component SERCA (fmol_per_sec).
C RATES(4) is d/dt q_Cai in component environment (fmol).
C RATES(2) is d/dt q_TRPN in component environment (fmol).
C RATES(3) is d/dt q_Ca_TRPN in component environment (fmol).
C RATES(5) is d/dt q_B in component environment (fmol).
C RATES(6) is d/dt q_U in component environment (fmol).
C RATES(7) is d/dt q_W in component environment (fmol).
C RATES(8) is d/dt q_S in component environment (fmol).
C RATES(9) is d/dt q_Cao in component environment (fmol).
C RATES(10) is d/dt q_Ki in component environment (fmol).
C RATES(11) is d/dt q_Ko in component environment (fmol).
C RATES(12) is d/dt q_PLB in component environment (fmol).
C RATES(13) is d/dt q_PKACI in component environment (fmol).
C RATES(14) is d/dt q_PLB_PKACI in component environment (fmol).
C RATES(15) is d/dt q_PP1 in component environment (fmol).
C RATES(16) is d/dt q_PLBp_PP1 in component environment (fmol).
C RATES(17) is d/dt q_PLBp in component environment (fmol).
C RATES(18) is d/dt q_Ip in component environment (fmol).
C RATES(19) is d/dt q_Ip_PP1 in component environment (fmol).
C RATES(20) is d/dt q_Ca_SR in component environment (fmol).
C RATES(21) is d/dt q_H in component environment (fmol).
C RATES(23) is d/dt q_MgADP in component environment (fmol).
C RATES(24) is d/dt q_MgATP in component environment (fmol).
C RATES(22) is d/dt q_P in component environment (fmol).
C RATES(1) is d/dt q_membrane in component environment (fC).
C RATES(42) is d/dt G_w in component crossbridge_TRPN (metre).
C RATES(41) is d/dt SL in component crossbridge_TRPN (metre).
C RATES(43) is d/dt G_s in component crossbridge_TRPN (metre).
C RATES(44) is d/dt Cdd in component crossbridge_TRPN (metre).
C RATES(25) is d/dt q_000_LCC in component environment (fmol).
C RATES(26) is d/dt q_010_LCC in component environment (fmol).
C RATES(27) is d/dt q_020_LCC in component environment (fmol).
C RATES(28) is d/dt q_100_LCC in component environment (fmol).
C RATES(29) is d/dt q_110_LCC in component environment (fmol).
C RATES(30) is d/dt q_120_LCC in component environment (fmol).
C RATES(31) is d/dt q_001_LCC in component environment (fmol).
C RATES(32) is d/dt q_011_LCC in component environment (fmol).
C RATES(33) is d/dt q_021_LCC in component environment (fmol).
C RATES(34) is d/dt q_101_LCC in component environment (fmol).
C RATES(35) is d/dt q_111_LCC in component environment (fmol).
C RATES(36) is d/dt q_121_LCC in component environment (fmol).
C RATES(40) is d/dt q_O_RyR in component environment (fmol).
C RATES(37) is d/dt q_C_RyR in component environment (fmol).
C RATES(38) is d/dt q_CI_RyR in component environment (fmol).
C RATES(39) is d/dt q_I_RyR in component environment (fmol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 52.0453
      CONSTS(2) = 2.90875
      CONSTS(3) = 0.332501
      CONSTS(4) = 0.00782356
      CONSTS(5) = 140824
      CONSTS(6) = 0.756518
      CONSTS(7) = 0.996063
      CONSTS(8) = 0.0247377
      CONSTS(9) = 0.0325707
      CONSTS(10) = 329.349
      CONSTS(11) = 16.6321
      CONSTS(12) = 21.8985
      CONSTS(13) = 0.141736
      CONSTS(14) = 0.00715766
      CONSTS(15) = 0.00942408
      CONSTS(16) = 31.7244
      CONSTS(17) = 157.535
      CONSTS(18) = 0.0136527
      CONSTS(19) = 0.0677958
      CONSTS(20) = 2.06197
      CONSTS(21) = 10.2392
      CONSTS(22) = 0.000887375
      CONSTS(23) = 0.00440648
      CONSTS(24) = 4498.81
      CONSTS(25) = 22339.9
      CONSTS(26) = 1.93608
      CONSTS(27) = 9.61406
      CONSTS(28) = 48530.3
      CONSTS(29) = 240989
      CONSTS(30) = 2450.77
      CONSTS(31) = 12169.9
      CONSTS(32) = 3226.79
      CONSTS(33) = 16023.4
      CONSTS(34) = 45.5263
      CONSTS(35) = 6.55904
      CONSTS(36) = 0.386674
      CONSTS(37) = 1.21269
      CONSTS(38) = 431.435
      CONSTS(39) = 2.85549e+06
      CONSTS(40) = 8.83262
      CONSTS(41) = 0.12618
      CONSTS(42) = 883.262
      CONSTS(43) = 73.6052
      CONSTS(44) = 1.57289e-05
      CONSTS(45) = 175018
      CONSTS(46) = 2.44679e+06
      CONSTS(47) = 175018
      CONSTS(48) = 3487.81
      CONSTS(49) = 4.39557e+06
      CONSTS(50) = 1.01676e+07
      CONSTS(51) = 1.01676e+07
      CONSTS(52) = 0.00441447
      CONSTS(53) = 1.27424
      CONSTS(54) = 1.27424
      CONSTS(55) = 0.11171
      CONSTS(56) = 6.06164
      CONSTS(57) = 2.27312
      CONSTS(58) = 14.8627
      CONSTS(59) = 3.71567e-06
      CONSTS(60) = 0.0343245
      CONSTS(61) = 0.227859
      CONSTS(62) = 0.000825154
      CONSTS(63) = 0.00547767
      CONSTS(64) = 1.47795
      CONSTS(65) = 29.2665
      CONSTS(66) = 22.2281
      CONSTS(67) = 0.297629
      CONSTS(68) = 5.89367
      CONSTS(69) = 4.47629
      CONSTS(70) = 3434.28
      CONSTS(71) = 68005.7
      CONSTS(72) = 51650.9
      CONSTS(73) = 691.593
      CONSTS(74) = 13695
      CONSTS(75) = 10401.4
      CONSTS(76) = 0.00235741
      CONSTS(77) = 0.395191
      CONSTS(78) = 0.638527
      CONSTS(79) = 0.361989
      CONSTS(80) = 0.203757
      CONSTS(81) = 0.326014
      CONSTS(82) = 5.41093
      CONSTS(83) = 0.0673793
      CONSTS(84) = 0.980701
      CONSTS(85) = 3962.58
      CONSTS(86) = 39.6258
      CONSTS(87) = 0.06793
      CONSTS(88) = 6.793
      CONSTS(89) = 38162.5
      CONSTS(90) = 128177
      CONSTS(91) = 10585.3
      CONSTS(92) = 147985
      CONSTS(93) = 0.728251
      CONSTS(94) = 29565.2
      CONSTS(95) = 4226.28
      CONSTS(96) = 2547.31
      CONSTS(97) = 4643.82
      CONSTS(98) = 212.798
      CONSTS(99) = 1254.32
      CONSTS(100) = 1.89443e-05
      CONSTS(101) = 0.0257106
      CONSTS(102) = 2
      CONSTS(103) = 2.1404
      CONSTS(104) = -1.1495
      CONSTS(105) = 0.72162
      CONSTS(106) = 4.2933
      CONSTS(107) = -2.1404
      CONSTS(108) = 1.8993
      CONSTS(109) = -0.52288
      CONSTS(110) = 0
      CONSTS(111) = 34.4
      CONSTS(112) = 500
      CONSTS(113) = 1.381e5
      CONSTS(114) = 96485
      STATES(1) = -13039
      STATES(2) = 2.57
      STATES(3) = 1e-18
      STATES(4) = 6.82e-3
      STATES(5) = 1
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 6.84
      STATES(10) = 5.51E+03
      STATES(11) = 2.05E+02
      STATES(12) = 4.028E+00
      STATES(13) = 2.234E-03
      STATES(14) = 1e-18
      STATES(15) = 3.382E-02
      STATES(16) = 1e-18
      STATES(17) = 1e-18
      STATES(18) = 1.999E-03
      STATES(19) = 1e-18
      STATES(20) = 0.641
      STATES(21) = 0.004028
      STATES(22) = 570
      STATES(23) = 1.3794
      STATES(24) = 3.8
      STATES(25) = 1.3798e-05
      STATES(26) = 1.3798e-05
      STATES(27) = 1.3798e-05
      STATES(28) = 1.3798e-05
      STATES(29) = 1.3798e-05
      STATES(30) = 1.3798e-05
      STATES(31) = 1.3798e-05
      STATES(32) = 1.3798e-05
      STATES(33) = 1.3798e-05
      STATES(34) = 1.3798e-05
      STATES(35) = 1.3798e-05
      STATES(36) = 1.3798e-05
      STATES(37) = 4.13941e-05
      STATES(38) = 4.13941e-05
      STATES(39) = 4.13941e-05
      STATES(40) = 4.13941e-05
      CONSTS(115) = 2.83036e-14
      CONSTS(116) = 2.83036e-14
      CONSTS(117) = 2.83036e-14
      CONSTS(118) = 2.83036e-14
      CONSTS(119) = 2.83036e-14
      CONSTS(120) = 2.83036e-14
      CONSTS(121) = 2.83036e-14
      CONSTS(122) = 2.83036e-14
      CONSTS(123) = 2.83036e-14
      STATES(41) = 1e-6
      CONSTS(124) = 4e-2
      CONSTS(125) = 4.1e-2
      CONSTS(126) = 1e8
      CONSTS(127) = 0
      CONSTS(128) = 2
      CONSTS(129) = 1
      CONSTS(130) = 8.31
      CONSTS(131) = 310
      CONSTS(132) = 1
      CONSTS(133) = -1
      CONSTS(134) = 1.1
      CONSTS(135) = -0.0118
      CONSTS(136) = 1
      CONSTS(137) = 2e-6
      CONSTS(138) = 1e-6
      CONSTS(139) = 0.25
      CONSTS(140) = 0.5
      CONSTS(141) = 25
      CONSTS(142) = 2.23
      CONSTS(143) = 26
      CONSTS(144) = 4
      STATES(42) = 1e-6
      STATES(43) = 1e-6
      CONSTS(145) = 0
      CONSTS(146) = 1
      CONSTS(147) = 40.5
      STATES(44) = 0
      CONSTS(148) = 200e3
      CONSTS(149) = 200e3
      CONSTS(150) = 7
      CONSTS(151) = 2.1
      CONSTS(152) = 1
      CONSTS(153) = 2
      CONSTS(154) = 2
      CONSTS(155) = 2
      CONSTS(156) = 2
      CONSTS(157) = 34.0
      CONSTS(158) = ( CONSTS(142)*CONSTS(143)*(1.00000 - CONSTS(140)))/CONSTS(140)
      CONSTS(159) = ( CONSTS(142)*CONSTS(144)*(1.00000 - CONSTS(139))*CONSTS(140))/CONSTS(139)
      CONSTS(160) = ( CONSTS(141)*CONSTS(139))/( (1.00000 - CONSTS(139))*CONSTS(140)+CONSTS(139))
      CONSTS(161) =  CONSTS(135)*CONSTS(130)*CONSTS(131)
      CONSTS(162) =  CONSTS(154)*CONSTS(118)+ CONSTS(154)*CONSTS(119)+ CONSTS(155)*CONSTS(120)
      CONSTS(163) =  CONSTS(157)*0.0350000
      CONSTS(185) = CONSTS(145)
      CONSTS(164) = CONSTS(160)
      CONSTS(165) =  - CONSTS(160)*CONSTS(145)*CONSTS(146)
      CONSTS(166) =  CONSTS(130)*CONSTS(131)*log( CONSTS(89)*CONSTS(115))
      CONSTS(167) = CONSTS(157)+CONSTS(163)
      CONSTS(168) =  - CONSTS(164)*CONSTS(145)*CONSTS(146)
      CONSTS(169) =  CONSTS(130)*CONSTS(131)*log( CONSTS(90)*CONSTS(116))
      CONSTS(170) =  CONSTS(130)*CONSTS(131)*log( CONSTS(91)*CONSTS(117))
      CONSTS(171) =  CONSTS(130)*CONSTS(131)*log( CONSTS(92)*CONSTS(118))
      CONSTS(172) =  CONSTS(130)*CONSTS(131)*log( CONSTS(93)*CONSTS(119))
      CONSTS(173) =  CONSTS(130)*CONSTS(131)*log( CONSTS(97)*CONSTS(123))
      CONSTS(174) =  CONSTS(130)*CONSTS(131)*log( CONSTS(94)*CONSTS(120))
      CONSTS(175) =  CONSTS(130)*CONSTS(131)*log( CONSTS(95)*CONSTS(121))
      CONSTS(176) =  CONSTS(130)*CONSTS(131)*log( CONSTS(96)*CONSTS(122))
      CONSTS(177) = CONSTS(173)
      CONSTS(178) = CONSTS(169)
      CONSTS(179) = CONSTS(172)
      CONSTS(180) = CONSTS(170)
      CONSTS(181) = CONSTS(171)
      CONSTS(182) = CONSTS(171)
      CONSTS(183) = CONSTS(174)
      CONSTS(184) = CONSTS(176)
      CONSTS(185) = CONSTS(176)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(41) = CONSTS(185)
      ALGBRC(8) =  CONSTS(158)*STATES(42)*CONSTS(146)
      ALGBRC(23) = - CONSTS(165) - ALGBRC(8)
      ALGBRC(27) = ALGBRC(23)/CONSTS(146)
      RATES(42) = ALGBRC(27)
      ALGBRC(9) =  CONSTS(159)*STATES(43)*CONSTS(146)
      ALGBRC(24) = - CONSTS(168) - ALGBRC(9)
      ALGBRC(28) = ALGBRC(24)/CONSTS(146)
      RATES(43) = ALGBRC(28)
      ALGBRC(31) =  CONSTS(130)*CONSTS(131)*log( CONSTS(55)*STATES(3))
      ALGBRC(32) =  CONSTS(130)*CONSTS(131)*log( CONSTS(56)*STATES(5))
      ALGBRC(35) =  CONSTS(130)*CONSTS(131)*log( CONSTS(57)*STATES(6))
      ALGBRC(44) =  CONSTS(2)*(EXP((ALGBRC(32)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(35)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))))
      RATES(5) = - ALGBRC(44)
      ALGBRC(2) =  CONSTS(130)*CONSTS(131)*log( CONSTS(76)*STATES(12))
      ALGBRC(22) =  CONSTS(130)*CONSTS(131)*log( CONSTS(77)*STATES(13))
      ALGBRC(26) =  CONSTS(130)*CONSTS(131)*log( CONSTS(78)*STATES(14))
      ALGBRC(45) =  CONSTS(34)*(EXP((ALGBRC(2)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(36) =  CONSTS(130)*CONSTS(131)*log( CONSTS(81)*STATES(17))
      ALGBRC(48) =  CONSTS(35)*(EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(36)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))))
      RATES(13) = ALGBRC(48) - ALGBRC(45)
      RATES(14) = ALGBRC(45) - ALGBRC(48)
      ALGBRC(38) =  CONSTS(130)*CONSTS(131)*log( CONSTS(58)*STATES(7))
      ALGBRC(47) =  CONSTS(3)*(EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(41) =  CONSTS(130)*CONSTS(131)*log( CONSTS(59)*STATES(8))
      ALGBRC(50) =  CONSTS(4)*(EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))))
      RATES(7) = ALGBRC(47) - ALGBRC(50)
      ALGBRC(30) =  CONSTS(130)*CONSTS(131)*log( CONSTS(79)*STATES(15))
      ALGBRC(33) =  CONSTS(130)*CONSTS(131)*log( CONSTS(80)*STATES(16))
      ALGBRC(51) =  CONSTS(36)*(EXP((ALGBRC(36)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))))
      RATES(17) = ALGBRC(48) - ALGBRC(51)
      ALGBRC(53) =  CONSTS(5)*(EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))))
      RATES(6) = (ALGBRC(44) - ALGBRC(47))+ALGBRC(53)
      RATES(8) = ALGBRC(50) - ALGBRC(53)
      ALGBRC(54) =  CONSTS(37)*(EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(2)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))))
      RATES(12) = ALGBRC(54) - ALGBRC(45)
      RATES(16) = ALGBRC(51) - ALGBRC(54)
      ALGBRC(39) =  CONSTS(130)*CONSTS(131)*log( CONSTS(82)*STATES(18))
      ALGBRC(42) =  CONSTS(130)*CONSTS(131)*log( CONSTS(83)*STATES(19))
      ALGBRC(56) =  CONSTS(38)*(EXP((ALGBRC(30)+ALGBRC(39))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(42)/( CONSTS(130)*CONSTS(131))))
      RATES(15) = (ALGBRC(54) - ALGBRC(51)) - ALGBRC(56)
      RATES(18) = - ALGBRC(56)
      RATES(19) = ALGBRC(56)
      ALGBRC(58) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))
      ALGBRC(88) = ALGBRC(58)
      ALGBRC(52) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(138) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(139) =  CONSTS(28)*(EXP(ALGBRC(88)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(138)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(46) = STATES(1)/CONSTS(113)
      ALGBRC(70) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(102) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(103) =  CONSTS(10)*(EXP(ALGBRC(70)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(102)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(76) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(114) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(115) =  CONSTS(16)*(EXP(ALGBRC(76)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(114)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(80) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(122) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(123) =  CONSTS(20)*(EXP(ALGBRC(80)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(122)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(141) = ((- ALGBRC(103) - ALGBRC(115)) - ALGBRC(123)) - ALGBRC(139)
      RATES(25) = ALGBRC(141)
      ALGBRC(71) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(104) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(105) =  CONSTS(13)*(EXP(ALGBRC(71)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(104)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(77) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(116) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(117) =  CONSTS(18)*(EXP(ALGBRC(77)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(116)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(81) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(124) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(125) =  CONSTS(22)*(EXP(ALGBRC(81)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(124)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(142) = ((- ALGBRC(105) - ALGBRC(117)) - ALGBRC(125)) - - ALGBRC(139)
      RATES(31) = ALGBRC(142)
      ALGBRC(59) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))
      ALGBRC(89) = ALGBRC(59)
      ALGBRC(140) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(143) =  CONSTS(30)*(EXP(ALGBRC(89)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(140)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(72) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(106) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(107) =  CONSTS(11)*(EXP(ALGBRC(72)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(106)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(84) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(130) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(131) =  CONSTS(24)*(EXP(ALGBRC(84)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(130)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(145) = ((- ALGBRC(107)+ALGBRC(115)) - ALGBRC(143)) - ALGBRC(131)
      RATES(26) = ALGBRC(145)
      ALGBRC(73) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(108) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(109) =  CONSTS(14)*(EXP(ALGBRC(73)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(108)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(85) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(132) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(133) =  CONSTS(26)*(EXP(ALGBRC(85)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(132)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(146) = ((- ALGBRC(109)+ALGBRC(117)) - - ALGBRC(143)) - ALGBRC(133)
      RATES(32) = ALGBRC(146)
      ALGBRC(60) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))
      ALGBRC(90) = ALGBRC(60)
      ALGBRC(144) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(147) =  CONSTS(32)*(EXP(ALGBRC(90)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(144)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(74) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(110) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(111) =  CONSTS(12)*(EXP(ALGBRC(74)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(110)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(149) = ((- ALGBRC(111)+ALGBRC(123)) - ALGBRC(147))+ALGBRC(131)
      RATES(27) = ALGBRC(149)
      ALGBRC(75) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(112) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(113) =  CONSTS(15)*(EXP(ALGBRC(75)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(112)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(150) = ((- ALGBRC(113)+ALGBRC(125)) - - ALGBRC(147))+ALGBRC(133)
      RATES(33) = ALGBRC(150)
      ALGBRC(61) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))
      ALGBRC(91) = ALGBRC(61)
      ALGBRC(148) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(151) =  CONSTS(29)*(EXP(ALGBRC(91)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(148)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(78) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(118) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(119) =  CONSTS(17)*(EXP(ALGBRC(78)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(118)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(82) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(126) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(127) =  CONSTS(21)*(EXP(ALGBRC(82)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(126)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(153) = ((- ALGBRC(127)+ALGBRC(103)) - ALGBRC(119)) - ALGBRC(151)
      RATES(28) = ALGBRC(153)
      ALGBRC(79) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(120) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(121) =  CONSTS(19)*(EXP(ALGBRC(79)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(120)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(83) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(128) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(129) =  CONSTS(23)*(EXP(ALGBRC(83)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(128)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(154) = ((- ALGBRC(129)+ALGBRC(105)) - ALGBRC(121)) - - ALGBRC(151)
      RATES(34) = ALGBRC(154)
      ALGBRC(62) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))
      ALGBRC(92) = ALGBRC(62)
      ALGBRC(152) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(155) =  CONSTS(31)*(EXP(ALGBRC(92)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(152)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(86) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(134) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(135) =  CONSTS(25)*(EXP(ALGBRC(86)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(134)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(157) = ((ALGBRC(119)+ALGBRC(107)) - ALGBRC(155)) - ALGBRC(135)
      RATES(29) = ALGBRC(157)
      ALGBRC(64) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))
      ALGBRC(93) = ALGBRC(64)
      ALGBRC(156) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(158) =  CONSTS(33)*(EXP(ALGBRC(93)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(156)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(160) = ((ALGBRC(127)+ALGBRC(111)) - ALGBRC(158))+ALGBRC(135)
      RATES(30) = ALGBRC(160)
      ALGBRC(63) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))
      ALGBRC(66) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
      ALGBRC(49) =  CONSTS(130)*CONSTS(131)*log( CONSTS(61)*STATES(9))
      ALGBRC(98) = ALGBRC(49)+ALGBRC(63)
      ALGBRC(94) =  CONSTS(128)*CONSTS(114)*ALGBRC(46)
      ALGBRC(159) = TERNRY(ALGBRC(94).EQ.0.00000,  CONSTS(6)*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(6)*ALGBRC(94))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(94)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(65) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))
      ALGBRC(67) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
      ALGBRC(99) = ALGBRC(49)+ALGBRC(65)
      ALGBRC(95) =  CONSTS(128)*CONSTS(114)*ALGBRC(46)
      ALGBRC(161) = TERNRY(ALGBRC(95).EQ.0.00000,  CONSTS(7)*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(7)*ALGBRC(95))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(95)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(162) = ALGBRC(161)+ALGBRC(159)
      RATES(9) = ALGBRC(162)
      ALGBRC(57) =  CONSTS(130)*CONSTS(131)*log( CONSTS(62)*STATES(10))
      ALGBRC(68) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
      ALGBRC(55) =  CONSTS(130)*CONSTS(131)*log( CONSTS(63)*STATES(11))
      ALGBRC(100) = ALGBRC(55)+ALGBRC(63)
      ALGBRC(96) =  CONSTS(129)*CONSTS(114)*ALGBRC(46)
      ALGBRC(163) = TERNRY(ALGBRC(96).EQ.0.00000,  CONSTS(8)*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(8)*ALGBRC(96))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(96)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(87) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(136) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(137) =  CONSTS(27)*(EXP(ALGBRC(87)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(136)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(167) = ((ALGBRC(121)+ALGBRC(109)) - ((((- ALGBRC(155)+ALGBRC(159)) - ALGBRC(159))+ALGBRC(163)) - ALGBRC(163))) - ALGBRC(137)
      RATES(35) = ALGBRC(167)
      ALGBRC(69) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
      ALGBRC(101) = ALGBRC(55)+ALGBRC(65)
      ALGBRC(97) =  CONSTS(129)*CONSTS(114)*ALGBRC(46)
      ALGBRC(165) = TERNRY(ALGBRC(97).EQ.0.00000,  CONSTS(9)*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(9)*ALGBRC(97))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(97)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(168) = - ALGBRC(163) - ALGBRC(165)
      RATES(10) = ALGBRC(168)
      ALGBRC(169) = ALGBRC(165)+ALGBRC(163)
      RATES(11) = ALGBRC(169)
      ALGBRC(172) = ((ALGBRC(129)+ALGBRC(113)) - ((((- ALGBRC(158)+ALGBRC(161)) - ALGBRC(161))+ALGBRC(165)) - ALGBRC(165)))+ALGBRC(137)
      RATES(36) = ALGBRC(172)
      ALGBRC(170) =  CONSTS(114)*((((((((((((((((((((((((((((((((((((( CONSTS(109)*ALGBRC(125) -  CONSTS(105)*ALGBRC(125))+ CONSTS(110)*ALGBRC(133)) -  CONSTS(106)*ALGBRC(133))+ CONSTS(108)*ALGBRC(121)) -  CONSTS(104)*ALGBRC(121))+ CONSTS(109)*ALGBRC(129)) -  CONSTS(105)*ALGBRC(129))+ CONSTS(110)*ALGBRC(137)) -  CONSTS(106)*ALGBRC(137))+ CONSTS(108)*ALGBRC(117)) -  CONSTS(104)*ALGBRC(117))+ CONSTS(107)*ALGBRC(113)) -  CONSTS(103)*ALGBRC(113))+ CONSTS(107)*ALGBRC(109)) -  CONSTS(103)*ALGBRC(109))+ CONSTS(107)*ALGBRC(105)) -  CONSTS(103)*ALGBRC(105))+ CONSTS(110)*ALGBRC(135)) -  CONSTS(106)*ALGBRC(135))+ CONSTS(109)*ALGBRC(127)) -  CONSTS(105)*ALGBRC(127))+ CONSTS(108)*ALGBRC(119)) -  CONSTS(104)*ALGBRC(119))+ CONSTS(110)*ALGBRC(131)) -  CONSTS(106)*ALGBRC(131))+ CONSTS(109)*ALGBRC(123)) -  CONSTS(105)*ALGBRC(123))+ CONSTS(108)*ALGBRC(115)) -  CONSTS(104)*ALGBRC(115))+ CONSTS(107)*ALGBRC(111)) -  CONSTS(103)*ALGBRC(111))+ CONSTS(107)*ALGBRC(107)) -  CONSTS(103)*ALGBRC(107))+ CONSTS(107)*ALGBRC(103)) -  CONSTS(103)*ALGBRC(103)) - ( CONSTS(128)*ALGBRC(159)+ CONSTS(128)*ALGBRC(161))) - ( CONSTS(129)*ALGBRC(163)+ CONSTS(129)*ALGBRC(165)))
      ALGBRC(10) = TERNRY(VOI.GT.CONSTS(124).AND.VOI.LT.CONSTS(125), CONSTS(126), CONSTS(127))
      ALGBRC(173) = - ALGBRC(170)+ALGBRC(10)
      RATES(1) = ALGBRC(173)
      ALGBRC(171) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(174) =  CONSTS(130)*CONSTS(131)*log( CONSTS(85)*STATES(37))
      ALGBRC(177) =  CONSTS(130)*CONSTS(131)*log( CONSTS(88)*STATES(40))
      ALGBRC(179) =  CONSTS(40)*(EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(174)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(175) =  CONSTS(130)*CONSTS(131)*log( CONSTS(86)*STATES(38))
      ALGBRC(180) =  CONSTS(41)*(EXP((ALGBRC(174)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(175)/( CONSTS(130)*CONSTS(131))))
      RATES(37) = ALGBRC(179) - ALGBRC(180)
      ALGBRC(176) =  CONSTS(130)*CONSTS(131)*log( CONSTS(87)*STATES(39))
      ALGBRC(181) =  CONSTS(42)*(EXP((ALGBRC(175)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))))
      RATES(38) = ALGBRC(180) - ALGBRC(181)
      ALGBRC(182) =  CONSTS(43)*(EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(177)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
      RATES(40) = ALGBRC(182) - ALGBRC(179)
      RATES(39) = ALGBRC(181) - ALGBRC(182)
      ALGBRC(189) =  CONSTS(130)*CONSTS(131)*log( CONSTS(101)*STATES(22))
      ALGBRC(190) = CONSTS(166)+ALGBRC(189)
      ALGBRC(199) =  CONSTS(52)*(EXP(CONSTS(177)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(190)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(200) = ALGBRC(199)
      RATES(22) = ALGBRC(200)
      ALGBRC(188) =  CONSTS(130)*CONSTS(131)*log( CONSTS(99)*STATES(24))
      ALGBRC(191) = CONSTS(166)+ALGBRC(188)
      ALGBRC(201) =  CONSTS(44)*(EXP(ALGBRC(191)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(178)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(202) = - ALGBRC(201)
      RATES(24) = ALGBRC(202)
      ALGBRC(187) =  CONSTS(130)*CONSTS(131)*log( CONSTS(100)*STATES(23))
      ALGBRC(192) = ALGBRC(187)+CONSTS(174)
      ALGBRC(203) =  CONSTS(48)*(EXP(CONSTS(179)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(192)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(205) = ALGBRC(203)
      RATES(23) = ALGBRC(205)
      ALGBRC(166) =  CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
      ALGBRC(178) =  CONSTS(39)*EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131)))*(EXP(ALGBRC(166)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(171)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(185) =  CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
      ALGBRC(196) = CONSTS(175)+ CONSTS(155)*ALGBRC(185)
      ALGBRC(212) =  CONSTS(49)*(EXP(CONSTS(183)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(196)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(220) =  CONSTS(155)*ALGBRC(212)
      RATES(20) = - ALGBRC(178)+ALGBRC(220)
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      ALGBRC(219) = ALGBRC(214)/CONSTS(146)
      RATES(44) = ALGBRC(219)
      ALGBRC(186) =  CONSTS(130)*CONSTS(131)*log( CONSTS(98)*STATES(21))
      ALGBRC(193) = CONSTS(169)+ALGBRC(186)
      ALGBRC(206) =  CONSTS(46)*(EXP(ALGBRC(193)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(180)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(195) = CONSTS(172)+ CONSTS(156)*ALGBRC(186)
      ALGBRC(209) =  CONSTS(47)*(EXP(CONSTS(182)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(195)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(197) = CONSTS(175)+ CONSTS(156)*ALGBRC(186)
      ALGBRC(221) =  CONSTS(50)*(EXP(ALGBRC(197)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(184)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(198) = ALGBRC(186)+CONSTS(173)
      ALGBRC(224) =  CONSTS(51)*(EXP(CONSTS(185)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(198)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(227) = (( CONSTS(156)*ALGBRC(209) - ALGBRC(206)) -  CONSTS(156)*ALGBRC(221))+ALGBRC(224)
      RATES(21) = ALGBRC(227)
      ALGBRC(25) =  CONSTS(130)*CONSTS(131)*log( CONSTS(53)*STATES(2))
      ALGBRC(29) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(34) =  (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(8)*(STATES(43)+CONSTS(137))
      ALGBRC(37) =  (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(7)*STATES(42)
      ALGBRC(40) = ALGBRC(34)+ALGBRC(37)
      ALGBRC(43) = ALGBRC(40)
      ALGBRC(223) =  (( CONSTS(151)*ALGBRC(213))/CONSTS(137))*ALGBRC(219)
      ALGBRC(226) = ALGBRC(223)
      ALGBRC(230) = ALGBRC(43)+ALGBRC(226)
      ALGBRC(231) =  ALGBRC(230)*1000.00
      ALGBRC(232) =  CONSTS(136)*CONSTS(161)*ALGBRC(231)
      ALGBRC(233) =  CONSTS(1)*(EXP((ALGBRC(29)+ALGBRC(25))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(31)+ALGBRC(232))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(164) = (( (CONSTS(102))*ALGBRC(151)+ CONSTS(102)*ALGBRC(155)+ CONSTS(102)*ALGBRC(147)+ CONSTS(102)*ALGBRC(143)+ CONSTS(102)*ALGBRC(139)+ CONSTS(102)*ALGBRC(158)) - ALGBRC(159)) - ALGBRC(161)
      ALGBRC(183) = (( CONSTS(153)*ALGBRC(179) -  CONSTS(152)*ALGBRC(180)) -  CONSTS(153)*ALGBRC(181))+ CONSTS(152)*ALGBRC(182)
      ALGBRC(184) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(194) = CONSTS(169)+ CONSTS(154)*ALGBRC(184)
      ALGBRC(207) =  CONSTS(45)*(EXP(ALGBRC(194)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(181)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(210) =  - CONSTS(154)*ALGBRC(207)
      RATES(4) = - ALGBRC(233)+ALGBRC(164)+ALGBRC(178)+ALGBRC(183)+ALGBRC(210)
      RATES(2) = - ALGBRC(233)
      RATES(3) = ALGBRC(233)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(8) =  CONSTS(158)*STATES(42)*CONSTS(146)
      ALGBRC(23) = - CONSTS(165) - ALGBRC(8)
      ALGBRC(27) = ALGBRC(23)/CONSTS(146)
      ALGBRC(9) =  CONSTS(159)*STATES(43)*CONSTS(146)
      ALGBRC(24) = - CONSTS(168) - ALGBRC(9)
      ALGBRC(28) = ALGBRC(24)/CONSTS(146)
      ALGBRC(31) =  CONSTS(130)*CONSTS(131)*log( CONSTS(55)*STATES(3))
      ALGBRC(32) =  CONSTS(130)*CONSTS(131)*log( CONSTS(56)*STATES(5))
      ALGBRC(35) =  CONSTS(130)*CONSTS(131)*log( CONSTS(57)*STATES(6))
      ALGBRC(44) =  CONSTS(2)*(EXP((ALGBRC(32)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(35)+ CONSTS(134)*ALGBRC(31))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(2) =  CONSTS(130)*CONSTS(131)*log( CONSTS(76)*STATES(12))
      ALGBRC(22) =  CONSTS(130)*CONSTS(131)*log( CONSTS(77)*STATES(13))
      ALGBRC(26) =  CONSTS(130)*CONSTS(131)*log( CONSTS(78)*STATES(14))
      ALGBRC(45) =  CONSTS(34)*(EXP((ALGBRC(2)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(36) =  CONSTS(130)*CONSTS(131)*log( CONSTS(81)*STATES(17))
      ALGBRC(48) =  CONSTS(35)*(EXP(ALGBRC(26)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(36)+ALGBRC(22))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(38) =  CONSTS(130)*CONSTS(131)*log( CONSTS(58)*STATES(7))
      ALGBRC(47) =  CONSTS(3)*(EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(41) =  CONSTS(130)*CONSTS(131)*log( CONSTS(59)*STATES(8))
      ALGBRC(50) =  CONSTS(4)*(EXP(ALGBRC(38)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(30) =  CONSTS(130)*CONSTS(131)*log( CONSTS(79)*STATES(15))
      ALGBRC(33) =  CONSTS(130)*CONSTS(131)*log( CONSTS(80)*STATES(16))
      ALGBRC(51) =  CONSTS(36)*(EXP((ALGBRC(36)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(53) =  CONSTS(5)*(EXP(ALGBRC(41)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(35)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(54) =  CONSTS(37)*(EXP(ALGBRC(33)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(2)+ALGBRC(30))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(39) =  CONSTS(130)*CONSTS(131)*log( CONSTS(82)*STATES(18))
      ALGBRC(42) =  CONSTS(130)*CONSTS(131)*log( CONSTS(83)*STATES(19))
      ALGBRC(56) =  CONSTS(38)*(EXP((ALGBRC(30)+ALGBRC(39))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(42)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(58) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))
      ALGBRC(88) = ALGBRC(58)
      ALGBRC(52) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(138) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(139) =  CONSTS(28)*(EXP(ALGBRC(88)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(138)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(46) = STATES(1)/CONSTS(113)
      ALGBRC(70) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(102) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(103) =  CONSTS(10)*(EXP(ALGBRC(70)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(102)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(76) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(114) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(115) =  CONSTS(16)*(EXP(ALGBRC(76)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(114)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(80) =  CONSTS(130)*CONSTS(131)*log( CONSTS(64)*STATES(25))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(122) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(123) =  CONSTS(20)*(EXP(ALGBRC(80)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(122)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(141) = ((- ALGBRC(103) - ALGBRC(115)) - ALGBRC(123)) - ALGBRC(139)
      ALGBRC(71) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(104) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(105) =  CONSTS(13)*(EXP(ALGBRC(71)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(104)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(77) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(116) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(117) =  CONSTS(18)*(EXP(ALGBRC(77)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(116)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(81) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(124) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(125) =  CONSTS(22)*(EXP(ALGBRC(81)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(124)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(142) = ((- ALGBRC(105) - ALGBRC(117)) - ALGBRC(125)) - - ALGBRC(139)
      ALGBRC(59) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))
      ALGBRC(89) = ALGBRC(59)
      ALGBRC(140) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(143) =  CONSTS(30)*(EXP(ALGBRC(89)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(140)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(72) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(106) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(107) =  CONSTS(11)*(EXP(ALGBRC(72)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(106)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(84) =  CONSTS(130)*CONSTS(131)*log( CONSTS(65)*STATES(26))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(130) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(131) =  CONSTS(24)*(EXP(ALGBRC(84)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(130)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(145) = ((- ALGBRC(107)+ALGBRC(115)) - ALGBRC(143)) - ALGBRC(131)
      ALGBRC(73) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(108) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(109) =  CONSTS(14)*(EXP(ALGBRC(73)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(108)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(85) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(132) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(133) =  CONSTS(26)*(EXP(ALGBRC(85)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(132)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(146) = ((- ALGBRC(109)+ALGBRC(117)) - - ALGBRC(143)) - ALGBRC(133)
      ALGBRC(60) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))
      ALGBRC(90) = ALGBRC(60)
      ALGBRC(144) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(147) =  CONSTS(32)*(EXP(ALGBRC(90)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(144)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(74) =  CONSTS(130)*CONSTS(131)*log( CONSTS(66)*STATES(27))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(110) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(111) =  CONSTS(12)*(EXP(ALGBRC(74)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(110)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(149) = ((- ALGBRC(111)+ALGBRC(123)) - ALGBRC(147))+ALGBRC(131)
      ALGBRC(75) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))+ CONSTS(103)*CONSTS(114)*ALGBRC(46)
      ALGBRC(112) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(107)*CONSTS(114)*ALGBRC(46)
      ALGBRC(113) =  CONSTS(15)*(EXP(ALGBRC(75)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(112)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(150) = ((- ALGBRC(113)+ALGBRC(125)) - - ALGBRC(147))+ALGBRC(133)
      ALGBRC(61) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))
      ALGBRC(91) = ALGBRC(61)
      ALGBRC(148) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(151) =  CONSTS(29)*(EXP(ALGBRC(91)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(148)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(78) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(118) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(119) =  CONSTS(17)*(EXP(ALGBRC(78)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(118)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(82) =  CONSTS(130)*CONSTS(131)*log( CONSTS(67)*STATES(28))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(126) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(127) =  CONSTS(21)*(EXP(ALGBRC(82)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(126)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(153) = ((- ALGBRC(127)+ALGBRC(103)) - ALGBRC(119)) - ALGBRC(151)
      ALGBRC(79) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(104)*CONSTS(114)*ALGBRC(46)
      ALGBRC(120) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(108)*CONSTS(114)*ALGBRC(46)
      ALGBRC(121) =  CONSTS(19)*(EXP(ALGBRC(79)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(120)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(83) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))+ CONSTS(105)*CONSTS(114)*ALGBRC(46)
      ALGBRC(128) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(109)*CONSTS(114)*ALGBRC(46)
      ALGBRC(129) =  CONSTS(23)*(EXP(ALGBRC(83)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(128)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(154) = ((- ALGBRC(129)+ALGBRC(105)) - ALGBRC(121)) - - ALGBRC(151)
      ALGBRC(62) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))
      ALGBRC(92) = ALGBRC(62)
      ALGBRC(152) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(155) =  CONSTS(31)*(EXP(ALGBRC(92)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(152)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(86) =  CONSTS(130)*CONSTS(131)*log( CONSTS(68)*STATES(29))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(134) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(135) =  CONSTS(25)*(EXP(ALGBRC(86)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(134)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(157) = ((ALGBRC(119)+ALGBRC(107)) - ALGBRC(155)) - ALGBRC(135)
      ALGBRC(64) =  CONSTS(130)*CONSTS(131)*log( CONSTS(69)*STATES(30))
      ALGBRC(93) = ALGBRC(64)
      ALGBRC(156) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(102)*ALGBRC(52)
      ALGBRC(158) =  CONSTS(33)*(EXP(ALGBRC(93)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(156)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(160) = ((ALGBRC(127)+ALGBRC(111)) - ALGBRC(158))+ALGBRC(135)
      ALGBRC(63) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))
      ALGBRC(66) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
      ALGBRC(49) =  CONSTS(130)*CONSTS(131)*log( CONSTS(61)*STATES(9))
      ALGBRC(98) = ALGBRC(49)+ALGBRC(63)
      ALGBRC(94) =  CONSTS(128)*CONSTS(114)*ALGBRC(46)
      ALGBRC(159) = TERNRY(ALGBRC(94).EQ.0.00000,  CONSTS(6)*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(6)*ALGBRC(94))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(94)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(66)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(98)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(65) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))
      ALGBRC(67) = ALGBRC(52)+ CONSTS(128)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
      ALGBRC(99) = ALGBRC(49)+ALGBRC(65)
      ALGBRC(95) =  CONSTS(128)*CONSTS(114)*ALGBRC(46)
      ALGBRC(161) = TERNRY(ALGBRC(95).EQ.0.00000,  CONSTS(7)*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(7)*ALGBRC(95))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(95)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(67)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(99)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(162) = ALGBRC(161)+ALGBRC(159)
      ALGBRC(57) =  CONSTS(130)*CONSTS(131)*log( CONSTS(62)*STATES(10))
      ALGBRC(68) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(63)
      ALGBRC(55) =  CONSTS(130)*CONSTS(131)*log( CONSTS(63)*STATES(11))
      ALGBRC(100) = ALGBRC(55)+ALGBRC(63)
      ALGBRC(96) =  CONSTS(129)*CONSTS(114)*ALGBRC(46)
      ALGBRC(163) = TERNRY(ALGBRC(96).EQ.0.00000,  CONSTS(8)*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(8)*ALGBRC(96))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(96)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(68)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(100)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(87) =  CONSTS(130)*CONSTS(131)*log( CONSTS(74)*STATES(35))+ CONSTS(106)*CONSTS(114)*ALGBRC(46)
      ALGBRC(136) =  CONSTS(130)*CONSTS(131)*log( CONSTS(75)*STATES(36))+ CONSTS(110)*CONSTS(114)*ALGBRC(46)
      ALGBRC(137) =  CONSTS(27)*(EXP(ALGBRC(87)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(136)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(167) = ((ALGBRC(121)+ALGBRC(109)) - ((((- ALGBRC(155)+ALGBRC(159)) - ALGBRC(159))+ALGBRC(163)) - ALGBRC(163))) - ALGBRC(137)
      ALGBRC(69) = ALGBRC(57)+ CONSTS(129)*CONSTS(114)*ALGBRC(46)+ALGBRC(65)
      ALGBRC(101) = ALGBRC(55)+ALGBRC(65)
      ALGBRC(97) =  CONSTS(129)*CONSTS(114)*ALGBRC(46)
      ALGBRC(165) = TERNRY(ALGBRC(97).EQ.0.00000,  CONSTS(9)*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))),  ((( CONSTS(9)*ALGBRC(97))/( CONSTS(130)*CONSTS(131)))/(EXP(ALGBRC(97)/( CONSTS(130)*CONSTS(131))) - 1.00000))*(EXP(ALGBRC(69)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(101)/( CONSTS(130)*CONSTS(131)))))
      ALGBRC(168) = - ALGBRC(163) - ALGBRC(165)
      ALGBRC(169) = ALGBRC(165)+ALGBRC(163)
      ALGBRC(172) = ((ALGBRC(129)+ALGBRC(113)) - ((((- ALGBRC(158)+ALGBRC(161)) - ALGBRC(161))+ALGBRC(165)) - ALGBRC(165)))+ALGBRC(137)
      ALGBRC(170) =  CONSTS(114)*((((((((((((((((((((((((((((((((((((( CONSTS(109)*ALGBRC(125) -  CONSTS(105)*ALGBRC(125))+ CONSTS(110)*ALGBRC(133)) -  CONSTS(106)*ALGBRC(133))+ CONSTS(108)*ALGBRC(121)) -  CONSTS(104)*ALGBRC(121))+ CONSTS(109)*ALGBRC(129)) -  CONSTS(105)*ALGBRC(129))+ CONSTS(110)*ALGBRC(137)) -  CONSTS(106)*ALGBRC(137))+ CONSTS(108)*ALGBRC(117)) -  CONSTS(104)*ALGBRC(117))+ CONSTS(107)*ALGBRC(113)) -  CONSTS(103)*ALGBRC(113))+ CONSTS(107)*ALGBRC(109)) -  CONSTS(103)*ALGBRC(109))+ CONSTS(107)*ALGBRC(105)) -  CONSTS(103)*ALGBRC(105))+ CONSTS(110)*ALGBRC(135)) -  CONSTS(106)*ALGBRC(135))+ CONSTS(109)*ALGBRC(127)) -  CONSTS(105)*ALGBRC(127))+ CONSTS(108)*ALGBRC(119)) -  CONSTS(104)*ALGBRC(119))+ CONSTS(110)*ALGBRC(131)) -  CONSTS(106)*ALGBRC(131))+ CONSTS(109)*ALGBRC(123)) -  CONSTS(105)*ALGBRC(123))+ CONSTS(108)*ALGBRC(115)) -  CONSTS(104)*ALGBRC(115))+ CONSTS(107)*ALGBRC(111)) -  CONSTS(103)*ALGBRC(111))+ CONSTS(107)*ALGBRC(107)) -  CONSTS(103)*ALGBRC(107))+ CONSTS(107)*ALGBRC(103)) -  CONSTS(103)*ALGBRC(103)) - ( CONSTS(128)*ALGBRC(159)+ CONSTS(128)*ALGBRC(161))) - ( CONSTS(129)*ALGBRC(163)+ CONSTS(129)*ALGBRC(165)))
      ALGBRC(10) = TERNRY(VOI.GT.CONSTS(124).AND.VOI.LT.CONSTS(125), CONSTS(126), CONSTS(127))
      ALGBRC(173) = - ALGBRC(170)+ALGBRC(10)
      ALGBRC(171) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(174) =  CONSTS(130)*CONSTS(131)*log( CONSTS(85)*STATES(37))
      ALGBRC(177) =  CONSTS(130)*CONSTS(131)*log( CONSTS(88)*STATES(40))
      ALGBRC(179) =  CONSTS(40)*(EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(174)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(175) =  CONSTS(130)*CONSTS(131)*log( CONSTS(86)*STATES(38))
      ALGBRC(180) =  CONSTS(41)*(EXP((ALGBRC(174)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(175)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(176) =  CONSTS(130)*CONSTS(131)*log( CONSTS(87)*STATES(39))
      ALGBRC(181) =  CONSTS(42)*(EXP((ALGBRC(175)+ CONSTS(153)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(182) =  CONSTS(43)*(EXP(ALGBRC(176)/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(177)+ CONSTS(152)*ALGBRC(171))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(189) =  CONSTS(130)*CONSTS(131)*log( CONSTS(101)*STATES(22))
      ALGBRC(190) = CONSTS(166)+ALGBRC(189)
      ALGBRC(199) =  CONSTS(52)*(EXP(CONSTS(177)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(190)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(200) = ALGBRC(199)
      ALGBRC(188) =  CONSTS(130)*CONSTS(131)*log( CONSTS(99)*STATES(24))
      ALGBRC(191) = CONSTS(166)+ALGBRC(188)
      ALGBRC(201) =  CONSTS(44)*(EXP(ALGBRC(191)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(178)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(202) = - ALGBRC(201)
      ALGBRC(187) =  CONSTS(130)*CONSTS(131)*log( CONSTS(100)*STATES(23))
      ALGBRC(192) = ALGBRC(187)+CONSTS(174)
      ALGBRC(203) =  CONSTS(48)*(EXP(CONSTS(179)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(192)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(205) = ALGBRC(203)
      ALGBRC(166) =  CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
      ALGBRC(178) =  CONSTS(39)*EXP(ALGBRC(177)/( CONSTS(130)*CONSTS(131)))*(EXP(ALGBRC(166)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(171)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(185) =  CONSTS(130)*CONSTS(131)*log( CONSTS(84)*STATES(20))
      ALGBRC(196) = CONSTS(175)+ CONSTS(155)*ALGBRC(185)
      ALGBRC(212) =  CONSTS(49)*(EXP(CONSTS(183)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(196)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(220) =  CONSTS(155)*ALGBRC(212)
      ALGBRC(219) = ALGBRC(214)/CONSTS(146)
      ALGBRC(186) =  CONSTS(130)*CONSTS(131)*log( CONSTS(98)*STATES(21))
      ALGBRC(193) = CONSTS(169)+ALGBRC(186)
      ALGBRC(206) =  CONSTS(46)*(EXP(ALGBRC(193)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(180)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(195) = CONSTS(172)+ CONSTS(156)*ALGBRC(186)
      ALGBRC(209) =  CONSTS(47)*(EXP(CONSTS(182)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(195)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(197) = CONSTS(175)+ CONSTS(156)*ALGBRC(186)
      ALGBRC(221) =  CONSTS(50)*(EXP(ALGBRC(197)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(184)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(198) = ALGBRC(186)+CONSTS(173)
      ALGBRC(224) =  CONSTS(51)*(EXP(CONSTS(185)/( CONSTS(130)*CONSTS(131))) - EXP(ALGBRC(198)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(227) = (( CONSTS(156)*ALGBRC(209) - ALGBRC(206)) -  CONSTS(156)*ALGBRC(221))+ALGBRC(224)
      ALGBRC(25) =  CONSTS(130)*CONSTS(131)*log( CONSTS(53)*STATES(2))
      ALGBRC(29) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(34) =  (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(8)*(STATES(43)+CONSTS(137))
      ALGBRC(37) =  (CONSTS(147)/( CONSTS(137)*CONSTS(138)*CONSTS(139)))*STATES(7)*STATES(42)
      ALGBRC(40) = ALGBRC(34)+ALGBRC(37)
      ALGBRC(43) = ALGBRC(40)
      ALGBRC(223) =  (( CONSTS(151)*ALGBRC(213))/CONSTS(137))*ALGBRC(219)
      ALGBRC(226) = ALGBRC(223)
      ALGBRC(230) = ALGBRC(43)+ALGBRC(226)
      ALGBRC(231) =  ALGBRC(230)*1000.00
      ALGBRC(232) =  CONSTS(136)*CONSTS(161)*ALGBRC(231)
      ALGBRC(233) =  CONSTS(1)*(EXP((ALGBRC(29)+ALGBRC(25))/( CONSTS(130)*CONSTS(131))) - EXP((ALGBRC(31)+ALGBRC(232))/( CONSTS(130)*CONSTS(131))))
      ALGBRC(164) = (( (CONSTS(102))*ALGBRC(151)+ CONSTS(102)*ALGBRC(155)+ CONSTS(102)*ALGBRC(147)+ CONSTS(102)*ALGBRC(143)+ CONSTS(102)*ALGBRC(139)+ CONSTS(102)*ALGBRC(158)) - ALGBRC(159)) - ALGBRC(161)
      ALGBRC(183) = (( CONSTS(153)*ALGBRC(179) -  CONSTS(152)*ALGBRC(180)) -  CONSTS(153)*ALGBRC(181))+ CONSTS(152)*ALGBRC(182)
      ALGBRC(184) =  CONSTS(130)*CONSTS(131)*log( CONSTS(60)*STATES(4))
      ALGBRC(194) = CONSTS(169)+ CONSTS(154)*ALGBRC(184)
      ALGBRC(207) =  CONSTS(45)*(EXP(ALGBRC(194)/( CONSTS(130)*CONSTS(131))) - EXP(CONSTS(181)/( CONSTS(130)*CONSTS(131))))
      ALGBRC(210) =  - CONSTS(154)*ALGBRC(207)
      ALGBRC(1) = STATES(1)/CONSTS(113)
      ALGBRC(3) = STATES(12)+STATES(14)+STATES(17)+STATES(16)
      ALGBRC(4) = STATES(13)+STATES(14)
      ALGBRC(5) = STATES(18)+STATES(19)
      ALGBRC(6) = STATES(10)
      ALGBRC(7) = STATES(11)
      ALGBRC(11) =  CONSTS(130)*CONSTS(131)*log( CONSTS(70)*STATES(31))
      ALGBRC(12) =  CONSTS(130)*CONSTS(131)*log( CONSTS(71)*STATES(32))
      ALGBRC(13) =  CONSTS(130)*CONSTS(131)*log( CONSTS(72)*STATES(33))
      ALGBRC(14) =  CONSTS(130)*CONSTS(131)*log( CONSTS(73)*STATES(34))
      ALGBRC(15) =  CONSTS(153)*STATES(40)+ CONSTS(152)*STATES(38)+ (CONSTS(152)+CONSTS(153))*STATES(39)
      ALGBRC(16) = STATES(4)/CONSTS(157)
      ALGBRC(17) = STATES(20)/CONSTS(163)
      ALGBRC(18) = STATES(21)/CONSTS(167)
      ALGBRC(19) = STATES(23)/CONSTS(157)
      ALGBRC(20) = STATES(24)/CONSTS(157)
      ALGBRC(21) = STATES(22)/CONSTS(157)
      ALGBRC(204) = ALGBRC(199) - ALGBRC(201)
      ALGBRC(208) = ALGBRC(206)
      ALGBRC(211) = ALGBRC(201) - ALGBRC(207)
      ALGBRC(217) = ALGBRC(207) - ALGBRC(209)
      ALGBRC(218) = ALGBRC(209) - ALGBRC(203)
      ALGBRC(222) = ALGBRC(203) - ALGBRC(212)
      ALGBRC(225) = ALGBRC(212) - ALGBRC(221)
      ALGBRC(228) = ALGBRC(221) - ALGBRC(224)
      ALGBRC(229) = ALGBRC(224) - ALGBRC(199)
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(213) - (TERNRY(ALGBRC(214).GT.0.00000, CONSTS(148), CONSTS(149))))
      minfunc_0[2] = abs(ALGBRC(215) -  (CONSTS(150)/ALGBRC(213))*STATES(44)*CONSTS(146))
      minfunc_0[3] = abs(ALGBRC(216) -  (CONSTS(150)/ALGBRC(213))*(STATES(41) - CONSTS(137))*CONSTS(146))
      minfunc_0[4] = abs(ALGBRC(214) - (- ALGBRC(215)+ALGBRC(216)))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END