/* There are a total of 295 entries in the algebraic variable array. There are a total of 70 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[14] is v_leak in component Ca_leak (fmol_per_sec). * ALGEBRAIC[22] is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[23] is v_R_BU in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[24] is v_R_UW in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[25] is v_R_WS in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[26] is v_R_SU in component crossbridge_TRPN (fmol_per_sec). * STATES[62] is T_total in component crossbridge_TRPN (kPa). * STATES[58] is SL in component crossbridge_TRPN (metre). * ALGEBRAIC[40] is v_diadDf in component diadCaDiffusion (fmol_per_sec). * ALGEBRAIC[262] is v_Ca_D_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[261] is v_Ca_o_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[260] is I_mem_LCC in component LCC (fA). * ALGEBRAIC[281] is v_Ca_o_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[282] is v_Ca_i_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[283] is v_Na_o_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[284] is v_Na_i_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[293] is I_mem_NCX in component NCX (fA). * ALGEBRAIC[140] is v_pCa_R1 in component pCa (fmol_per_sec). * ALGEBRAIC[141] is v_pCa_R2 in component pCa (fmol_per_sec). * ALGEBRAIC[142] is I_mem_pCa in component pCa (fA). * ALGEBRAIC[151] is v_PLBph1 in component PLB (fmol_per_sec). * ALGEBRAIC[152] is v_PLBph2 in component PLB (fmol_per_sec). * ALGEBRAIC[153] is v_PLBd1 in component PLB (fmol_per_sec). * ALGEBRAIC[154] is v_PLBd2 in component PLB (fmol_per_sec). * ALGEBRAIC[155] is v_Inh in component PLB (fmol_per_sec). * ALGEBRAIC[291] is v_RyR in component RyR (fmol_per_sec). * ALGEBRAIC[292] is v_RyRgate_Ca_D in component RyR (fmol_per_sec). * ALGEBRAIC[219] is v_Ca_i_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[220] is v_Ca_SR_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[221] is v_H_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[222] is v_MgADP_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[223] is v_MgATP_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[224] is v_Pi_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[0] is V_m in component environment (volt). * ALGEBRAIC[2] is I_stim in component environment (fA). * CONSTANTS[135] is stimPeriod in component environment (second). * CONSTANTS[136] is stimDuration in component environment (second). * ALGEBRAIC[1] is tPeriod in component environment (second). * ALGEBRAIC[3] is v_Ca_stim in component environment (fmol_per_sec). * ALGEBRAIC[294] is sum_I in component environment (fA). * ALGEBRAIC[4] 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[5] is cCai in component environment (mM). * ALGEBRAIC[6] is cCao in component environment (mM). * ALGEBRAIC[8] is cCaSR in component environment (mM). * ALGEBRAIC[7] is cCaD in component environment (mM). * ALGEBRAIC[9] is cNai in component environment (mM). * ALGEBRAIC[10] is pHi in component environment (dimensionless). * ALGEBRAIC[11] 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[12] is mu_Ca_SR in component Ca_leak (J_per_mol). * ALGEBRAIC[13] is mu_Ca_i in component Ca_leak (J_per_mol). * ALGEBRAIC[15] is mu_TRPN in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[16] is mu_Ca_i in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[17] is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[18] is mu_B_XB in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[19] is mu_U_XB in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[20] is mu_W_XB in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[21] 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[32] is nu_1 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[28] is nu_2 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[29] is nu_3 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[33] is nu_4 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[30] is nu_5 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[31] is nu_6 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[234] is v_1 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[27] is v_2 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[235] is v_4 in component crossbridge_TRPN (m_per_s). * CONSTANTS[156] is v_to_mu in component crossbridge_TRPN (Js_per_m2). * ALGEBRAIC[36] is T_active in component crossbridge_TRPN (kPa). * CONSTANTS[157] is T_ref in component crossbridge_TRPN (kPa). * ALGEBRAIC[34] is F_T_S in component crossbridge_TRPN (kPa). * ALGEBRAIC[35] is F_T_W in component crossbridge_TRPN (kPa). * STATES[61] is Cdd in component crossbridge_TRPN (metre). * STATES[63] 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). * STATES[64] is v_Cdd in component crossbridge_TRPN (m_per_s). * STATES[65] is nu_d in component crossbridge_TRPN (m_per_s). * STATES[66] is nu_k in component crossbridge_TRPN (m_per_s). * STATES[67] is nu_f in component crossbridge_TRPN (m_per_s). * STATES[68] is F_d in component crossbridge_TRPN (kPa). * ALGEBRAIC[37] is F_1 in component crossbridge_TRPN (kPa). * STATES[69] is T_passive in component crossbridge_TRPN (kPa). * CONSTANTS[162] is b in component crossbridge_TRPN (dimensionless). * ALGEBRAIC[38] is mu_Ca_D in component diadCaDiffusion (J_per_mol). * ALGEBRAIC[39] is mu_Ca_i in component diadCaDiffusion (J_per_mol). * ALGEBRAIC[42] is mu_Ca_o in component LCC (J_per_mol). * ALGEBRAIC[43] is mu_Ca_D in component LCC (J_per_mol). * ALGEBRAIC[41] is V_mem in component LCC (volt). * ALGEBRAIC[56] is Af_LCC_Ca1 in component LCC (J_per_mol). * ALGEBRAIC[84] is Ar_LCC_Ca1 in component LCC (J_per_mol). * ALGEBRAIC[82] is Am_LCC_Ca1 in component LCC (J_per_mol). * ALGEBRAIC[110] is v_LCC_Ca1 in component LCC (fmol_per_sec). * ALGEBRAIC[57] is Af_LCC_Ca2 in component LCC (J_per_mol). * ALGEBRAIC[85] is Ar_LCC_Ca2 in component LCC (J_per_mol). * ALGEBRAIC[83] is Am_LCC_Ca2 in component LCC (J_per_mol). * ALGEBRAIC[111] is v_LCC_Ca2 in component LCC (fmol_per_sec). * ALGEBRAIC[76] is Af_fCa00 in component LCC (J_per_mol). * ALGEBRAIC[104] is Ar_fCa00 in component LCC (J_per_mol). * ALGEBRAIC[254] is v_fCa00 in component LCC (fmol_per_sec). * ALGEBRAIC[77] is Af_fCa01 in component LCC (J_per_mol). * ALGEBRAIC[105] is Ar_fCa01 in component LCC (J_per_mol). * ALGEBRAIC[255] is v_fCa01 in component LCC (fmol_per_sec). * ALGEBRAIC[78] is Af_fCa02 in component LCC (J_per_mol). * ALGEBRAIC[106] is Ar_fCa02 in component LCC (J_per_mol). * ALGEBRAIC[256] is v_fCa02 in component LCC (fmol_per_sec). * ALGEBRAIC[79] is Af_fCa10 in component LCC (J_per_mol). * ALGEBRAIC[107] is Ar_fCa10 in component LCC (J_per_mol). * ALGEBRAIC[257] is v_fCa10 in component LCC (fmol_per_sec). * ALGEBRAIC[80] is Af_fCa11 in component LCC (J_per_mol). * ALGEBRAIC[108] is Ar_fCa11 in component LCC (J_per_mol). * ALGEBRAIC[258] is v_fCa11 in component LCC (fmol_per_sec). * ALGEBRAIC[81] is Af_fCa12 in component LCC (J_per_mol). * ALGEBRAIC[109] is Ar_fCa12 in component LCC (J_per_mol). * ALGEBRAIC[259] is v_fCa12 in component LCC (fmol_per_sec). * ALGEBRAIC[44] is mu_S000_LCC in component LCC (J_per_mol). * ALGEBRAIC[263] is v_S000_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[46] is mu_S010_LCC in component LCC (J_per_mol). * ALGEBRAIC[265] is v_S010_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[48] is mu_S020_LCC in component LCC (J_per_mol). * ALGEBRAIC[267] is v_S020_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[50] is mu_S100_LCC in component LCC (J_per_mol). * ALGEBRAIC[269] is v_S100_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[52] is mu_S110_LCC in component LCC (J_per_mol). * ALGEBRAIC[271] is v_S110_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[54] is mu_S120_LCC in component LCC (J_per_mol). * ALGEBRAIC[273] is v_S120_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[45] is mu_S001_LCC in component LCC (J_per_mol). * ALGEBRAIC[264] is v_S001_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[47] is mu_S011_LCC in component LCC (J_per_mol). * ALGEBRAIC[266] is v_S011_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[49] is mu_S021_LCC in component LCC (J_per_mol). * ALGEBRAIC[268] is v_S021_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[51] is mu_S101_LCC in component LCC (J_per_mol). * ALGEBRAIC[270] is v_S101_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[53] is mu_S111_LCC in component LCC (J_per_mol). * ALGEBRAIC[272] is v_S111_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[55] is mu_S121_LCC in component LCC (J_per_mol). * ALGEBRAIC[274] is v_S121_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[58] is Af_d000 in component LCC (J_per_mol). * ALGEBRAIC[86] is Ar_d000 in component LCC (J_per_mol). * ALGEBRAIC[236] is v_d000 in component LCC (fmol_per_sec). * ALGEBRAIC[60] is Af_d010 in component LCC (J_per_mol). * ALGEBRAIC[88] is Ar_d010 in component LCC (J_per_mol). * ALGEBRAIC[238] is v_d010 in component LCC (fmol_per_sec). * ALGEBRAIC[62] is Af_d020 in component LCC (J_per_mol). * ALGEBRAIC[90] is Ar_d020 in component LCC (J_per_mol). * ALGEBRAIC[240] is v_d020 in component LCC (fmol_per_sec). * ALGEBRAIC[64] is Af_f1_000 in component LCC (J_per_mol). * ALGEBRAIC[92] is Ar_f1_000 in component LCC (J_per_mol). * ALGEBRAIC[242] is v_f1_000 in component LCC (fmol_per_sec). * ALGEBRAIC[68] is Af_f2_000 in component LCC (J_per_mol). * ALGEBRAIC[96] is Ar_f2_000 in component LCC (J_per_mol). * ALGEBRAIC[246] is v_f2_000 in component LCC (fmol_per_sec). * ALGEBRAIC[72] is Af_f3_010 in component LCC (J_per_mol). * ALGEBRAIC[100] is Ar_f3_010 in component LCC (J_per_mol). * ALGEBRAIC[250] is v_f3_010 in component LCC (fmol_per_sec). * ALGEBRAIC[66] is Af_f1_100 in component LCC (J_per_mol). * ALGEBRAIC[94] is Ar_f1_100 in component LCC (J_per_mol). * ALGEBRAIC[244] is v_f1_100 in component LCC (fmol_per_sec). * ALGEBRAIC[70] is Af_f2_100 in component LCC (J_per_mol). * ALGEBRAIC[98] is Ar_f2_100 in component LCC (J_per_mol). * ALGEBRAIC[248] is v_f2_100 in component LCC (fmol_per_sec). * ALGEBRAIC[74] is Af_f3_110 in component LCC (J_per_mol). * ALGEBRAIC[102] is Ar_f3_110 in component LCC (J_per_mol). * ALGEBRAIC[252] is v_f3_110 in component LCC (fmol_per_sec). * ALGEBRAIC[59] is Af_d001 in component LCC (J_per_mol). * ALGEBRAIC[87] is Ar_d001 in component LCC (J_per_mol). * ALGEBRAIC[237] is v_d001 in component LCC (fmol_per_sec). * ALGEBRAIC[61] is Af_d011 in component LCC (J_per_mol). * ALGEBRAIC[89] is Ar_d011 in component LCC (J_per_mol). * ALGEBRAIC[239] is v_d011 in component LCC (fmol_per_sec). * ALGEBRAIC[63] is Af_d021 in component LCC (J_per_mol). * ALGEBRAIC[91] is Ar_d021 in component LCC (J_per_mol). * ALGEBRAIC[241] is v_d021 in component LCC (fmol_per_sec). * ALGEBRAIC[65] is Af_f1_001 in component LCC (J_per_mol). * ALGEBRAIC[93] is Ar_f1_001 in component LCC (J_per_mol). * ALGEBRAIC[243] is v_f1_001 in component LCC (fmol_per_sec). * ALGEBRAIC[69] is Af_f2_001 in component LCC (J_per_mol). * ALGEBRAIC[97] is Ar_f2_001 in component LCC (J_per_mol). * ALGEBRAIC[247] is v_f2_001 in component LCC (fmol_per_sec). * ALGEBRAIC[73] is Af_f3_011 in component LCC (J_per_mol). * ALGEBRAIC[101] is Ar_f3_011 in component LCC (J_per_mol). * ALGEBRAIC[251] is v_f3_011 in component LCC (fmol_per_sec). * ALGEBRAIC[67] is Af_f1_101 in component LCC (J_per_mol). * ALGEBRAIC[95] is Ar_f1_101 in component LCC (J_per_mol). * ALGEBRAIC[245] is v_f1_101 in component LCC (fmol_per_sec). * ALGEBRAIC[71] is Af_f2_101 in component LCC (J_per_mol). * ALGEBRAIC[99] is Ar_f2_101 in component LCC (J_per_mol). * ALGEBRAIC[249] is v_f2_101 in component LCC (fmol_per_sec). * ALGEBRAIC[75] is Af_f3_111 in component LCC (J_per_mol). * ALGEBRAIC[103] is Ar_f3_111 in component LCC (J_per_mol). * ALGEBRAIC[253] is v_f3_111 in component LCC (fmol_per_sec). * ALGEBRAIC[112] is Ca_tot in component LCC (fmol). * ALGEBRAIC[114] is mu_Ca_o in component NCX (J_per_mol). * ALGEBRAIC[115] is mu_Ca_i in component NCX (J_per_mol). * ALGEBRAIC[116] is mu_Na_o in component NCX (J_per_mol). * ALGEBRAIC[117] is mu_Na_i in component NCX (J_per_mol). * ALGEBRAIC[118] is mu_P1_NCX in component NCX (J_per_mol). * ALGEBRAIC[119] is mu_P2_NCX in component NCX (J_per_mol). * ALGEBRAIC[120] is mu_P3_NCX in component NCX (J_per_mol). * ALGEBRAIC[121] is mu_P4_NCX in component NCX (J_per_mol). * ALGEBRAIC[122] is mu_P5_NCX in component NCX (J_per_mol). * ALGEBRAIC[123] is mu_P6_NCX in component NCX (J_per_mol). * ALGEBRAIC[285] is v_P1_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[286] is v_P2_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[287] is v_P3_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[288] is v_P4_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[289] is v_P5_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[290] is v_P6_NCX in component NCX (fmol_per_sec). * ALGEBRAIC[124] is Af_r1 in component NCX (J_per_mol). * ALGEBRAIC[130] is Ar_r1 in component NCX (J_per_mol). * ALGEBRAIC[275] is v_r1 in component NCX (fmol_per_sec). * ALGEBRAIC[125] is Af_r2 in component NCX (J_per_mol). * ALGEBRAIC[131] is Ar_r2 in component NCX (J_per_mol). * ALGEBRAIC[276] is v_r2 in component NCX (fmol_per_sec). * ALGEBRAIC[126] is Af_r3 in component NCX (J_per_mol). * ALGEBRAIC[132] is Ar_r3 in component NCX (J_per_mol). * ALGEBRAIC[277] is v_r3 in component NCX (fmol_per_sec). * ALGEBRAIC[127] is Af_r4 in component NCX (J_per_mol). * ALGEBRAIC[133] is Ar_r4 in component NCX (J_per_mol). * ALGEBRAIC[278] is v_r4 in component NCX (fmol_per_sec). * ALGEBRAIC[128] is Af_r5 in component NCX (J_per_mol). * ALGEBRAIC[134] is Ar_r5 in component NCX (J_per_mol). * ALGEBRAIC[279] is v_r5 in component NCX (fmol_per_sec). * ALGEBRAIC[129] is Af_r6 in component NCX (J_per_mol). * ALGEBRAIC[135] is Ar_r6 in component NCX (J_per_mol). * ALGEBRAIC[280] is v_r6 in component NCX (fmol_per_sec). * ALGEBRAIC[113] is V_mem in component NCX (volt). * ALGEBRAIC[136] is mu_Ca_i in component pCa (J_per_mol). * ALGEBRAIC[137] is mu_Ca_o in component pCa (J_per_mol). * ALGEBRAIC[138] is mu_pCa in component pCa (J_per_mol). * ALGEBRAIC[139] is mu_pCa_Ca in component pCa (J_per_mol). * ALGEBRAIC[143] is mu_PLB in component PLB (J_per_mol). * ALGEBRAIC[144] is mu_PKACI in component PLB (J_per_mol). * ALGEBRAIC[145] is mu_PLB_PKACI in component PLB (J_per_mol). * ALGEBRAIC[146] is mu_PP1 in component PLB (J_per_mol). * ALGEBRAIC[147] is mu_PLBp_PP1 in component PLB (J_per_mol). * ALGEBRAIC[148] is mu_PLBp in component PLB (J_per_mol). * ALGEBRAIC[149] is mu_Ip in component PLB (J_per_mol). * ALGEBRAIC[150] 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[156] is mu_Ca_SR in component RyR (J_per_mol). * ALGEBRAIC[157] is mu_Ca_D in component RyR (J_per_mol). * ALGEBRAIC[158] is q_Ca_gate_complexes in component RyR (fmol). * ALGEBRAIC[159] is mu_C_RyR in component RyR (J_per_mol). * ALGEBRAIC[160] is mu_CI_RyR in component RyR (J_per_mol). * ALGEBRAIC[161] is mu_I_RyR in component RyR (J_per_mol). * ALGEBRAIC[162] is mu_O_RyR in component RyR (J_per_mol). * ALGEBRAIC[163] is v_OC in component RyR (fmol_per_sec). * ALGEBRAIC[164] is v_CCI in component RyR (fmol_per_sec). * ALGEBRAIC[165] is v_CII in component RyR (fmol_per_sec). * ALGEBRAIC[166] 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[167] is q_SERCA_Ca_complexes in component SERCA (fmol). * ALGEBRAIC[168] is c_Ca_i in component SERCA (mM). * ALGEBRAIC[169] is c_Ca_SR in component SERCA (mM). * ALGEBRAIC[170] is c_H in component SERCA (mM). * ALGEBRAIC[171] is c_MgADP in component SERCA (mM). * ALGEBRAIC[172] is c_MgATP in component SERCA (mM). * ALGEBRAIC[173] 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[174] is mu_Ca_i in component SERCA (J_per_mol). * ALGEBRAIC[175] is mu_Ca_SR in component SERCA (J_per_mol). * ALGEBRAIC[176] is mu_H in component SERCA (J_per_mol). * ALGEBRAIC[177] is mu_MgADP in component SERCA (J_per_mol). * ALGEBRAIC[178] is mu_MgATP in component SERCA (J_per_mol). * ALGEBRAIC[180] is mu_P1 in component SERCA (J_per_mol). * ALGEBRAIC[225] is v_P1 in component SERCA (fmol_per_sec). * ALGEBRAIC[179] is mu_Pi in component SERCA (J_per_mol). * ALGEBRAIC[181] is mu_P2 in component SERCA (J_per_mol). * ALGEBRAIC[226] is v_P2 in component SERCA (fmol_per_sec). * ALGEBRAIC[182] is mu_P2a in component SERCA (J_per_mol). * ALGEBRAIC[227] is v_P2a in component SERCA (fmol_per_sec). * ALGEBRAIC[183] is mu_P4 in component SERCA (J_per_mol). * ALGEBRAIC[228] is v_P4 in component SERCA (fmol_per_sec). * ALGEBRAIC[184] is mu_P5 in component SERCA (J_per_mol). * ALGEBRAIC[229] is v_P5 in component SERCA (fmol_per_sec). * ALGEBRAIC[186] is mu_P6 in component SERCA (J_per_mol). * ALGEBRAIC[231] is v_P6 in component SERCA (fmol_per_sec). * ALGEBRAIC[187] is mu_P8 in component SERCA (J_per_mol). * ALGEBRAIC[232] is v_P8 in component SERCA (fmol_per_sec). * ALGEBRAIC[188] is mu_P9 in component SERCA (J_per_mol). * ALGEBRAIC[233] is v_P9 in component SERCA (fmol_per_sec). * ALGEBRAIC[185] is mu_P10 in component SERCA (J_per_mol). * ALGEBRAIC[230] is v_P10 in component SERCA (fmol_per_sec). * ALGEBRAIC[189] is mu_PLBp_baseline in component SERCA (J_per_mol). * CONSTANTS[173] is mu_PLBp in component SERCA (J_per_mol). * ALGEBRAIC[190] is mu_PLBp_eff in component SERCA (J_per_mol). * ALGEBRAIC[193] is Af_R1_2 in component SERCA (J_per_mol). * ALGEBRAIC[194] is Ar_R1_2 in component SERCA (J_per_mol). * ALGEBRAIC[211] is v_SERCA_R1_2 in component SERCA (fmol_per_sec). * ALGEBRAIC[195] is Af_R5_6 in component SERCA (J_per_mol). * ALGEBRAIC[196] is Ar_R5_6 in component SERCA (J_per_mol). * ALGEBRAIC[212] is v_SERCA_R5_6 in component SERCA (fmol_per_sec). * ALGEBRAIC[197] is Af_R2_2a in component SERCA (J_per_mol). * ALGEBRAIC[198] is Ar_R2_2a in component SERCA (J_per_mol). * ALGEBRAIC[213] is v_SERCA_R2_2a in component SERCA (fmol_per_sec). * ALGEBRAIC[199] is Af_R2_4 in component SERCA (J_per_mol). * ALGEBRAIC[200] is Ar_R2_4 in component SERCA (J_per_mol). * ALGEBRAIC[214] is v_SERCA_R2_4 in component SERCA (fmol_per_sec). * ALGEBRAIC[201] is Af_R4_5 in component SERCA (J_per_mol). * ALGEBRAIC[202] is Ar_R4_5 in component SERCA (J_per_mol). * ALGEBRAIC[215] is v_SERCA_R4_5 in component SERCA (fmol_per_sec). * ALGEBRAIC[203] is Af_R6_8 in component SERCA (J_per_mol). * ALGEBRAIC[204] is Ar_R6_8 in component SERCA (J_per_mol). * ALGEBRAIC[216] is v_SERCA_R6_8 in component SERCA (fmol_per_sec). * ALGEBRAIC[205] is Af_R8_9 in component SERCA (J_per_mol). * ALGEBRAIC[206] is Ar_R8_9 in component SERCA (J_per_mol). * ALGEBRAIC[217] is v_SERCA_R8_9 in component SERCA (fmol_per_sec). * ALGEBRAIC[207] is Af_R9_10 in component SERCA (J_per_mol). * ALGEBRAIC[208] is Ar_R9_10 in component SERCA (J_per_mol). * ALGEBRAIC[218] is v_SERCA_R9_10 in component SERCA (fmol_per_sec). * ALGEBRAIC[191] is Af_R10_1 in component SERCA (J_per_mol). * ALGEBRAIC[192] is Ar_R10_1 in component SERCA (J_per_mol). * ALGEBRAIC[209] is v_SERCA_R10_1 in component SERCA (fmol_per_sec). * ALGEBRAIC[210] 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). * There are a total of 8 condition variables. */ 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]; STATES[62] = 0.1001; STATES[63] = 0.1001; STATES[64] = 0.1001; STATES[65] = 0.1001; STATES[66] = 0.1001; STATES[67] = 0.1001; STATES[68] = 0.1001; STATES[69] = 0.1001; RATES[20] = 0.1001; RATES[4] = 0.1001; RATES[2] = 0.1001; RATES[3] = 0.1001; RATES[7] = 0.1001; RATES[8] = 0.1001; RATES[9] = 0.1001; RATES[10] = 0.1001; RATES[1] = 0.1001; RATES[11] = 0.1001; RATES[5] = 0.1001; RATES[6] = 0.1001; RATES[21] = 0.1001; RATES[22] = 0.1001; RATES[12] = 0.1001; RATES[13] = 0.1001; RATES[14] = 0.1001; RATES[15] = 0.1001; RATES[16] = 0.1001; RATES[17] = 0.1001; RATES[18] = 0.1001; RATES[19] = 0.1001; RATES[23] = 0.1001; RATES[25] = 0.1001; RATES[26] = 0.1001; RATES[24] = 0.1001; RATES[0] = 0.1001; RATES[59] = 0.1001; RATES[58] = 0.1001; RATES[60] = 0.1001; RATES[61] = 0.1001; RATES[27] = 0.1001; RATES[28] = 0.1001; RATES[29] = 0.1001; RATES[30] = 0.1001; RATES[31] = 0.1001; RATES[32] = 0.1001; RATES[33] = 0.1001; RATES[34] = 0.1001; RATES[35] = 0.1001; RATES[36] = 0.1001; RATES[37] = 0.1001; RATES[38] = 0.1001; RATES[39] = 0.1001; RATES[40] = 0.1001; RATES[41] = 0.1001; RATES[42] = 0.1001; RATES[43] = 0.1001; RATES[44] = 0.1001; RATES[48] = 0.1001; RATES[45] = 0.1001; RATES[46] = 0.1001; RATES[47] = 0.1001; RATES[49] = 0.1001; RATES[50] = 0.1001; RATES[51] = 0.1001; RATES[52] = 0.1001; RATES[53] = 0.1001; RATES[57] = 0.1001; RATES[54] = 0.1001; RATES[55] = 0.1001; RATES[56] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[20] - (- ALGEBRAIC[14] - ALGEBRAIC[291])+ALGEBRAIC[220]; resid[1] = RATES[4] - (( - ALGEBRAIC[22]*CONSTANTS[131]+ALGEBRAIC[219]+ALGEBRAIC[14]+ALGEBRAIC[40]+ALGEBRAIC[282]) - ALGEBRAIC[140])+ALGEBRAIC[3]; resid[2] = RATES[2] - - ALGEBRAIC[22]; resid[3] = RATES[3] - ALGEBRAIC[22]; resid[4] = RATES[7] - - ALGEBRAIC[23]; resid[5] = RATES[8] - (ALGEBRAIC[23] - ALGEBRAIC[24])+ALGEBRAIC[26]; resid[6] = RATES[9] - ALGEBRAIC[24] - ALGEBRAIC[25]; resid[7] = RATES[10] - ALGEBRAIC[25] - ALGEBRAIC[26]; resid[8] = RATES[1] - - ALGEBRAIC[40]+ALGEBRAIC[262]+ALGEBRAIC[291]+ALGEBRAIC[292]; resid[9] = RATES[11] - ALGEBRAIC[261]+ALGEBRAIC[281]+ALGEBRAIC[141]; resid[10] = RATES[5] - ALGEBRAIC[284]; resid[11] = RATES[6] - ALGEBRAIC[283]; resid[12] = RATES[21] - - ALGEBRAIC[140]+ALGEBRAIC[141]; resid[13] = RATES[22] - ALGEBRAIC[140] - ALGEBRAIC[141]; resid[14] = RATES[12] - ALGEBRAIC[154] - ALGEBRAIC[151]; resid[15] = RATES[13] - ALGEBRAIC[152] - ALGEBRAIC[151]; resid[16] = RATES[14] - ALGEBRAIC[151] - ALGEBRAIC[152]; resid[17] = RATES[15] - ALGEBRAIC[151] - ALGEBRAIC[152]; resid[18] = RATES[16] - ALGEBRAIC[153] - ALGEBRAIC[154]; resid[19] = RATES[17] - ALGEBRAIC[152] - ALGEBRAIC[153]; resid[20] = RATES[18] - - ALGEBRAIC[155]; resid[21] = RATES[19] - ALGEBRAIC[155]; resid[22] = RATES[23] - ALGEBRAIC[221]; resid[23] = RATES[25] - ALGEBRAIC[222]; resid[24] = RATES[26] - ALGEBRAIC[223]; resid[25] = RATES[24] - ALGEBRAIC[224]; resid[26] = RATES[0] - ALGEBRAIC[294]; resid[27] = RATES[59] - ALGEBRAIC[234]; resid[28] = RATES[58] - ALGEBRAIC[27]; resid[29] = RATES[60] - ALGEBRAIC[235]; resid[30] = STATES[63] - (CONDVAR[7]>0.00000 ? CONSTANTS[158] : CONSTANTS[159]); resid[31] = RATES[61] - STATES[64]; resid[32] = STATES[64] - STATES[65]; resid[33] = STATES[66] - (CONSTANTS[160]/STATES[63])*STATES[61]; resid[34] = STATES[67] - (CONSTANTS[160]/STATES[63])*(STATES[58] - CONSTANTS[148]); resid[35] = STATES[65] - - STATES[66]+STATES[67]; resid[36] = STATES[68] - (( CONSTANTS[161]*STATES[63])/CONSTANTS[148])*STATES[64]; resid[37] = STATES[69] - ALGEBRAIC[37]+STATES[68]; resid[38] = STATES[62] - ALGEBRAIC[36]+STATES[69]; resid[39] = RATES[27] - ALGEBRAIC[263]; resid[40] = RATES[28] - ALGEBRAIC[265]; resid[41] = RATES[29] - ALGEBRAIC[267]; resid[42] = RATES[30] - ALGEBRAIC[269]; resid[43] = RATES[31] - ALGEBRAIC[271]; resid[44] = RATES[32] - ALGEBRAIC[273]; resid[45] = RATES[33] - ALGEBRAIC[264]; resid[46] = RATES[34] - ALGEBRAIC[266]; resid[47] = RATES[35] - ALGEBRAIC[268]; resid[48] = RATES[36] - ALGEBRAIC[270]; resid[49] = RATES[37] - ALGEBRAIC[272]; resid[50] = RATES[38] - ALGEBRAIC[274]; resid[51] = RATES[39] - ALGEBRAIC[285]; resid[52] = RATES[40] - ALGEBRAIC[286]; resid[53] = RATES[41] - ALGEBRAIC[287]; resid[54] = RATES[42] - ALGEBRAIC[288]; resid[55] = RATES[43] - ALGEBRAIC[289]; resid[56] = RATES[44] - ALGEBRAIC[290]; resid[57] = RATES[48] - ALGEBRAIC[166] - ALGEBRAIC[163]; resid[58] = RATES[45] - ALGEBRAIC[163] - ALGEBRAIC[164]; resid[59] = RATES[46] - ALGEBRAIC[164] - ALGEBRAIC[165]; resid[60] = RATES[47] - ALGEBRAIC[165] - ALGEBRAIC[166]; resid[61] = RATES[49] - ALGEBRAIC[225]; resid[62] = RATES[50] - ALGEBRAIC[226]; resid[63] = RATES[51] - ALGEBRAIC[227]; resid[64] = RATES[52] - ALGEBRAIC[228]; resid[65] = RATES[53] - ALGEBRAIC[229]; resid[66] = RATES[57] - ALGEBRAIC[230]; resid[67] = RATES[54] - ALGEBRAIC[231]; resid[68] = RATES[55] - ALGEBRAIC[232]; resid[69] = RATES[56] - ALGEBRAIC[233]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = STATES[0]/CONSTANTS[133]; ALGEBRAIC[5] = STATES[4]/CONSTANTS[132]; ALGEBRAIC[6] = STATES[11]/CONSTANTS[139]; ALGEBRAIC[7] = STATES[1]/CONSTANTS[140]; ALGEBRAIC[8] = STATES[20]/CONSTANTS[170]; ALGEBRAIC[9] = STATES[5]/CONSTANTS[132]; ALGEBRAIC[10] = - arbitrary_log(STATES[23]/CONSTANTS[132], 10); ALGEBRAIC[11] = 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[45] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33]); ALGEBRAIC[47] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34]); ALGEBRAIC[49] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35]); ALGEBRAIC[51] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36]); ALGEBRAIC[112] = STATES[11]+STATES[1]+STATES[33]+STATES[34]+STATES[35]+STATES[36]+STATES[37]+STATES[38]; ALGEBRAIC[158] = CONSTANTS[164]*STATES[48]+ CONSTANTS[163]*STATES[46]+ (CONSTANTS[163]+CONSTANTS[164])*STATES[47]; ALGEBRAIC[167] = CONSTANTS[165]*STATES[52]+ CONSTANTS[165]*STATES[53]+ CONSTANTS[166]*STATES[54]; ALGEBRAIC[168] = STATES[4]/CONSTANTS[168]; ALGEBRAIC[169] = STATES[20]/CONSTANTS[174]; ALGEBRAIC[170] = STATES[23]/CONSTANTS[176]; ALGEBRAIC[171] = STATES[25]/CONSTANTS[168]; ALGEBRAIC[172] = STATES[26]/CONSTANTS[168]; ALGEBRAIC[173] = STATES[24]/CONSTANTS[168]; ALGEBRAIC[210] = CONSTANTS[59]*(exp(ALGEBRAIC[185]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[180]+ALGEBRAIC[179])/( CONSTANTS[141]*CONSTANTS[142]))); } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[1] = VOI - floor(VOI/CONSTANTS[135])*CONSTANTS[135]; ALGEBRAIC[3] = (CONDVAR[2]>=0.00000&&CONDVAR[3]<=0.00000 ? 0.00000*8000.00 : 0.00000); ALGEBRAIC[12] = ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[60]*STATES[20]); ALGEBRAIC[13] = ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[61]*STATES[4]); ALGEBRAIC[14] = CONSTANTS[0]*(exp(ALGEBRAIC[12]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[13]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[15] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[62]*STATES[2]); ALGEBRAIC[16] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]); ALGEBRAIC[17] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[63]*STATES[3]); ALGEBRAIC[22] = CONSTANTS[1]*(exp(( CONSTANTS[131]*ALGEBRAIC[16]+ALGEBRAIC[15])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[17]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[18] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[64]*STATES[7]); ALGEBRAIC[19] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[65]*STATES[8]); ALGEBRAIC[23] = CONSTANTS[2]*(exp((ALGEBRAIC[18]+ CONSTANTS[130]*ALGEBRAIC[17])/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[19]+ CONSTANTS[130]*ALGEBRAIC[17])/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[20] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[66]*STATES[9]); ALGEBRAIC[24] = CONSTANTS[3]*(exp(ALGEBRAIC[19]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[20]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[21] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[67]*STATES[10]); ALGEBRAIC[25] = CONSTANTS[4]*(exp(ALGEBRAIC[20]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[21]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[26] = CONSTANTS[5]*(exp(ALGEBRAIC[21]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[19]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[4] = (CONDVAR[4]<0.00000 ? 1.00000*CONSTANTS[137] : CONDVAR[5]>0.00000&&CONDVAR[6]<0.00000 ? - 0.100000*CONSTANTS[137] : 0.00000*CONSTANTS[137]); ALGEBRAIC[27] = ALGEBRAIC[4]; ALGEBRAIC[34] = (CONSTANTS[157]/( CONSTANTS[148]*CONSTANTS[149]*CONSTANTS[150]))*STATES[10]*(STATES[60]+CONSTANTS[148]); ALGEBRAIC[35] = (CONSTANTS[157]/( CONSTANTS[148]*CONSTANTS[149]*CONSTANTS[150]))*STATES[9]*STATES[59]; ALGEBRAIC[36] = ALGEBRAIC[34]+ALGEBRAIC[35]; ALGEBRAIC[37] = CONSTANTS[161]*(exp(( CONSTANTS[162]*STATES[58])/CONSTANTS[148]) - 1.00000); ALGEBRAIC[38] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]); ALGEBRAIC[39] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]); ALGEBRAIC[40] = CONSTANTS[6]*(exp(ALGEBRAIC[38]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[39]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[136] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]); ALGEBRAIC[138] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[90]*STATES[21]); ALGEBRAIC[139] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[91]*STATES[22]); ALGEBRAIC[140] = CONSTANTS[39]*(exp((ALGEBRAIC[136]+ALGEBRAIC[138])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[139]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[137] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]); ALGEBRAIC[141] = CONSTANTS[40]*(exp(ALGEBRAIC[139]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[137]+ALGEBRAIC[138])/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[143] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[92]*STATES[12]); ALGEBRAIC[144] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[93]*STATES[13]); ALGEBRAIC[145] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[94]*STATES[14]); ALGEBRAIC[151] = CONSTANTS[41]*(exp((ALGEBRAIC[143]+ALGEBRAIC[144])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[145]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[148] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]); ALGEBRAIC[152] = CONSTANTS[42]*(exp(ALGEBRAIC[145]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[148]+ALGEBRAIC[144])/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[146] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[95]*STATES[15]); ALGEBRAIC[147] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[96]*STATES[16]); ALGEBRAIC[153] = CONSTANTS[43]*(exp((ALGEBRAIC[148]+ALGEBRAIC[146])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[147]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[154] = CONSTANTS[44]*(exp(ALGEBRAIC[147]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[143]+ALGEBRAIC[146])/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[149] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[98]*STATES[18]); ALGEBRAIC[150] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[99]*STATES[19]); ALGEBRAIC[155] = CONSTANTS[45]*(exp((ALGEBRAIC[146]+ALGEBRAIC[149])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[150]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[157] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]); ALGEBRAIC[159] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[100]*STATES[45]); ALGEBRAIC[162] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[103]*STATES[48]); ALGEBRAIC[163] = CONSTANTS[47]*(exp(ALGEBRAIC[162]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[159]+ CONSTANTS[164]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[160] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[101]*STATES[46]); ALGEBRAIC[164] = CONSTANTS[48]*(exp((ALGEBRAIC[159]+ CONSTANTS[163]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[160]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[161] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[102]*STATES[47]); ALGEBRAIC[165] = CONSTANTS[49]*(exp((ALGEBRAIC[160]+ CONSTANTS[164]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[161]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[166] = CONSTANTS[50]*(exp(ALGEBRAIC[161]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[162]+ CONSTANTS[163]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[174] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]); ALGEBRAIC[181] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[105]*STATES[50]); ALGEBRAIC[189] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]); ALGEBRAIC[190] = ALGEBRAIC[189] - CONSTANTS[173]; ALGEBRAIC[199] = ALGEBRAIC[181]+ CONSTANTS[165]*ALGEBRAIC[174]+ALGEBRAIC[190]; ALGEBRAIC[183] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[107]*STATES[52]); ALGEBRAIC[200] = ALGEBRAIC[183]+ALGEBRAIC[190]; ALGEBRAIC[214] = CONSTANTS[52]*(exp(ALGEBRAIC[199]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[200]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[219] = - CONSTANTS[165]*ALGEBRAIC[214]; ALGEBRAIC[186] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[109]*STATES[54]); ALGEBRAIC[203] = ALGEBRAIC[186]+ALGEBRAIC[190]; ALGEBRAIC[175] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]); ALGEBRAIC[187] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[110]*STATES[55]); ALGEBRAIC[204] = ALGEBRAIC[187]+ CONSTANTS[166]*ALGEBRAIC[175]+ALGEBRAIC[190]; ALGEBRAIC[216] = CONSTANTS[56]*(exp(ALGEBRAIC[203]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[204]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[220] = CONSTANTS[166]*ALGEBRAIC[216]; ALGEBRAIC[176] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[113]*STATES[23]); ALGEBRAIC[197] = ALGEBRAIC[181]+ALGEBRAIC[176]+ALGEBRAIC[190]; ALGEBRAIC[182] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[106]*STATES[51]); ALGEBRAIC[198] = ALGEBRAIC[182]+ALGEBRAIC[190]; ALGEBRAIC[213] = CONSTANTS[53]*(exp(ALGEBRAIC[197]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[198]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[201] = ALGEBRAIC[183]+ALGEBRAIC[190]; ALGEBRAIC[184] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[108]*STATES[53]); ALGEBRAIC[202] = ALGEBRAIC[184]+ CONSTANTS[167]*ALGEBRAIC[176]+ALGEBRAIC[190]; ALGEBRAIC[215] = CONSTANTS[54]*(exp(ALGEBRAIC[201]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[202]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[205] = ALGEBRAIC[187]+ CONSTANTS[167]*ALGEBRAIC[176]+ALGEBRAIC[190]; ALGEBRAIC[188] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[111]*STATES[56]); ALGEBRAIC[206] = ALGEBRAIC[188]+ALGEBRAIC[190]; ALGEBRAIC[217] = CONSTANTS[57]*(exp(ALGEBRAIC[205]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[206]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[207] = ALGEBRAIC[188]+ALGEBRAIC[190]; ALGEBRAIC[185] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[112]*STATES[57]); ALGEBRAIC[208] = ALGEBRAIC[176]+ALGEBRAIC[185]+ALGEBRAIC[190]; ALGEBRAIC[218] = CONSTANTS[58]*(exp(ALGEBRAIC[207]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[208]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[221] = (( CONSTANTS[167]*ALGEBRAIC[215] - ALGEBRAIC[213]) - CONSTANTS[167]*ALGEBRAIC[217])+ALGEBRAIC[218]; ALGEBRAIC[195] = ALGEBRAIC[184]+ALGEBRAIC[190]; ALGEBRAIC[177] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[115]*STATES[25]); ALGEBRAIC[196] = ALGEBRAIC[177]+ALGEBRAIC[186]+ALGEBRAIC[190]; ALGEBRAIC[212] = CONSTANTS[55]*(exp(ALGEBRAIC[195]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[196]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[222] = ALGEBRAIC[212]; ALGEBRAIC[178] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[114]*STATES[26]); ALGEBRAIC[180] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[104]*STATES[49]); ALGEBRAIC[193] = ALGEBRAIC[180]+ALGEBRAIC[178]+ALGEBRAIC[190]; ALGEBRAIC[194] = ALGEBRAIC[181]+ALGEBRAIC[190]; ALGEBRAIC[211] = CONSTANTS[51]*(exp(ALGEBRAIC[193]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[194]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[223] = - ALGEBRAIC[211]; ALGEBRAIC[191] = ALGEBRAIC[185]+ALGEBRAIC[190]; ALGEBRAIC[179] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[116]*STATES[24]); ALGEBRAIC[192] = ALGEBRAIC[180]+ALGEBRAIC[179]+ALGEBRAIC[190]; ALGEBRAIC[209] = CONSTANTS[59]*(exp(ALGEBRAIC[191]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[192]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[224] = ALGEBRAIC[209]; ALGEBRAIC[225] = ALGEBRAIC[209] - ALGEBRAIC[211]; ALGEBRAIC[226] = ALGEBRAIC[211] - ALGEBRAIC[214]; ALGEBRAIC[227] = ALGEBRAIC[213]; ALGEBRAIC[228] = ALGEBRAIC[214] - ALGEBRAIC[215]; ALGEBRAIC[229] = ALGEBRAIC[215] - ALGEBRAIC[212]; ALGEBRAIC[230] = ALGEBRAIC[218] - ALGEBRAIC[209]; ALGEBRAIC[231] = ALGEBRAIC[212] - ALGEBRAIC[216]; ALGEBRAIC[232] = ALGEBRAIC[216] - ALGEBRAIC[217]; ALGEBRAIC[233] = ALGEBRAIC[217] - ALGEBRAIC[218]; ALGEBRAIC[28] = - CONSTANTS[169]*ALGEBRAIC[27]; ALGEBRAIC[29] = CONSTANTS[171]*STATES[59]; ALGEBRAIC[32] = - ALGEBRAIC[28] - ALGEBRAIC[29]; ALGEBRAIC[234] = ALGEBRAIC[32]; ALGEBRAIC[30] = - CONSTANTS[175]*ALGEBRAIC[27]; ALGEBRAIC[31] = CONSTANTS[172]*STATES[60]; ALGEBRAIC[33] = - ALGEBRAIC[30] - ALGEBRAIC[31]; ALGEBRAIC[235] = ALGEBRAIC[33]; ALGEBRAIC[43] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]); ALGEBRAIC[41] = STATES[0]/CONSTANTS[133]; ALGEBRAIC[53] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37]); ALGEBRAIC[56] = ALGEBRAIC[43]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41]+ALGEBRAIC[53]; ALGEBRAIC[42] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]); ALGEBRAIC[84] = ALGEBRAIC[42]+ALGEBRAIC[53]; ALGEBRAIC[82] = CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[110] = (ALGEBRAIC[82]==0.00000 ? CONSTANTS[7]*(exp(ALGEBRAIC[56]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[84]/( CONSTANTS[141]*CONSTANTS[142]))) : ((( CONSTANTS[7]*ALGEBRAIC[82])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[82]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[56]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[84]/( CONSTANTS[141]*CONSTANTS[142])))); ALGEBRAIC[55] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38]); ALGEBRAIC[57] = ALGEBRAIC[43]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41]+ALGEBRAIC[55]; ALGEBRAIC[85] = ALGEBRAIC[42]+ALGEBRAIC[55]; ALGEBRAIC[83] = CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[111] = (ALGEBRAIC[83]==0.00000 ? CONSTANTS[8]*(exp(ALGEBRAIC[57]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[85]/( CONSTANTS[141]*CONSTANTS[142]))) : ((( CONSTANTS[8]*ALGEBRAIC[83])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[83]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[57]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[85]/( CONSTANTS[141]*CONSTANTS[142])))); ALGEBRAIC[261] = ALGEBRAIC[111]+ALGEBRAIC[110]; ALGEBRAIC[44] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27]); ALGEBRAIC[76] = ALGEBRAIC[44]; ALGEBRAIC[104] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[117]*ALGEBRAIC[43]; ALGEBRAIC[254] = CONSTANTS[27]*(exp(ALGEBRAIC[76]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[104]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[46] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28]); ALGEBRAIC[77] = ALGEBRAIC[46]; ALGEBRAIC[105] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[117]*ALGEBRAIC[43]; ALGEBRAIC[255] = CONSTANTS[29]*(exp(ALGEBRAIC[77]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[105]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[48] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29]); ALGEBRAIC[78] = ALGEBRAIC[48]; ALGEBRAIC[106] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[117]*ALGEBRAIC[43]; ALGEBRAIC[256] = CONSTANTS[31]*(exp(ALGEBRAIC[78]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[106]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[50] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30]); ALGEBRAIC[79] = ALGEBRAIC[50]; ALGEBRAIC[107] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[117]*ALGEBRAIC[43]; ALGEBRAIC[257] = CONSTANTS[28]*(exp(ALGEBRAIC[79]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[107]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[52] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31]); ALGEBRAIC[80] = ALGEBRAIC[52]; ALGEBRAIC[108] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[117]*ALGEBRAIC[43]; ALGEBRAIC[258] = CONSTANTS[30]*(exp(ALGEBRAIC[80]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[108]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[54] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32]); ALGEBRAIC[81] = ALGEBRAIC[54]; ALGEBRAIC[109] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[117]*ALGEBRAIC[43]; ALGEBRAIC[259] = CONSTANTS[32]*(exp(ALGEBRAIC[81]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[109]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[262] = (( (CONSTANTS[117])*ALGEBRAIC[257]+ CONSTANTS[117]*ALGEBRAIC[258]+ CONSTANTS[117]*ALGEBRAIC[256]+ CONSTANTS[117]*ALGEBRAIC[255]+ CONSTANTS[117]*ALGEBRAIC[254]+ CONSTANTS[117]*ALGEBRAIC[259]) - ALGEBRAIC[110]) - ALGEBRAIC[111]; ALGEBRAIC[58] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[86] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[236] = CONSTANTS[9]*(exp(ALGEBRAIC[58]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[86]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[64] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[92] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[242] = CONSTANTS[15]*(exp(ALGEBRAIC[64]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[92]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[68] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[96] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[246] = CONSTANTS[19]*(exp(ALGEBRAIC[68]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[96]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[263] = ((- ALGEBRAIC[236] - ALGEBRAIC[242]) - ALGEBRAIC[246]) - ALGEBRAIC[254]; ALGEBRAIC[59] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[87] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[237] = CONSTANTS[12]*(exp(ALGEBRAIC[59]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[87]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[65] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[93] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[243] = CONSTANTS[17]*(exp(ALGEBRAIC[65]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[93]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[69] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[97] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[247] = CONSTANTS[21]*(exp(ALGEBRAIC[69]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[97]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[264] = ((- ALGEBRAIC[237] - ALGEBRAIC[243]) - ALGEBRAIC[247]) - - ALGEBRAIC[254]; ALGEBRAIC[60] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[88] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[238] = CONSTANTS[10]*(exp(ALGEBRAIC[60]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[88]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[72] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[100] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[250] = CONSTANTS[23]*(exp(ALGEBRAIC[72]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[100]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[265] = ((- ALGEBRAIC[238]+ALGEBRAIC[242]) - ALGEBRAIC[255]) - ALGEBRAIC[250]; ALGEBRAIC[61] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[89] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[239] = CONSTANTS[13]*(exp(ALGEBRAIC[61]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[89]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[73] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[101] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[251] = CONSTANTS[25]*(exp(ALGEBRAIC[73]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[101]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[266] = ((- ALGEBRAIC[239]+ALGEBRAIC[243]) - - ALGEBRAIC[255]) - ALGEBRAIC[251]; ALGEBRAIC[62] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[90] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[240] = CONSTANTS[11]*(exp(ALGEBRAIC[62]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[90]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[267] = ((- ALGEBRAIC[240]+ALGEBRAIC[246]) - ALGEBRAIC[256])+ALGEBRAIC[250]; ALGEBRAIC[63] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[91] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[241] = CONSTANTS[14]*(exp(ALGEBRAIC[63]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[91]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[268] = ((- ALGEBRAIC[241]+ALGEBRAIC[247]) - - ALGEBRAIC[256])+ALGEBRAIC[251]; ALGEBRAIC[66] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[94] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[244] = CONSTANTS[16]*(exp(ALGEBRAIC[66]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[94]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[70] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[98] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[248] = CONSTANTS[20]*(exp(ALGEBRAIC[70]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[98]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[269] = ((- ALGEBRAIC[248]+ALGEBRAIC[236]) - ALGEBRAIC[244]) - ALGEBRAIC[257]; ALGEBRAIC[67] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[95] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[245] = CONSTANTS[18]*(exp(ALGEBRAIC[67]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[95]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[71] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[99] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[249] = CONSTANTS[22]*(exp(ALGEBRAIC[71]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[99]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[270] = ((- ALGEBRAIC[249]+ALGEBRAIC[237]) - ALGEBRAIC[245]) - - ALGEBRAIC[257]; ALGEBRAIC[74] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[102] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[252] = CONSTANTS[24]*(exp(ALGEBRAIC[74]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[102]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[271] = ((ALGEBRAIC[244]+ALGEBRAIC[238]) - ALGEBRAIC[258]) - ALGEBRAIC[252]; ALGEBRAIC[75] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[103] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41]; ALGEBRAIC[253] = CONSTANTS[26]*(exp(ALGEBRAIC[75]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[103]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[272] = ((ALGEBRAIC[245]+ALGEBRAIC[239]) - ((- ALGEBRAIC[258]+ALGEBRAIC[110]) - ALGEBRAIC[110])) - ALGEBRAIC[253]; ALGEBRAIC[273] = ((ALGEBRAIC[248]+ALGEBRAIC[240]) - ALGEBRAIC[259])+ALGEBRAIC[252]; ALGEBRAIC[274] = ((ALGEBRAIC[249]+ALGEBRAIC[241]) - ((- ALGEBRAIC[259]+ALGEBRAIC[111]) - ALGEBRAIC[111]))+ALGEBRAIC[253]; ALGEBRAIC[121] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[85]*STATES[42]); ALGEBRAIC[127] = ALGEBRAIC[121]; ALGEBRAIC[114] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]); ALGEBRAIC[122] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[86]*STATES[43]); ALGEBRAIC[133] = ALGEBRAIC[122]+ALGEBRAIC[114]; ALGEBRAIC[278] = CONSTANTS[36]*(exp(ALGEBRAIC[127]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[133]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[281] = ALGEBRAIC[278]; ALGEBRAIC[115] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]); ALGEBRAIC[119] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[83]*STATES[40]); ALGEBRAIC[125] = ALGEBRAIC[119]+ALGEBRAIC[115]; ALGEBRAIC[120] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[84]*STATES[41]); ALGEBRAIC[131] = ALGEBRAIC[120]; ALGEBRAIC[276] = CONSTANTS[34]*(exp(ALGEBRAIC[125]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[131]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[282] = - ALGEBRAIC[276]; ALGEBRAIC[116] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[89]*STATES[6]); ALGEBRAIC[128] = ALGEBRAIC[122]+ CONSTANTS[126]*ALGEBRAIC[116]; ALGEBRAIC[123] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44]); ALGEBRAIC[134] = ALGEBRAIC[123]; ALGEBRAIC[279] = CONSTANTS[37]*(exp(ALGEBRAIC[128]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[134]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[283] = - CONSTANTS[126]*ALGEBRAIC[279]; ALGEBRAIC[118] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39]); ALGEBRAIC[124] = ALGEBRAIC[118]; ALGEBRAIC[117] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[88]*STATES[5]); ALGEBRAIC[130] = CONSTANTS[127]*ALGEBRAIC[117]+ALGEBRAIC[119]; ALGEBRAIC[275] = CONSTANTS[33]*(exp(ALGEBRAIC[124]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[130]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[284] = CONSTANTS[127]*ALGEBRAIC[275]; ALGEBRAIC[113] = STATES[0]/CONSTANTS[133]; ALGEBRAIC[129] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44])+ CONSTANTS[128]*CONSTANTS[143]*ALGEBRAIC[113]; ALGEBRAIC[135] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39])+ CONSTANTS[129]*CONSTANTS[143]*ALGEBRAIC[113]; ALGEBRAIC[280] = CONSTANTS[38]*(exp(ALGEBRAIC[129]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[135]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[285] = ALGEBRAIC[280] - ALGEBRAIC[275]; ALGEBRAIC[286] = ALGEBRAIC[275] - ALGEBRAIC[276]; ALGEBRAIC[126] = ALGEBRAIC[120]; ALGEBRAIC[132] = ALGEBRAIC[121]; ALGEBRAIC[277] = CONSTANTS[35]*(exp(ALGEBRAIC[126]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[132]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[287] = ALGEBRAIC[276] - ALGEBRAIC[277]; ALGEBRAIC[288] = - ALGEBRAIC[278]+ALGEBRAIC[277]; ALGEBRAIC[289] = - ALGEBRAIC[279]+ALGEBRAIC[278]; ALGEBRAIC[290] = - ALGEBRAIC[280]+ALGEBRAIC[279]; ALGEBRAIC[156] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]); ALGEBRAIC[291] = CONSTANTS[46]*exp(ALGEBRAIC[162]/( CONSTANTS[141]*CONSTANTS[142]))*(exp(ALGEBRAIC[156]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[157]/( CONSTANTS[141]*CONSTANTS[142]))); ALGEBRAIC[292] = (( CONSTANTS[164]*ALGEBRAIC[163] - CONSTANTS[163]*ALGEBRAIC[164]) - CONSTANTS[164]*ALGEBRAIC[165])+ CONSTANTS[163]*ALGEBRAIC[166]; ALGEBRAIC[260] = CONSTANTS[143]*(((((((((((((((((((((((((((((((((((( CONSTANTS[124]*ALGEBRAIC[247] - CONSTANTS[120]*ALGEBRAIC[247])+ CONSTANTS[125]*ALGEBRAIC[251]) - CONSTANTS[121]*ALGEBRAIC[251])+ CONSTANTS[123]*ALGEBRAIC[245]) - CONSTANTS[119]*ALGEBRAIC[245])+ CONSTANTS[124]*ALGEBRAIC[249]) - CONSTANTS[120]*ALGEBRAIC[249])+ CONSTANTS[125]*ALGEBRAIC[253]) - CONSTANTS[121]*ALGEBRAIC[253])+ CONSTANTS[123]*ALGEBRAIC[243]) - CONSTANTS[119]*ALGEBRAIC[243])+ CONSTANTS[122]*ALGEBRAIC[241]) - CONSTANTS[118]*ALGEBRAIC[241])+ CONSTANTS[122]*ALGEBRAIC[239]) - CONSTANTS[118]*ALGEBRAIC[239])+ CONSTANTS[122]*ALGEBRAIC[237]) - CONSTANTS[118]*ALGEBRAIC[237])+ CONSTANTS[125]*ALGEBRAIC[252]) - CONSTANTS[121]*ALGEBRAIC[252])+ CONSTANTS[124]*ALGEBRAIC[248]) - CONSTANTS[120]*ALGEBRAIC[248])+ CONSTANTS[123]*ALGEBRAIC[244]) - CONSTANTS[119]*ALGEBRAIC[244])+ CONSTANTS[125]*ALGEBRAIC[250]) - CONSTANTS[121]*ALGEBRAIC[250])+ CONSTANTS[124]*ALGEBRAIC[246]) - CONSTANTS[120]*ALGEBRAIC[246])+ CONSTANTS[123]*ALGEBRAIC[242]) - CONSTANTS[119]*ALGEBRAIC[242])+ CONSTANTS[122]*ALGEBRAIC[240]) - CONSTANTS[118]*ALGEBRAIC[240])+ CONSTANTS[122]*ALGEBRAIC[238]) - CONSTANTS[118]*ALGEBRAIC[238])+ CONSTANTS[122]*ALGEBRAIC[236]) - CONSTANTS[118]*ALGEBRAIC[236]) - ( CONSTANTS[147]*ALGEBRAIC[110]+ CONSTANTS[147]*ALGEBRAIC[111])); ALGEBRAIC[293] = CONSTANTS[143]*( - CONSTANTS[128]*ALGEBRAIC[280]+ CONSTANTS[129]*ALGEBRAIC[280]); ALGEBRAIC[142] = - CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[141]; ALGEBRAIC[2] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<=0.00000 ? ( 0.00000*0.0300000*CONSTANTS[133])/CONSTANTS[136] : 0.00000); ALGEBRAIC[294] = ALGEBRAIC[260]+ALGEBRAIC[293]+ALGEBRAIC[142]+ALGEBRAIC[2]; } void getStateInformation(double* SI) { SI[0] = 1.0; SI[1] = 1.0; SI[2] = 1.0; SI[3] = 1.0; SI[4] = 1.0; SI[5] = 1.0; SI[6] = 1.0; SI[7] = 1.0; SI[8] = 1.0; SI[9] = 1.0; SI[10] = 1.0; SI[11] = 1.0; SI[12] = 1.0; SI[13] = 1.0; SI[14] = 1.0; SI[15] = 1.0; SI[16] = 1.0; SI[17] = 1.0; SI[18] = 1.0; SI[19] = 1.0; SI[20] = 1.0; SI[21] = 1.0; SI[22] = 1.0; SI[23] = 1.0; SI[24] = 1.0; SI[25] = 1.0; SI[26] = 1.0; SI[27] = 1.0; SI[28] = 1.0; SI[29] = 1.0; SI[30] = 1.0; SI[31] = 1.0; SI[32] = 1.0; SI[33] = 1.0; SI[34] = 1.0; SI[35] = 1.0; SI[36] = 1.0; SI[37] = 1.0; SI[38] = 1.0; SI[39] = 1.0; SI[40] = 1.0; SI[41] = 1.0; SI[42] = 1.0; SI[43] = 1.0; SI[44] = 1.0; SI[45] = 1.0; SI[46] = 1.0; SI[47] = 1.0; SI[48] = 1.0; SI[49] = 1.0; SI[50] = 1.0; SI[51] = 1.0; SI[52] = 1.0; SI[53] = 1.0; SI[54] = 1.0; SI[55] = 1.0; SI[56] = 1.0; SI[57] = 1.0; SI[62] = 0.0; SI[58] = 1.0; SI[59] = 1.0; SI[60] = 1.0; SI[61] = 1.0; SI[63] = 0.0; SI[64] = 0.0; SI[65] = 0.0; SI[66] = 0.0; SI[67] = 0.0; SI[68] = 0.0; SI[69] = 0.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = ALGEBRAIC[1] - 0.300000; CONDVAR[1] = ALGEBRAIC[1] - (0.300000+CONSTANTS[136]); CONDVAR[2] = ALGEBRAIC[1] - 0.300000; CONDVAR[3] = ALGEBRAIC[1] - (0.300000+CONSTANTS[136]); CONDVAR[4] = VOI - 0.0600000; CONDVAR[5] = VOI - 0.500000; CONDVAR[6] = VOI - 1.10000; CONDVAR[7] = STATES[64] - 0.00000; }