Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
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;
}