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 295 entries in the algebraic variable array.
There are a total of 70 entries in each of the rate and state variable arrays.
There are a total of 177 entries in the constant variable array.
*/
/*
* CONSTANTS[0] is kappa_leak in component BG_parameters (fmol_per_sec).
* CONSTANTS[1] is kappa_R_TRPNCa in component BG_parameters (fmol_per_sec).
* CONSTANTS[2] is kappa_R_BU in component BG_parameters (fmol_per_sec).
* CONSTANTS[3] is kappa_R_UW in component BG_parameters (fmol_per_sec).
* CONSTANTS[4] is kappa_R_WS in component BG_parameters (fmol_per_sec).
* CONSTANTS[5] is kappa_R_SU in component BG_parameters (fmol_per_sec).
* CONSTANTS[6] is kappa_diadDf in component BG_parameters (fmol_per_sec).
* CONSTANTS[7] is kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[8] is kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[9] is kappa_d000 in component BG_parameters (fmol_per_sec).
* CONSTANTS[10] is kappa_d010 in component BG_parameters (fmol_per_sec).
* CONSTANTS[11] is kappa_d020 in component BG_parameters (fmol_per_sec).
* CONSTANTS[12] is kappa_d001 in component BG_parameters (fmol_per_sec).
* CONSTANTS[13] is kappa_d011 in component BG_parameters (fmol_per_sec).
* CONSTANTS[14] is kappa_d021 in component BG_parameters (fmol_per_sec).
* CONSTANTS[15] is kappa_f1_000 in component BG_parameters (fmol_per_sec).
* CONSTANTS[16] is kappa_f1_100 in component BG_parameters (fmol_per_sec).
* CONSTANTS[17] is kappa_f1_001 in component BG_parameters (fmol_per_sec).
* CONSTANTS[18] is kappa_f1_101 in component BG_parameters (fmol_per_sec).
* CONSTANTS[19] is kappa_f2_000 in component BG_parameters (fmol_per_sec).
* CONSTANTS[20] is kappa_f2_100 in component BG_parameters (fmol_per_sec).
* CONSTANTS[21] is kappa_f2_001 in component BG_parameters (fmol_per_sec).
* CONSTANTS[22] is kappa_f2_101 in component BG_parameters (fmol_per_sec).
* CONSTANTS[23] is kappa_f3_010 in component BG_parameters (fmol_per_sec).
* CONSTANTS[24] is kappa_f3_110 in component BG_parameters (fmol_per_sec).
* CONSTANTS[25] is kappa_f3_011 in component BG_parameters (fmol_per_sec).
* CONSTANTS[26] is kappa_f3_111 in component BG_parameters (fmol_per_sec).
* CONSTANTS[27] is kappa_fCa000 in component BG_parameters (fmol_per_sec).
* CONSTANTS[28] is kappa_fCa100 in component BG_parameters (fmol_per_sec).
* CONSTANTS[29] is kappa_fCa010 in component BG_parameters (fmol_per_sec).
* CONSTANTS[30] is kappa_fCa110 in component BG_parameters (fmol_per_sec).
* CONSTANTS[31] is kappa_fCa020 in component BG_parameters (fmol_per_sec).
* CONSTANTS[32] is kappa_fCa120 in component BG_parameters (fmol_per_sec).
* CONSTANTS[33] is kappa_1_NCX in component BG_parameters (fmol_per_sec).
* CONSTANTS[34] is kappa_2_NCX in component BG_parameters (fmol_per_sec).
* CONSTANTS[35] is kappa_3_NCX in component BG_parameters (fmol_per_sec).
* CONSTANTS[36] is kappa_4_NCX in component BG_parameters (fmol_per_sec).
* CONSTANTS[37] is kappa_5_NCX in component BG_parameters (fmol_per_sec).
* CONSTANTS[38] is kappa_6_NCX in component BG_parameters (fmol_per_sec).
* CONSTANTS[39] is kappa_pCa_R1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[40] is kappa_pCa_R2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[41] is kappa_PLBph1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[42] is kappa_PLBph2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[43] is kappa_PLBd1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[44] is kappa_PLBd2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[45] is kappa_Inh in component BG_parameters (fmol_per_sec).
* CONSTANTS[46] is kappa_RyR in component BG_parameters (fmol_per_sec).
* CONSTANTS[47] is kappa_OC in component BG_parameters (fmol_per_sec).
* CONSTANTS[48] is kappa_CCI in component BG_parameters (fmol_per_sec).
* CONSTANTS[49] is kappa_CII in component BG_parameters (fmol_per_sec).
* CONSTANTS[50] is kappa_IO in component BG_parameters (fmol_per_sec).
* CONSTANTS[51] is kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[52] is kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec).
* CONSTANTS[53] is kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec).
* CONSTANTS[54] is kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec).
* CONSTANTS[55] is kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec).
* CONSTANTS[56] is kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec).
* CONSTANTS[57] is kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec).
* CONSTANTS[58] is kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec).
* CONSTANTS[59] is kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[60] is K_Ca_SR in component BG_parameters (per_fmol).
* CONSTANTS[61] is K_Ca_i in component BG_parameters (per_fmol).
* CONSTANTS[62] is K_TRPN in component BG_parameters (per_fmol).
* CONSTANTS[63] is K_Ca_TRPN in component BG_parameters (per_fmol).
* CONSTANTS[64] is K_B_XB in component BG_parameters (per_fmol).
* CONSTANTS[65] is K_U_XB in component BG_parameters (per_fmol).
* CONSTANTS[66] is K_W_XB in component BG_parameters (per_fmol).
* CONSTANTS[67] is K_S_XB in component BG_parameters (per_fmol).
* CONSTANTS[68] is K_Ca_D in component BG_parameters (per_fmol).
* CONSTANTS[69] is K_Ca_o in component BG_parameters (per_fmol).
* CONSTANTS[70] is K_000_LCC in component BG_parameters (per_fmol).
* CONSTANTS[71] is K_010_LCC in component BG_parameters (per_fmol).
* CONSTANTS[72] is K_020_LCC in component BG_parameters (per_fmol).
* CONSTANTS[73] is K_100_LCC in component BG_parameters (per_fmol).
* CONSTANTS[74] is K_110_LCC in component BG_parameters (per_fmol).
* CONSTANTS[75] is K_120_LCC in component BG_parameters (per_fmol).
* CONSTANTS[76] is K_001_LCC in component BG_parameters (per_fmol).
* CONSTANTS[77] is K_011_LCC in component BG_parameters (per_fmol).
* CONSTANTS[78] is K_021_LCC in component BG_parameters (per_fmol).
* CONSTANTS[79] is K_101_LCC in component BG_parameters (per_fmol).
* CONSTANTS[80] is K_111_LCC in component BG_parameters (per_fmol).
* CONSTANTS[81] is K_121_LCC in component BG_parameters (per_fmol).
* CONSTANTS[82] is K_1_NCX in component BG_parameters (per_fmol).
* CONSTANTS[83] is K_2_NCX in component BG_parameters (per_fmol).
* CONSTANTS[84] is K_3_NCX in component BG_parameters (per_fmol).
* CONSTANTS[85] is K_4_NCX in component BG_parameters (per_fmol).
* CONSTANTS[86] is K_5_NCX in component BG_parameters (per_fmol).
* CONSTANTS[87] is K_6_NCX in component BG_parameters (per_fmol).
* CONSTANTS[88] is K_Na_i in component BG_parameters (per_fmol).
* CONSTANTS[89] is K_Na_o in component BG_parameters (per_fmol).
* CONSTANTS[90] is K_pCa in component BG_parameters (per_fmol).
* CONSTANTS[91] is K_pCa_Ca in component BG_parameters (per_fmol).
* CONSTANTS[92] is K_PLB in component BG_parameters (per_fmol).
* CONSTANTS[93] is K_PKACI in component BG_parameters (per_fmol).
* CONSTANTS[94] is K_PLB_PKACI in component BG_parameters (per_fmol).
* CONSTANTS[95] is K_PP1 in component BG_parameters (per_fmol).
* CONSTANTS[96] is K_PLBp_PP1 in component BG_parameters (per_fmol).
* CONSTANTS[97] is K_PLBp in component BG_parameters (per_fmol).
* CONSTANTS[98] is K_Ip in component BG_parameters (per_fmol).
* CONSTANTS[99] is K_Ip_PP1 in component BG_parameters (per_fmol).
* CONSTANTS[100] is K_C_RyR in component BG_parameters (per_fmol).
* CONSTANTS[101] is K_CI_RyR in component BG_parameters (per_fmol).
* CONSTANTS[102] is K_I_RyR in component BG_parameters (per_fmol).
* CONSTANTS[103] is K_O_RyR in component BG_parameters (per_fmol).
* CONSTANTS[104] is K_P1_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[105] is K_P2_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[106] is K_P2a_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[107] is K_P4_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[108] is K_P5_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[109] is K_P6_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[110] is K_P8_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[111] is K_P9_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[112] is K_P10_SERCA in component BG_parameters (per_fmol).
* CONSTANTS[113] is K_H in component BG_parameters (per_fmol).
* CONSTANTS[114] is K_MgATP in component BG_parameters (per_fmol).
* CONSTANTS[115] is K_MgADP in component BG_parameters (per_fmol).
* CONSTANTS[116] is K_Pi in component BG_parameters (per_fmol).
* CONSTANTS[117] is z_rCa in component BG_parameters (dimensionless).
* CONSTANTS[118] is z_fd in component BG_parameters (dimensionless).
* CONSTANTS[119] is z_ff1 in component BG_parameters (dimensionless).
* CONSTANTS[120] is z_ff2 in component BG_parameters (dimensionless).
* CONSTANTS[121] is z_ff3 in component BG_parameters (dimensionless).
* CONSTANTS[122] is z_rd in component BG_parameters (dimensionless).
* CONSTANTS[123] is z_rf1 in component BG_parameters (dimensionless).
* CONSTANTS[124] is z_rf2 in component BG_parameters (dimensionless).
* CONSTANTS[125] is z_rf3 in component BG_parameters (dimensionless).
* CONSTANTS[126] is nNa_o_NCX in component BG_parameters (dimensionless).
* CONSTANTS[127] is nNa_i_NCX in component BG_parameters (dimensionless).
* CONSTANTS[128] is zf_NCX in component BG_parameters (dimensionless).
* CONSTANTS[129] is zr_NCX in component BG_parameters (dimensionless).
* CONSTANTS[130] is n_Tm in component BG_parameters (dimensionless).
* CONSTANTS[131] is n_CaTRPN in component BG_parameters (dimensionless).
* VOI is time in component environment (second).
* CONSTANTS[132] is vol_myo in component environment (pL).
* CONSTANTS[133] is C_m in component environment (fF).
* STATES[0] is q_membrane in component environment (fC).
* STATES[1] is q_Ca_D in component environment (fmol).
* STATES[2] is q_TRPN in component environment (fmol).
* STATES[3] is q_Ca_TRPN in component environment (fmol).
* STATES[4] is q_Ca_i in component environment (fmol).
* STATES[5] is q_Na_i in component environment (fmol).
* STATES[6] is q_Na_o in component environment (fmol).
* STATES[7] is q_B_XB in component environment (fmol).
* STATES[8] is q_U_XB in component environment (fmol).
* STATES[9] is q_W_XB in component environment (fmol).
* STATES[10] is q_S_XB in component environment (fmol).
* STATES[11] is q_Ca_o in component environment (fmol).
* STATES[12] is q_PLB in component environment (fmol).
* STATES[13] is q_PKACI in component environment (fmol).
* STATES[14] is q_PLB_PKACI in component environment (fmol).
* STATES[15] is q_PP1 in component environment (fmol).
* STATES[16] is q_PLBp_PP1 in component environment (fmol).
* STATES[17] is q_PLBp in component environment (fmol).
* CONSTANTS[134] is q_PLBp_baseline in component environment (fmol).
* STATES[18] is q_Ip in component environment (fmol).
* STATES[19] is q_Ip_PP1 in component environment (fmol).
* STATES[20] is q_Ca_SR in component environment (fmol).
* STATES[21] is q_pCa in component environment (fmol).
* STATES[22] is q_pCa_Ca in component environment (fmol).
* STATES[23] is q_H in component environment (fmol).
* STATES[24] is q_Pi in component environment (fmol).
* STATES[25] is q_MgADP in component environment (fmol).
* STATES[26] is q_MgATP in component environment (fmol).
* STATES[27] is q_000_LCC in component environment (fmol).
* STATES[28] is q_010_LCC in component environment (fmol).
* STATES[29] is q_020_LCC in component environment (fmol).
* STATES[30] is q_100_LCC in component environment (fmol).
* STATES[31] is q_110_LCC in component environment (fmol).
* STATES[32] is q_120_LCC in component environment (fmol).
* STATES[33] is q_001_LCC in component environment (fmol).
* STATES[34] is q_011_LCC in component environment (fmol).
* STATES[35] is q_021_LCC in component environment (fmol).
* STATES[36] is q_101_LCC in component environment (fmol).
* STATES[37] is q_111_LCC in component environment (fmol).
* STATES[38] is q_121_LCC in component environment (fmol).
* STATES[39] is q_P1_NCX in component environment (fmol).
* STATES[40] is q_P2_NCX in component environment (fmol).
* STATES[41] is q_P3_NCX in component environment (fmol).
* STATES[42] is q_P4_NCX in component environment (fmol).
* STATES[43] is q_P5_NCX in component environment (fmol).
* STATES[44] is q_P6_NCX in component environment (fmol).
* STATES[45] is q_C_RyR in component environment (fmol).
* STATES[46] is q_CI_RyR in component environment (fmol).
* STATES[47] is q_I_RyR in component environment (fmol).
* STATES[48] is q_O_RyR in component environment (fmol).
* STATES[49] is q_P1_SERCA in component environment (fmol).
* STATES[50] is q_P2_SERCA in component environment (fmol).
* STATES[51] is q_P2a_SERCA in component environment (fmol).
* STATES[52] is q_P4_SERCA in component environment (fmol).
* STATES[53] is q_P5_SERCA in component environment (fmol).
* STATES[54] is q_P6_SERCA in component environment (fmol).
* STATES[55] is q_P8_SERCA in component environment (fmol).
* STATES[56] is q_P9_SERCA in component environment (fmol).
* STATES[57] is q_P10_SERCA in component environment (fmol).
* ALGEBRAIC[14] is v_leak in component Ca_leak (fmol_per_sec).
* ALGEBRAIC[22] is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec).
* ALGEBRAIC[23] is v_R_BU in component crossbridge_TRPN (fmol_per_sec).
* ALGEBRAIC[24] is v_R_UW in component crossbridge_TRPN (fmol_per_sec).
* ALGEBRAIC[25] is v_R_WS in component crossbridge_TRPN (fmol_per_sec).
* ALGEBRAIC[26] is v_R_SU in component crossbridge_TRPN (fmol_per_sec).
* STATES[62] is T_total in component crossbridge_TRPN (kPa).
* STATES[58] is SL in component crossbridge_TRPN (metre).
* ALGEBRAIC[40] is v_diadDf in component diadCaDiffusion (fmol_per_sec).
* ALGEBRAIC[262] is v_Ca_D_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[261] is v_Ca_o_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[260] is I_mem_LCC in component LCC (fA).
* ALGEBRAIC[281] is v_Ca_o_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[282] is v_Ca_i_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[283] is v_Na_o_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[284] is v_Na_i_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[293] is I_mem_NCX in component NCX (fA).
* ALGEBRAIC[140] is v_pCa_R1 in component pCa (fmol_per_sec).
* ALGEBRAIC[141] is v_pCa_R2 in component pCa (fmol_per_sec).
* ALGEBRAIC[142] is I_mem_pCa in component pCa (fA).
* ALGEBRAIC[151] is v_PLBph1 in component PLB (fmol_per_sec).
* ALGEBRAIC[152] is v_PLBph2 in component PLB (fmol_per_sec).
* ALGEBRAIC[153] is v_PLBd1 in component PLB (fmol_per_sec).
* ALGEBRAIC[154] is v_PLBd2 in component PLB (fmol_per_sec).
* ALGEBRAIC[155] is v_Inh in component PLB (fmol_per_sec).
* ALGEBRAIC[291] is v_RyR in component RyR (fmol_per_sec).
* ALGEBRAIC[292] is v_RyRgate_Ca_D in component RyR (fmol_per_sec).
* ALGEBRAIC[219] is v_Ca_i_SERCA in component SERCA (fmol_per_sec).
* ALGEBRAIC[220] is v_Ca_SR_SERCA in component SERCA (fmol_per_sec).
* ALGEBRAIC[221] is v_H_SERCA in component SERCA (fmol_per_sec).
* ALGEBRAIC[222] is v_MgADP_SERCA in component SERCA (fmol_per_sec).
* ALGEBRAIC[223] is v_MgATP_SERCA in component SERCA (fmol_per_sec).
* ALGEBRAIC[224] is v_Pi_SERCA in component SERCA (fmol_per_sec).
* ALGEBRAIC[0] is V_m in component environment (volt).
* ALGEBRAIC[2] is I_stim in component environment (fA).
* CONSTANTS[135] is stimPeriod in component environment (second).
* CONSTANTS[136] is stimDuration in component environment (second).
* ALGEBRAIC[1] is tPeriod in component environment (second).
* ALGEBRAIC[3] is v_Ca_stim in component environment (fmol_per_sec).
* ALGEBRAIC[294] is sum_I in component environment (fA).
* ALGEBRAIC[4] is dSL in component environment (m_per_s).
* CONSTANTS[137] is V_max in component environment (m_per_s).
* CONSTANTS[138] is vol_cell in component environment (pL).
* CONSTANTS[170] is vol_SR in component environment (pL).
* CONSTANTS[139] is vol_o in component environment (pL).
* CONSTANTS[140] is vol_D in component environment (pL).
* ALGEBRAIC[5] is cCai in component environment (mM).
* ALGEBRAIC[6] is cCao in component environment (mM).
* ALGEBRAIC[8] is cCaSR in component environment (mM).
* ALGEBRAIC[7] is cCaD in component environment (mM).
* ALGEBRAIC[9] is cNai in component environment (mM).
* ALGEBRAIC[10] is pHi in component environment (dimensionless).
* ALGEBRAIC[11] is Ca_tot in component environment (fmol).
* CONSTANTS[141] is R in component constants (J_per_K_per_mol).
* CONSTANTS[142] is T in component constants (kelvin).
* CONSTANTS[143] is F in component constants (C_per_mol).
* CONSTANTS[144] is zNa in component ion_valences (dimensionless).
* CONSTANTS[145] is zCl in component ion_valences (dimensionless).
* CONSTANTS[146] is zK in component ion_valences (dimensionless).
* CONSTANTS[147] is zCa in component ion_valences (dimensionless).
* ALGEBRAIC[12] is mu_Ca_SR in component Ca_leak (J_per_mol).
* ALGEBRAIC[13] is mu_Ca_i in component Ca_leak (J_per_mol).
* ALGEBRAIC[15] is mu_TRPN in component crossbridge_TRPN (J_per_mol).
* ALGEBRAIC[16] is mu_Ca_i in component crossbridge_TRPN (J_per_mol).
* ALGEBRAIC[17] is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol).
* ALGEBRAIC[18] is mu_B_XB in component crossbridge_TRPN (J_per_mol).
* ALGEBRAIC[19] is mu_U_XB in component crossbridge_TRPN (J_per_mol).
* ALGEBRAIC[20] is mu_W_XB in component crossbridge_TRPN (J_per_mol).
* ALGEBRAIC[21] is mu_S_XB in component crossbridge_TRPN (J_per_mol).
* CONSTANTS[148] is SL_0 in component crossbridge_TRPN (metre).
* CONSTANTS[149] is q_MS in component crossbridge_TRPN (fmol).
* CONSTANTS[150] is r_s in component crossbridge_TRPN (dimensionless).
* CONSTANTS[151] is r_w in component crossbridge_TRPN (dimensionless).
* CONSTANTS[152] is A_eff in component crossbridge_TRPN (dimensionless).
* CONSTANTS[153] is phi in component crossbridge_TRPN (dimensionless).
* CONSTANTS[154] is k_uw in component crossbridge_TRPN (per_sec).
* CONSTANTS[155] is k_ws in component crossbridge_TRPN (per_sec).
* STATES[59] is G_w in component crossbridge_TRPN (metre).
* STATES[60] is G_s in component crossbridge_TRPN (metre).
* CONSTANTS[171] is c_w in component crossbridge_TRPN (per_sec).
* CONSTANTS[172] is c_s in component crossbridge_TRPN (per_sec).
* CONSTANTS[169] is A_w in component crossbridge_TRPN (dimensionless).
* CONSTANTS[175] is A_s in component crossbridge_TRPN (dimensionless).
* ALGEBRAIC[32] is nu_1 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[28] is nu_2 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[29] is nu_3 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[33] is nu_4 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[30] is nu_5 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[31] is nu_6 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[234] is v_1 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[27] is v_2 in component crossbridge_TRPN (m_per_s).
* ALGEBRAIC[235] is v_4 in component crossbridge_TRPN (m_per_s).
* CONSTANTS[156] is v_to_mu in component crossbridge_TRPN (Js_per_m2).
* ALGEBRAIC[36] is T_active in component crossbridge_TRPN (kPa).
* CONSTANTS[157] is T_ref in component crossbridge_TRPN (kPa).
* ALGEBRAIC[34] is F_T_S in component crossbridge_TRPN (kPa).
* ALGEBRAIC[35] is F_T_W in component crossbridge_TRPN (kPa).
* STATES[61] is Cdd in component crossbridge_TRPN (metre).
* STATES[63] is eta in component crossbridge_TRPN (second).
* CONSTANTS[158] is eta_l in component crossbridge_TRPN (second).
* CONSTANTS[159] is eta_s in component crossbridge_TRPN (second).
* CONSTANTS[160] is k in component crossbridge_TRPN (dimensionless).
* CONSTANTS[161] is alpha in component crossbridge_TRPN (kPa).
* STATES[64] is v_Cdd in component crossbridge_TRPN (m_per_s).
* STATES[65] is nu_d in component crossbridge_TRPN (m_per_s).
* STATES[66] is nu_k in component crossbridge_TRPN (m_per_s).
* STATES[67] is nu_f in component crossbridge_TRPN (m_per_s).
* STATES[68] is F_d in component crossbridge_TRPN (kPa).
* ALGEBRAIC[37] is F_1 in component crossbridge_TRPN (kPa).
* STATES[69] is T_passive in component crossbridge_TRPN (kPa).
* CONSTANTS[162] is b in component crossbridge_TRPN (dimensionless).
* ALGEBRAIC[38] is mu_Ca_D in component diadCaDiffusion (J_per_mol).
* ALGEBRAIC[39] is mu_Ca_i in component diadCaDiffusion (J_per_mol).
* ALGEBRAIC[42] is mu_Ca_o in component LCC (J_per_mol).
* ALGEBRAIC[43] is mu_Ca_D in component LCC (J_per_mol).
* ALGEBRAIC[41] is V_mem in component LCC (volt).
* ALGEBRAIC[56] is Af_LCC_Ca1 in component LCC (J_per_mol).
* ALGEBRAIC[84] is Ar_LCC_Ca1 in component LCC (J_per_mol).
* ALGEBRAIC[82] is Am_LCC_Ca1 in component LCC (J_per_mol).
* ALGEBRAIC[110] is v_LCC_Ca1 in component LCC (fmol_per_sec).
* ALGEBRAIC[57] is Af_LCC_Ca2 in component LCC (J_per_mol).
* ALGEBRAIC[85] is Ar_LCC_Ca2 in component LCC (J_per_mol).
* ALGEBRAIC[83] is Am_LCC_Ca2 in component LCC (J_per_mol).
* ALGEBRAIC[111] is v_LCC_Ca2 in component LCC (fmol_per_sec).
* ALGEBRAIC[76] is Af_fCa00 in component LCC (J_per_mol).
* ALGEBRAIC[104] is Ar_fCa00 in component LCC (J_per_mol).
* ALGEBRAIC[254] is v_fCa00 in component LCC (fmol_per_sec).
* ALGEBRAIC[77] is Af_fCa01 in component LCC (J_per_mol).
* ALGEBRAIC[105] is Ar_fCa01 in component LCC (J_per_mol).
* ALGEBRAIC[255] is v_fCa01 in component LCC (fmol_per_sec).
* ALGEBRAIC[78] is Af_fCa02 in component LCC (J_per_mol).
* ALGEBRAIC[106] is Ar_fCa02 in component LCC (J_per_mol).
* ALGEBRAIC[256] is v_fCa02 in component LCC (fmol_per_sec).
* ALGEBRAIC[79] is Af_fCa10 in component LCC (J_per_mol).
* ALGEBRAIC[107] is Ar_fCa10 in component LCC (J_per_mol).
* ALGEBRAIC[257] is v_fCa10 in component LCC (fmol_per_sec).
* ALGEBRAIC[80] is Af_fCa11 in component LCC (J_per_mol).
* ALGEBRAIC[108] is Ar_fCa11 in component LCC (J_per_mol).
* ALGEBRAIC[258] is v_fCa11 in component LCC (fmol_per_sec).
* ALGEBRAIC[81] is Af_fCa12 in component LCC (J_per_mol).
* ALGEBRAIC[109] is Ar_fCa12 in component LCC (J_per_mol).
* ALGEBRAIC[259] is v_fCa12 in component LCC (fmol_per_sec).
* ALGEBRAIC[44] is mu_S000_LCC in component LCC (J_per_mol).
* ALGEBRAIC[263] is v_S000_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[46] is mu_S010_LCC in component LCC (J_per_mol).
* ALGEBRAIC[265] is v_S010_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[48] is mu_S020_LCC in component LCC (J_per_mol).
* ALGEBRAIC[267] is v_S020_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[50] is mu_S100_LCC in component LCC (J_per_mol).
* ALGEBRAIC[269] is v_S100_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[52] is mu_S110_LCC in component LCC (J_per_mol).
* ALGEBRAIC[271] is v_S110_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[54] is mu_S120_LCC in component LCC (J_per_mol).
* ALGEBRAIC[273] is v_S120_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[45] is mu_S001_LCC in component LCC (J_per_mol).
* ALGEBRAIC[264] is v_S001_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[47] is mu_S011_LCC in component LCC (J_per_mol).
* ALGEBRAIC[266] is v_S011_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[49] is mu_S021_LCC in component LCC (J_per_mol).
* ALGEBRAIC[268] is v_S021_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[51] is mu_S101_LCC in component LCC (J_per_mol).
* ALGEBRAIC[270] is v_S101_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[53] is mu_S111_LCC in component LCC (J_per_mol).
* ALGEBRAIC[272] is v_S111_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[55] is mu_S121_LCC in component LCC (J_per_mol).
* ALGEBRAIC[274] is v_S121_LCC in component LCC (fmol_per_sec).
* ALGEBRAIC[58] is Af_d000 in component LCC (J_per_mol).
* ALGEBRAIC[86] is Ar_d000 in component LCC (J_per_mol).
* ALGEBRAIC[236] is v_d000 in component LCC (fmol_per_sec).
* ALGEBRAIC[60] is Af_d010 in component LCC (J_per_mol).
* ALGEBRAIC[88] is Ar_d010 in component LCC (J_per_mol).
* ALGEBRAIC[238] is v_d010 in component LCC (fmol_per_sec).
* ALGEBRAIC[62] is Af_d020 in component LCC (J_per_mol).
* ALGEBRAIC[90] is Ar_d020 in component LCC (J_per_mol).
* ALGEBRAIC[240] is v_d020 in component LCC (fmol_per_sec).
* ALGEBRAIC[64] is Af_f1_000 in component LCC (J_per_mol).
* ALGEBRAIC[92] is Ar_f1_000 in component LCC (J_per_mol).
* ALGEBRAIC[242] is v_f1_000 in component LCC (fmol_per_sec).
* ALGEBRAIC[68] is Af_f2_000 in component LCC (J_per_mol).
* ALGEBRAIC[96] is Ar_f2_000 in component LCC (J_per_mol).
* ALGEBRAIC[246] is v_f2_000 in component LCC (fmol_per_sec).
* ALGEBRAIC[72] is Af_f3_010 in component LCC (J_per_mol).
* ALGEBRAIC[100] is Ar_f3_010 in component LCC (J_per_mol).
* ALGEBRAIC[250] is v_f3_010 in component LCC (fmol_per_sec).
* ALGEBRAIC[66] is Af_f1_100 in component LCC (J_per_mol).
* ALGEBRAIC[94] is Ar_f1_100 in component LCC (J_per_mol).
* ALGEBRAIC[244] is v_f1_100 in component LCC (fmol_per_sec).
* ALGEBRAIC[70] is Af_f2_100 in component LCC (J_per_mol).
* ALGEBRAIC[98] is Ar_f2_100 in component LCC (J_per_mol).
* ALGEBRAIC[248] is v_f2_100 in component LCC (fmol_per_sec).
* ALGEBRAIC[74] is Af_f3_110 in component LCC (J_per_mol).
* ALGEBRAIC[102] is Ar_f3_110 in component LCC (J_per_mol).
* ALGEBRAIC[252] is v_f3_110 in component LCC (fmol_per_sec).
* ALGEBRAIC[59] is Af_d001 in component LCC (J_per_mol).
* ALGEBRAIC[87] is Ar_d001 in component LCC (J_per_mol).
* ALGEBRAIC[237] is v_d001 in component LCC (fmol_per_sec).
* ALGEBRAIC[61] is Af_d011 in component LCC (J_per_mol).
* ALGEBRAIC[89] is Ar_d011 in component LCC (J_per_mol).
* ALGEBRAIC[239] is v_d011 in component LCC (fmol_per_sec).
* ALGEBRAIC[63] is Af_d021 in component LCC (J_per_mol).
* ALGEBRAIC[91] is Ar_d021 in component LCC (J_per_mol).
* ALGEBRAIC[241] is v_d021 in component LCC (fmol_per_sec).
* ALGEBRAIC[65] is Af_f1_001 in component LCC (J_per_mol).
* ALGEBRAIC[93] is Ar_f1_001 in component LCC (J_per_mol).
* ALGEBRAIC[243] is v_f1_001 in component LCC (fmol_per_sec).
* ALGEBRAIC[69] is Af_f2_001 in component LCC (J_per_mol).
* ALGEBRAIC[97] is Ar_f2_001 in component LCC (J_per_mol).
* ALGEBRAIC[247] is v_f2_001 in component LCC (fmol_per_sec).
* ALGEBRAIC[73] is Af_f3_011 in component LCC (J_per_mol).
* ALGEBRAIC[101] is Ar_f3_011 in component LCC (J_per_mol).
* ALGEBRAIC[251] is v_f3_011 in component LCC (fmol_per_sec).
* ALGEBRAIC[67] is Af_f1_101 in component LCC (J_per_mol).
* ALGEBRAIC[95] is Ar_f1_101 in component LCC (J_per_mol).
* ALGEBRAIC[245] is v_f1_101 in component LCC (fmol_per_sec).
* ALGEBRAIC[71] is Af_f2_101 in component LCC (J_per_mol).
* ALGEBRAIC[99] is Ar_f2_101 in component LCC (J_per_mol).
* ALGEBRAIC[249] is v_f2_101 in component LCC (fmol_per_sec).
* ALGEBRAIC[75] is Af_f3_111 in component LCC (J_per_mol).
* ALGEBRAIC[103] is Ar_f3_111 in component LCC (J_per_mol).
* ALGEBRAIC[253] is v_f3_111 in component LCC (fmol_per_sec).
* ALGEBRAIC[112] is Ca_tot in component LCC (fmol).
* ALGEBRAIC[114] is mu_Ca_o in component NCX (J_per_mol).
* ALGEBRAIC[115] is mu_Ca_i in component NCX (J_per_mol).
* ALGEBRAIC[116] is mu_Na_o in component NCX (J_per_mol).
* ALGEBRAIC[117] is mu_Na_i in component NCX (J_per_mol).
* ALGEBRAIC[118] is mu_P1_NCX in component NCX (J_per_mol).
* ALGEBRAIC[119] is mu_P2_NCX in component NCX (J_per_mol).
* ALGEBRAIC[120] is mu_P3_NCX in component NCX (J_per_mol).
* ALGEBRAIC[121] is mu_P4_NCX in component NCX (J_per_mol).
* ALGEBRAIC[122] is mu_P5_NCX in component NCX (J_per_mol).
* ALGEBRAIC[123] is mu_P6_NCX in component NCX (J_per_mol).
* ALGEBRAIC[285] is v_P1_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[286] is v_P2_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[287] is v_P3_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[288] is v_P4_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[289] is v_P5_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[290] is v_P6_NCX in component NCX (fmol_per_sec).
* ALGEBRAIC[124] is Af_r1 in component NCX (J_per_mol).
* ALGEBRAIC[130] is Ar_r1 in component NCX (J_per_mol).
* ALGEBRAIC[275] is v_r1 in component NCX (fmol_per_sec).
* ALGEBRAIC[125] is Af_r2 in component NCX (J_per_mol).
* ALGEBRAIC[131] is Ar_r2 in component NCX (J_per_mol).
* ALGEBRAIC[276] is v_r2 in component NCX (fmol_per_sec).
* ALGEBRAIC[126] is Af_r3 in component NCX (J_per_mol).
* ALGEBRAIC[132] is Ar_r3 in component NCX (J_per_mol).
* ALGEBRAIC[277] is v_r3 in component NCX (fmol_per_sec).
* ALGEBRAIC[127] is Af_r4 in component NCX (J_per_mol).
* ALGEBRAIC[133] is Ar_r4 in component NCX (J_per_mol).
* ALGEBRAIC[278] is v_r4 in component NCX (fmol_per_sec).
* ALGEBRAIC[128] is Af_r5 in component NCX (J_per_mol).
* ALGEBRAIC[134] is Ar_r5 in component NCX (J_per_mol).
* ALGEBRAIC[279] is v_r5 in component NCX (fmol_per_sec).
* ALGEBRAIC[129] is Af_r6 in component NCX (J_per_mol).
* ALGEBRAIC[135] is Ar_r6 in component NCX (J_per_mol).
* ALGEBRAIC[280] is v_r6 in component NCX (fmol_per_sec).
* ALGEBRAIC[113] is V_mem in component NCX (volt).
* ALGEBRAIC[136] is mu_Ca_i in component pCa (J_per_mol).
* ALGEBRAIC[137] is mu_Ca_o in component pCa (J_per_mol).
* ALGEBRAIC[138] is mu_pCa in component pCa (J_per_mol).
* ALGEBRAIC[139] is mu_pCa_Ca in component pCa (J_per_mol).
* ALGEBRAIC[143] is mu_PLB in component PLB (J_per_mol).
* ALGEBRAIC[144] is mu_PKACI in component PLB (J_per_mol).
* ALGEBRAIC[145] is mu_PLB_PKACI in component PLB (J_per_mol).
* ALGEBRAIC[146] is mu_PP1 in component PLB (J_per_mol).
* ALGEBRAIC[147] is mu_PLBp_PP1 in component PLB (J_per_mol).
* ALGEBRAIC[148] is mu_PLBp in component PLB (J_per_mol).
* ALGEBRAIC[149] is mu_Ip in component PLB (J_per_mol).
* ALGEBRAIC[150] is mu_Ip_PP1 in component PLB (J_per_mol).
* CONSTANTS[163] is nCa_1 in component RyR (dimensionless).
* CONSTANTS[164] is nCa_2 in component RyR (dimensionless).
* ALGEBRAIC[156] is mu_Ca_SR in component RyR (J_per_mol).
* ALGEBRAIC[157] is mu_Ca_D in component RyR (J_per_mol).
* ALGEBRAIC[158] is q_Ca_gate_complexes in component RyR (fmol).
* ALGEBRAIC[159] is mu_C_RyR in component RyR (J_per_mol).
* ALGEBRAIC[160] is mu_CI_RyR in component RyR (J_per_mol).
* ALGEBRAIC[161] is mu_I_RyR in component RyR (J_per_mol).
* ALGEBRAIC[162] is mu_O_RyR in component RyR (J_per_mol).
* ALGEBRAIC[163] is v_OC in component RyR (fmol_per_sec).
* ALGEBRAIC[164] is v_CCI in component RyR (fmol_per_sec).
* ALGEBRAIC[165] is v_CII in component RyR (fmol_per_sec).
* ALGEBRAIC[166] is v_IO in component RyR (fmol_per_sec).
* CONSTANTS[165] is n_Ca_i in component SERCA (dimensionless).
* CONSTANTS[166] is n_Ca_SR in component SERCA (dimensionless).
* CONSTANTS[167] is n_H in component SERCA (dimensionless).
* ALGEBRAIC[167] is q_SERCA_Ca_complexes in component SERCA (fmol).
* ALGEBRAIC[168] is c_Ca_i in component SERCA (mM).
* ALGEBRAIC[169] is c_Ca_SR in component SERCA (mM).
* ALGEBRAIC[170] is c_H in component SERCA (mM).
* ALGEBRAIC[171] is c_MgADP in component SERCA (mM).
* ALGEBRAIC[172] is c_MgATP in component SERCA (mM).
* ALGEBRAIC[173] is c_Pi in component SERCA (mM).
* CONSTANTS[168] is vol_i in component SERCA (pL).
* CONSTANTS[174] is vol_sr in component SERCA (pL).
* CONSTANTS[176] is vol_isr in component SERCA (pL).
* ALGEBRAIC[174] is mu_Ca_i in component SERCA (J_per_mol).
* ALGEBRAIC[175] is mu_Ca_SR in component SERCA (J_per_mol).
* ALGEBRAIC[176] is mu_H in component SERCA (J_per_mol).
* ALGEBRAIC[177] is mu_MgADP in component SERCA (J_per_mol).
* ALGEBRAIC[178] is mu_MgATP in component SERCA (J_per_mol).
* ALGEBRAIC[180] is mu_P1 in component SERCA (J_per_mol).
* ALGEBRAIC[225] is v_P1 in component SERCA (fmol_per_sec).
* ALGEBRAIC[179] is mu_Pi in component SERCA (J_per_mol).
* ALGEBRAIC[181] is mu_P2 in component SERCA (J_per_mol).
* ALGEBRAIC[226] is v_P2 in component SERCA (fmol_per_sec).
* ALGEBRAIC[182] is mu_P2a in component SERCA (J_per_mol).
* ALGEBRAIC[227] is v_P2a in component SERCA (fmol_per_sec).
* ALGEBRAIC[183] is mu_P4 in component SERCA (J_per_mol).
* ALGEBRAIC[228] is v_P4 in component SERCA (fmol_per_sec).
* ALGEBRAIC[184] is mu_P5 in component SERCA (J_per_mol).
* ALGEBRAIC[229] is v_P5 in component SERCA (fmol_per_sec).
* ALGEBRAIC[186] is mu_P6 in component SERCA (J_per_mol).
* ALGEBRAIC[231] is v_P6 in component SERCA (fmol_per_sec).
* ALGEBRAIC[187] is mu_P8 in component SERCA (J_per_mol).
* ALGEBRAIC[232] is v_P8 in component SERCA (fmol_per_sec).
* ALGEBRAIC[188] is mu_P9 in component SERCA (J_per_mol).
* ALGEBRAIC[233] is v_P9 in component SERCA (fmol_per_sec).
* ALGEBRAIC[185] is mu_P10 in component SERCA (J_per_mol).
* ALGEBRAIC[230] is v_P10 in component SERCA (fmol_per_sec).
* ALGEBRAIC[189] is mu_PLBp_baseline in component SERCA (J_per_mol).
* CONSTANTS[173] is mu_PLBp in component SERCA (J_per_mol).
* ALGEBRAIC[190] is mu_PLBp_eff in component SERCA (J_per_mol).
* ALGEBRAIC[193] is Af_R1_2 in component SERCA (J_per_mol).
* ALGEBRAIC[194] is Ar_R1_2 in component SERCA (J_per_mol).
* ALGEBRAIC[211] is v_SERCA_R1_2 in component SERCA (fmol_per_sec).
* ALGEBRAIC[195] is Af_R5_6 in component SERCA (J_per_mol).
* ALGEBRAIC[196] is Ar_R5_6 in component SERCA (J_per_mol).
* ALGEBRAIC[212] is v_SERCA_R5_6 in component SERCA (fmol_per_sec).
* ALGEBRAIC[197] is Af_R2_2a in component SERCA (J_per_mol).
* ALGEBRAIC[198] is Ar_R2_2a in component SERCA (J_per_mol).
* ALGEBRAIC[213] is v_SERCA_R2_2a in component SERCA (fmol_per_sec).
* ALGEBRAIC[199] is Af_R2_4 in component SERCA (J_per_mol).
* ALGEBRAIC[200] is Ar_R2_4 in component SERCA (J_per_mol).
* ALGEBRAIC[214] is v_SERCA_R2_4 in component SERCA (fmol_per_sec).
* ALGEBRAIC[201] is Af_R4_5 in component SERCA (J_per_mol).
* ALGEBRAIC[202] is Ar_R4_5 in component SERCA (J_per_mol).
* ALGEBRAIC[215] is v_SERCA_R4_5 in component SERCA (fmol_per_sec).
* ALGEBRAIC[203] is Af_R6_8 in component SERCA (J_per_mol).
* ALGEBRAIC[204] is Ar_R6_8 in component SERCA (J_per_mol).
* ALGEBRAIC[216] is v_SERCA_R6_8 in component SERCA (fmol_per_sec).
* ALGEBRAIC[205] is Af_R8_9 in component SERCA (J_per_mol).
* ALGEBRAIC[206] is Ar_R8_9 in component SERCA (J_per_mol).
* ALGEBRAIC[217] is v_SERCA_R8_9 in component SERCA (fmol_per_sec).
* ALGEBRAIC[207] is Af_R9_10 in component SERCA (J_per_mol).
* ALGEBRAIC[208] is Ar_R9_10 in component SERCA (J_per_mol).
* ALGEBRAIC[218] is v_SERCA_R9_10 in component SERCA (fmol_per_sec).
* ALGEBRAIC[191] is Af_R10_1 in component SERCA (J_per_mol).
* ALGEBRAIC[192] is Ar_R10_1 in component SERCA (J_per_mol).
* ALGEBRAIC[209] is v_SERCA_R10_1 in component SERCA (fmol_per_sec).
* ALGEBRAIC[210] is v_SERCA_R10_1_noPLBp in component SERCA (fmol_per_sec).
* RATES[20] is d/dt q_Ca_SR in component environment (fmol).
* RATES[4] is d/dt q_Ca_i in component environment (fmol).
* RATES[2] is d/dt q_TRPN in component environment (fmol).
* RATES[3] is d/dt q_Ca_TRPN in component environment (fmol).
* RATES[7] is d/dt q_B_XB in component environment (fmol).
* RATES[8] is d/dt q_U_XB in component environment (fmol).
* RATES[9] is d/dt q_W_XB in component environment (fmol).
* RATES[10] is d/dt q_S_XB in component environment (fmol).
* RATES[1] is d/dt q_Ca_D in component environment (fmol).
* RATES[11] is d/dt q_Ca_o in component environment (fmol).
* RATES[5] is d/dt q_Na_i in component environment (fmol).
* RATES[6] is d/dt q_Na_o in component environment (fmol).
* RATES[21] is d/dt q_pCa in component environment (fmol).
* RATES[22] is d/dt q_pCa_Ca in component environment (fmol).
* RATES[12] is d/dt q_PLB in component environment (fmol).
* RATES[13] is d/dt q_PKACI in component environment (fmol).
* RATES[14] is d/dt q_PLB_PKACI in component environment (fmol).
* RATES[15] is d/dt q_PP1 in component environment (fmol).
* RATES[16] is d/dt q_PLBp_PP1 in component environment (fmol).
* RATES[17] is d/dt q_PLBp in component environment (fmol).
* RATES[18] is d/dt q_Ip in component environment (fmol).
* RATES[19] is d/dt q_Ip_PP1 in component environment (fmol).
* RATES[23] is d/dt q_H in component environment (fmol).
* RATES[25] is d/dt q_MgADP in component environment (fmol).
* RATES[26] is d/dt q_MgATP in component environment (fmol).
* RATES[24] is d/dt q_Pi in component environment (fmol).
* RATES[0] is d/dt q_membrane in component environment (fC).
* RATES[59] is d/dt G_w in component crossbridge_TRPN (metre).
* RATES[58] is d/dt SL in component crossbridge_TRPN (metre).
* RATES[60] is d/dt G_s in component crossbridge_TRPN (metre).
* RATES[61] is d/dt Cdd in component crossbridge_TRPN (metre).
* RATES[27] is d/dt q_000_LCC in component environment (fmol).
* RATES[28] is d/dt q_010_LCC in component environment (fmol).
* RATES[29] is d/dt q_020_LCC in component environment (fmol).
* RATES[30] is d/dt q_100_LCC in component environment (fmol).
* RATES[31] is d/dt q_110_LCC in component environment (fmol).
* RATES[32] is d/dt q_120_LCC in component environment (fmol).
* RATES[33] is d/dt q_001_LCC in component environment (fmol).
* RATES[34] is d/dt q_011_LCC in component environment (fmol).
* RATES[35] is d/dt q_021_LCC in component environment (fmol).
* RATES[36] is d/dt q_101_LCC in component environment (fmol).
* RATES[37] is d/dt q_111_LCC in component environment (fmol).
* RATES[38] is d/dt q_121_LCC in component environment (fmol).
* RATES[39] is d/dt q_P1_NCX in component environment (fmol).
* RATES[40] is d/dt q_P2_NCX in component environment (fmol).
* RATES[41] is d/dt q_P3_NCX in component environment (fmol).
* RATES[42] is d/dt q_P4_NCX in component environment (fmol).
* RATES[43] is d/dt q_P5_NCX in component environment (fmol).
* RATES[44] is d/dt q_P6_NCX in component environment (fmol).
* RATES[48] is d/dt q_O_RyR in component environment (fmol).
* RATES[45] is d/dt q_C_RyR in component environment (fmol).
* RATES[46] is d/dt q_CI_RyR in component environment (fmol).
* RATES[47] is d/dt q_I_RyR in component environment (fmol).
* RATES[49] is d/dt q_P1_SERCA in component environment (fmol).
* RATES[50] is d/dt q_P2_SERCA in component environment (fmol).
* RATES[51] is d/dt q_P2a_SERCA in component environment (fmol).
* RATES[52] is d/dt q_P4_SERCA in component environment (fmol).
* RATES[53] is d/dt q_P5_SERCA in component environment (fmol).
* RATES[57] is d/dt q_P10_SERCA in component environment (fmol).
* RATES[54] is d/dt q_P6_SERCA in component environment (fmol).
* RATES[55] is d/dt q_P8_SERCA in component environment (fmol).
* RATES[56] is d/dt q_P9_SERCA in component environment (fmol).
* There are a total of 8 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1.19278e-07;
CONSTANTS[1] = 61.975;
CONSTANTS[2] = 6.70563;
CONSTANTS[3] = 0.295103;
CONSTANTS[4] = 0.00694361;
CONSTANTS[5] = 124985;
CONSTANTS[6] = 0.0319966;
CONSTANTS[7] = 1.45746;
CONSTANTS[8] = 1.91895;
CONSTANTS[9] = 306.532;
CONSTANTS[10] = 15.4798;
CONSTANTS[11] = 20.3814;
CONSTANTS[12] = 0.0998404;
CONSTANTS[13] = 0.00504193;
CONSTANTS[14] = 0.00663841;
CONSTANTS[15] = 29.5266;
CONSTANTS[16] = 146.621;
CONSTANTS[17] = 0.00961708;
CONSTANTS[18] = 0.047756;
CONSTANTS[19] = 1.91912;
CONSTANTS[20] = 9.52987;
CONSTANTS[21] = 0.000625075;
CONSTANTS[22] = 0.00310397;
CONSTANTS[23] = 4187.14;
CONSTANTS[24] = 20792.3;
CONSTANTS[25] = 1.36379;
CONSTANTS[26] = 6.77224;
CONSTANTS[27] = 34185.2;
CONSTANTS[28] = 169755;
CONSTANTS[29] = 1726.35;
CONSTANTS[30] = 8572.61;
CONSTANTS[31] = 2272.98;
CONSTANTS[32] = 11287.1;
CONSTANTS[33] = 429216;
CONSTANTS[34] = 1.06923e+08;
CONSTANTS[35] = 3622.96;
CONSTANTS[36] = 1.06923e+08;
CONSTANTS[37] = 429216;
CONSTANTS[38] = 0.380974;
CONSTANTS[39] = 10726.4;
CONSTANTS[40] = 1.40741e-05;
CONSTANTS[41] = 12.9276;
CONSTANTS[42] = 1.86249;
CONSTANTS[43] = 0.0640148;
CONSTANTS[44] = 0.200763;
CONSTANTS[45] = 1269.27;
CONSTANTS[46] = 1.13157e+06;
CONSTANTS[47] = 8.50173;
CONSTANTS[48] = 0.164843;
CONSTANTS[49] = 626.393;
CONSTANTS[50] = 52.1994;
CONSTANTS[51] = 3.59662e-07;
CONSTANTS[52] = 6.73317;
CONSTANTS[53] = 32063.4;
CONSTANTS[54] = 6.73317;
CONSTANTS[55] = 0.408566;
CONSTANTS[56] = 313.54;
CONSTANTS[57] = 6675;
CONSTANTS[58] = 6675;
CONSTANTS[59] = 5.05704e-06;
CONSTANTS[60] = 4.69975;
CONSTANTS[61] = 1.11404;
CONSTANTS[62] = 12.0566;
CONSTANTS[63] = 0.062444;
CONSTANTS[64] = 9.09233;
CONSTANTS[65] = 3.40963;
CONSTANTS[66] = 22.2937;
CONSTANTS[67] = 5.57343e-06;
CONSTANTS[68] = 66.2656;
CONSTANTS[69] = 0.0265011;
CONSTANTS[70] = 1.58796;
CONSTANTS[71] = 31.4449;
CONSTANTS[72] = 23.8827;
CONSTANTS[73] = 0.319783;
CONSTANTS[74] = 6.33236;
CONSTANTS[75] = 4.80948;
CONSTANTS[76] = 4875.4;
CONSTANTS[77] = 96542.9;
CONSTANTS[78] = 73325.1;
CONSTANTS[79] = 981.805;
CONSTANTS[80] = 19441.8;
CONSTANTS[81] = 14766.2;
CONSTANTS[82] = 2.85254e+08;
CONSTANTS[83] = 16.9249;
CONSTANTS[84] = 34379.9;
CONSTANTS[85] = 313010;
CONSTANTS[86] = 6.01862e+07;
CONSTANTS[87] = 929.686;
CONSTANTS[88] = 0.012842;
CONSTANTS[89] = 0.0107799;
CONSTANTS[90] = 0.0126813;
CONSTANTS[91] = 0.0178289;
CONSTANTS[92] = 0.482706;
CONSTANTS[93] = 0.0120458;
CONSTANTS[94] = 2.99358;
CONSTANTS[95] = 0.0189254;
CONSTANTS[96] = 1.63848;
CONSTANTS[97] = 66.755;
CONSTANTS[98] = 62.347;
CONSTANTS[99] = 0.0304897;
CONSTANTS[100] = 2234.81;
CONSTANTS[101] = 30.332;
CONSTANTS[102] = 0.0957866;
CONSTANTS[103] = 7.05739;
CONSTANTS[104] = 62.115;
CONSTANTS[105] = 1.88391;
CONSTANTS[106] = 0.271481;
CONSTANTS[107] = 1292.79;
CONSTANTS[108] = 0.00208939;
CONSTANTS[109] = 139.3;
CONSTANTS[110] = 0.710561;
CONSTANTS[111] = 1.30406;
CONSTANTS[112] = 1.3624;
CONSTANTS[113] = 486.59;
CONSTANTS[114] = 15.0787;
CONSTANTS[115] = 1.53572e-05;
CONSTANTS[116] = 0.00616948;
CONSTANTS[117] = 2;
CONSTANTS[118] = 2.1404;
CONSTANTS[119] = -1.1495;
CONSTANTS[120] = 0.72162;
CONSTANTS[121] = 4.2933;
CONSTANTS[122] = -2.1404;
CONSTANTS[123] = 1.8993;
CONSTANTS[124] = -0.52288;
CONSTANTS[125] = 0;
CONSTANTS[126] = 3;
CONSTANTS[127] = 3;
CONSTANTS[128] = -0.55679;
CONSTANTS[129] = 0.44321;
CONSTANTS[130] = 1.1;
CONSTANTS[131] = 2;
CONSTANTS[132] = 25.84;
CONSTANTS[133] = 153400;
STATES[0] = -13039;
STATES[1] = 3.68876E-07;
STATES[2] = 2.57;
STATES[3] = 1e-18;
STATES[4] = 0.005663912;
STATES[5] = 481.46;
STATES[6] = 725.48;
STATES[7] = 1;
STATES[8] = 0;
STATES[9] = 0;
STATES[10] = 0;
STATES[11] = 6.84;
STATES[12] = 4.028E+00;
STATES[13] = 2.234E-03;
STATES[14] = 1e-18;
STATES[15] = 3.382E-02;
STATES[16] = 1e-18;
STATES[17] = 1e-18;
CONSTANTS[134] = 1e-18;
STATES[18] = 1.999E-03;
STATES[19] = 1e-18;
STATES[20] = 0.641;
STATES[21] = 0.0032;
STATES[22] = 1e-9;
STATES[23] = 2.584e-6;
STATES[24] = 570;
STATES[25] = 1.3794;
STATES[26] = 3.8;
STATES[27] = 1.3798e-05;
STATES[28] = 1.3798e-05;
STATES[29] = 1.3798e-05;
STATES[30] = 1.3798e-05;
STATES[31] = 1.3798e-05;
STATES[32] = 1.3798e-05;
STATES[33] = 1.3798e-05;
STATES[34] = 1.3798e-05;
STATES[35] = 1.3798e-05;
STATES[36] = 1.3798e-05;
STATES[37] = 1.3798e-05;
STATES[38] = 1.3798e-05;
STATES[39] = 4.24554e-07;
STATES[40] = 4.24554e-07;
STATES[41] = 4.24554e-07;
STATES[42] = 4.24554e-07;
STATES[43] = 4.24554e-07;
STATES[44] = 4.24554e-07;
STATES[45] = 4.13941e-05;
STATES[46] = 4.13941e-05;
STATES[47] = 4.13941e-05;
STATES[48] = 4.13941e-05;
STATES[49] = 2.83036e-08;
STATES[50] = 2.83036e-08;
STATES[51] = 2.83036e-08;
STATES[52] = 2.83036e-08;
STATES[53] = 2.83036e-08;
STATES[54] = 2.83036e-08;
STATES[55] = 2.83036e-08;
STATES[56] = 2.83036e-08;
STATES[57] = 2.83036e-08;
STATES[58] = 2e-6;
CONSTANTS[135] = 1;
CONSTANTS[136] = 0.001;
CONSTANTS[137] = 19e-6;
CONSTANTS[138] = 38;
CONSTANTS[139] = 5.182;
CONSTANTS[140] = 0.0025375;
CONSTANTS[141] = 8.31;
CONSTANTS[142] = 310;
CONSTANTS[143] = 96485;
CONSTANTS[144] = 1;
CONSTANTS[145] = -1;
CONSTANTS[146] = 1;
CONSTANTS[147] = 2;
CONSTANTS[148] = 2e-6;
CONSTANTS[149] = 1e-6;
CONSTANTS[150] = 0.25;
CONSTANTS[151] = 0.5;
CONSTANTS[152] = 25;
CONSTANTS[153] = 2.23;
CONSTANTS[154] = 26;
CONSTANTS[155] = 4;
STATES[59] = 1e-6;
STATES[60] = 1e-6;
CONSTANTS[156] = 1;
CONSTANTS[157] = 40.5;
STATES[61] = 0;
CONSTANTS[158] = 200e3;
CONSTANTS[159] = 20e3;
CONSTANTS[160] = 7;
CONSTANTS[161] = 2.1;
CONSTANTS[162] = 9.1;
CONSTANTS[163] = 1;
CONSTANTS[164] = 2;
CONSTANTS[165] = 2;
CONSTANTS[166] = 2;
CONSTANTS[167] = 2;
CONSTANTS[168] = 34.0;
CONSTANTS[169] = ( CONSTANTS[152]*CONSTANTS[150])/( (1.00000 - CONSTANTS[150])*CONSTANTS[151]+CONSTANTS[150]);
CONSTANTS[170] = CONSTANTS[138]*0.0350000;
CONSTANTS[171] = ( CONSTANTS[153]*CONSTANTS[154]*(1.00000 - CONSTANTS[151]))/CONSTANTS[151];
CONSTANTS[172] = ( CONSTANTS[153]*CONSTANTS[155]*(1.00000 - CONSTANTS[150])*CONSTANTS[151])/CONSTANTS[150];
CONSTANTS[173] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*CONSTANTS[134]);
CONSTANTS[174] = CONSTANTS[168]*0.0350000;
CONSTANTS[175] = CONSTANTS[169];
CONSTANTS[176] = CONSTANTS[168]+CONSTANTS[174];
STATES[62] = 0.1001;
STATES[63] = 0.1001;
STATES[64] = 0.1001;
STATES[65] = 0.1001;
STATES[66] = 0.1001;
STATES[67] = 0.1001;
STATES[68] = 0.1001;
STATES[69] = 0.1001;
RATES[20] = 0.1001;
RATES[4] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[1] = 0.1001;
RATES[11] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[21] = 0.1001;
RATES[22] = 0.1001;
RATES[12] = 0.1001;
RATES[13] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[16] = 0.1001;
RATES[17] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[23] = 0.1001;
RATES[25] = 0.1001;
RATES[26] = 0.1001;
RATES[24] = 0.1001;
RATES[0] = 0.1001;
RATES[59] = 0.1001;
RATES[58] = 0.1001;
RATES[60] = 0.1001;
RATES[61] = 0.1001;
RATES[27] = 0.1001;
RATES[28] = 0.1001;
RATES[29] = 0.1001;
RATES[30] = 0.1001;
RATES[31] = 0.1001;
RATES[32] = 0.1001;
RATES[33] = 0.1001;
RATES[34] = 0.1001;
RATES[35] = 0.1001;
RATES[36] = 0.1001;
RATES[37] = 0.1001;
RATES[38] = 0.1001;
RATES[39] = 0.1001;
RATES[40] = 0.1001;
RATES[41] = 0.1001;
RATES[42] = 0.1001;
RATES[43] = 0.1001;
RATES[44] = 0.1001;
RATES[48] = 0.1001;
RATES[45] = 0.1001;
RATES[46] = 0.1001;
RATES[47] = 0.1001;
RATES[49] = 0.1001;
RATES[50] = 0.1001;
RATES[51] = 0.1001;
RATES[52] = 0.1001;
RATES[53] = 0.1001;
RATES[57] = 0.1001;
RATES[54] = 0.1001;
RATES[55] = 0.1001;
RATES[56] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[20] - (- ALGEBRAIC[14] - ALGEBRAIC[291])+ALGEBRAIC[220];
resid[1] = RATES[4] - (( - ALGEBRAIC[22]*CONSTANTS[131]+ALGEBRAIC[219]+ALGEBRAIC[14]+ALGEBRAIC[40]+ALGEBRAIC[282]) - ALGEBRAIC[140])+ALGEBRAIC[3];
resid[2] = RATES[2] - - ALGEBRAIC[22];
resid[3] = RATES[3] - ALGEBRAIC[22];
resid[4] = RATES[7] - - ALGEBRAIC[23];
resid[5] = RATES[8] - (ALGEBRAIC[23] - ALGEBRAIC[24])+ALGEBRAIC[26];
resid[6] = RATES[9] - ALGEBRAIC[24] - ALGEBRAIC[25];
resid[7] = RATES[10] - ALGEBRAIC[25] - ALGEBRAIC[26];
resid[8] = RATES[1] - - ALGEBRAIC[40]+ALGEBRAIC[262]+ALGEBRAIC[291]+ALGEBRAIC[292];
resid[9] = RATES[11] - ALGEBRAIC[261]+ALGEBRAIC[281]+ALGEBRAIC[141];
resid[10] = RATES[5] - ALGEBRAIC[284];
resid[11] = RATES[6] - ALGEBRAIC[283];
resid[12] = RATES[21] - - ALGEBRAIC[140]+ALGEBRAIC[141];
resid[13] = RATES[22] - ALGEBRAIC[140] - ALGEBRAIC[141];
resid[14] = RATES[12] - ALGEBRAIC[154] - ALGEBRAIC[151];
resid[15] = RATES[13] - ALGEBRAIC[152] - ALGEBRAIC[151];
resid[16] = RATES[14] - ALGEBRAIC[151] - ALGEBRAIC[152];
resid[17] = RATES[15] - ALGEBRAIC[151] - ALGEBRAIC[152];
resid[18] = RATES[16] - ALGEBRAIC[153] - ALGEBRAIC[154];
resid[19] = RATES[17] - ALGEBRAIC[152] - ALGEBRAIC[153];
resid[20] = RATES[18] - - ALGEBRAIC[155];
resid[21] = RATES[19] - ALGEBRAIC[155];
resid[22] = RATES[23] - ALGEBRAIC[221];
resid[23] = RATES[25] - ALGEBRAIC[222];
resid[24] = RATES[26] - ALGEBRAIC[223];
resid[25] = RATES[24] - ALGEBRAIC[224];
resid[26] = RATES[0] - ALGEBRAIC[294];
resid[27] = RATES[59] - ALGEBRAIC[234];
resid[28] = RATES[58] - ALGEBRAIC[27];
resid[29] = RATES[60] - ALGEBRAIC[235];
resid[30] = STATES[63] - (CONDVAR[7]>0.00000 ? CONSTANTS[158] : CONSTANTS[159]);
resid[31] = RATES[61] - STATES[64];
resid[32] = STATES[64] - STATES[65];
resid[33] = STATES[66] - (CONSTANTS[160]/STATES[63])*STATES[61];
resid[34] = STATES[67] - (CONSTANTS[160]/STATES[63])*(STATES[58] - CONSTANTS[148]);
resid[35] = STATES[65] - - STATES[66]+STATES[67];
resid[36] = STATES[68] - (( CONSTANTS[161]*STATES[63])/CONSTANTS[148])*STATES[64];
resid[37] = STATES[69] - ALGEBRAIC[37]+STATES[68];
resid[38] = STATES[62] - ALGEBRAIC[36]+STATES[69];
resid[39] = RATES[27] - ALGEBRAIC[263];
resid[40] = RATES[28] - ALGEBRAIC[265];
resid[41] = RATES[29] - ALGEBRAIC[267];
resid[42] = RATES[30] - ALGEBRAIC[269];
resid[43] = RATES[31] - ALGEBRAIC[271];
resid[44] = RATES[32] - ALGEBRAIC[273];
resid[45] = RATES[33] - ALGEBRAIC[264];
resid[46] = RATES[34] - ALGEBRAIC[266];
resid[47] = RATES[35] - ALGEBRAIC[268];
resid[48] = RATES[36] - ALGEBRAIC[270];
resid[49] = RATES[37] - ALGEBRAIC[272];
resid[50] = RATES[38] - ALGEBRAIC[274];
resid[51] = RATES[39] - ALGEBRAIC[285];
resid[52] = RATES[40] - ALGEBRAIC[286];
resid[53] = RATES[41] - ALGEBRAIC[287];
resid[54] = RATES[42] - ALGEBRAIC[288];
resid[55] = RATES[43] - ALGEBRAIC[289];
resid[56] = RATES[44] - ALGEBRAIC[290];
resid[57] = RATES[48] - ALGEBRAIC[166] - ALGEBRAIC[163];
resid[58] = RATES[45] - ALGEBRAIC[163] - ALGEBRAIC[164];
resid[59] = RATES[46] - ALGEBRAIC[164] - ALGEBRAIC[165];
resid[60] = RATES[47] - ALGEBRAIC[165] - ALGEBRAIC[166];
resid[61] = RATES[49] - ALGEBRAIC[225];
resid[62] = RATES[50] - ALGEBRAIC[226];
resid[63] = RATES[51] - ALGEBRAIC[227];
resid[64] = RATES[52] - ALGEBRAIC[228];
resid[65] = RATES[53] - ALGEBRAIC[229];
resid[66] = RATES[57] - ALGEBRAIC[230];
resid[67] = RATES[54] - ALGEBRAIC[231];
resid[68] = RATES[55] - ALGEBRAIC[232];
resid[69] = RATES[56] - ALGEBRAIC[233];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[5] = STATES[4]/CONSTANTS[132];
ALGEBRAIC[6] = STATES[11]/CONSTANTS[139];
ALGEBRAIC[7] = STATES[1]/CONSTANTS[140];
ALGEBRAIC[8] = STATES[20]/CONSTANTS[170];
ALGEBRAIC[9] = STATES[5]/CONSTANTS[132];
ALGEBRAIC[10] = - arbitrary_log(STATES[23]/CONSTANTS[132], 10);
ALGEBRAIC[11] = STATES[4]+STATES[11]+STATES[3]+STATES[20]+STATES[1]+STATES[22]+ 2.00000*(STATES[52]+STATES[53]+STATES[54])+STATES[41]+STATES[42];
ALGEBRAIC[45] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33]);
ALGEBRAIC[47] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34]);
ALGEBRAIC[49] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35]);
ALGEBRAIC[51] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36]);
ALGEBRAIC[112] = STATES[11]+STATES[1]+STATES[33]+STATES[34]+STATES[35]+STATES[36]+STATES[37]+STATES[38];
ALGEBRAIC[158] = CONSTANTS[164]*STATES[48]+ CONSTANTS[163]*STATES[46]+ (CONSTANTS[163]+CONSTANTS[164])*STATES[47];
ALGEBRAIC[167] = CONSTANTS[165]*STATES[52]+ CONSTANTS[165]*STATES[53]+ CONSTANTS[166]*STATES[54];
ALGEBRAIC[168] = STATES[4]/CONSTANTS[168];
ALGEBRAIC[169] = STATES[20]/CONSTANTS[174];
ALGEBRAIC[170] = STATES[23]/CONSTANTS[176];
ALGEBRAIC[171] = STATES[25]/CONSTANTS[168];
ALGEBRAIC[172] = STATES[26]/CONSTANTS[168];
ALGEBRAIC[173] = STATES[24]/CONSTANTS[168];
ALGEBRAIC[210] = CONSTANTS[59]*(exp(ALGEBRAIC[185]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[180]+ALGEBRAIC[179])/( CONSTANTS[141]*CONSTANTS[142])));
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = VOI - floor(VOI/CONSTANTS[135])*CONSTANTS[135];
ALGEBRAIC[3] = (CONDVAR[2]>=0.00000&&CONDVAR[3]<=0.00000 ? 0.00000*8000.00 : 0.00000);
ALGEBRAIC[12] = ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[13] = ( CONSTANTS[141]*CONSTANTS[142])*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[14] = CONSTANTS[0]*(exp(ALGEBRAIC[12]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[13]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[15] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[62]*STATES[2]);
ALGEBRAIC[16] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[17] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[63]*STATES[3]);
ALGEBRAIC[22] = CONSTANTS[1]*(exp(( CONSTANTS[131]*ALGEBRAIC[16]+ALGEBRAIC[15])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[17]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[18] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[64]*STATES[7]);
ALGEBRAIC[19] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[65]*STATES[8]);
ALGEBRAIC[23] = CONSTANTS[2]*(exp((ALGEBRAIC[18]+ CONSTANTS[130]*ALGEBRAIC[17])/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[19]+ CONSTANTS[130]*ALGEBRAIC[17])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[20] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[66]*STATES[9]);
ALGEBRAIC[24] = CONSTANTS[3]*(exp(ALGEBRAIC[19]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[20]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[21] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[67]*STATES[10]);
ALGEBRAIC[25] = CONSTANTS[4]*(exp(ALGEBRAIC[20]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[21]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[26] = CONSTANTS[5]*(exp(ALGEBRAIC[21]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[19]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[4] = (CONDVAR[4]<0.00000 ? 1.00000*CONSTANTS[137] : CONDVAR[5]>0.00000&&CONDVAR[6]<0.00000 ? - 0.100000*CONSTANTS[137] : 0.00000*CONSTANTS[137]);
ALGEBRAIC[27] = ALGEBRAIC[4];
ALGEBRAIC[34] = (CONSTANTS[157]/( CONSTANTS[148]*CONSTANTS[149]*CONSTANTS[150]))*STATES[10]*(STATES[60]+CONSTANTS[148]);
ALGEBRAIC[35] = (CONSTANTS[157]/( CONSTANTS[148]*CONSTANTS[149]*CONSTANTS[150]))*STATES[9]*STATES[59];
ALGEBRAIC[36] = ALGEBRAIC[34]+ALGEBRAIC[35];
ALGEBRAIC[37] = CONSTANTS[161]*(exp(( CONSTANTS[162]*STATES[58])/CONSTANTS[148]) - 1.00000);
ALGEBRAIC[38] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[39] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[40] = CONSTANTS[6]*(exp(ALGEBRAIC[38]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[39]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[136] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[138] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[90]*STATES[21]);
ALGEBRAIC[139] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[91]*STATES[22]);
ALGEBRAIC[140] = CONSTANTS[39]*(exp((ALGEBRAIC[136]+ALGEBRAIC[138])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[139]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[137] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[141] = CONSTANTS[40]*(exp(ALGEBRAIC[139]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[137]+ALGEBRAIC[138])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[143] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[92]*STATES[12]);
ALGEBRAIC[144] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[93]*STATES[13]);
ALGEBRAIC[145] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[94]*STATES[14]);
ALGEBRAIC[151] = CONSTANTS[41]*(exp((ALGEBRAIC[143]+ALGEBRAIC[144])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[145]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[148] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]);
ALGEBRAIC[152] = CONSTANTS[42]*(exp(ALGEBRAIC[145]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[148]+ALGEBRAIC[144])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[146] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[95]*STATES[15]);
ALGEBRAIC[147] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[96]*STATES[16]);
ALGEBRAIC[153] = CONSTANTS[43]*(exp((ALGEBRAIC[148]+ALGEBRAIC[146])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[147]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[154] = CONSTANTS[44]*(exp(ALGEBRAIC[147]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[143]+ALGEBRAIC[146])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[149] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[98]*STATES[18]);
ALGEBRAIC[150] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[99]*STATES[19]);
ALGEBRAIC[155] = CONSTANTS[45]*(exp((ALGEBRAIC[146]+ALGEBRAIC[149])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[150]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[157] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[159] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[100]*STATES[45]);
ALGEBRAIC[162] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[103]*STATES[48]);
ALGEBRAIC[163] = CONSTANTS[47]*(exp(ALGEBRAIC[162]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[159]+ CONSTANTS[164]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[160] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[101]*STATES[46]);
ALGEBRAIC[164] = CONSTANTS[48]*(exp((ALGEBRAIC[159]+ CONSTANTS[163]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[160]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[161] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[102]*STATES[47]);
ALGEBRAIC[165] = CONSTANTS[49]*(exp((ALGEBRAIC[160]+ CONSTANTS[164]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[161]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[166] = CONSTANTS[50]*(exp(ALGEBRAIC[161]/( CONSTANTS[141]*CONSTANTS[142])) - exp((ALGEBRAIC[162]+ CONSTANTS[163]*ALGEBRAIC[157])/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[174] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[181] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[105]*STATES[50]);
ALGEBRAIC[189] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[97]*STATES[17]);
ALGEBRAIC[190] = ALGEBRAIC[189] - CONSTANTS[173];
ALGEBRAIC[199] = ALGEBRAIC[181]+ CONSTANTS[165]*ALGEBRAIC[174]+ALGEBRAIC[190];
ALGEBRAIC[183] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[107]*STATES[52]);
ALGEBRAIC[200] = ALGEBRAIC[183]+ALGEBRAIC[190];
ALGEBRAIC[214] = CONSTANTS[52]*(exp(ALGEBRAIC[199]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[200]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[219] = - CONSTANTS[165]*ALGEBRAIC[214];
ALGEBRAIC[186] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[109]*STATES[54]);
ALGEBRAIC[203] = ALGEBRAIC[186]+ALGEBRAIC[190];
ALGEBRAIC[175] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[187] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[110]*STATES[55]);
ALGEBRAIC[204] = ALGEBRAIC[187]+ CONSTANTS[166]*ALGEBRAIC[175]+ALGEBRAIC[190];
ALGEBRAIC[216] = CONSTANTS[56]*(exp(ALGEBRAIC[203]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[204]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[220] = CONSTANTS[166]*ALGEBRAIC[216];
ALGEBRAIC[176] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[113]*STATES[23]);
ALGEBRAIC[197] = ALGEBRAIC[181]+ALGEBRAIC[176]+ALGEBRAIC[190];
ALGEBRAIC[182] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[106]*STATES[51]);
ALGEBRAIC[198] = ALGEBRAIC[182]+ALGEBRAIC[190];
ALGEBRAIC[213] = CONSTANTS[53]*(exp(ALGEBRAIC[197]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[198]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[201] = ALGEBRAIC[183]+ALGEBRAIC[190];
ALGEBRAIC[184] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[108]*STATES[53]);
ALGEBRAIC[202] = ALGEBRAIC[184]+ CONSTANTS[167]*ALGEBRAIC[176]+ALGEBRAIC[190];
ALGEBRAIC[215] = CONSTANTS[54]*(exp(ALGEBRAIC[201]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[202]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[205] = ALGEBRAIC[187]+ CONSTANTS[167]*ALGEBRAIC[176]+ALGEBRAIC[190];
ALGEBRAIC[188] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[111]*STATES[56]);
ALGEBRAIC[206] = ALGEBRAIC[188]+ALGEBRAIC[190];
ALGEBRAIC[217] = CONSTANTS[57]*(exp(ALGEBRAIC[205]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[206]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[207] = ALGEBRAIC[188]+ALGEBRAIC[190];
ALGEBRAIC[185] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[112]*STATES[57]);
ALGEBRAIC[208] = ALGEBRAIC[176]+ALGEBRAIC[185]+ALGEBRAIC[190];
ALGEBRAIC[218] = CONSTANTS[58]*(exp(ALGEBRAIC[207]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[208]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[221] = (( CONSTANTS[167]*ALGEBRAIC[215] - ALGEBRAIC[213]) - CONSTANTS[167]*ALGEBRAIC[217])+ALGEBRAIC[218];
ALGEBRAIC[195] = ALGEBRAIC[184]+ALGEBRAIC[190];
ALGEBRAIC[177] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[115]*STATES[25]);
ALGEBRAIC[196] = ALGEBRAIC[177]+ALGEBRAIC[186]+ALGEBRAIC[190];
ALGEBRAIC[212] = CONSTANTS[55]*(exp(ALGEBRAIC[195]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[196]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[222] = ALGEBRAIC[212];
ALGEBRAIC[178] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[114]*STATES[26]);
ALGEBRAIC[180] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[104]*STATES[49]);
ALGEBRAIC[193] = ALGEBRAIC[180]+ALGEBRAIC[178]+ALGEBRAIC[190];
ALGEBRAIC[194] = ALGEBRAIC[181]+ALGEBRAIC[190];
ALGEBRAIC[211] = CONSTANTS[51]*(exp(ALGEBRAIC[193]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[194]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[223] = - ALGEBRAIC[211];
ALGEBRAIC[191] = ALGEBRAIC[185]+ALGEBRAIC[190];
ALGEBRAIC[179] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[116]*STATES[24]);
ALGEBRAIC[192] = ALGEBRAIC[180]+ALGEBRAIC[179]+ALGEBRAIC[190];
ALGEBRAIC[209] = CONSTANTS[59]*(exp(ALGEBRAIC[191]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[192]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[224] = ALGEBRAIC[209];
ALGEBRAIC[225] = ALGEBRAIC[209] - ALGEBRAIC[211];
ALGEBRAIC[226] = ALGEBRAIC[211] - ALGEBRAIC[214];
ALGEBRAIC[227] = ALGEBRAIC[213];
ALGEBRAIC[228] = ALGEBRAIC[214] - ALGEBRAIC[215];
ALGEBRAIC[229] = ALGEBRAIC[215] - ALGEBRAIC[212];
ALGEBRAIC[230] = ALGEBRAIC[218] - ALGEBRAIC[209];
ALGEBRAIC[231] = ALGEBRAIC[212] - ALGEBRAIC[216];
ALGEBRAIC[232] = ALGEBRAIC[216] - ALGEBRAIC[217];
ALGEBRAIC[233] = ALGEBRAIC[217] - ALGEBRAIC[218];
ALGEBRAIC[28] = - CONSTANTS[169]*ALGEBRAIC[27];
ALGEBRAIC[29] = CONSTANTS[171]*STATES[59];
ALGEBRAIC[32] = - ALGEBRAIC[28] - ALGEBRAIC[29];
ALGEBRAIC[234] = ALGEBRAIC[32];
ALGEBRAIC[30] = - CONSTANTS[175]*ALGEBRAIC[27];
ALGEBRAIC[31] = CONSTANTS[172]*STATES[60];
ALGEBRAIC[33] = - ALGEBRAIC[30] - ALGEBRAIC[31];
ALGEBRAIC[235] = ALGEBRAIC[33];
ALGEBRAIC[43] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[68]*STATES[1]);
ALGEBRAIC[41] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[53] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37]);
ALGEBRAIC[56] = ALGEBRAIC[43]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41]+ALGEBRAIC[53];
ALGEBRAIC[42] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[84] = ALGEBRAIC[42]+ALGEBRAIC[53];
ALGEBRAIC[82] = CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[110] = (ALGEBRAIC[82]==0.00000 ? CONSTANTS[7]*(exp(ALGEBRAIC[56]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[84]/( CONSTANTS[141]*CONSTANTS[142]))) : ((( CONSTANTS[7]*ALGEBRAIC[82])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[82]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[56]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[84]/( CONSTANTS[141]*CONSTANTS[142]))));
ALGEBRAIC[55] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38]);
ALGEBRAIC[57] = ALGEBRAIC[43]+ CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41]+ALGEBRAIC[55];
ALGEBRAIC[85] = ALGEBRAIC[42]+ALGEBRAIC[55];
ALGEBRAIC[83] = CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[111] = (ALGEBRAIC[83]==0.00000 ? CONSTANTS[8]*(exp(ALGEBRAIC[57]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[85]/( CONSTANTS[141]*CONSTANTS[142]))) : ((( CONSTANTS[8]*ALGEBRAIC[83])/( CONSTANTS[141]*CONSTANTS[142]))/(exp(ALGEBRAIC[83]/( CONSTANTS[141]*CONSTANTS[142])) - 1.00000))*(exp(ALGEBRAIC[57]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[85]/( CONSTANTS[141]*CONSTANTS[142]))));
ALGEBRAIC[261] = ALGEBRAIC[111]+ALGEBRAIC[110];
ALGEBRAIC[44] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27]);
ALGEBRAIC[76] = ALGEBRAIC[44];
ALGEBRAIC[104] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[117]*ALGEBRAIC[43];
ALGEBRAIC[254] = CONSTANTS[27]*(exp(ALGEBRAIC[76]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[104]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[46] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28]);
ALGEBRAIC[77] = ALGEBRAIC[46];
ALGEBRAIC[105] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[117]*ALGEBRAIC[43];
ALGEBRAIC[255] = CONSTANTS[29]*(exp(ALGEBRAIC[77]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[105]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[48] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29]);
ALGEBRAIC[78] = ALGEBRAIC[48];
ALGEBRAIC[106] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[117]*ALGEBRAIC[43];
ALGEBRAIC[256] = CONSTANTS[31]*(exp(ALGEBRAIC[78]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[106]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[50] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30]);
ALGEBRAIC[79] = ALGEBRAIC[50];
ALGEBRAIC[107] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[117]*ALGEBRAIC[43];
ALGEBRAIC[257] = CONSTANTS[28]*(exp(ALGEBRAIC[79]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[107]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[52] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31]);
ALGEBRAIC[80] = ALGEBRAIC[52];
ALGEBRAIC[108] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[117]*ALGEBRAIC[43];
ALGEBRAIC[258] = CONSTANTS[30]*(exp(ALGEBRAIC[80]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[108]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[54] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32]);
ALGEBRAIC[81] = ALGEBRAIC[54];
ALGEBRAIC[109] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[117]*ALGEBRAIC[43];
ALGEBRAIC[259] = CONSTANTS[32]*(exp(ALGEBRAIC[81]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[109]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[262] = (( (CONSTANTS[117])*ALGEBRAIC[257]+ CONSTANTS[117]*ALGEBRAIC[258]+ CONSTANTS[117]*ALGEBRAIC[256]+ CONSTANTS[117]*ALGEBRAIC[255]+ CONSTANTS[117]*ALGEBRAIC[254]+ CONSTANTS[117]*ALGEBRAIC[259]) - ALGEBRAIC[110]) - ALGEBRAIC[111];
ALGEBRAIC[58] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[86] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[236] = CONSTANTS[9]*(exp(ALGEBRAIC[58]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[86]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[64] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[92] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[242] = CONSTANTS[15]*(exp(ALGEBRAIC[64]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[92]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[68] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[70]*STATES[27])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[96] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[246] = CONSTANTS[19]*(exp(ALGEBRAIC[68]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[96]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[263] = ((- ALGEBRAIC[236] - ALGEBRAIC[242]) - ALGEBRAIC[246]) - ALGEBRAIC[254];
ALGEBRAIC[59] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[87] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[237] = CONSTANTS[12]*(exp(ALGEBRAIC[59]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[87]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[65] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[93] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[243] = CONSTANTS[17]*(exp(ALGEBRAIC[65]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[93]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[69] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[76]*STATES[33])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[97] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[247] = CONSTANTS[21]*(exp(ALGEBRAIC[69]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[97]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[264] = ((- ALGEBRAIC[237] - ALGEBRAIC[243]) - ALGEBRAIC[247]) - - ALGEBRAIC[254];
ALGEBRAIC[60] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[88] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[238] = CONSTANTS[10]*(exp(ALGEBRAIC[60]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[88]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[72] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[71]*STATES[28])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[100] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[250] = CONSTANTS[23]*(exp(ALGEBRAIC[72]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[100]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[265] = ((- ALGEBRAIC[238]+ALGEBRAIC[242]) - ALGEBRAIC[255]) - ALGEBRAIC[250];
ALGEBRAIC[61] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[89] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[239] = CONSTANTS[13]*(exp(ALGEBRAIC[61]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[89]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[73] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[77]*STATES[34])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[101] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[251] = CONSTANTS[25]*(exp(ALGEBRAIC[73]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[101]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[266] = ((- ALGEBRAIC[239]+ALGEBRAIC[243]) - - ALGEBRAIC[255]) - ALGEBRAIC[251];
ALGEBRAIC[62] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[72]*STATES[29])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[90] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[240] = CONSTANTS[11]*(exp(ALGEBRAIC[62]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[90]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[267] = ((- ALGEBRAIC[240]+ALGEBRAIC[246]) - ALGEBRAIC[256])+ALGEBRAIC[250];
ALGEBRAIC[63] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[78]*STATES[35])+ CONSTANTS[118]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[91] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[122]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[241] = CONSTANTS[14]*(exp(ALGEBRAIC[63]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[91]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[268] = ((- ALGEBRAIC[241]+ALGEBRAIC[247]) - - ALGEBRAIC[256])+ALGEBRAIC[251];
ALGEBRAIC[66] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[94] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[244] = CONSTANTS[16]*(exp(ALGEBRAIC[66]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[94]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[70] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[73]*STATES[30])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[98] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[248] = CONSTANTS[20]*(exp(ALGEBRAIC[70]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[98]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[269] = ((- ALGEBRAIC[248]+ALGEBRAIC[236]) - ALGEBRAIC[244]) - ALGEBRAIC[257];
ALGEBRAIC[67] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[119]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[95] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[123]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[245] = CONSTANTS[18]*(exp(ALGEBRAIC[67]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[95]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[71] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[79]*STATES[36])+ CONSTANTS[120]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[99] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[124]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[249] = CONSTANTS[22]*(exp(ALGEBRAIC[71]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[99]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[270] = ((- ALGEBRAIC[249]+ALGEBRAIC[237]) - ALGEBRAIC[245]) - - ALGEBRAIC[257];
ALGEBRAIC[74] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[74]*STATES[31])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[102] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[75]*STATES[32])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[252] = CONSTANTS[24]*(exp(ALGEBRAIC[74]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[102]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[271] = ((ALGEBRAIC[244]+ALGEBRAIC[238]) - ALGEBRAIC[258]) - ALGEBRAIC[252];
ALGEBRAIC[75] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[80]*STATES[37])+ CONSTANTS[121]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[103] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[81]*STATES[38])+ CONSTANTS[125]*CONSTANTS[143]*ALGEBRAIC[41];
ALGEBRAIC[253] = CONSTANTS[26]*(exp(ALGEBRAIC[75]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[103]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[272] = ((ALGEBRAIC[245]+ALGEBRAIC[239]) - ((- ALGEBRAIC[258]+ALGEBRAIC[110]) - ALGEBRAIC[110])) - ALGEBRAIC[253];
ALGEBRAIC[273] = ((ALGEBRAIC[248]+ALGEBRAIC[240]) - ALGEBRAIC[259])+ALGEBRAIC[252];
ALGEBRAIC[274] = ((ALGEBRAIC[249]+ALGEBRAIC[241]) - ((- ALGEBRAIC[259]+ALGEBRAIC[111]) - ALGEBRAIC[111]))+ALGEBRAIC[253];
ALGEBRAIC[121] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[85]*STATES[42]);
ALGEBRAIC[127] = ALGEBRAIC[121];
ALGEBRAIC[114] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[69]*STATES[11]);
ALGEBRAIC[122] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[86]*STATES[43]);
ALGEBRAIC[133] = ALGEBRAIC[122]+ALGEBRAIC[114];
ALGEBRAIC[278] = CONSTANTS[36]*(exp(ALGEBRAIC[127]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[133]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[281] = ALGEBRAIC[278];
ALGEBRAIC[115] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[61]*STATES[4]);
ALGEBRAIC[119] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[83]*STATES[40]);
ALGEBRAIC[125] = ALGEBRAIC[119]+ALGEBRAIC[115];
ALGEBRAIC[120] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[84]*STATES[41]);
ALGEBRAIC[131] = ALGEBRAIC[120];
ALGEBRAIC[276] = CONSTANTS[34]*(exp(ALGEBRAIC[125]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[131]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[282] = - ALGEBRAIC[276];
ALGEBRAIC[116] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[89]*STATES[6]);
ALGEBRAIC[128] = ALGEBRAIC[122]+ CONSTANTS[126]*ALGEBRAIC[116];
ALGEBRAIC[123] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44]);
ALGEBRAIC[134] = ALGEBRAIC[123];
ALGEBRAIC[279] = CONSTANTS[37]*(exp(ALGEBRAIC[128]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[134]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[283] = - CONSTANTS[126]*ALGEBRAIC[279];
ALGEBRAIC[118] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39]);
ALGEBRAIC[124] = ALGEBRAIC[118];
ALGEBRAIC[117] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[88]*STATES[5]);
ALGEBRAIC[130] = CONSTANTS[127]*ALGEBRAIC[117]+ALGEBRAIC[119];
ALGEBRAIC[275] = CONSTANTS[33]*(exp(ALGEBRAIC[124]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[130]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[284] = CONSTANTS[127]*ALGEBRAIC[275];
ALGEBRAIC[113] = STATES[0]/CONSTANTS[133];
ALGEBRAIC[129] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[87]*STATES[44])+ CONSTANTS[128]*CONSTANTS[143]*ALGEBRAIC[113];
ALGEBRAIC[135] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[82]*STATES[39])+ CONSTANTS[129]*CONSTANTS[143]*ALGEBRAIC[113];
ALGEBRAIC[280] = CONSTANTS[38]*(exp(ALGEBRAIC[129]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[135]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[285] = ALGEBRAIC[280] - ALGEBRAIC[275];
ALGEBRAIC[286] = ALGEBRAIC[275] - ALGEBRAIC[276];
ALGEBRAIC[126] = ALGEBRAIC[120];
ALGEBRAIC[132] = ALGEBRAIC[121];
ALGEBRAIC[277] = CONSTANTS[35]*(exp(ALGEBRAIC[126]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[132]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[287] = ALGEBRAIC[276] - ALGEBRAIC[277];
ALGEBRAIC[288] = - ALGEBRAIC[278]+ALGEBRAIC[277];
ALGEBRAIC[289] = - ALGEBRAIC[279]+ALGEBRAIC[278];
ALGEBRAIC[290] = - ALGEBRAIC[280]+ALGEBRAIC[279];
ALGEBRAIC[156] = CONSTANTS[141]*CONSTANTS[142]*log( CONSTANTS[60]*STATES[20]);
ALGEBRAIC[291] = CONSTANTS[46]*exp(ALGEBRAIC[162]/( CONSTANTS[141]*CONSTANTS[142]))*(exp(ALGEBRAIC[156]/( CONSTANTS[141]*CONSTANTS[142])) - exp(ALGEBRAIC[157]/( CONSTANTS[141]*CONSTANTS[142])));
ALGEBRAIC[292] = (( CONSTANTS[164]*ALGEBRAIC[163] - CONSTANTS[163]*ALGEBRAIC[164]) - CONSTANTS[164]*ALGEBRAIC[165])+ CONSTANTS[163]*ALGEBRAIC[166];
ALGEBRAIC[260] = CONSTANTS[143]*(((((((((((((((((((((((((((((((((((( CONSTANTS[124]*ALGEBRAIC[247] - CONSTANTS[120]*ALGEBRAIC[247])+ CONSTANTS[125]*ALGEBRAIC[251]) - CONSTANTS[121]*ALGEBRAIC[251])+ CONSTANTS[123]*ALGEBRAIC[245]) - CONSTANTS[119]*ALGEBRAIC[245])+ CONSTANTS[124]*ALGEBRAIC[249]) - CONSTANTS[120]*ALGEBRAIC[249])+ CONSTANTS[125]*ALGEBRAIC[253]) - CONSTANTS[121]*ALGEBRAIC[253])+ CONSTANTS[123]*ALGEBRAIC[243]) - CONSTANTS[119]*ALGEBRAIC[243])+ CONSTANTS[122]*ALGEBRAIC[241]) - CONSTANTS[118]*ALGEBRAIC[241])+ CONSTANTS[122]*ALGEBRAIC[239]) - CONSTANTS[118]*ALGEBRAIC[239])+ CONSTANTS[122]*ALGEBRAIC[237]) - CONSTANTS[118]*ALGEBRAIC[237])+ CONSTANTS[125]*ALGEBRAIC[252]) - CONSTANTS[121]*ALGEBRAIC[252])+ CONSTANTS[124]*ALGEBRAIC[248]) - CONSTANTS[120]*ALGEBRAIC[248])+ CONSTANTS[123]*ALGEBRAIC[244]) - CONSTANTS[119]*ALGEBRAIC[244])+ CONSTANTS[125]*ALGEBRAIC[250]) - CONSTANTS[121]*ALGEBRAIC[250])+ CONSTANTS[124]*ALGEBRAIC[246]) - CONSTANTS[120]*ALGEBRAIC[246])+ CONSTANTS[123]*ALGEBRAIC[242]) - CONSTANTS[119]*ALGEBRAIC[242])+ CONSTANTS[122]*ALGEBRAIC[240]) - CONSTANTS[118]*ALGEBRAIC[240])+ CONSTANTS[122]*ALGEBRAIC[238]) - CONSTANTS[118]*ALGEBRAIC[238])+ CONSTANTS[122]*ALGEBRAIC[236]) - CONSTANTS[118]*ALGEBRAIC[236]) - ( CONSTANTS[147]*ALGEBRAIC[110]+ CONSTANTS[147]*ALGEBRAIC[111]));
ALGEBRAIC[293] = CONSTANTS[143]*( - CONSTANTS[128]*ALGEBRAIC[280]+ CONSTANTS[129]*ALGEBRAIC[280]);
ALGEBRAIC[142] = - CONSTANTS[147]*CONSTANTS[143]*ALGEBRAIC[141];
ALGEBRAIC[2] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<=0.00000 ? ( 0.00000*0.0300000*CONSTANTS[133])/CONSTANTS[136] : 0.00000);
ALGEBRAIC[294] = ALGEBRAIC[260]+ALGEBRAIC[293]+ALGEBRAIC[142]+ALGEBRAIC[2];
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[16] = 1.0;
SI[17] = 1.0;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[21] = 1.0;
SI[22] = 1.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[25] = 1.0;
SI[26] = 1.0;
SI[27] = 1.0;
SI[28] = 1.0;
SI[29] = 1.0;
SI[30] = 1.0;
SI[31] = 1.0;
SI[32] = 1.0;
SI[33] = 1.0;
SI[34] = 1.0;
SI[35] = 1.0;
SI[36] = 1.0;
SI[37] = 1.0;
SI[38] = 1.0;
SI[39] = 1.0;
SI[40] = 1.0;
SI[41] = 1.0;
SI[42] = 1.0;
SI[43] = 1.0;
SI[44] = 1.0;
SI[45] = 1.0;
SI[46] = 1.0;
SI[47] = 1.0;
SI[48] = 1.0;
SI[49] = 1.0;
SI[50] = 1.0;
SI[51] = 1.0;
SI[52] = 1.0;
SI[53] = 1.0;
SI[54] = 1.0;
SI[55] = 1.0;
SI[56] = 1.0;
SI[57] = 1.0;
SI[62] = 0.0;
SI[58] = 1.0;
SI[59] = 1.0;
SI[60] = 1.0;
SI[61] = 1.0;
SI[63] = 0.0;
SI[64] = 0.0;
SI[65] = 0.0;
SI[66] = 0.0;
SI[67] = 0.0;
SI[68] = 0.0;
SI[69] = 0.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = ALGEBRAIC[1] - 0.300000;
CONDVAR[1] = ALGEBRAIC[1] - (0.300000+CONSTANTS[136]);
CONDVAR[2] = ALGEBRAIC[1] - 0.300000;
CONDVAR[3] = ALGEBRAIC[1] - (0.300000+CONSTANTS[136]);
CONDVAR[4] = VOI - 0.0600000;
CONDVAR[5] = VOI - 0.500000;
CONDVAR[6] = VOI - 1.10000;
CONDVAR[7] = STATES[64] - 0.00000;
}