Generated Code

/*
   There are a total of 303 entries in the algebraic variable array.
   There are a total of 62 entries in each of the rate and state variable arrays.
   There are a total of 177 entries in the constant variable array.
 */
/*
 * CONSTANTS[0] is kappa_leak in component BG_parameters (fmol_per_sec).
 * CONSTANTS[1] is kappa_R_TRPNCa in component BG_parameters (fmol_per_sec).
 * CONSTANTS[2] is kappa_R_BU in component BG_parameters (fmol_per_sec).
 * CONSTANTS[3] is kappa_R_UW in component BG_parameters (fmol_per_sec).
 * CONSTANTS[4] is kappa_R_WS in component BG_parameters (fmol_per_sec).
 * CONSTANTS[5] is kappa_R_SU in component BG_parameters (fmol_per_sec).
 * CONSTANTS[6] is kappa_diadDf in component BG_parameters (fmol_per_sec).
 * CONSTANTS[7] is kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[8] is kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[9] is kappa_d000 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[10] is kappa_d010 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[11] is kappa_d020 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[12] is kappa_d001 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[13] is kappa_d011 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[14] is kappa_d021 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[15] is kappa_f1_000 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[16] is kappa_f1_100 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[17] is kappa_f1_001 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[18] is kappa_f1_101 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[19] is kappa_f2_000 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[20] is kappa_f2_100 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[21] is kappa_f2_001 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[22] is kappa_f2_101 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[23] is kappa_f3_010 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[24] is kappa_f3_110 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[25] is kappa_f3_011 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[26] is kappa_f3_111 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[27] is kappa_fCa000 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[28] is kappa_fCa100 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[29] is kappa_fCa010 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[30] is kappa_fCa110 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[31] is kappa_fCa020 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[32] is kappa_fCa120 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[33] is kappa_1_NCX in component BG_parameters (fmol_per_sec).
 * CONSTANTS[34] is kappa_2_NCX in component BG_parameters (fmol_per_sec).
 * CONSTANTS[35] is kappa_3_NCX in component BG_parameters (fmol_per_sec).
 * CONSTANTS[36] is kappa_4_NCX in component BG_parameters (fmol_per_sec).
 * CONSTANTS[37] is kappa_5_NCX in component BG_parameters (fmol_per_sec).
 * CONSTANTS[38] is kappa_6_NCX in component BG_parameters (fmol_per_sec).
 * CONSTANTS[39] is kappa_pCa_R1 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[40] is kappa_pCa_R2 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[41] is kappa_PLBph1 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[42] is kappa_PLBph2 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[43] is kappa_PLBd1 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[44] is kappa_PLBd2 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[45] is kappa_Inh in component BG_parameters (fmol_per_sec).
 * CONSTANTS[46] is kappa_RyR in component BG_parameters (fmol_per_sec).
 * CONSTANTS[47] is kappa_OC in component BG_parameters (fmol_per_sec).
 * CONSTANTS[48] is kappa_CCI in component BG_parameters (fmol_per_sec).
 * CONSTANTS[49] is kappa_CII in component BG_parameters (fmol_per_sec).
 * CONSTANTS[50] is kappa_IO in component BG_parameters (fmol_per_sec).
 * CONSTANTS[51] is kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[52] is kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[53] is kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec).
 * CONSTANTS[54] is kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[55] is kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[56] is kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[57] is kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[58] is kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[59] is kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec).
 * CONSTANTS[60] is K_Ca_SR in component BG_parameters (per_fmol).
 * CONSTANTS[61] is K_Ca_i in component BG_parameters (per_fmol).
 * CONSTANTS[62] is K_TRPN in component BG_parameters (per_fmol).
 * CONSTANTS[63] is K_Ca_TRPN in component BG_parameters (per_fmol).
 * CONSTANTS[64] is K_B_XB in component BG_parameters (per_fmol).
 * CONSTANTS[65] is K_U_XB in component BG_parameters (per_fmol).
 * CONSTANTS[66] is K_W_XB in component BG_parameters (per_fmol).
 * CONSTANTS[67] is K_S_XB in component BG_parameters (per_fmol).
 * CONSTANTS[68] is K_Ca_D in component BG_parameters (per_fmol).
 * CONSTANTS[69] is K_Ca_o in component BG_parameters (per_fmol).
 * CONSTANTS[70] is K_000_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[71] is K_010_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[72] is K_020_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[73] is K_100_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[74] is K_110_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[75] is K_120_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[76] is K_001_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[77] is K_011_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[78] is K_021_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[79] is K_101_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[80] is K_111_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[81] is K_121_LCC in component BG_parameters (per_fmol).
 * CONSTANTS[82] is K_1_NCX in component BG_parameters (per_fmol).
 * CONSTANTS[83] is K_2_NCX in component BG_parameters (per_fmol).
 * CONSTANTS[84] is K_3_NCX in component BG_parameters (per_fmol).
 * CONSTANTS[85] is K_4_NCX in component BG_parameters (per_fmol).
 * CONSTANTS[86] is K_5_NCX in component BG_parameters (per_fmol).
 * CONSTANTS[87] is K_6_NCX in component BG_parameters (per_fmol).
 * CONSTANTS[88] is K_Na_i in component BG_parameters (per_fmol).
 * CONSTANTS[89] is K_Na_o in component BG_parameters (per_fmol).
 * CONSTANTS[90] is K_pCa in component BG_parameters (per_fmol).
 * CONSTANTS[91] is K_pCa_Ca in component BG_parameters (per_fmol).
 * CONSTANTS[92] is K_PLB in component BG_parameters (per_fmol).
 * CONSTANTS[93] is K_PKACI in component BG_parameters (per_fmol).
 * CONSTANTS[94] is K_PLB_PKACI in component BG_parameters (per_fmol).
 * CONSTANTS[95] is K_PP1 in component BG_parameters (per_fmol).
 * CONSTANTS[96] is K_PLBp_PP1 in component BG_parameters (per_fmol).
 * CONSTANTS[97] is K_PLBp in component BG_parameters (per_fmol).
 * CONSTANTS[98] is K_Ip in component BG_parameters (per_fmol).
 * CONSTANTS[99] is K_Ip_PP1 in component BG_parameters (per_fmol).
 * CONSTANTS[100] is K_C_RyR in component BG_parameters (per_fmol).
 * CONSTANTS[101] is K_CI_RyR in component BG_parameters (per_fmol).
 * CONSTANTS[102] is K_I_RyR in component BG_parameters (per_fmol).
 * CONSTANTS[103] is K_O_RyR in component BG_parameters (per_fmol).
 * CONSTANTS[104] is K_P1_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[105] is K_P2_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[106] is K_P2a_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[107] is K_P4_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[108] is K_P5_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[109] is K_P6_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[110] is K_P8_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[111] is K_P9_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[112] is K_P10_SERCA in component BG_parameters (per_fmol).
 * CONSTANTS[113] is K_H in component BG_parameters (per_fmol).
 * CONSTANTS[114] is K_MgATP in component BG_parameters (per_fmol).
 * CONSTANTS[115] is K_MgADP in component BG_parameters (per_fmol).
 * CONSTANTS[116] is K_Pi in component BG_parameters (per_fmol).
 * CONSTANTS[117] is z_rCa in component BG_parameters (dimensionless).
 * CONSTANTS[118] is z_fd in component BG_parameters (dimensionless).
 * CONSTANTS[119] is z_ff1 in component BG_parameters (dimensionless).
 * CONSTANTS[120] is z_ff2 in component BG_parameters (dimensionless).
 * CONSTANTS[121] is z_ff3 in component BG_parameters (dimensionless).
 * CONSTANTS[122] is z_rd in component BG_parameters (dimensionless).
 * CONSTANTS[123] is z_rf1 in component BG_parameters (dimensionless).
 * CONSTANTS[124] is z_rf2 in component BG_parameters (dimensionless).
 * CONSTANTS[125] is z_rf3 in component BG_parameters (dimensionless).
 * CONSTANTS[126] is nNa_o_NCX in component BG_parameters (dimensionless).
 * CONSTANTS[127] is nNa_i_NCX in component BG_parameters (dimensionless).
 * CONSTANTS[128] is zf_NCX in component BG_parameters (dimensionless).
 * CONSTANTS[129] is zr_NCX in component BG_parameters (dimensionless).
 * CONSTANTS[130] is n_Tm in component BG_parameters (dimensionless).
 * CONSTANTS[131] is n_CaTRPN in component BG_parameters (dimensionless).
 * VOI is time in component environment (second).
 * CONSTANTS[132] is vol_myo in component environment (pL).
 * CONSTANTS[133] is C_m in component environment (fF).
 * STATES[0] is q_membrane in component environment (fC).
 * STATES[1] is q_Ca_D in component environment (fmol).
 * STATES[2] is q_TRPN in component environment (fmol).
 * STATES[3] is q_Ca_TRPN in component environment (fmol).
 * STATES[4] is q_Ca_i in component environment (fmol).
 * STATES[5] is q_Na_i in component environment (fmol).
 * STATES[6] is q_Na_o in component environment (fmol).
 * STATES[7] is q_B_XB in component environment (fmol).
 * STATES[8] is q_U_XB in component environment (fmol).
 * STATES[9] is q_W_XB in component environment (fmol).
 * STATES[10] is q_S_XB in component environment (fmol).
 * STATES[11] is q_Ca_o in component environment (fmol).
 * STATES[12] is q_PLB in component environment (fmol).
 * STATES[13] is q_PKACI in component environment (fmol).
 * STATES[14] is q_PLB_PKACI in component environment (fmol).
 * STATES[15] is q_PP1 in component environment (fmol).
 * STATES[16] is q_PLBp_PP1 in component environment (fmol).
 * STATES[17] is q_PLBp in component environment (fmol).
 * CONSTANTS[134] is q_PLBp_baseline in component environment (fmol).
 * STATES[18] is q_Ip in component environment (fmol).
 * STATES[19] is q_Ip_PP1 in component environment (fmol).
 * STATES[20] is q_Ca_SR in component environment (fmol).
 * STATES[21] is q_pCa in component environment (fmol).
 * STATES[22] is q_pCa_Ca in component environment (fmol).
 * STATES[23] is q_H in component environment (fmol).
 * STATES[24] is q_Pi in component environment (fmol).
 * STATES[25] is q_MgADP in component environment (fmol).
 * STATES[26] is q_MgATP in component environment (fmol).
 * STATES[27] is q_000_LCC in component environment (fmol).
 * STATES[28] is q_010_LCC in component environment (fmol).
 * STATES[29] is q_020_LCC in component environment (fmol).
 * STATES[30] is q_100_LCC in component environment (fmol).
 * STATES[31] is q_110_LCC in component environment (fmol).
 * STATES[32] is q_120_LCC in component environment (fmol).
 * STATES[33] is q_001_LCC in component environment (fmol).
 * STATES[34] is q_011_LCC in component environment (fmol).
 * STATES[35] is q_021_LCC in component environment (fmol).
 * STATES[36] is q_101_LCC in component environment (fmol).
 * STATES[37] is q_111_LCC in component environment (fmol).
 * STATES[38] is q_121_LCC in component environment (fmol).
 * STATES[39] is q_P1_NCX in component environment (fmol).
 * STATES[40] is q_P2_NCX in component environment (fmol).
 * STATES[41] is q_P3_NCX in component environment (fmol).
 * STATES[42] is q_P4_NCX in component environment (fmol).
 * STATES[43] is q_P5_NCX in component environment (fmol).
 * STATES[44] is q_P6_NCX in component environment (fmol).
 * STATES[45] is q_C_RyR in component environment (fmol).
 * STATES[46] is q_CI_RyR in component environment (fmol).
 * STATES[47] is q_I_RyR in component environment (fmol).
 * STATES[48] is q_O_RyR in component environment (fmol).
 * STATES[49] is q_P1_SERCA in component environment (fmol).
 * STATES[50] is q_P2_SERCA in component environment (fmol).
 * STATES[51] is q_P2a_SERCA in component environment (fmol).
 * STATES[52] is q_P4_SERCA in component environment (fmol).
 * STATES[53] is q_P5_SERCA in component environment (fmol).
 * STATES[54] is q_P6_SERCA in component environment (fmol).
 * STATES[55] is q_P8_SERCA in component environment (fmol).
 * STATES[56] is q_P9_SERCA in component environment (fmol).
 * STATES[57] is q_P10_SERCA in component environment (fmol).
 * ALGEBRAIC[54] is v_leak in component Ca_leak (fmol_per_sec).
 * ALGEBRAIC[66] is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec).
 * ALGEBRAIC[73] is v_R_BU in component crossbridge_TRPN (fmol_per_sec).
 * ALGEBRAIC[75] is v_R_UW in component crossbridge_TRPN (fmol_per_sec).
 * ALGEBRAIC[77] is v_R_WS in component crossbridge_TRPN (fmol_per_sec).
 * ALGEBRAIC[79] is v_R_SU in component crossbridge_TRPN (fmol_per_sec).
 * ALGEBRAIC[59] is T_total in component crossbridge_TRPN (kPa).
 * STATES[58] is SL in component crossbridge_TRPN (metre).
 * ALGEBRAIC[72] is v_diadDf in component diadCaDiffusion (fmol_per_sec).
 * ALGEBRAIC[182] is v_Ca_D_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[181] is v_Ca_o_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[180] is I_mem_LCC in component LCC (fA).
 * ALGEBRAIC[221] is v_Ca_o_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[222] is v_Ca_i_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[214] is v_Na_o_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[203] is v_Na_i_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[220] is I_mem_NCX in component NCX (fA).
 * ALGEBRAIC[227] is v_pCa_R1 in component pCa (fmol_per_sec).
 * ALGEBRAIC[228] is v_pCa_R2 in component pCa (fmol_per_sec).
 * ALGEBRAIC[230] is I_mem_pCa in component pCa (fA).
 * ALGEBRAIC[55] is v_PLBph1 in component PLB (fmol_per_sec).
 * ALGEBRAIC[60] is v_PLBph2 in component PLB (fmol_per_sec).
 * ALGEBRAIC[62] is v_PLBd1 in component PLB (fmol_per_sec).
 * ALGEBRAIC[65] is v_PLBd2 in component PLB (fmol_per_sec).
 * ALGEBRAIC[56] is v_Inh in component PLB (fmol_per_sec).
 * ALGEBRAIC[237] is v_RyR in component RyR (fmol_per_sec).
 * ALGEBRAIC[246] is v_RyRgate_Ca_D in component RyR (fmol_per_sec).
 * ALGEBRAIC[289] is v_Ca_i_SERCA in component SERCA (fmol_per_sec).
 * ALGEBRAIC[295] is v_Ca_SR_SERCA in component SERCA (fmol_per_sec).
 * ALGEBRAIC[300] is v_H_SERCA in component SERCA (fmol_per_sec).
 * ALGEBRAIC[286] is v_MgADP_SERCA in component SERCA (fmol_per_sec).
 * ALGEBRAIC[283] is v_MgATP_SERCA in component SERCA (fmol_per_sec).
 * ALGEBRAIC[281] is v_Pi_SERCA in component SERCA (fmol_per_sec).
 * ALGEBRAIC[1] is V_m in component environment (volt).
 * ALGEBRAIC[28] is I_stim in component environment (fA).
 * CONSTANTS[135] is stimPeriod in component environment (second).
 * CONSTANTS[136] is stimDuration in component environment (second).
 * ALGEBRAIC[11] is tPeriod in component environment (second).
 * ALGEBRAIC[33] is v_Ca_stim in component environment (fmol_per_sec).
 * ALGEBRAIC[232] is sum_I in component environment (fA).
 * ALGEBRAIC[2] is dSL in component environment (m_per_s).
 * CONSTANTS[137] is V_max in component environment (m_per_s).
 * CONSTANTS[138] is vol_cell in component environment (pL).
 * CONSTANTS[170] is vol_SR in component environment (pL).
 * CONSTANTS[139] is vol_o in component environment (pL).
 * CONSTANTS[140] is vol_D in component environment (pL).
 * ALGEBRAIC[3] is cCai in component environment (mM).
 * ALGEBRAIC[4] is cCao in component environment (mM).
 * ALGEBRAIC[5] is cCaSR in component environment (mM).
 * ALGEBRAIC[6] is cCaD in component environment (mM).
 * ALGEBRAIC[7] is cNai in component environment (mM).
 * ALGEBRAIC[8] is pHi in component environment (dimensionless).
 * ALGEBRAIC[9] is Ca_tot in component environment (fmol).
 * CONSTANTS[141] is R in component constants (J_per_K_per_mol).
 * CONSTANTS[142] is T in component constants (kelvin).
 * CONSTANTS[143] is F in component constants (C_per_mol).
 * CONSTANTS[144] is zNa in component ion_valences (dimensionless).
 * CONSTANTS[145] is zCl in component ion_valences (dimensionless).
 * CONSTANTS[146] is zK in component ion_valences (dimensionless).
 * CONSTANTS[147] is zCa in component ion_valences (dimensionless).
 * ALGEBRAIC[38] is mu_Ca_SR in component Ca_leak (J_per_mol).
 * ALGEBRAIC[48] is mu_Ca_i in component Ca_leak (J_per_mol).
 * ALGEBRAIC[58] is mu_TRPN in component crossbridge_TRPN (J_per_mol).
 * ALGEBRAIC[61] is mu_Ca_i in component crossbridge_TRPN (J_per_mol).
 * ALGEBRAIC[63] is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol).
 * ALGEBRAIC[64] is mu_B_XB in component crossbridge_TRPN (J_per_mol).
 * ALGEBRAIC[67] is mu_U_XB in component crossbridge_TRPN (J_per_mol).
 * ALGEBRAIC[69] is mu_W_XB in component crossbridge_TRPN (J_per_mol).
 * ALGEBRAIC[71] is mu_S_XB in component crossbridge_TRPN (J_per_mol).
 * CONSTANTS[148] is SL_0 in component crossbridge_TRPN (metre).
 * CONSTANTS[149] is q_MS in component crossbridge_TRPN (fmol).
 * CONSTANTS[150] is r_s in component crossbridge_TRPN (dimensionless).
 * CONSTANTS[151] is r_w in component crossbridge_TRPN (dimensionless).
 * CONSTANTS[152] is A_eff in component crossbridge_TRPN (dimensionless).
 * CONSTANTS[153] is phi in component crossbridge_TRPN (dimensionless).
 * CONSTANTS[154] is k_uw in component crossbridge_TRPN (per_sec).
 * CONSTANTS[155] is k_ws in component crossbridge_TRPN (per_sec).
 * STATES[59] is G_w in component crossbridge_TRPN (metre).
 * STATES[60] is G_s in component crossbridge_TRPN (metre).
 * CONSTANTS[171] is c_w in component crossbridge_TRPN (per_sec).
 * CONSTANTS[172] is c_s in component crossbridge_TRPN (per_sec).
 * CONSTANTS[169] is A_w in component crossbridge_TRPN (dimensionless).
 * CONSTANTS[175] is A_s in component crossbridge_TRPN (dimensionless).
 * ALGEBRAIC[41] is nu_1 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[30] is nu_2 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[35] is nu_3 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[42] is nu_4 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[31] is nu_5 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[36] is nu_6 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[51] is v_1 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[26] is v_2 in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[52] is v_4 in component crossbridge_TRPN (m_per_s).
 * CONSTANTS[156] is v_to_mu in component crossbridge_TRPN (Js_per_m2).
 * ALGEBRAIC[32] is T_active in component crossbridge_TRPN (kPa).
 * CONSTANTS[157] is T_ref in component crossbridge_TRPN (kPa).
 * ALGEBRAIC[10] is F_T_S in component crossbridge_TRPN (kPa).
 * ALGEBRAIC[27] is F_T_W in component crossbridge_TRPN (kPa).
 * STATES[61] is Cdd in component crossbridge_TRPN (metre).
 * ALGEBRAIC[43] is eta in component crossbridge_TRPN (second).
 * CONSTANTS[158] is eta_l in component crossbridge_TRPN (second).
 * CONSTANTS[159] is eta_s in component crossbridge_TRPN (second).
 * CONSTANTS[160] is k in component crossbridge_TRPN (dimensionless).
 * CONSTANTS[161] is alpha in component crossbridge_TRPN (kPa).
 * ALGEBRAIC[44] is v_Cdd in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[45] is nu_d in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[46] is nu_k in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[47] is nu_f in component crossbridge_TRPN (m_per_s).
 * ALGEBRAIC[53] is F_d in component crossbridge_TRPN (kPa).
 * ALGEBRAIC[37] is F_1 in component crossbridge_TRPN (kPa).
 * ALGEBRAIC[57] is T_passive in component crossbridge_TRPN (kPa).
 * CONSTANTS[162] is b in component crossbridge_TRPN (dimensionless).
 * ALGEBRAIC[68] is mu_Ca_D in component diadCaDiffusion (J_per_mol).
 * ALGEBRAIC[70] is mu_Ca_i in component diadCaDiffusion (J_per_mol).
 * ALGEBRAIC[76] is mu_Ca_o in component LCC (J_per_mol).
 * ALGEBRAIC[78] is mu_Ca_D in component LCC (J_per_mol).
 * ALGEBRAIC[74] is V_mem in component LCC (volt).
 * ALGEBRAIC[88] is Af_LCC_Ca1 in component LCC (J_per_mol).
 * ALGEBRAIC[116] is Ar_LCC_Ca1 in component LCC (J_per_mol).
 * ALGEBRAIC[114] is Am_LCC_Ca1 in component LCC (J_per_mol).
 * ALGEBRAIC[175] is v_LCC_Ca1 in component LCC (fmol_per_sec).
 * ALGEBRAIC[89] is Af_LCC_Ca2 in component LCC (J_per_mol).
 * ALGEBRAIC[117] is Ar_LCC_Ca2 in component LCC (J_per_mol).
 * ALGEBRAIC[115] is Am_LCC_Ca2 in component LCC (J_per_mol).
 * ALGEBRAIC[177] is v_LCC_Ca2 in component LCC (fmol_per_sec).
 * ALGEBRAIC[108] is Af_fCa00 in component LCC (J_per_mol).
 * ALGEBRAIC[154] is Ar_fCa00 in component LCC (J_per_mol).
 * ALGEBRAIC[155] is v_fCa00 in component LCC (fmol_per_sec).
 * ALGEBRAIC[109] is Af_fCa01 in component LCC (J_per_mol).
 * ALGEBRAIC[156] is Ar_fCa01 in component LCC (J_per_mol).
 * ALGEBRAIC[159] is v_fCa01 in component LCC (fmol_per_sec).
 * ALGEBRAIC[110] is Af_fCa02 in component LCC (J_per_mol).
 * ALGEBRAIC[160] is Ar_fCa02 in component LCC (J_per_mol).
 * ALGEBRAIC[163] is v_fCa02 in component LCC (fmol_per_sec).
 * ALGEBRAIC[111] is Af_fCa10 in component LCC (J_per_mol).
 * ALGEBRAIC[164] is Ar_fCa10 in component LCC (J_per_mol).
 * ALGEBRAIC[167] is v_fCa10 in component LCC (fmol_per_sec).
 * ALGEBRAIC[112] is Af_fCa11 in component LCC (J_per_mol).
 * ALGEBRAIC[168] is Ar_fCa11 in component LCC (J_per_mol).
 * ALGEBRAIC[171] is v_fCa11 in component LCC (fmol_per_sec).
 * ALGEBRAIC[113] is Af_fCa12 in component LCC (J_per_mol).
 * ALGEBRAIC[172] is Ar_fCa12 in component LCC (J_per_mol).
 * ALGEBRAIC[174] is v_fCa12 in component LCC (fmol_per_sec).
 * ALGEBRAIC[80] is mu_S000_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[157] is v_S000_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[81] is mu_S010_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[161] is v_S010_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[82] is mu_S020_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[165] is v_S020_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[83] is mu_S100_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[169] is v_S100_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[84] is mu_S110_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[173] is v_S110_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[86] is mu_S120_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[176] is v_S120_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[12] is mu_S001_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[158] is v_S001_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[13] is mu_S011_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[162] is v_S011_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[14] is mu_S021_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[166] is v_S021_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[15] is mu_S101_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[170] is v_S101_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[85] is mu_S111_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[178] is v_S111_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[87] is mu_S121_LCC in component LCC (J_per_mol).
 * ALGEBRAIC[179] is v_S121_LCC in component LCC (fmol_per_sec).
 * ALGEBRAIC[90] is Af_d000 in component LCC (J_per_mol).
 * ALGEBRAIC[118] is Ar_d000 in component LCC (J_per_mol).
 * ALGEBRAIC[119] is v_d000 in component LCC (fmol_per_sec).
 * ALGEBRAIC[92] is Af_d010 in component LCC (J_per_mol).
 * ALGEBRAIC[122] is Ar_d010 in component LCC (J_per_mol).
 * ALGEBRAIC[123] is v_d010 in component LCC (fmol_per_sec).
 * ALGEBRAIC[94] is Af_d020 in component LCC (J_per_mol).
 * ALGEBRAIC[126] is Ar_d020 in component LCC (J_per_mol).
 * ALGEBRAIC[127] is v_d020 in component LCC (fmol_per_sec).
 * ALGEBRAIC[96] is Af_f1_000 in component LCC (J_per_mol).
 * ALGEBRAIC[130] is Ar_f1_000 in component LCC (J_per_mol).
 * ALGEBRAIC[131] is v_f1_000 in component LCC (fmol_per_sec).
 * ALGEBRAIC[100] is Af_f2_000 in component LCC (J_per_mol).
 * ALGEBRAIC[138] is Ar_f2_000 in component LCC (J_per_mol).
 * ALGEBRAIC[139] is v_f2_000 in component LCC (fmol_per_sec).
 * ALGEBRAIC[104] is Af_f3_010 in component LCC (J_per_mol).
 * ALGEBRAIC[146] is Ar_f3_010 in component LCC (J_per_mol).
 * ALGEBRAIC[147] is v_f3_010 in component LCC (fmol_per_sec).
 * ALGEBRAIC[98] is Af_f1_100 in component LCC (J_per_mol).
 * ALGEBRAIC[134] is Ar_f1_100 in component LCC (J_per_mol).
 * ALGEBRAIC[135] is v_f1_100 in component LCC (fmol_per_sec).
 * ALGEBRAIC[102] is Af_f2_100 in component LCC (J_per_mol).
 * ALGEBRAIC[142] is Ar_f2_100 in component LCC (J_per_mol).
 * ALGEBRAIC[143] is v_f2_100 in component LCC (fmol_per_sec).
 * ALGEBRAIC[106] is Af_f3_110 in component LCC (J_per_mol).
 * ALGEBRAIC[150] is Ar_f3_110 in component LCC (J_per_mol).
 * ALGEBRAIC[151] is v_f3_110 in component LCC (fmol_per_sec).
 * ALGEBRAIC[91] is Af_d001 in component LCC (J_per_mol).
 * ALGEBRAIC[120] is Ar_d001 in component LCC (J_per_mol).
 * ALGEBRAIC[121] is v_d001 in component LCC (fmol_per_sec).
 * ALGEBRAIC[93] is Af_d011 in component LCC (J_per_mol).
 * ALGEBRAIC[124] is Ar_d011 in component LCC (J_per_mol).
 * ALGEBRAIC[125] is v_d011 in component LCC (fmol_per_sec).
 * ALGEBRAIC[95] is Af_d021 in component LCC (J_per_mol).
 * ALGEBRAIC[128] is Ar_d021 in component LCC (J_per_mol).
 * ALGEBRAIC[129] is v_d021 in component LCC (fmol_per_sec).
 * ALGEBRAIC[97] is Af_f1_001 in component LCC (J_per_mol).
 * ALGEBRAIC[132] is Ar_f1_001 in component LCC (J_per_mol).
 * ALGEBRAIC[133] is v_f1_001 in component LCC (fmol_per_sec).
 * ALGEBRAIC[101] is Af_f2_001 in component LCC (J_per_mol).
 * ALGEBRAIC[140] is Ar_f2_001 in component LCC (J_per_mol).
 * ALGEBRAIC[141] is v_f2_001 in component LCC (fmol_per_sec).
 * ALGEBRAIC[105] is Af_f3_011 in component LCC (J_per_mol).
 * ALGEBRAIC[148] is Ar_f3_011 in component LCC (J_per_mol).
 * ALGEBRAIC[149] is v_f3_011 in component LCC (fmol_per_sec).
 * ALGEBRAIC[99] is Af_f1_101 in component LCC (J_per_mol).
 * ALGEBRAIC[136] is Ar_f1_101 in component LCC (J_per_mol).
 * ALGEBRAIC[137] is v_f1_101 in component LCC (fmol_per_sec).
 * ALGEBRAIC[103] is Af_f2_101 in component LCC (J_per_mol).
 * ALGEBRAIC[144] is Ar_f2_101 in component LCC (J_per_mol).
 * ALGEBRAIC[145] is v_f2_101 in component LCC (fmol_per_sec).
 * ALGEBRAIC[107] is Af_f3_111 in component LCC (J_per_mol).
 * ALGEBRAIC[152] is Ar_f3_111 in component LCC (J_per_mol).
 * ALGEBRAIC[153] is v_f3_111 in component LCC (fmol_per_sec).
 * ALGEBRAIC[16] is Ca_tot in component LCC (fmol).
 * ALGEBRAIC[184] is mu_Ca_o in component NCX (J_per_mol).
 * ALGEBRAIC[185] is mu_Ca_i in component NCX (J_per_mol).
 * ALGEBRAIC[186] is mu_Na_o in component NCX (J_per_mol).
 * ALGEBRAIC[187] is mu_Na_i in component NCX (J_per_mol).
 * ALGEBRAIC[188] is mu_P1_NCX in component NCX (J_per_mol).
 * ALGEBRAIC[189] is mu_P2_NCX in component NCX (J_per_mol).
 * ALGEBRAIC[190] is mu_P3_NCX in component NCX (J_per_mol).
 * ALGEBRAIC[191] is mu_P4_NCX in component NCX (J_per_mol).
 * ALGEBRAIC[192] is mu_P5_NCX in component NCX (J_per_mol).
 * ALGEBRAIC[193] is mu_P6_NCX in component NCX (J_per_mol).
 * ALGEBRAIC[218] is v_P1_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[206] is v_P2_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[209] is v_P3_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[212] is v_P4_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[215] is v_P5_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[219] is v_P6_NCX in component NCX (fmol_per_sec).
 * ALGEBRAIC[194] is Af_r1 in component NCX (J_per_mol).
 * ALGEBRAIC[200] is Ar_r1 in component NCX (J_per_mol).
 * ALGEBRAIC[201] is v_r1 in component NCX (fmol_per_sec).
 * ALGEBRAIC[195] is Af_r2 in component NCX (J_per_mol).
 * ALGEBRAIC[202] is Ar_r2 in component NCX (J_per_mol).
 * ALGEBRAIC[204] is v_r2 in component NCX (fmol_per_sec).
 * ALGEBRAIC[196] is Af_r3 in component NCX (J_per_mol).
 * ALGEBRAIC[205] is Ar_r3 in component NCX (J_per_mol).
 * ALGEBRAIC[207] is v_r3 in component NCX (fmol_per_sec).
 * ALGEBRAIC[197] is Af_r4 in component NCX (J_per_mol).
 * ALGEBRAIC[208] is Ar_r4 in component NCX (J_per_mol).
 * ALGEBRAIC[210] is v_r4 in component NCX (fmol_per_sec).
 * ALGEBRAIC[198] is Af_r5 in component NCX (J_per_mol).
 * ALGEBRAIC[211] is Ar_r5 in component NCX (J_per_mol).
 * ALGEBRAIC[213] is v_r5 in component NCX (fmol_per_sec).
 * ALGEBRAIC[199] is Af_r6 in component NCX (J_per_mol).
 * ALGEBRAIC[216] is Ar_r6 in component NCX (J_per_mol).
 * ALGEBRAIC[217] is v_r6 in component NCX (fmol_per_sec).
 * ALGEBRAIC[183] is V_mem in component NCX (volt).
 * ALGEBRAIC[223] is mu_Ca_i in component pCa (J_per_mol).
 * ALGEBRAIC[224] is mu_Ca_o in component pCa (J_per_mol).
 * ALGEBRAIC[225] is mu_pCa in component pCa (J_per_mol).
 * ALGEBRAIC[226] is mu_pCa_Ca in component pCa (J_per_mol).
 * ALGEBRAIC[0] is mu_PLB in component PLB (J_per_mol).
 * ALGEBRAIC[25] is mu_PKACI in component PLB (J_per_mol).
 * ALGEBRAIC[29] is mu_PLB_PKACI in component PLB (J_per_mol).
 * ALGEBRAIC[34] is mu_PP1 in component PLB (J_per_mol).
 * ALGEBRAIC[39] is mu_PLBp_PP1 in component PLB (J_per_mol).
 * ALGEBRAIC[49] is mu_PLBp in component PLB (J_per_mol).
 * ALGEBRAIC[40] is mu_Ip in component PLB (J_per_mol).
 * ALGEBRAIC[50] is mu_Ip_PP1 in component PLB (J_per_mol).
 * CONSTANTS[163] is nCa_1 in component RyR (dimensionless).
 * CONSTANTS[164] is nCa_2 in component RyR (dimensionless).
 * ALGEBRAIC[229] is mu_Ca_SR in component RyR (J_per_mol).
 * ALGEBRAIC[231] is mu_Ca_D in component RyR (J_per_mol).
 * ALGEBRAIC[17] is q_Ca_gate_complexes in component RyR (fmol).
 * ALGEBRAIC[233] is mu_C_RyR in component RyR (J_per_mol).
 * ALGEBRAIC[234] is mu_CI_RyR in component RyR (J_per_mol).
 * ALGEBRAIC[235] is mu_I_RyR in component RyR (J_per_mol).
 * ALGEBRAIC[236] is mu_O_RyR in component RyR (J_per_mol).
 * ALGEBRAIC[238] is v_OC in component RyR (fmol_per_sec).
 * ALGEBRAIC[240] is v_CCI in component RyR (fmol_per_sec).
 * ALGEBRAIC[242] is v_CII in component RyR (fmol_per_sec).
 * ALGEBRAIC[244] is v_IO in component RyR (fmol_per_sec).
 * CONSTANTS[165] is n_Ca_i in component SERCA (dimensionless).
 * CONSTANTS[166] is n_Ca_SR in component SERCA (dimensionless).
 * CONSTANTS[167] is n_H in component SERCA (dimensionless).
 * ALGEBRAIC[18] is q_SERCA_Ca_complexes in component SERCA (fmol).
 * ALGEBRAIC[19] is c_Ca_i in component SERCA (mM).
 * ALGEBRAIC[20] is c_Ca_SR in component SERCA (mM).
 * ALGEBRAIC[21] is c_H in component SERCA (mM).
 * ALGEBRAIC[22] is c_MgADP in component SERCA (mM).
 * ALGEBRAIC[23] is c_MgATP in component SERCA (mM).
 * ALGEBRAIC[24] is c_Pi in component SERCA (mM).
 * CONSTANTS[168] is vol_i in component SERCA (pL).
 * CONSTANTS[174] is vol_sr in component SERCA (pL).
 * CONSTANTS[176] is vol_isr in component SERCA (pL).
 * ALGEBRAIC[239] is mu_Ca_i in component SERCA (J_per_mol).
 * ALGEBRAIC[241] is mu_Ca_SR in component SERCA (J_per_mol).
 * ALGEBRAIC[243] is mu_H in component SERCA (J_per_mol).
 * ALGEBRAIC[245] is mu_MgADP in component SERCA (J_per_mol).
 * ALGEBRAIC[247] is mu_MgATP in component SERCA (J_per_mol).
 * ALGEBRAIC[249] is mu_P1 in component SERCA (J_per_mol).
 * ALGEBRAIC[284] is v_P1 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[248] is mu_Pi in component SERCA (J_per_mol).
 * ALGEBRAIC[250] is mu_P2 in component SERCA (J_per_mol).
 * ALGEBRAIC[291] is v_P2 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[251] is mu_P2a in component SERCA (J_per_mol).
 * ALGEBRAIC[288] is v_P2a in component SERCA (fmol_per_sec).
 * ALGEBRAIC[252] is mu_P4 in component SERCA (J_per_mol).
 * ALGEBRAIC[293] is v_P4 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[253] is mu_P5 in component SERCA (J_per_mol).
 * ALGEBRAIC[294] is v_P5 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[255] is mu_P6 in component SERCA (J_per_mol).
 * ALGEBRAIC[297] is v_P6 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[257] is mu_P8 in component SERCA (J_per_mol).
 * ALGEBRAIC[299] is v_P8 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[258] is mu_P9 in component SERCA (J_per_mol).
 * ALGEBRAIC[301] is v_P9 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[254] is mu_P10 in component SERCA (J_per_mol).
 * ALGEBRAIC[302] is v_P10 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[259] is mu_PLBp_baseline in component SERCA (J_per_mol).
 * CONSTANTS[173] is mu_PLBp in component SERCA (J_per_mol).
 * ALGEBRAIC[260] is mu_PLBp_eff in component SERCA (J_per_mol).
 * ALGEBRAIC[263] is Af_R1_2 in component SERCA (J_per_mol).
 * ALGEBRAIC[264] is Ar_R1_2 in component SERCA (J_per_mol).
 * ALGEBRAIC[280] is v_SERCA_R1_2 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[265] is Af_R5_6 in component SERCA (J_per_mol).
 * ALGEBRAIC[266] is Ar_R5_6 in component SERCA (J_per_mol).
 * ALGEBRAIC[282] is v_SERCA_R5_6 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[267] is Af_R2_2a in component SERCA (J_per_mol).
 * ALGEBRAIC[268] is Ar_R2_2a in component SERCA (J_per_mol).
 * ALGEBRAIC[285] is v_SERCA_R2_2a in component SERCA (fmol_per_sec).
 * ALGEBRAIC[269] is Af_R2_4 in component SERCA (J_per_mol).
 * ALGEBRAIC[270] is Ar_R2_4 in component SERCA (J_per_mol).
 * ALGEBRAIC[287] is v_SERCA_R2_4 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[271] is Af_R4_5 in component SERCA (J_per_mol).
 * ALGEBRAIC[272] is Ar_R4_5 in component SERCA (J_per_mol).
 * ALGEBRAIC[290] is v_SERCA_R4_5 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[273] is Af_R6_8 in component SERCA (J_per_mol).
 * ALGEBRAIC[274] is Ar_R6_8 in component SERCA (J_per_mol).
 * ALGEBRAIC[292] is v_SERCA_R6_8 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[275] is Af_R8_9 in component SERCA (J_per_mol).
 * ALGEBRAIC[276] is Ar_R8_9 in component SERCA (J_per_mol).
 * ALGEBRAIC[296] is v_SERCA_R8_9 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[277] is Af_R9_10 in component SERCA (J_per_mol).
 * ALGEBRAIC[278] is Ar_R9_10 in component SERCA (J_per_mol).
 * ALGEBRAIC[298] is v_SERCA_R9_10 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[261] is Af_R10_1 in component SERCA (J_per_mol).
 * ALGEBRAIC[262] is Ar_R10_1 in component SERCA (J_per_mol).
 * ALGEBRAIC[279] is v_SERCA_R10_1 in component SERCA (fmol_per_sec).
 * ALGEBRAIC[256] is v_SERCA_R10_1_noPLBp in component SERCA (fmol_per_sec).
 * RATES[20] is d/dt q_Ca_SR in component environment (fmol).
 * RATES[4] is d/dt q_Ca_i in component environment (fmol).
 * RATES[2] is d/dt q_TRPN in component environment (fmol).
 * RATES[3] is d/dt q_Ca_TRPN in component environment (fmol).
 * RATES[7] is d/dt q_B_XB in component environment (fmol).
 * RATES[8] is d/dt q_U_XB in component environment (fmol).
 * RATES[9] is d/dt q_W_XB in component environment (fmol).
 * RATES[10] is d/dt q_S_XB in component environment (fmol).
 * RATES[1] is d/dt q_Ca_D in component environment (fmol).
 * RATES[11] is d/dt q_Ca_o in component environment (fmol).
 * RATES[5] is d/dt q_Na_i in component environment (fmol).
 * RATES[6] is d/dt q_Na_o in component environment (fmol).
 * RATES[21] is d/dt q_pCa in component environment (fmol).
 * RATES[22] is d/dt q_pCa_Ca in component environment (fmol).
 * RATES[12] is d/dt q_PLB in component environment (fmol).
 * RATES[13] is d/dt q_PKACI in component environment (fmol).
 * RATES[14] is d/dt q_PLB_PKACI in component environment (fmol).
 * RATES[15] is d/dt q_PP1 in component environment (fmol).
 * RATES[16] is d/dt q_PLBp_PP1 in component environment (fmol).
 * RATES[17] is d/dt q_PLBp in component environment (fmol).
 * RATES[18] is d/dt q_Ip in component environment (fmol).
 * RATES[19] is d/dt q_Ip_PP1 in component environment (fmol).
 * RATES[23] is d/dt q_H in component environment (fmol).
 * RATES[25] is d/dt q_MgADP in component environment (fmol).
 * RATES[26] is d/dt q_MgATP in component environment (fmol).
 * RATES[24] is d/dt q_Pi in component environment (fmol).
 * RATES[0] is d/dt q_membrane in component environment (fC).
 * RATES[59] is d/dt G_w in component crossbridge_TRPN (metre).
 * RATES[58] is d/dt SL in component crossbridge_TRPN (metre).
 * RATES[60] is d/dt G_s in component crossbridge_TRPN (metre).
 * RATES[61] is d/dt Cdd in component crossbridge_TRPN (metre).
 * RATES[27] is d/dt q_000_LCC in component environment (fmol).
 * RATES[28] is d/dt q_010_LCC in component environment (fmol).
 * RATES[29] is d/dt q_020_LCC in component environment (fmol).
 * RATES[30] is d/dt q_100_LCC in component environment (fmol).
 * RATES[31] is d/dt q_110_LCC in component environment (fmol).
 * RATES[32] is d/dt q_120_LCC in component environment (fmol).
 * RATES[33] is d/dt q_001_LCC in component environment (fmol).
 * RATES[34] is d/dt q_011_LCC in component environment (fmol).
 * RATES[35] is d/dt q_021_LCC in component environment (fmol).
 * RATES[36] is d/dt q_101_LCC in component environment (fmol).
 * RATES[37] is d/dt q_111_LCC in component environment (fmol).
 * RATES[38] is d/dt q_121_LCC in component environment (fmol).
 * RATES[39] is d/dt q_P1_NCX in component environment (fmol).
 * RATES[40] is d/dt q_P2_NCX in component environment (fmol).
 * RATES[41] is d/dt q_P3_NCX in component environment (fmol).
 * RATES[42] is d/dt q_P4_NCX in component environment (fmol).
 * RATES[43] is d/dt q_P5_NCX in component environment (fmol).
 * RATES[44] is d/dt q_P6_NCX in component environment (fmol).
 * RATES[48] is d/dt q_O_RyR in component environment (fmol).
 * RATES[45] is d/dt q_C_RyR in component environment (fmol).
 * RATES[46] is d/dt q_CI_RyR in component environment (fmol).
 * RATES[47] is d/dt q_I_RyR in component environment (fmol).
 * RATES[49] is d/dt q_P1_SERCA in component environment (fmol).
 * RATES[50] is d/dt q_P2_SERCA in component environment (fmol).
 * RATES[51] is d/dt q_P2a_SERCA in component environment (fmol).
 * RATES[52] is d/dt q_P4_SERCA in component environment (fmol).
 * RATES[53] is d/dt q_P5_SERCA in component environment (fmol).
 * RATES[57] is d/dt q_P10_SERCA in component environment (fmol).
 * RATES[54] is d/dt q_P6_SERCA in component environment (fmol).
 * RATES[55] is d/dt q_P8_SERCA in component environment (fmol).
 * RATES[56] is d/dt q_P9_SERCA in component environment (fmol).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1.19278e-07;
CONSTANTS[1] = 61.975;
CONSTANTS[2] = 6.70563;
CONSTANTS[3] = 0.295103;
CONSTANTS[4] = 0.00694361;
CONSTANTS[5] = 124985;
CONSTANTS[6] = 0.0319966;
CONSTANTS[7] = 1.45746;
CONSTANTS[8] = 1.91895;
CONSTANTS[9] = 306.532;
CONSTANTS[10] = 15.4798;
CONSTANTS[11] = 20.3814;
CONSTANTS[12] = 0.0998404;
CONSTANTS[13] = 0.00504193;
CONSTANTS[14] = 0.00663841;
CONSTANTS[15] = 29.5266;
CONSTANTS[16] = 146.621;
CONSTANTS[17] = 0.00961708;
CONSTANTS[18] = 0.047756;
CONSTANTS[19] = 1.91912;
CONSTANTS[20] = 9.52987;
CONSTANTS[21] = 0.000625075;
CONSTANTS[22] = 0.00310397;
CONSTANTS[23] = 4187.14;
CONSTANTS[24] = 20792.3;
CONSTANTS[25] = 1.36379;
CONSTANTS[26] = 6.77224;
CONSTANTS[27] = 34185.2;
CONSTANTS[28] = 169755;
CONSTANTS[29] = 1726.35;
CONSTANTS[30] = 8572.61;
CONSTANTS[31] = 2272.98;
CONSTANTS[32] = 11287.1;
CONSTANTS[33] = 429216;
CONSTANTS[34] = 1.06923e+08;
CONSTANTS[35] = 3622.96;
CONSTANTS[36] = 1.06923e+08;
CONSTANTS[37] = 429216;
CONSTANTS[38] = 0.380974;
CONSTANTS[39] = 10726.4;
CONSTANTS[40] = 1.40741e-05;
CONSTANTS[41] = 12.9276;
CONSTANTS[42] = 1.86249;
CONSTANTS[43] = 0.0640148;
CONSTANTS[44] = 0.200763;
CONSTANTS[45] = 1269.27;
CONSTANTS[46] = 1.13157e+06;
CONSTANTS[47] = 8.50173;
CONSTANTS[48] = 0.164843;
CONSTANTS[49] = 626.393;
CONSTANTS[50] = 52.1994;
CONSTANTS[51] = 3.59662e-07;
CONSTANTS[52] = 6.73317;
CONSTANTS[53] = 32063.4;
CONSTANTS[54] = 6.73317;
CONSTANTS[55] = 0.408566;
CONSTANTS[56] = 313.54;
CONSTANTS[57] = 6675;
CONSTANTS[58] = 6675;
CONSTANTS[59] = 5.05704e-06;
CONSTANTS[60] = 4.69975;
CONSTANTS[61] = 1.11404;
CONSTANTS[62] = 12.0566;
CONSTANTS[63] = 0.062444;
CONSTANTS[64] = 9.09233;
CONSTANTS[65] = 3.40963;
CONSTANTS[66] = 22.2937;
CONSTANTS[67] = 5.57343e-06;
CONSTANTS[68] = 66.2656;
CONSTANTS[69] = 0.0265011;
CONSTANTS[70] = 1.58796;
CONSTANTS[71] = 31.4449;
CONSTANTS[72] = 23.8827;
CONSTANTS[73] = 0.319783;
CONSTANTS[74] = 6.33236;
CONSTANTS[75] = 4.80948;
CONSTANTS[76] = 4875.4;
CONSTANTS[77] = 96542.9;
CONSTANTS[78] = 73325.1;
CONSTANTS[79] = 981.805;
CONSTANTS[80] = 19441.8;
CONSTANTS[81] = 14766.2;
CONSTANTS[82] = 2.85254e+08;
CONSTANTS[83] = 16.9249;
CONSTANTS[84] = 34379.9;
CONSTANTS[85] = 313010;
CONSTANTS[86] = 6.01862e+07;
CONSTANTS[87] = 929.686;
CONSTANTS[88] = 0.012842;
CONSTANTS[89] = 0.0107799;
CONSTANTS[90] = 0.0126813;
CONSTANTS[91] = 0.0178289;
CONSTANTS[92] = 0.482706;
CONSTANTS[93] = 0.0120458;
CONSTANTS[94] = 2.99358;
CONSTANTS[95] = 0.0189254;
CONSTANTS[96] = 1.63848;
CONSTANTS[97] = 66.755;
CONSTANTS[98] = 62.347;
CONSTANTS[99] = 0.0304897;
CONSTANTS[100] = 2234.81;
CONSTANTS[101] = 30.332;
CONSTANTS[102] = 0.0957866;
CONSTANTS[103] = 7.05739;
CONSTANTS[104] = 62.115;
CONSTANTS[105] = 1.88391;
CONSTANTS[106] = 0.271481;
CONSTANTS[107] = 1292.79;
CONSTANTS[108] = 0.00208939;
CONSTANTS[109] = 139.3;
CONSTANTS[110] = 0.710561;
CONSTANTS[111] = 1.30406;
CONSTANTS[112] = 1.3624;
CONSTANTS[113] = 486.59;
CONSTANTS[114] = 15.0787;
CONSTANTS[115] = 1.53572e-05;
CONSTANTS[116] = 0.00616948;
CONSTANTS[117] = 2;
CONSTANTS[118] = 2.1404;
CONSTANTS[119] = -1.1495;
CONSTANTS[120] = 0.72162;
CONSTANTS[121] = 4.2933;
CONSTANTS[122] = -2.1404;
CONSTANTS[123] = 1.8993;
CONSTANTS[124] = -0.52288;
CONSTANTS[125] = 0;
CONSTANTS[126] = 3;
CONSTANTS[127] = 3;
CONSTANTS[128] = -0.55679;
CONSTANTS[129] = 0.44321;
CONSTANTS[130] = 1.1;
CONSTANTS[131] = 2;
CONSTANTS[132] = 25.84;
CONSTANTS[133] = 153400;
STATES[0] = -13039;
STATES[1] = 3.68876E-07;
STATES[2] = 2.57;
STATES[3] = 1e-18;
STATES[4] = 0.005663912;
STATES[5] = 481.46;
STATES[6] = 725.48;
STATES[7] = 1;
STATES[8] = 0;
STATES[9] = 0;
STATES[10] = 0;
STATES[11] = 6.84;
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;
CONSTANTS[134] = 1e-18;
STATES[18] = 1.999E-03;
STATES[19] = 1e-18;
STATES[20] = 0.641;
STATES[21] = 0.0032;
STATES[22] = 1e-9;
STATES[23] = 2.584e-6;
STATES[24] = 570;
STATES[25] = 1.3794;
STATES[26] = 3.8;
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] = 1.3798e-05;
STATES[38] = 1.3798e-05;
STATES[39] = 4.24554e-07;
STATES[40] = 4.24554e-07;
STATES[41] = 4.24554e-07;
STATES[42] = 4.24554e-07;
STATES[43] = 4.24554e-07;
STATES[44] = 4.24554e-07;
STATES[45] = 4.13941e-05;
STATES[46] = 4.13941e-05;
STATES[47] = 4.13941e-05;
STATES[48] = 4.13941e-05;
STATES[49] = 2.83036e-08;
STATES[50] = 2.83036e-08;
STATES[51] = 2.83036e-08;
STATES[52] = 2.83036e-08;
STATES[53] = 2.83036e-08;
STATES[54] = 2.83036e-08;
STATES[55] = 2.83036e-08;
STATES[56] = 2.83036e-08;
STATES[57] = 2.83036e-08;
STATES[58] = 2e-6;
CONSTANTS[135] = 1;
CONSTANTS[136] = 0.001;
CONSTANTS[137] = 19e-6;
CONSTANTS[138] = 38;
CONSTANTS[139] = 5.182;
CONSTANTS[140] = 0.0025375;
CONSTANTS[141] = 8.31;
CONSTANTS[142] = 310;
CONSTANTS[143] = 96485;
CONSTANTS[144] = 1;
CONSTANTS[145] = -1;
CONSTANTS[146] = 1;
CONSTANTS[147] = 2;
CONSTANTS[148] = 2e-6;
CONSTANTS[149] = 1e-6;
CONSTANTS[150] = 0.25;
CONSTANTS[151] = 0.5;
CONSTANTS[152] = 25;
CONSTANTS[153] = 2.23;
CONSTANTS[154] = 26;
CONSTANTS[155] = 4;
STATES[59] = 1e-6;
STATES[60] = 1e-6;
CONSTANTS[156] = 1;
CONSTANTS[157] = 40.5;
STATES[61] = 0;
CONSTANTS[158] = 200e3;
CONSTANTS[159] = 20e3;
CONSTANTS[160] = 7;
CONSTANTS[161] = 2.1;
CONSTANTS[162] = 9.1;
CONSTANTS[163] = 1;
CONSTANTS[164] = 2;
CONSTANTS[165] = 2;
CONSTANTS[166] = 2;
CONSTANTS[167] = 2;
CONSTANTS[168] = 34.0;
CONSTANTS[169] = ( CONSTANTS[152]*CONSTANTS[150])/( (1.00000 - CONSTANTS[150])*CONSTANTS[151]+CONSTANTS[150]);
CONSTANTS[170] =  CONSTANTS[138]*0.0350000;
CONSTANTS[171] = ( CONSTANTS[153]*CONSTANTS[154]*(1.00000 - CONSTANTS[151]))/CONSTANTS[151];
CONSTANTS[172] = ( CONSTANTS[153]*CONSTANTS[155]*(1.00000 - CONSTANTS[150])*CONSTANTS[151])/CONSTANTS[150];
CONSTANTS[173] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*CONSTANTS[134]);
CONSTANTS[174] =  CONSTANTS[168]*0.0350000;
CONSTANTS[175] = CONSTANTS[169];
CONSTANTS[176] = CONSTANTS[168]+CONSTANTS[174];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = (VOI<0.0600000 ?  1.00000*CONSTANTS[137] : VOI>0.500000&&VOI<1.10000 ?  - 0.100000*CONSTANTS[137] :  0.00000*CONSTANTS[137]);
ALGEBRAIC[26] = ALGEBRAIC[2];
RATES[58] = ALGEBRAIC[26];
rootfind_0(VOI, CONSTANTS, RATES, STATES, ALGEBRAIC, pret);
RATES[61] = ALGEBRAIC[44];
ALGEBRAIC[30] =  - CONSTANTS[169]*ALGEBRAIC[26];
ALGEBRAIC[35] =  CONSTANTS[171]*STATES[59];
ALGEBRAIC[41] = - ALGEBRAIC[30] - ALGEBRAIC[35];
ALGEBRAIC[51] = ALGEBRAIC[41];
RATES[59] = ALGEBRAIC[51];
ALGEBRAIC[31] =  - CONSTANTS[175]*ALGEBRAIC[26];
ALGEBRAIC[36] =  CONSTANTS[172]*STATES[60];
ALGEBRAIC[42] = - ALGEBRAIC[31] - ALGEBRAIC[36];
ALGEBRAIC[52] = ALGEBRAIC[42];
RATES[60] = ALGEBRAIC[52];
ALGEBRAIC[34] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[95]*STATES[15]);
ALGEBRAIC[40] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[98]*STATES[18]);
ALGEBRAIC[50] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[99]*STATES[19]);
ALGEBRAIC[56] =  CONSTANTS[45]*(exp((ALGEBRAIC[34]+ALGEBRAIC[40])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[50]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[18] = - ALGEBRAIC[56];
RATES[19] = ALGEBRAIC[56];
ALGEBRAIC[0] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[92]*STATES[12]);
ALGEBRAIC[25] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[93]*STATES[13]);
ALGEBRAIC[29] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[94]*STATES[14]);
ALGEBRAIC[55] =  CONSTANTS[41]*(exp((ALGEBRAIC[0]+ALGEBRAIC[25])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[29]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[49] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]);
ALGEBRAIC[60] =  CONSTANTS[42]*(exp(ALGEBRAIC[29]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[49]+ALGEBRAIC[25])/( CONSTANTS[141]*CONSTANTS[142])));
RATES[13] = ALGEBRAIC[60] - ALGEBRAIC[55];
RATES[14] = ALGEBRAIC[55] - ALGEBRAIC[60];
RATES[15] = ALGEBRAIC[55] - ALGEBRAIC[60];
ALGEBRAIC[39] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[96]*STATES[16]);
ALGEBRAIC[62] =  CONSTANTS[43]*(exp((ALGEBRAIC[49]+ALGEBRAIC[34])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[39]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[17] = ALGEBRAIC[60] - ALGEBRAIC[62];
ALGEBRAIC[58] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[62]*STATES[2]);
ALGEBRAIC[61] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[63] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[63]*STATES[3]);
ALGEBRAIC[66] =  CONSTANTS[1]*(exp(( CONSTANTS[131]*ALGEBRAIC[61]+ALGEBRAIC[58])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[63]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[2] = - ALGEBRAIC[66];
RATES[3] = ALGEBRAIC[66];
ALGEBRAIC[65] =  CONSTANTS[44]*(exp(ALGEBRAIC[39]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[0]+ALGEBRAIC[34])/( CONSTANTS[141]*CONSTANTS[142])));
RATES[12] = ALGEBRAIC[65] - ALGEBRAIC[55];
RATES[16] = ALGEBRAIC[62] - ALGEBRAIC[65];
ALGEBRAIC[64] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[64]*STATES[7]);
ALGEBRAIC[67] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[65]*STATES[8]);
ALGEBRAIC[73] =  CONSTANTS[2]*(exp((ALGEBRAIC[64]+ CONSTANTS[130]*ALGEBRAIC[63])/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[67]+ CONSTANTS[130]*ALGEBRAIC[63])/( CONSTANTS[141]*CONSTANTS[142])));
RATES[7] = - ALGEBRAIC[73];
ALGEBRAIC[69] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[66]*STATES[9]);
ALGEBRAIC[75] =  CONSTANTS[3]*(exp(ALGEBRAIC[67]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[69]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[71] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[67]*STATES[10]);
ALGEBRAIC[77] =  CONSTANTS[4]*(exp(ALGEBRAIC[69]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[71]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[9] = ALGEBRAIC[75] - ALGEBRAIC[77];
ALGEBRAIC[79] =  CONSTANTS[5]*(exp(ALGEBRAIC[71]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[67]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[8] = (ALGEBRAIC[73] - ALGEBRAIC[75])+ALGEBRAIC[79];
RATES[10] = ALGEBRAIC[77] - ALGEBRAIC[79];
ALGEBRAIC[80] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27]);
ALGEBRAIC[108] = ALGEBRAIC[80];
ALGEBRAIC[78] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[154] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[155] =  CONSTANTS[27]*(exp(ALGEBRAIC[108]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[154]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[74] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[90] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[118] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[119] =  CONSTANTS[9]*(exp(ALGEBRAIC[90]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[118]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[96] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[130] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[131] =  CONSTANTS[15]*(exp(ALGEBRAIC[96]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[130]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[100] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[138] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[139] =  CONSTANTS[19]*(exp(ALGEBRAIC[100]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[138]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[157] = ((- ALGEBRAIC[119] - ALGEBRAIC[131]) - ALGEBRAIC[139]) - ALGEBRAIC[155];
RATES[27] = ALGEBRAIC[157];
ALGEBRAIC[91] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[120] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[121] =  CONSTANTS[12]*(exp(ALGEBRAIC[91]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[120]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[97] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[132] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[133] =  CONSTANTS[17]*(exp(ALGEBRAIC[97]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[132]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[101] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[140] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[141] =  CONSTANTS[21]*(exp(ALGEBRAIC[101]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[140]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[158] = ((- ALGEBRAIC[121] - ALGEBRAIC[133]) - ALGEBRAIC[141]) - - ALGEBRAIC[155];
RATES[33] = ALGEBRAIC[158];
ALGEBRAIC[81] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28]);
ALGEBRAIC[109] = ALGEBRAIC[81];
ALGEBRAIC[156] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[159] =  CONSTANTS[29]*(exp(ALGEBRAIC[109]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[156]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[92] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[122] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[123] =  CONSTANTS[10]*(exp(ALGEBRAIC[92]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[122]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[104] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[146] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[147] =  CONSTANTS[23]*(exp(ALGEBRAIC[104]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[146]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[161] = ((- ALGEBRAIC[123]+ALGEBRAIC[131]) - ALGEBRAIC[159]) - ALGEBRAIC[147];
RATES[28] = ALGEBRAIC[161];
ALGEBRAIC[93] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[124] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[125] =  CONSTANTS[13]*(exp(ALGEBRAIC[93]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[124]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[105] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[148] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[149] =  CONSTANTS[25]*(exp(ALGEBRAIC[105]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[148]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[162] = ((- ALGEBRAIC[125]+ALGEBRAIC[133]) - - ALGEBRAIC[159]) - ALGEBRAIC[149];
RATES[34] = ALGEBRAIC[162];
ALGEBRAIC[82] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29]);
ALGEBRAIC[110] = ALGEBRAIC[82];
ALGEBRAIC[160] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[163] =  CONSTANTS[31]*(exp(ALGEBRAIC[110]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[160]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[94] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[126] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[127] =  CONSTANTS[11]*(exp(ALGEBRAIC[94]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[126]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[165] = ((- ALGEBRAIC[127]+ALGEBRAIC[139]) - ALGEBRAIC[163])+ALGEBRAIC[147];
RATES[29] = ALGEBRAIC[165];
ALGEBRAIC[95] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[128] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[129] =  CONSTANTS[14]*(exp(ALGEBRAIC[95]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[128]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[166] = ((- ALGEBRAIC[129]+ALGEBRAIC[141]) - - ALGEBRAIC[163])+ALGEBRAIC[149];
RATES[35] = ALGEBRAIC[166];
ALGEBRAIC[83] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30]);
ALGEBRAIC[111] = ALGEBRAIC[83];
ALGEBRAIC[164] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[167] =  CONSTANTS[28]*(exp(ALGEBRAIC[111]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[164]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[98] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[134] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[135] =  CONSTANTS[16]*(exp(ALGEBRAIC[98]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[134]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[102] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[142] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[143] =  CONSTANTS[20]*(exp(ALGEBRAIC[102]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[142]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[169] = ((- ALGEBRAIC[143]+ALGEBRAIC[119]) - ALGEBRAIC[135]) - ALGEBRAIC[167];
RATES[30] = ALGEBRAIC[169];
ALGEBRAIC[99] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[136] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[137] =  CONSTANTS[18]*(exp(ALGEBRAIC[99]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[136]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[103] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[144] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[145] =  CONSTANTS[22]*(exp(ALGEBRAIC[103]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[144]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[170] = ((- ALGEBRAIC[145]+ALGEBRAIC[121]) - ALGEBRAIC[137]) - - ALGEBRAIC[167];
RATES[36] = ALGEBRAIC[170];
ALGEBRAIC[84] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31]);
ALGEBRAIC[112] = ALGEBRAIC[84];
ALGEBRAIC[168] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[171] =  CONSTANTS[30]*(exp(ALGEBRAIC[112]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[168]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[106] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[150] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[151] =  CONSTANTS[24]*(exp(ALGEBRAIC[106]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[150]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[173] = ((ALGEBRAIC[135]+ALGEBRAIC[123]) - ALGEBRAIC[171]) - ALGEBRAIC[151];
RATES[31] = ALGEBRAIC[173];
ALGEBRAIC[86] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32]);
ALGEBRAIC[113] = ALGEBRAIC[86];
ALGEBRAIC[172] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[174] =  CONSTANTS[32]*(exp(ALGEBRAIC[113]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[172]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[176] = ((ALGEBRAIC[143]+ALGEBRAIC[127]) - ALGEBRAIC[174])+ALGEBRAIC[151];
RATES[32] = ALGEBRAIC[176];
ALGEBRAIC[85] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37]);
ALGEBRAIC[88] = ALGEBRAIC[78]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74]+ALGEBRAIC[85];
ALGEBRAIC[76] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[116] = ALGEBRAIC[76]+ALGEBRAIC[85];
ALGEBRAIC[114] =  CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[175] = (ALGEBRAIC[114]==0.00000 ?  CONSTANTS[7]*(exp(ALGEBRAIC[88]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[116]/( CONSTANTS[141]*CONSTANTS[142]))) :  ((( CONSTANTS[7]*ALGEBRAIC[114])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[114]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[88]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[116]/( CONSTANTS[141]*CONSTANTS[142]))));
ALGEBRAIC[107] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[152] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[153] =  CONSTANTS[26]*(exp(ALGEBRAIC[107]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[152]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[178] = ((ALGEBRAIC[137]+ALGEBRAIC[125]) - ((- ALGEBRAIC[171]+ALGEBRAIC[175]) - ALGEBRAIC[175])) - ALGEBRAIC[153];
RATES[37] = ALGEBRAIC[178];
ALGEBRAIC[87] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38]);
ALGEBRAIC[89] = ALGEBRAIC[78]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74]+ALGEBRAIC[87];
ALGEBRAIC[117] = ALGEBRAIC[76]+ALGEBRAIC[87];
ALGEBRAIC[115] =  CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[177] = (ALGEBRAIC[115]==0.00000 ?  CONSTANTS[8]*(exp(ALGEBRAIC[89]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[117]/( CONSTANTS[141]*CONSTANTS[142]))) :  ((( CONSTANTS[8]*ALGEBRAIC[115])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[115]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[89]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[117]/( CONSTANTS[141]*CONSTANTS[142]))));
ALGEBRAIC[179] = ((ALGEBRAIC[145]+ALGEBRAIC[129]) - ((- ALGEBRAIC[174]+ALGEBRAIC[177]) - ALGEBRAIC[177]))+ALGEBRAIC[153];
RATES[38] = ALGEBRAIC[179];
ALGEBRAIC[188] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39]);
ALGEBRAIC[194] = ALGEBRAIC[188];
ALGEBRAIC[187] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[88]*STATES[5]);
ALGEBRAIC[189] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[83]*STATES[40]);
ALGEBRAIC[200] =  CONSTANTS[127]*ALGEBRAIC[187]+ALGEBRAIC[189];
ALGEBRAIC[201] =  CONSTANTS[33]*(exp(ALGEBRAIC[194]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[200]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[203] =  CONSTANTS[127]*ALGEBRAIC[201];
RATES[5] = ALGEBRAIC[203];
ALGEBRAIC[185] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[195] = ALGEBRAIC[189]+ALGEBRAIC[185];
ALGEBRAIC[190] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[84]*STATES[41]);
ALGEBRAIC[202] = ALGEBRAIC[190];
ALGEBRAIC[204] =  CONSTANTS[34]*(exp(ALGEBRAIC[195]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[202]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[206] = ALGEBRAIC[201] - ALGEBRAIC[204];
RATES[40] = ALGEBRAIC[206];
ALGEBRAIC[196] = ALGEBRAIC[190];
ALGEBRAIC[191] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[85]*STATES[42]);
ALGEBRAIC[205] = ALGEBRAIC[191];
ALGEBRAIC[207] =  CONSTANTS[35]*(exp(ALGEBRAIC[196]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[205]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[209] = ALGEBRAIC[204] - ALGEBRAIC[207];
RATES[41] = ALGEBRAIC[209];
ALGEBRAIC[197] = ALGEBRAIC[191];
ALGEBRAIC[184] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[192] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[86]*STATES[43]);
ALGEBRAIC[208] = ALGEBRAIC[192]+ALGEBRAIC[184];
ALGEBRAIC[210] =  CONSTANTS[36]*(exp(ALGEBRAIC[197]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[208]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[212] = - ALGEBRAIC[210]+ALGEBRAIC[207];
RATES[42] = ALGEBRAIC[212];
ALGEBRAIC[186] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[89]*STATES[6]);
ALGEBRAIC[198] = ALGEBRAIC[192]+ CONSTANTS[126]*ALGEBRAIC[186];
ALGEBRAIC[193] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44]);
ALGEBRAIC[211] = ALGEBRAIC[193];
ALGEBRAIC[213] =  CONSTANTS[37]*(exp(ALGEBRAIC[198]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[211]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[214] =  - CONSTANTS[126]*ALGEBRAIC[213];
RATES[6] = ALGEBRAIC[214];
ALGEBRAIC[215] = - ALGEBRAIC[213]+ALGEBRAIC[210];
RATES[43] = ALGEBRAIC[215];
ALGEBRAIC[183] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[199] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44])+ CONSTANTS[128]*CONSTANTS[143]*ALGEBRAIC[183];
ALGEBRAIC[216] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39])+ CONSTANTS[129]*CONSTANTS[143]*ALGEBRAIC[183];
ALGEBRAIC[217] =  CONSTANTS[38]*(exp(ALGEBRAIC[199]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[216]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[218] = ALGEBRAIC[217] - ALGEBRAIC[201];
RATES[39] = ALGEBRAIC[218];
ALGEBRAIC[219] = - ALGEBRAIC[217]+ALGEBRAIC[213];
RATES[44] = ALGEBRAIC[219];
ALGEBRAIC[181] = ALGEBRAIC[177]+ALGEBRAIC[175];
ALGEBRAIC[221] = ALGEBRAIC[210];
ALGEBRAIC[224] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[225] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[90]*STATES[21]);
ALGEBRAIC[226] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[91]*STATES[22]);
ALGEBRAIC[228] =  CONSTANTS[40]*(exp(ALGEBRAIC[226]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[224]+ALGEBRAIC[225])/( CONSTANTS[141]*CONSTANTS[142])));
RATES[11] = ALGEBRAIC[181]+ALGEBRAIC[221]+ALGEBRAIC[228];
ALGEBRAIC[223] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[227] =  CONSTANTS[39]*(exp((ALGEBRAIC[223]+ALGEBRAIC[225])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[226]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[21] = - ALGEBRAIC[227]+ALGEBRAIC[228];
RATES[22] = ALGEBRAIC[227] - ALGEBRAIC[228];
ALGEBRAIC[180] =  CONSTANTS[143]*(((((((((((((((((((((((((((((((((((( CONSTANTS[124]*ALGEBRAIC[141] -  CONSTANTS[120]*ALGEBRAIC[141])+ CONSTANTS[125]*ALGEBRAIC[149]) -  CONSTANTS[121]*ALGEBRAIC[149])+ CONSTANTS[123]*ALGEBRAIC[137]) -  CONSTANTS[119]*ALGEBRAIC[137])+ CONSTANTS[124]*ALGEBRAIC[145]) -  CONSTANTS[120]*ALGEBRAIC[145])+ CONSTANTS[125]*ALGEBRAIC[153]) -  CONSTANTS[121]*ALGEBRAIC[153])+ CONSTANTS[123]*ALGEBRAIC[133]) -  CONSTANTS[119]*ALGEBRAIC[133])+ CONSTANTS[122]*ALGEBRAIC[129]) -  CONSTANTS[118]*ALGEBRAIC[129])+ CONSTANTS[122]*ALGEBRAIC[125]) -  CONSTANTS[118]*ALGEBRAIC[125])+ CONSTANTS[122]*ALGEBRAIC[121]) -  CONSTANTS[118]*ALGEBRAIC[121])+ CONSTANTS[125]*ALGEBRAIC[151]) -  CONSTANTS[121]*ALGEBRAIC[151])+ CONSTANTS[124]*ALGEBRAIC[143]) -  CONSTANTS[120]*ALGEBRAIC[143])+ CONSTANTS[123]*ALGEBRAIC[135]) -  CONSTANTS[119]*ALGEBRAIC[135])+ CONSTANTS[125]*ALGEBRAIC[147]) -  CONSTANTS[121]*ALGEBRAIC[147])+ CONSTANTS[124]*ALGEBRAIC[139]) -  CONSTANTS[120]*ALGEBRAIC[139])+ CONSTANTS[123]*ALGEBRAIC[131]) -  CONSTANTS[119]*ALGEBRAIC[131])+ CONSTANTS[122]*ALGEBRAIC[127]) -  CONSTANTS[118]*ALGEBRAIC[127])+ CONSTANTS[122]*ALGEBRAIC[123]) -  CONSTANTS[118]*ALGEBRAIC[123])+ CONSTANTS[122]*ALGEBRAIC[119]) -  CONSTANTS[118]*ALGEBRAIC[119]) - ( CONSTANTS[147]*ALGEBRAIC[175]+ CONSTANTS[147]*ALGEBRAIC[177]));
ALGEBRAIC[220] =  CONSTANTS[143]*( - CONSTANTS[128]*ALGEBRAIC[217]+ CONSTANTS[129]*ALGEBRAIC[217]);
ALGEBRAIC[230] =  - CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[228];
ALGEBRAIC[11] = VOI -  floor(VOI/CONSTANTS[135])*CONSTANTS[135];
ALGEBRAIC[28] = (ALGEBRAIC[11]>=0.300000&&ALGEBRAIC[11]<=0.300000+CONSTANTS[136] ? ( 0.00000*0.0300000*CONSTANTS[133])/CONSTANTS[136] : 0.00000);
ALGEBRAIC[232] = ALGEBRAIC[180]+ALGEBRAIC[220]+ALGEBRAIC[230]+ALGEBRAIC[28];
RATES[0] = ALGEBRAIC[232];
ALGEBRAIC[231] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[233] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[100]*STATES[45]);
ALGEBRAIC[236] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[103]*STATES[48]);
ALGEBRAIC[238] =  CONSTANTS[47]*(exp(ALGEBRAIC[236]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[233]+ CONSTANTS[164]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[234] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[101]*STATES[46]);
ALGEBRAIC[240] =  CONSTANTS[48]*(exp((ALGEBRAIC[233]+ CONSTANTS[163]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[234]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[45] = ALGEBRAIC[238] - ALGEBRAIC[240];
ALGEBRAIC[235] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[102]*STATES[47]);
ALGEBRAIC[242] =  CONSTANTS[49]*(exp((ALGEBRAIC[234]+ CONSTANTS[164]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[235]/( CONSTANTS[141]*CONSTANTS[142])));
RATES[46] = ALGEBRAIC[240] - ALGEBRAIC[242];
ALGEBRAIC[244] =  CONSTANTS[50]*(exp(ALGEBRAIC[235]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[236]+ CONSTANTS[163]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])));
RATES[48] = ALGEBRAIC[244] - ALGEBRAIC[238];
RATES[47] = ALGEBRAIC[242] - ALGEBRAIC[244];
ALGEBRAIC[68] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[70] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[72] =  CONSTANTS[6]*(exp(ALGEBRAIC[68]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[70]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[182] = (( (CONSTANTS[117])*ALGEBRAIC[167]+ CONSTANTS[117]*ALGEBRAIC[171]+ CONSTANTS[117]*ALGEBRAIC[163]+ CONSTANTS[117]*ALGEBRAIC[159]+ CONSTANTS[117]*ALGEBRAIC[155]+ CONSTANTS[117]*ALGEBRAIC[174]) - ALGEBRAIC[175]) - ALGEBRAIC[177];
ALGEBRAIC[229] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[237] =  CONSTANTS[46]*exp(ALGEBRAIC[236]/( CONSTANTS[141]*CONSTANTS[142]))*(exp(ALGEBRAIC[229]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[231]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[246] = (( CONSTANTS[164]*ALGEBRAIC[238] -  CONSTANTS[163]*ALGEBRAIC[240]) -  CONSTANTS[164]*ALGEBRAIC[242])+ CONSTANTS[163]*ALGEBRAIC[244];
RATES[1] = - ALGEBRAIC[72]+ALGEBRAIC[182]+ALGEBRAIC[237]+ALGEBRAIC[246];
ALGEBRAIC[254] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[112]*STATES[57]);
ALGEBRAIC[259] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]);
ALGEBRAIC[260] = ALGEBRAIC[259] - CONSTANTS[173];
ALGEBRAIC[261] = ALGEBRAIC[254]+ALGEBRAIC[260];
ALGEBRAIC[249] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[104]*STATES[49]);
ALGEBRAIC[248] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[116]*STATES[24]);
ALGEBRAIC[262] = ALGEBRAIC[249]+ALGEBRAIC[248]+ALGEBRAIC[260];
ALGEBRAIC[279] =  CONSTANTS[59]*(exp(ALGEBRAIC[261]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[262]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[281] = ALGEBRAIC[279];
RATES[24] = ALGEBRAIC[281];
ALGEBRAIC[247] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[114]*STATES[26]);
ALGEBRAIC[263] = ALGEBRAIC[249]+ALGEBRAIC[247]+ALGEBRAIC[260];
ALGEBRAIC[250] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[105]*STATES[50]);
ALGEBRAIC[264] = ALGEBRAIC[250]+ALGEBRAIC[260];
ALGEBRAIC[280] =  CONSTANTS[51]*(exp(ALGEBRAIC[263]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[264]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[283] = - ALGEBRAIC[280];
RATES[26] = ALGEBRAIC[283];
ALGEBRAIC[284] = ALGEBRAIC[279] - ALGEBRAIC[280];
RATES[49] = ALGEBRAIC[284];
ALGEBRAIC[253] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[108]*STATES[53]);
ALGEBRAIC[265] = ALGEBRAIC[253]+ALGEBRAIC[260];
ALGEBRAIC[245] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[115]*STATES[25]);
ALGEBRAIC[255] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[109]*STATES[54]);
ALGEBRAIC[266] = ALGEBRAIC[245]+ALGEBRAIC[255]+ALGEBRAIC[260];
ALGEBRAIC[282] =  CONSTANTS[55]*(exp(ALGEBRAIC[265]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[266]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[286] = ALGEBRAIC[282];
RATES[25] = ALGEBRAIC[286];
ALGEBRAIC[243] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[113]*STATES[23]);
ALGEBRAIC[267] = ALGEBRAIC[250]+ALGEBRAIC[243]+ALGEBRAIC[260];
ALGEBRAIC[251] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[106]*STATES[51]);
ALGEBRAIC[268] = ALGEBRAIC[251]+ALGEBRAIC[260];
ALGEBRAIC[285] =  CONSTANTS[53]*(exp(ALGEBRAIC[267]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[268]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[288] = ALGEBRAIC[285];
RATES[51] = ALGEBRAIC[288];
ALGEBRAIC[38] =  ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[48] =  ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[54] =  CONSTANTS[0]*(exp(ALGEBRAIC[38]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[48]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[222] = - ALGEBRAIC[204];
ALGEBRAIC[239] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[269] = ALGEBRAIC[250]+ CONSTANTS[165]*ALGEBRAIC[239]+ALGEBRAIC[260];
ALGEBRAIC[252] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[107]*STATES[52]);
ALGEBRAIC[270] = ALGEBRAIC[252]+ALGEBRAIC[260];
ALGEBRAIC[287] =  CONSTANTS[52]*(exp(ALGEBRAIC[269]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[270]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[289] =  - CONSTANTS[165]*ALGEBRAIC[287];
ALGEBRAIC[33] = (ALGEBRAIC[11]>=0.300000&&ALGEBRAIC[11]<=0.300000+CONSTANTS[136] ?  0.00000*8000.00 : 0.00000);
RATES[4] = (( - ALGEBRAIC[66]*CONSTANTS[131]+ALGEBRAIC[289]+ALGEBRAIC[54]+ALGEBRAIC[72]+ALGEBRAIC[222]) - ALGEBRAIC[227])+ALGEBRAIC[33];
ALGEBRAIC[291] = ALGEBRAIC[280] - ALGEBRAIC[287];
RATES[50] = ALGEBRAIC[291];
ALGEBRAIC[271] = ALGEBRAIC[252]+ALGEBRAIC[260];
ALGEBRAIC[272] = ALGEBRAIC[253]+ CONSTANTS[167]*ALGEBRAIC[243]+ALGEBRAIC[260];
ALGEBRAIC[290] =  CONSTANTS[54]*(exp(ALGEBRAIC[271]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[272]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[293] = ALGEBRAIC[287] - ALGEBRAIC[290];
RATES[52] = ALGEBRAIC[293];
ALGEBRAIC[294] = ALGEBRAIC[290] - ALGEBRAIC[282];
RATES[53] = ALGEBRAIC[294];
ALGEBRAIC[273] = ALGEBRAIC[255]+ALGEBRAIC[260];
ALGEBRAIC[241] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[257] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[110]*STATES[55]);
ALGEBRAIC[274] = ALGEBRAIC[257]+ CONSTANTS[166]*ALGEBRAIC[241]+ALGEBRAIC[260];
ALGEBRAIC[292] =  CONSTANTS[56]*(exp(ALGEBRAIC[273]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[274]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[295] =  CONSTANTS[166]*ALGEBRAIC[292];
RATES[20] = (- ALGEBRAIC[54] - ALGEBRAIC[237])+ALGEBRAIC[295];
ALGEBRAIC[297] = ALGEBRAIC[282] - ALGEBRAIC[292];
RATES[54] = ALGEBRAIC[297];
ALGEBRAIC[275] = ALGEBRAIC[257]+ CONSTANTS[167]*ALGEBRAIC[243]+ALGEBRAIC[260];
ALGEBRAIC[258] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[111]*STATES[56]);
ALGEBRAIC[276] = ALGEBRAIC[258]+ALGEBRAIC[260];
ALGEBRAIC[296] =  CONSTANTS[57]*(exp(ALGEBRAIC[275]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[276]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[299] = ALGEBRAIC[292] - ALGEBRAIC[296];
RATES[55] = ALGEBRAIC[299];
ALGEBRAIC[277] = ALGEBRAIC[258]+ALGEBRAIC[260];
ALGEBRAIC[278] = ALGEBRAIC[243]+ALGEBRAIC[254]+ALGEBRAIC[260];
ALGEBRAIC[298] =  CONSTANTS[58]*(exp(ALGEBRAIC[277]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[278]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[300] = (( CONSTANTS[167]*ALGEBRAIC[290] - ALGEBRAIC[285]) -  CONSTANTS[167]*ALGEBRAIC[296])+ALGEBRAIC[298];
RATES[23] = ALGEBRAIC[300];
ALGEBRAIC[302] = ALGEBRAIC[298] - ALGEBRAIC[279];
RATES[57] = ALGEBRAIC[302];
ALGEBRAIC[301] = ALGEBRAIC[296] - ALGEBRAIC[298];
RATES[56] = ALGEBRAIC[301];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = (VOI<0.0600000 ?  1.00000*CONSTANTS[137] : VOI>0.500000&&VOI<1.10000 ?  - 0.100000*CONSTANTS[137] :  0.00000*CONSTANTS[137]);
ALGEBRAIC[26] = ALGEBRAIC[2];
ALGEBRAIC[30] =  - CONSTANTS[169]*ALGEBRAIC[26];
ALGEBRAIC[35] =  CONSTANTS[171]*STATES[59];
ALGEBRAIC[41] = - ALGEBRAIC[30] - ALGEBRAIC[35];
ALGEBRAIC[51] = ALGEBRAIC[41];
ALGEBRAIC[31] =  - CONSTANTS[175]*ALGEBRAIC[26];
ALGEBRAIC[36] =  CONSTANTS[172]*STATES[60];
ALGEBRAIC[42] = - ALGEBRAIC[31] - ALGEBRAIC[36];
ALGEBRAIC[52] = ALGEBRAIC[42];
ALGEBRAIC[34] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[95]*STATES[15]);
ALGEBRAIC[40] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[98]*STATES[18]);
ALGEBRAIC[50] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[99]*STATES[19]);
ALGEBRAIC[56] =  CONSTANTS[45]*(exp((ALGEBRAIC[34]+ALGEBRAIC[40])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[50]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[0] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[92]*STATES[12]);
ALGEBRAIC[25] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[93]*STATES[13]);
ALGEBRAIC[29] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[94]*STATES[14]);
ALGEBRAIC[55] =  CONSTANTS[41]*(exp((ALGEBRAIC[0]+ALGEBRAIC[25])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[29]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[49] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]);
ALGEBRAIC[60] =  CONSTANTS[42]*(exp(ALGEBRAIC[29]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[49]+ALGEBRAIC[25])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[39] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[96]*STATES[16]);
ALGEBRAIC[62] =  CONSTANTS[43]*(exp((ALGEBRAIC[49]+ALGEBRAIC[34])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[39]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[58] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[62]*STATES[2]);
ALGEBRAIC[61] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[63] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[63]*STATES[3]);
ALGEBRAIC[66] =  CONSTANTS[1]*(exp(( CONSTANTS[131]*ALGEBRAIC[61]+ALGEBRAIC[58])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[63]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[65] =  CONSTANTS[44]*(exp(ALGEBRAIC[39]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[0]+ALGEBRAIC[34])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[64] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[64]*STATES[7]);
ALGEBRAIC[67] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[65]*STATES[8]);
ALGEBRAIC[73] =  CONSTANTS[2]*(exp((ALGEBRAIC[64]+ CONSTANTS[130]*ALGEBRAIC[63])/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[67]+ CONSTANTS[130]*ALGEBRAIC[63])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[69] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[66]*STATES[9]);
ALGEBRAIC[75] =  CONSTANTS[3]*(exp(ALGEBRAIC[67]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[69]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[71] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[67]*STATES[10]);
ALGEBRAIC[77] =  CONSTANTS[4]*(exp(ALGEBRAIC[69]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[71]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[79] =  CONSTANTS[5]*(exp(ALGEBRAIC[71]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[67]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[80] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27]);
ALGEBRAIC[108] = ALGEBRAIC[80];
ALGEBRAIC[78] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[154] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[155] =  CONSTANTS[27]*(exp(ALGEBRAIC[108]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[154]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[74] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[90] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[118] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[119] =  CONSTANTS[9]*(exp(ALGEBRAIC[90]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[118]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[96] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[130] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[131] =  CONSTANTS[15]*(exp(ALGEBRAIC[96]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[130]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[100] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[138] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[139] =  CONSTANTS[19]*(exp(ALGEBRAIC[100]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[138]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[157] = ((- ALGEBRAIC[119] - ALGEBRAIC[131]) - ALGEBRAIC[139]) - ALGEBRAIC[155];
ALGEBRAIC[91] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[120] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[121] =  CONSTANTS[12]*(exp(ALGEBRAIC[91]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[120]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[97] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[132] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[133] =  CONSTANTS[17]*(exp(ALGEBRAIC[97]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[132]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[101] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[140] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[141] =  CONSTANTS[21]*(exp(ALGEBRAIC[101]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[140]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[158] = ((- ALGEBRAIC[121] - ALGEBRAIC[133]) - ALGEBRAIC[141]) - - ALGEBRAIC[155];
ALGEBRAIC[81] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28]);
ALGEBRAIC[109] = ALGEBRAIC[81];
ALGEBRAIC[156] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[159] =  CONSTANTS[29]*(exp(ALGEBRAIC[109]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[156]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[92] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[122] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[123] =  CONSTANTS[10]*(exp(ALGEBRAIC[92]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[122]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[104] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[146] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[147] =  CONSTANTS[23]*(exp(ALGEBRAIC[104]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[146]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[161] = ((- ALGEBRAIC[123]+ALGEBRAIC[131]) - ALGEBRAIC[159]) - ALGEBRAIC[147];
ALGEBRAIC[93] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[124] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[125] =  CONSTANTS[13]*(exp(ALGEBRAIC[93]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[124]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[105] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[148] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[149] =  CONSTANTS[25]*(exp(ALGEBRAIC[105]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[148]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[162] = ((- ALGEBRAIC[125]+ALGEBRAIC[133]) - - ALGEBRAIC[159]) - ALGEBRAIC[149];
ALGEBRAIC[82] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29]);
ALGEBRAIC[110] = ALGEBRAIC[82];
ALGEBRAIC[160] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[163] =  CONSTANTS[31]*(exp(ALGEBRAIC[110]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[160]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[94] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[126] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[127] =  CONSTANTS[11]*(exp(ALGEBRAIC[94]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[126]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[165] = ((- ALGEBRAIC[127]+ALGEBRAIC[139]) - ALGEBRAIC[163])+ALGEBRAIC[147];
ALGEBRAIC[95] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[128] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[129] =  CONSTANTS[14]*(exp(ALGEBRAIC[95]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[128]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[166] = ((- ALGEBRAIC[129]+ALGEBRAIC[141]) - - ALGEBRAIC[163])+ALGEBRAIC[149];
ALGEBRAIC[83] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30]);
ALGEBRAIC[111] = ALGEBRAIC[83];
ALGEBRAIC[164] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[167] =  CONSTANTS[28]*(exp(ALGEBRAIC[111]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[164]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[98] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[134] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[135] =  CONSTANTS[16]*(exp(ALGEBRAIC[98]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[134]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[102] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[142] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[143] =  CONSTANTS[20]*(exp(ALGEBRAIC[102]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[142]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[169] = ((- ALGEBRAIC[143]+ALGEBRAIC[119]) - ALGEBRAIC[135]) - ALGEBRAIC[167];
ALGEBRAIC[99] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[136] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[137] =  CONSTANTS[18]*(exp(ALGEBRAIC[99]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[136]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[103] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[144] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[145] =  CONSTANTS[22]*(exp(ALGEBRAIC[103]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[144]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[170] = ((- ALGEBRAIC[145]+ALGEBRAIC[121]) - ALGEBRAIC[137]) - - ALGEBRAIC[167];
ALGEBRAIC[84] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31]);
ALGEBRAIC[112] = ALGEBRAIC[84];
ALGEBRAIC[168] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[171] =  CONSTANTS[30]*(exp(ALGEBRAIC[112]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[168]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[106] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[150] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[151] =  CONSTANTS[24]*(exp(ALGEBRAIC[106]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[150]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[173] = ((ALGEBRAIC[135]+ALGEBRAIC[123]) - ALGEBRAIC[171]) - ALGEBRAIC[151];
ALGEBRAIC[86] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32]);
ALGEBRAIC[113] = ALGEBRAIC[86];
ALGEBRAIC[172] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[117]*ALGEBRAIC[78];
ALGEBRAIC[174] =  CONSTANTS[32]*(exp(ALGEBRAIC[113]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[172]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[176] = ((ALGEBRAIC[143]+ALGEBRAIC[127]) - ALGEBRAIC[174])+ALGEBRAIC[151];
ALGEBRAIC[85] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37]);
ALGEBRAIC[88] = ALGEBRAIC[78]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74]+ALGEBRAIC[85];
ALGEBRAIC[76] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[116] = ALGEBRAIC[76]+ALGEBRAIC[85];
ALGEBRAIC[114] =  CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[175] = (ALGEBRAIC[114]==0.00000 ?  CONSTANTS[7]*(exp(ALGEBRAIC[88]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[116]/( CONSTANTS[141]*CONSTANTS[142]))) :  ((( CONSTANTS[7]*ALGEBRAIC[114])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[114]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[88]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[116]/( CONSTANTS[141]*CONSTANTS[142]))));
ALGEBRAIC[107] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[152] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[153] =  CONSTANTS[26]*(exp(ALGEBRAIC[107]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[152]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[178] = ((ALGEBRAIC[137]+ALGEBRAIC[125]) - ((- ALGEBRAIC[171]+ALGEBRAIC[175]) - ALGEBRAIC[175])) - ALGEBRAIC[153];
ALGEBRAIC[87] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38]);
ALGEBRAIC[89] = ALGEBRAIC[78]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74]+ALGEBRAIC[87];
ALGEBRAIC[117] = ALGEBRAIC[76]+ALGEBRAIC[87];
ALGEBRAIC[115] =  CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[74];
ALGEBRAIC[177] = (ALGEBRAIC[115]==0.00000 ?  CONSTANTS[8]*(exp(ALGEBRAIC[89]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[117]/( CONSTANTS[141]*CONSTANTS[142]))) :  ((( CONSTANTS[8]*ALGEBRAIC[115])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[115]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[89]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[117]/( CONSTANTS[141]*CONSTANTS[142]))));
ALGEBRAIC[179] = ((ALGEBRAIC[145]+ALGEBRAIC[129]) - ((- ALGEBRAIC[174]+ALGEBRAIC[177]) - ALGEBRAIC[177]))+ALGEBRAIC[153];
ALGEBRAIC[188] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39]);
ALGEBRAIC[194] = ALGEBRAIC[188];
ALGEBRAIC[187] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[88]*STATES[5]);
ALGEBRAIC[189] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[83]*STATES[40]);
ALGEBRAIC[200] =  CONSTANTS[127]*ALGEBRAIC[187]+ALGEBRAIC[189];
ALGEBRAIC[201] =  CONSTANTS[33]*(exp(ALGEBRAIC[194]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[200]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[203] =  CONSTANTS[127]*ALGEBRAIC[201];
ALGEBRAIC[185] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[195] = ALGEBRAIC[189]+ALGEBRAIC[185];
ALGEBRAIC[190] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[84]*STATES[41]);
ALGEBRAIC[202] = ALGEBRAIC[190];
ALGEBRAIC[204] =  CONSTANTS[34]*(exp(ALGEBRAIC[195]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[202]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[206] = ALGEBRAIC[201] - ALGEBRAIC[204];
ALGEBRAIC[196] = ALGEBRAIC[190];
ALGEBRAIC[191] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[85]*STATES[42]);
ALGEBRAIC[205] = ALGEBRAIC[191];
ALGEBRAIC[207] =  CONSTANTS[35]*(exp(ALGEBRAIC[196]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[205]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[209] = ALGEBRAIC[204] - ALGEBRAIC[207];
ALGEBRAIC[197] = ALGEBRAIC[191];
ALGEBRAIC[184] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[192] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[86]*STATES[43]);
ALGEBRAIC[208] = ALGEBRAIC[192]+ALGEBRAIC[184];
ALGEBRAIC[210] =  CONSTANTS[36]*(exp(ALGEBRAIC[197]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[208]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[212] = - ALGEBRAIC[210]+ALGEBRAIC[207];
ALGEBRAIC[186] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[89]*STATES[6]);
ALGEBRAIC[198] = ALGEBRAIC[192]+ CONSTANTS[126]*ALGEBRAIC[186];
ALGEBRAIC[193] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44]);
ALGEBRAIC[211] = ALGEBRAIC[193];
ALGEBRAIC[213] =  CONSTANTS[37]*(exp(ALGEBRAIC[198]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[211]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[214] =  - CONSTANTS[126]*ALGEBRAIC[213];
ALGEBRAIC[215] = - ALGEBRAIC[213]+ALGEBRAIC[210];
ALGEBRAIC[183] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[199] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44])+ CONSTANTS[128]*CONSTANTS[143]*ALGEBRAIC[183];
ALGEBRAIC[216] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39])+ CONSTANTS[129]*CONSTANTS[143]*ALGEBRAIC[183];
ALGEBRAIC[217] =  CONSTANTS[38]*(exp(ALGEBRAIC[199]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[216]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[218] = ALGEBRAIC[217] - ALGEBRAIC[201];
ALGEBRAIC[219] = - ALGEBRAIC[217]+ALGEBRAIC[213];
ALGEBRAIC[181] = ALGEBRAIC[177]+ALGEBRAIC[175];
ALGEBRAIC[221] = ALGEBRAIC[210];
ALGEBRAIC[224] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[225] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[90]*STATES[21]);
ALGEBRAIC[226] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[91]*STATES[22]);
ALGEBRAIC[228] =  CONSTANTS[40]*(exp(ALGEBRAIC[226]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[224]+ALGEBRAIC[225])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[223] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[227] =  CONSTANTS[39]*(exp((ALGEBRAIC[223]+ALGEBRAIC[225])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[226]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[180] =  CONSTANTS[143]*(((((((((((((((((((((((((((((((((((( CONSTANTS[124]*ALGEBRAIC[141] -  CONSTANTS[120]*ALGEBRAIC[141])+ CONSTANTS[125]*ALGEBRAIC[149]) -  CONSTANTS[121]*ALGEBRAIC[149])+ CONSTANTS[123]*ALGEBRAIC[137]) -  CONSTANTS[119]*ALGEBRAIC[137])+ CONSTANTS[124]*ALGEBRAIC[145]) -  CONSTANTS[120]*ALGEBRAIC[145])+ CONSTANTS[125]*ALGEBRAIC[153]) -  CONSTANTS[121]*ALGEBRAIC[153])+ CONSTANTS[123]*ALGEBRAIC[133]) -  CONSTANTS[119]*ALGEBRAIC[133])+ CONSTANTS[122]*ALGEBRAIC[129]) -  CONSTANTS[118]*ALGEBRAIC[129])+ CONSTANTS[122]*ALGEBRAIC[125]) -  CONSTANTS[118]*ALGEBRAIC[125])+ CONSTANTS[122]*ALGEBRAIC[121]) -  CONSTANTS[118]*ALGEBRAIC[121])+ CONSTANTS[125]*ALGEBRAIC[151]) -  CONSTANTS[121]*ALGEBRAIC[151])+ CONSTANTS[124]*ALGEBRAIC[143]) -  CONSTANTS[120]*ALGEBRAIC[143])+ CONSTANTS[123]*ALGEBRAIC[135]) -  CONSTANTS[119]*ALGEBRAIC[135])+ CONSTANTS[125]*ALGEBRAIC[147]) -  CONSTANTS[121]*ALGEBRAIC[147])+ CONSTANTS[124]*ALGEBRAIC[139]) -  CONSTANTS[120]*ALGEBRAIC[139])+ CONSTANTS[123]*ALGEBRAIC[131]) -  CONSTANTS[119]*ALGEBRAIC[131])+ CONSTANTS[122]*ALGEBRAIC[127]) -  CONSTANTS[118]*ALGEBRAIC[127])+ CONSTANTS[122]*ALGEBRAIC[123]) -  CONSTANTS[118]*ALGEBRAIC[123])+ CONSTANTS[122]*ALGEBRAIC[119]) -  CONSTANTS[118]*ALGEBRAIC[119]) - ( CONSTANTS[147]*ALGEBRAIC[175]+ CONSTANTS[147]*ALGEBRAIC[177]));
ALGEBRAIC[220] =  CONSTANTS[143]*( - CONSTANTS[128]*ALGEBRAIC[217]+ CONSTANTS[129]*ALGEBRAIC[217]);
ALGEBRAIC[230] =  - CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[228];
ALGEBRAIC[11] = VOI -  floor(VOI/CONSTANTS[135])*CONSTANTS[135];
ALGEBRAIC[28] = (ALGEBRAIC[11]>=0.300000&&ALGEBRAIC[11]<=0.300000+CONSTANTS[136] ? ( 0.00000*0.0300000*CONSTANTS[133])/CONSTANTS[136] : 0.00000);
ALGEBRAIC[232] = ALGEBRAIC[180]+ALGEBRAIC[220]+ALGEBRAIC[230]+ALGEBRAIC[28];
ALGEBRAIC[231] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[233] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[100]*STATES[45]);
ALGEBRAIC[236] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[103]*STATES[48]);
ALGEBRAIC[238] =  CONSTANTS[47]*(exp(ALGEBRAIC[236]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[233]+ CONSTANTS[164]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[234] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[101]*STATES[46]);
ALGEBRAIC[240] =  CONSTANTS[48]*(exp((ALGEBRAIC[233]+ CONSTANTS[163]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[234]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[235] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[102]*STATES[47]);
ALGEBRAIC[242] =  CONSTANTS[49]*(exp((ALGEBRAIC[234]+ CONSTANTS[164]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[235]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[244] =  CONSTANTS[50]*(exp(ALGEBRAIC[235]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[236]+ CONSTANTS[163]*ALGEBRAIC[231])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[68] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[70] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[72] =  CONSTANTS[6]*(exp(ALGEBRAIC[68]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[70]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[182] = (( (CONSTANTS[117])*ALGEBRAIC[167]+ CONSTANTS[117]*ALGEBRAIC[171]+ CONSTANTS[117]*ALGEBRAIC[163]+ CONSTANTS[117]*ALGEBRAIC[159]+ CONSTANTS[117]*ALGEBRAIC[155]+ CONSTANTS[117]*ALGEBRAIC[174]) - ALGEBRAIC[175]) - ALGEBRAIC[177];
ALGEBRAIC[229] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[237] =  CONSTANTS[46]*exp(ALGEBRAIC[236]/( CONSTANTS[141]*CONSTANTS[142]))*(exp(ALGEBRAIC[229]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[231]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[246] = (( CONSTANTS[164]*ALGEBRAIC[238] -  CONSTANTS[163]*ALGEBRAIC[240]) -  CONSTANTS[164]*ALGEBRAIC[242])+ CONSTANTS[163]*ALGEBRAIC[244];
ALGEBRAIC[254] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[112]*STATES[57]);
ALGEBRAIC[259] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]);
ALGEBRAIC[260] = ALGEBRAIC[259] - CONSTANTS[173];
ALGEBRAIC[261] = ALGEBRAIC[254]+ALGEBRAIC[260];
ALGEBRAIC[249] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[104]*STATES[49]);
ALGEBRAIC[248] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[116]*STATES[24]);
ALGEBRAIC[262] = ALGEBRAIC[249]+ALGEBRAIC[248]+ALGEBRAIC[260];
ALGEBRAIC[279] =  CONSTANTS[59]*(exp(ALGEBRAIC[261]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[262]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[281] = ALGEBRAIC[279];
ALGEBRAIC[247] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[114]*STATES[26]);
ALGEBRAIC[263] = ALGEBRAIC[249]+ALGEBRAIC[247]+ALGEBRAIC[260];
ALGEBRAIC[250] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[105]*STATES[50]);
ALGEBRAIC[264] = ALGEBRAIC[250]+ALGEBRAIC[260];
ALGEBRAIC[280] =  CONSTANTS[51]*(exp(ALGEBRAIC[263]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[264]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[283] = - ALGEBRAIC[280];
ALGEBRAIC[284] = ALGEBRAIC[279] - ALGEBRAIC[280];
ALGEBRAIC[253] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[108]*STATES[53]);
ALGEBRAIC[265] = ALGEBRAIC[253]+ALGEBRAIC[260];
ALGEBRAIC[245] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[115]*STATES[25]);
ALGEBRAIC[255] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[109]*STATES[54]);
ALGEBRAIC[266] = ALGEBRAIC[245]+ALGEBRAIC[255]+ALGEBRAIC[260];
ALGEBRAIC[282] =  CONSTANTS[55]*(exp(ALGEBRAIC[265]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[266]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[286] = ALGEBRAIC[282];
ALGEBRAIC[243] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[113]*STATES[23]);
ALGEBRAIC[267] = ALGEBRAIC[250]+ALGEBRAIC[243]+ALGEBRAIC[260];
ALGEBRAIC[251] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[106]*STATES[51]);
ALGEBRAIC[268] = ALGEBRAIC[251]+ALGEBRAIC[260];
ALGEBRAIC[285] =  CONSTANTS[53]*(exp(ALGEBRAIC[267]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[268]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[288] = ALGEBRAIC[285];
ALGEBRAIC[38] =  ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[48] =  ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[54] =  CONSTANTS[0]*(exp(ALGEBRAIC[38]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[48]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[222] = - ALGEBRAIC[204];
ALGEBRAIC[239] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[269] = ALGEBRAIC[250]+ CONSTANTS[165]*ALGEBRAIC[239]+ALGEBRAIC[260];
ALGEBRAIC[252] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[107]*STATES[52]);
ALGEBRAIC[270] = ALGEBRAIC[252]+ALGEBRAIC[260];
ALGEBRAIC[287] =  CONSTANTS[52]*(exp(ALGEBRAIC[269]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[270]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[289] =  - CONSTANTS[165]*ALGEBRAIC[287];
ALGEBRAIC[33] = (ALGEBRAIC[11]>=0.300000&&ALGEBRAIC[11]<=0.300000+CONSTANTS[136] ?  0.00000*8000.00 : 0.00000);
ALGEBRAIC[291] = ALGEBRAIC[280] - ALGEBRAIC[287];
ALGEBRAIC[271] = ALGEBRAIC[252]+ALGEBRAIC[260];
ALGEBRAIC[272] = ALGEBRAIC[253]+ CONSTANTS[167]*ALGEBRAIC[243]+ALGEBRAIC[260];
ALGEBRAIC[290] =  CONSTANTS[54]*(exp(ALGEBRAIC[271]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[272]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[293] = ALGEBRAIC[287] - ALGEBRAIC[290];
ALGEBRAIC[294] = ALGEBRAIC[290] - ALGEBRAIC[282];
ALGEBRAIC[273] = ALGEBRAIC[255]+ALGEBRAIC[260];
ALGEBRAIC[241] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[257] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[110]*STATES[55]);
ALGEBRAIC[274] = ALGEBRAIC[257]+ CONSTANTS[166]*ALGEBRAIC[241]+ALGEBRAIC[260];
ALGEBRAIC[292] =  CONSTANTS[56]*(exp(ALGEBRAIC[273]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[274]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[295] =  CONSTANTS[166]*ALGEBRAIC[292];
ALGEBRAIC[297] = ALGEBRAIC[282] - ALGEBRAIC[292];
ALGEBRAIC[275] = ALGEBRAIC[257]+ CONSTANTS[167]*ALGEBRAIC[243]+ALGEBRAIC[260];
ALGEBRAIC[258] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[111]*STATES[56]);
ALGEBRAIC[276] = ALGEBRAIC[258]+ALGEBRAIC[260];
ALGEBRAIC[296] =  CONSTANTS[57]*(exp(ALGEBRAIC[275]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[276]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[299] = ALGEBRAIC[292] - ALGEBRAIC[296];
ALGEBRAIC[277] = ALGEBRAIC[258]+ALGEBRAIC[260];
ALGEBRAIC[278] = ALGEBRAIC[243]+ALGEBRAIC[254]+ALGEBRAIC[260];
ALGEBRAIC[298] =  CONSTANTS[58]*(exp(ALGEBRAIC[277]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[278]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[300] = (( CONSTANTS[167]*ALGEBRAIC[290] - ALGEBRAIC[285]) -  CONSTANTS[167]*ALGEBRAIC[296])+ALGEBRAIC[298];
ALGEBRAIC[302] = ALGEBRAIC[298] - ALGEBRAIC[279];
ALGEBRAIC[301] = ALGEBRAIC[296] - ALGEBRAIC[298];
ALGEBRAIC[1] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[3] = STATES[4]/CONSTANTS[132];
ALGEBRAIC[4] = STATES[11]/CONSTANTS[139];
ALGEBRAIC[5] = STATES[20]/CONSTANTS[170];
ALGEBRAIC[6] = STATES[1]/CONSTANTS[140];
ALGEBRAIC[7] = STATES[5]/CONSTANTS[132];
ALGEBRAIC[8] = - arbitrary_log(STATES[23]/CONSTANTS[132], 10);
ALGEBRAIC[9] = STATES[4]+STATES[11]+STATES[3]+STATES[20]+STATES[1]+STATES[22]+ 2.00000*(STATES[52]+STATES[53]+STATES[54])+STATES[41]+STATES[42];
ALGEBRAIC[10] =  (CONSTANTS[157]/( CONSTANTS[148]*CONSTANTS[149]*CONSTANTS[150]))*STATES[10]*(STATES[60]+CONSTANTS[148]);
ALGEBRAIC[12] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33]);
ALGEBRAIC[13] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34]);
ALGEBRAIC[14] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35]);
ALGEBRAIC[15] =  CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36]);
ALGEBRAIC[16] = STATES[11]+STATES[1]+STATES[33]+STATES[34]+STATES[35]+STATES[36]+STATES[37]+STATES[38];
ALGEBRAIC[17] =  CONSTANTS[164]*STATES[48]+ CONSTANTS[163]*STATES[46]+ (CONSTANTS[163]+CONSTANTS[164])*STATES[47];
ALGEBRAIC[18] =  CONSTANTS[165]*STATES[52]+ CONSTANTS[165]*STATES[53]+ CONSTANTS[166]*STATES[54];
ALGEBRAIC[19] = STATES[4]/CONSTANTS[168];
ALGEBRAIC[20] = STATES[20]/CONSTANTS[174];
ALGEBRAIC[21] = STATES[23]/CONSTANTS[176];
ALGEBRAIC[22] = STATES[25]/CONSTANTS[168];
ALGEBRAIC[23] = STATES[26]/CONSTANTS[168];
ALGEBRAIC[24] = STATES[24]/CONSTANTS[168];
ALGEBRAIC[27] =  (CONSTANTS[157]/( CONSTANTS[148]*CONSTANTS[149]*CONSTANTS[150]))*STATES[9]*STATES[59];
ALGEBRAIC[32] = ALGEBRAIC[10]+ALGEBRAIC[27];
ALGEBRAIC[37] =  CONSTANTS[161]*(exp(( CONSTANTS[162]*STATES[58])/CONSTANTS[148]) - 1.00000);
ALGEBRAIC[53] =  (( CONSTANTS[161]*ALGEBRAIC[43])/CONSTANTS[148])*ALGEBRAIC[44];
ALGEBRAIC[57] = ALGEBRAIC[37]+ALGEBRAIC[53];
ALGEBRAIC[59] = ALGEBRAIC[32]+ALGEBRAIC[57];
ALGEBRAIC[256] =  CONSTANTS[59]*(exp(ALGEBRAIC[254]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[249]+ALGEBRAIC[248])/( CONSTANTS[141]*CONSTANTS[142])));
}

void objfunc_0(double *p, double *hx, int m, int n, void *adata)
{
  struct rootfind_info* rfi = (struct rootfind_info*)adata;
#define VOI rfi->aVOI
#define CONSTANTS rfi->aCONSTANTS
#define RATES rfi->aRATES
#define STATES rfi->aSTATES
#define ALGEBRAIC rfi->aALGEBRAIC
#define pret rfi->aPRET
  ALGEBRAIC[43] = p[0];
  ALGEBRAIC[44] = p[1];
  ALGEBRAIC[45] = p[2];
  ALGEBRAIC[46] = p[3];
  ALGEBRAIC[47] = p[4];
  hx[0] = ALGEBRAIC[43] - (ALGEBRAIC[44]>0.00000 ? CONSTANTS[158] : CONSTANTS[159]);
  hx[1] = ALGEBRAIC[44] - ALGEBRAIC[45];
  hx[2] = ALGEBRAIC[46] -  (CONSTANTS[160]/ALGEBRAIC[43])*STATES[61];
  hx[3] = ALGEBRAIC[47] -  (CONSTANTS[160]/ALGEBRAIC[43])*(STATES[58] - CONSTANTS[148]);
  hx[4] = ALGEBRAIC[45] - (- ALGEBRAIC[46]+ALGEBRAIC[47]);
#undef VOI
#undef CONSTANTS
#undef RATES
#undef STATES
#undef ALGEBRAIC
#undef pret
}
void rootfind_0(double VOI, double* CONSTANTS, double* RATES,
double* STATES, double* ALGEBRAIC, int* pret)
{
  static double p[5] = {0.1,0.1,0.1,0.1,0.1};
  double bp[5], work[LM_DIF_WORKSZ(5, 5)];
  struct rootfind_info rfi;
  rfi.aVOI = VOI;
  rfi.aCONSTANTS = CONSTANTS;
  rfi.aRATES = RATES;
  rfi.aSTATES = STATES;
  rfi.aALGEBRAIC = ALGEBRAIC;
  rfi.aPRET = pret;
  do_levmar(objfunc_0, p, bp, work, pret, 5, &rfi);
  ALGEBRAIC[43] = p[0];
  ALGEBRAIC[44] = p[1];
  ALGEBRAIC[45] = p[2];
  ALGEBRAIC[46] = p[3];
  ALGEBRAIC[47] = p[4];
}