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