/* 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]; }