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