# Size of variable arrays: sizeAlgebraic = 303 sizeStates = 62 sizeConstants = 177 from math import * from numpy import * def createLegends(): legend_states = [""] * sizeStates legend_rates = [""] * sizeStates legend_algebraic = [""] * sizeAlgebraic legend_voi = "" legend_constants = [""] * sizeConstants legend_constants[0] = "kappa_leak in component BG_parameters (fmol_per_sec)" legend_constants[1] = "kappa_R_TRPNCa in component BG_parameters (fmol_per_sec)" legend_constants[2] = "kappa_R_BU in component BG_parameters (fmol_per_sec)" legend_constants[3] = "kappa_R_UW in component BG_parameters (fmol_per_sec)" legend_constants[4] = "kappa_R_WS in component BG_parameters (fmol_per_sec)" legend_constants[5] = "kappa_R_SU in component BG_parameters (fmol_per_sec)" legend_constants[6] = "kappa_diadDf in component BG_parameters (fmol_per_sec)" legend_constants[7] = "kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec)" legend_constants[8] = "kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec)" legend_constants[9] = "kappa_d000 in component BG_parameters (fmol_per_sec)" legend_constants[10] = "kappa_d010 in component BG_parameters (fmol_per_sec)" legend_constants[11] = "kappa_d020 in component BG_parameters (fmol_per_sec)" legend_constants[12] = "kappa_d001 in component BG_parameters (fmol_per_sec)" legend_constants[13] = "kappa_d011 in component BG_parameters (fmol_per_sec)" legend_constants[14] = "kappa_d021 in component BG_parameters (fmol_per_sec)" legend_constants[15] = "kappa_f1_000 in component BG_parameters (fmol_per_sec)" legend_constants[16] = "kappa_f1_100 in component BG_parameters (fmol_per_sec)" legend_constants[17] = "kappa_f1_001 in component BG_parameters (fmol_per_sec)" legend_constants[18] = "kappa_f1_101 in component BG_parameters (fmol_per_sec)" legend_constants[19] = "kappa_f2_000 in component BG_parameters (fmol_per_sec)" legend_constants[20] = "kappa_f2_100 in component BG_parameters (fmol_per_sec)" legend_constants[21] = "kappa_f2_001 in component BG_parameters (fmol_per_sec)" legend_constants[22] = "kappa_f2_101 in component BG_parameters (fmol_per_sec)" legend_constants[23] = "kappa_f3_010 in component BG_parameters (fmol_per_sec)" legend_constants[24] = "kappa_f3_110 in component BG_parameters (fmol_per_sec)" legend_constants[25] = "kappa_f3_011 in component BG_parameters (fmol_per_sec)" legend_constants[26] = "kappa_f3_111 in component BG_parameters (fmol_per_sec)" legend_constants[27] = "kappa_fCa000 in component BG_parameters (fmol_per_sec)" legend_constants[28] = "kappa_fCa100 in component BG_parameters (fmol_per_sec)" legend_constants[29] = "kappa_fCa010 in component BG_parameters (fmol_per_sec)" legend_constants[30] = "kappa_fCa110 in component BG_parameters (fmol_per_sec)" legend_constants[31] = "kappa_fCa020 in component BG_parameters (fmol_per_sec)" legend_constants[32] = "kappa_fCa120 in component BG_parameters (fmol_per_sec)" legend_constants[33] = "kappa_1_NCX in component BG_parameters (fmol_per_sec)" legend_constants[34] = "kappa_2_NCX in component BG_parameters (fmol_per_sec)" legend_constants[35] = "kappa_3_NCX in component BG_parameters (fmol_per_sec)" legend_constants[36] = "kappa_4_NCX in component BG_parameters (fmol_per_sec)" legend_constants[37] = "kappa_5_NCX in component BG_parameters (fmol_per_sec)" legend_constants[38] = "kappa_6_NCX in component BG_parameters (fmol_per_sec)" legend_constants[39] = "kappa_pCa_R1 in component BG_parameters (fmol_per_sec)" legend_constants[40] = "kappa_pCa_R2 in component BG_parameters (fmol_per_sec)" legend_constants[41] = "kappa_PLBph1 in component BG_parameters (fmol_per_sec)" legend_constants[42] = "kappa_PLBph2 in component BG_parameters (fmol_per_sec)" legend_constants[43] = "kappa_PLBd1 in component BG_parameters (fmol_per_sec)" legend_constants[44] = "kappa_PLBd2 in component BG_parameters (fmol_per_sec)" legend_constants[45] = "kappa_Inh in component BG_parameters (fmol_per_sec)" legend_constants[46] = "kappa_RyR in component BG_parameters (fmol_per_sec)" legend_constants[47] = "kappa_OC in component BG_parameters (fmol_per_sec)" legend_constants[48] = "kappa_CCI in component BG_parameters (fmol_per_sec)" legend_constants[49] = "kappa_CII in component BG_parameters (fmol_per_sec)" legend_constants[50] = "kappa_IO in component BG_parameters (fmol_per_sec)" legend_constants[51] = "kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec)" legend_constants[52] = "kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec)" legend_constants[53] = "kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec)" legend_constants[54] = "kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec)" legend_constants[55] = "kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec)" legend_constants[56] = "kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec)" legend_constants[57] = "kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec)" legend_constants[58] = "kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec)" legend_constants[59] = "kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec)" legend_constants[60] = "K_Ca_SR in component BG_parameters (per_fmol)" legend_constants[61] = "K_Ca_i in component BG_parameters (per_fmol)" legend_constants[62] = "K_TRPN in component BG_parameters (per_fmol)" legend_constants[63] = "K_Ca_TRPN in component BG_parameters (per_fmol)" legend_constants[64] = "K_B_XB in component BG_parameters (per_fmol)" legend_constants[65] = "K_U_XB in component BG_parameters (per_fmol)" legend_constants[66] = "K_W_XB in component BG_parameters (per_fmol)" legend_constants[67] = "K_S_XB in component BG_parameters (per_fmol)" legend_constants[68] = "K_Ca_D in component BG_parameters (per_fmol)" legend_constants[69] = "K_Ca_o in component BG_parameters (per_fmol)" legend_constants[70] = "K_000_LCC in component BG_parameters (per_fmol)" legend_constants[71] = "K_010_LCC in component BG_parameters (per_fmol)" legend_constants[72] = "K_020_LCC in component BG_parameters (per_fmol)" legend_constants[73] = "K_100_LCC in component BG_parameters (per_fmol)" legend_constants[74] = "K_110_LCC in component BG_parameters (per_fmol)" legend_constants[75] = "K_120_LCC in component BG_parameters (per_fmol)" legend_constants[76] = "K_001_LCC in component BG_parameters (per_fmol)" legend_constants[77] = "K_011_LCC in component BG_parameters (per_fmol)" legend_constants[78] = "K_021_LCC in component BG_parameters (per_fmol)" legend_constants[79] = "K_101_LCC in component BG_parameters (per_fmol)" legend_constants[80] = "K_111_LCC in component BG_parameters (per_fmol)" legend_constants[81] = "K_121_LCC in component BG_parameters (per_fmol)" legend_constants[82] = "K_1_NCX in component BG_parameters (per_fmol)" legend_constants[83] = "K_2_NCX in component BG_parameters (per_fmol)" legend_constants[84] = "K_3_NCX in component BG_parameters (per_fmol)" legend_constants[85] = "K_4_NCX in component BG_parameters (per_fmol)" legend_constants[86] = "K_5_NCX in component BG_parameters (per_fmol)" legend_constants[87] = "K_6_NCX in component BG_parameters (per_fmol)" legend_constants[88] = "K_Na_i in component BG_parameters (per_fmol)" legend_constants[89] = "K_Na_o in component BG_parameters (per_fmol)" legend_constants[90] = "K_pCa in component BG_parameters (per_fmol)" legend_constants[91] = "K_pCa_Ca in component BG_parameters (per_fmol)" legend_constants[92] = "K_PLB in component BG_parameters (per_fmol)" legend_constants[93] = "K_PKACI in component BG_parameters (per_fmol)" legend_constants[94] = "K_PLB_PKACI in component BG_parameters (per_fmol)" legend_constants[95] = "K_PP1 in component BG_parameters (per_fmol)" legend_constants[96] = "K_PLBp_PP1 in component BG_parameters (per_fmol)" legend_constants[97] = "K_PLBp in component BG_parameters (per_fmol)" legend_constants[98] = "K_Ip in component BG_parameters (per_fmol)" legend_constants[99] = "K_Ip_PP1 in component BG_parameters (per_fmol)" legend_constants[100] = "K_C_RyR in component BG_parameters (per_fmol)" legend_constants[101] = "K_CI_RyR in component BG_parameters (per_fmol)" legend_constants[102] = "K_I_RyR in component BG_parameters (per_fmol)" legend_constants[103] = "K_O_RyR in component BG_parameters (per_fmol)" legend_constants[104] = "K_P1_SERCA in component BG_parameters (per_fmol)" legend_constants[105] = "K_P2_SERCA in component BG_parameters (per_fmol)" legend_constants[106] = "K_P2a_SERCA in component BG_parameters (per_fmol)" legend_constants[107] = "K_P4_SERCA in component BG_parameters (per_fmol)" legend_constants[108] = "K_P5_SERCA in component BG_parameters (per_fmol)" legend_constants[109] = "K_P6_SERCA in component BG_parameters (per_fmol)" legend_constants[110] = "K_P8_SERCA in component BG_parameters (per_fmol)" legend_constants[111] = "K_P9_SERCA in component BG_parameters (per_fmol)" legend_constants[112] = "K_P10_SERCA in component BG_parameters (per_fmol)" legend_constants[113] = "K_H in component BG_parameters (per_fmol)" legend_constants[114] = "K_MgATP in component BG_parameters (per_fmol)" legend_constants[115] = "K_MgADP in component BG_parameters (per_fmol)" legend_constants[116] = "K_Pi in component BG_parameters (per_fmol)" legend_constants[117] = "z_rCa in component BG_parameters (dimensionless)" legend_constants[118] = "z_fd in component BG_parameters (dimensionless)" legend_constants[119] = "z_ff1 in component BG_parameters (dimensionless)" legend_constants[120] = "z_ff2 in component BG_parameters (dimensionless)" legend_constants[121] = "z_ff3 in component BG_parameters (dimensionless)" legend_constants[122] = "z_rd in component BG_parameters (dimensionless)" legend_constants[123] = "z_rf1 in component BG_parameters (dimensionless)" legend_constants[124] = "z_rf2 in component BG_parameters (dimensionless)" legend_constants[125] = "z_rf3 in component BG_parameters (dimensionless)" legend_constants[126] = "nNa_o_NCX in component BG_parameters (dimensionless)" legend_constants[127] = "nNa_i_NCX in component BG_parameters (dimensionless)" legend_constants[128] = "zf_NCX in component BG_parameters (dimensionless)" legend_constants[129] = "zr_NCX in component BG_parameters (dimensionless)" legend_constants[130] = "n_Tm in component BG_parameters (dimensionless)" legend_constants[131] = "n_CaTRPN in component BG_parameters (dimensionless)" legend_voi = "time in component environment (second)" legend_constants[132] = "vol_myo in component environment (pL)" legend_constants[133] = "C_m in component environment (fF)" legend_states[0] = "q_membrane in component environment (fC)" legend_states[1] = "q_Ca_D in component environment (fmol)" legend_states[2] = "q_TRPN in component environment (fmol)" legend_states[3] = "q_Ca_TRPN in component environment (fmol)" legend_states[4] = "q_Ca_i in component environment (fmol)" legend_states[5] = "q_Na_i in component environment (fmol)" legend_states[6] = "q_Na_o in component environment (fmol)" legend_states[7] = "q_B_XB in component environment (fmol)" legend_states[8] = "q_U_XB in component environment (fmol)" legend_states[9] = "q_W_XB in component environment (fmol)" legend_states[10] = "q_S_XB in component environment (fmol)" legend_states[11] = "q_Ca_o in component environment (fmol)" legend_states[12] = "q_PLB in component environment (fmol)" legend_states[13] = "q_PKACI in component environment (fmol)" legend_states[14] = "q_PLB_PKACI in component environment (fmol)" legend_states[15] = "q_PP1 in component environment (fmol)" legend_states[16] = "q_PLBp_PP1 in component environment (fmol)" legend_states[17] = "q_PLBp in component environment (fmol)" legend_constants[134] = "q_PLBp_baseline in component environment (fmol)" legend_states[18] = "q_Ip in component environment (fmol)" legend_states[19] = "q_Ip_PP1 in component environment (fmol)" legend_states[20] = "q_Ca_SR in component environment (fmol)" legend_states[21] = "q_pCa in component environment (fmol)" legend_states[22] = "q_pCa_Ca in component environment (fmol)" legend_states[23] = "q_H in component environment (fmol)" legend_states[24] = "q_Pi in component environment (fmol)" legend_states[25] = "q_MgADP in component environment (fmol)" legend_states[26] = "q_MgATP in component environment (fmol)" legend_states[27] = "q_000_LCC in component environment (fmol)" legend_states[28] = "q_010_LCC in component environment (fmol)" legend_states[29] = "q_020_LCC in component environment (fmol)" legend_states[30] = "q_100_LCC in component environment (fmol)" legend_states[31] = "q_110_LCC in component environment (fmol)" legend_states[32] = "q_120_LCC in component environment (fmol)" legend_states[33] = "q_001_LCC in component environment (fmol)" legend_states[34] = "q_011_LCC in component environment (fmol)" legend_states[35] = "q_021_LCC in component environment (fmol)" legend_states[36] = "q_101_LCC in component environment (fmol)" legend_states[37] = "q_111_LCC in component environment (fmol)" legend_states[38] = "q_121_LCC in component environment (fmol)" legend_states[39] = "q_P1_NCX in component environment (fmol)" legend_states[40] = "q_P2_NCX in component environment (fmol)" legend_states[41] = "q_P3_NCX in component environment (fmol)" legend_states[42] = "q_P4_NCX in component environment (fmol)" legend_states[43] = "q_P5_NCX in component environment (fmol)" legend_states[44] = "q_P6_NCX in component environment (fmol)" legend_states[45] = "q_C_RyR in component environment (fmol)" legend_states[46] = "q_CI_RyR in component environment (fmol)" legend_states[47] = "q_I_RyR in component environment (fmol)" legend_states[48] = "q_O_RyR in component environment (fmol)" legend_states[49] = "q_P1_SERCA in component environment (fmol)" legend_states[50] = "q_P2_SERCA in component environment (fmol)" legend_states[51] = "q_P2a_SERCA in component environment (fmol)" legend_states[52] = "q_P4_SERCA in component environment (fmol)" legend_states[53] = "q_P5_SERCA in component environment (fmol)" legend_states[54] = "q_P6_SERCA in component environment (fmol)" legend_states[55] = "q_P8_SERCA in component environment (fmol)" legend_states[56] = "q_P9_SERCA in component environment (fmol)" legend_states[57] = "q_P10_SERCA in component environment (fmol)" legend_algebraic[54] = "v_leak in component Ca_leak (fmol_per_sec)" legend_algebraic[66] = "v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec)" legend_algebraic[73] = "v_R_BU in component crossbridge_TRPN (fmol_per_sec)" legend_algebraic[75] = "v_R_UW in component crossbridge_TRPN (fmol_per_sec)" legend_algebraic[77] = "v_R_WS in component crossbridge_TRPN (fmol_per_sec)" legend_algebraic[79] = "v_R_SU in component crossbridge_TRPN (fmol_per_sec)" legend_algebraic[59] = "T_total in component crossbridge_TRPN (kPa)" legend_states[58] = "SL in component crossbridge_TRPN (metre)" legend_algebraic[72] = "v_diadDf in component diadCaDiffusion (fmol_per_sec)" legend_algebraic[182] = "v_Ca_D_LCC in component LCC (fmol_per_sec)" legend_algebraic[181] = "v_Ca_o_LCC in component LCC (fmol_per_sec)" legend_algebraic[180] = "I_mem_LCC in component LCC (fA)" legend_algebraic[221] = "v_Ca_o_NCX in component NCX (fmol_per_sec)" legend_algebraic[222] = "v_Ca_i_NCX in component NCX (fmol_per_sec)" legend_algebraic[214] = "v_Na_o_NCX in component NCX (fmol_per_sec)" legend_algebraic[203] = "v_Na_i_NCX in component NCX (fmol_per_sec)" legend_algebraic[220] = "I_mem_NCX in component NCX (fA)" legend_algebraic[227] = "v_pCa_R1 in component pCa (fmol_per_sec)" legend_algebraic[228] = "v_pCa_R2 in component pCa (fmol_per_sec)" legend_algebraic[230] = "I_mem_pCa in component pCa (fA)" legend_algebraic[55] = "v_PLBph1 in component PLB (fmol_per_sec)" legend_algebraic[60] = "v_PLBph2 in component PLB (fmol_per_sec)" legend_algebraic[62] = "v_PLBd1 in component PLB (fmol_per_sec)" legend_algebraic[65] = "v_PLBd2 in component PLB (fmol_per_sec)" legend_algebraic[56] = "v_Inh in component PLB (fmol_per_sec)" legend_algebraic[237] = "v_RyR in component RyR (fmol_per_sec)" legend_algebraic[246] = "v_RyRgate_Ca_D in component RyR (fmol_per_sec)" legend_algebraic[289] = "v_Ca_i_SERCA in component SERCA (fmol_per_sec)" legend_algebraic[295] = "v_Ca_SR_SERCA in component SERCA (fmol_per_sec)" legend_algebraic[300] = "v_H_SERCA in component SERCA (fmol_per_sec)" legend_algebraic[286] = "v_MgADP_SERCA in component SERCA (fmol_per_sec)" legend_algebraic[283] = "v_MgATP_SERCA in component SERCA (fmol_per_sec)" legend_algebraic[281] = "v_Pi_SERCA in component SERCA (fmol_per_sec)" legend_algebraic[1] = "V_m in component environment (volt)" legend_algebraic[28] = "I_stim in component environment (fA)" legend_constants[135] = "stimPeriod in component environment (second)" legend_constants[136] = "stimDuration in component environment (second)" legend_algebraic[11] = "tPeriod in component environment (second)" legend_algebraic[33] = "v_Ca_stim in component environment (fmol_per_sec)" legend_algebraic[232] = "sum_I in component environment (fA)" legend_algebraic[2] = "dSL in component environment (m_per_s)" legend_constants[137] = "V_max in component environment (m_per_s)" legend_constants[138] = "vol_cell in component environment (pL)" legend_constants[170] = "vol_SR in component environment (pL)" legend_constants[139] = "vol_o in component environment (pL)" legend_constants[140] = "vol_D in component environment (pL)" legend_algebraic[3] = "cCai in component environment (mM)" legend_algebraic[4] = "cCao in component environment (mM)" legend_algebraic[5] = "cCaSR in component environment (mM)" legend_algebraic[6] = "cCaD in component environment (mM)" legend_algebraic[7] = "cNai in component environment (mM)" legend_algebraic[8] = "pHi in component environment (dimensionless)" legend_algebraic[9] = "Ca_tot in component environment (fmol)" legend_constants[141] = "R in component constants (J_per_K_per_mol)" legend_constants[142] = "T in component constants (kelvin)" legend_constants[143] = "F in component constants (C_per_mol)" legend_constants[144] = "zNa in component ion_valences (dimensionless)" legend_constants[145] = "zCl in component ion_valences (dimensionless)" legend_constants[146] = "zK in component ion_valences (dimensionless)" legend_constants[147] = "zCa in component ion_valences (dimensionless)" legend_algebraic[38] = "mu_Ca_SR in component Ca_leak (J_per_mol)" legend_algebraic[48] = "mu_Ca_i in component Ca_leak (J_per_mol)" legend_algebraic[58] = "mu_TRPN in component crossbridge_TRPN (J_per_mol)" legend_algebraic[61] = "mu_Ca_i in component crossbridge_TRPN (J_per_mol)" legend_algebraic[63] = "mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol)" legend_algebraic[64] = "mu_B_XB in component crossbridge_TRPN (J_per_mol)" legend_algebraic[67] = "mu_U_XB in component crossbridge_TRPN (J_per_mol)" legend_algebraic[69] = "mu_W_XB in component crossbridge_TRPN (J_per_mol)" legend_algebraic[71] = "mu_S_XB in component crossbridge_TRPN (J_per_mol)" legend_constants[148] = "SL_0 in component crossbridge_TRPN (metre)" legend_constants[149] = "q_MS in component crossbridge_TRPN (fmol)" legend_constants[150] = "r_s in component crossbridge_TRPN (dimensionless)" legend_constants[151] = "r_w in component crossbridge_TRPN (dimensionless)" legend_constants[152] = "A_eff in component crossbridge_TRPN (dimensionless)" legend_constants[153] = "phi in component crossbridge_TRPN (dimensionless)" legend_constants[154] = "k_uw in component crossbridge_TRPN (per_sec)" legend_constants[155] = "k_ws in component crossbridge_TRPN (per_sec)" legend_states[59] = "G_w in component crossbridge_TRPN (metre)" legend_states[60] = "G_s in component crossbridge_TRPN (metre)" legend_constants[171] = "c_w in component crossbridge_TRPN (per_sec)" legend_constants[172] = "c_s in component crossbridge_TRPN (per_sec)" legend_constants[169] = "A_w in component crossbridge_TRPN (dimensionless)" legend_constants[175] = "A_s in component crossbridge_TRPN (dimensionless)" legend_algebraic[41] = "nu_1 in component crossbridge_TRPN (m_per_s)" legend_algebraic[30] = "nu_2 in component crossbridge_TRPN (m_per_s)" legend_algebraic[35] = "nu_3 in component crossbridge_TRPN (m_per_s)" legend_algebraic[42] = "nu_4 in component crossbridge_TRPN (m_per_s)" legend_algebraic[31] = "nu_5 in component crossbridge_TRPN (m_per_s)" legend_algebraic[36] = "nu_6 in component crossbridge_TRPN (m_per_s)" legend_algebraic[51] = "v_1 in component crossbridge_TRPN (m_per_s)" legend_algebraic[26] = "v_2 in component crossbridge_TRPN (m_per_s)" legend_algebraic[52] = "v_4 in component crossbridge_TRPN (m_per_s)" legend_constants[156] = "v_to_mu in component crossbridge_TRPN (Js_per_m2)" legend_algebraic[32] = "T_active in component crossbridge_TRPN (kPa)" legend_constants[157] = "T_ref in component crossbridge_TRPN (kPa)" legend_algebraic[10] = "F_T_S in component crossbridge_TRPN (kPa)" legend_algebraic[27] = "F_T_W in component crossbridge_TRPN (kPa)" legend_states[61] = "Cdd in component crossbridge_TRPN (metre)" legend_algebraic[43] = "eta in component crossbridge_TRPN (second)" legend_constants[158] = "eta_l in component crossbridge_TRPN (second)" legend_constants[159] = "eta_s in component crossbridge_TRPN (second)" legend_constants[160] = "k in component crossbridge_TRPN (dimensionless)" legend_constants[161] = "alpha in component crossbridge_TRPN (kPa)" legend_algebraic[44] = "v_Cdd in component crossbridge_TRPN (m_per_s)" legend_algebraic[45] = "nu_d in component crossbridge_TRPN (m_per_s)" legend_algebraic[46] = "nu_k in component crossbridge_TRPN (m_per_s)" legend_algebraic[47] = "nu_f in component crossbridge_TRPN (m_per_s)" legend_algebraic[53] = "F_d in component crossbridge_TRPN (kPa)" legend_algebraic[37] = "F_1 in component crossbridge_TRPN (kPa)" legend_algebraic[57] = "T_passive in component crossbridge_TRPN (kPa)" legend_constants[162] = "b in component crossbridge_TRPN (dimensionless)" legend_algebraic[68] = "mu_Ca_D in component diadCaDiffusion (J_per_mol)" legend_algebraic[70] = "mu_Ca_i in component diadCaDiffusion (J_per_mol)" legend_algebraic[76] = "mu_Ca_o in component LCC (J_per_mol)" legend_algebraic[78] = "mu_Ca_D in component LCC (J_per_mol)" legend_algebraic[74] = "V_mem in component LCC (volt)" legend_algebraic[88] = "Af_LCC_Ca1 in component LCC (J_per_mol)" legend_algebraic[116] = "Ar_LCC_Ca1 in component LCC (J_per_mol)" legend_algebraic[114] = "Am_LCC_Ca1 in component LCC (J_per_mol)" legend_algebraic[175] = "v_LCC_Ca1 in component LCC (fmol_per_sec)" legend_algebraic[89] = "Af_LCC_Ca2 in component LCC (J_per_mol)" legend_algebraic[117] = "Ar_LCC_Ca2 in component LCC (J_per_mol)" legend_algebraic[115] = "Am_LCC_Ca2 in component LCC (J_per_mol)" legend_algebraic[177] = "v_LCC_Ca2 in component LCC (fmol_per_sec)" legend_algebraic[108] = "Af_fCa00 in component LCC (J_per_mol)" legend_algebraic[154] = "Ar_fCa00 in component LCC (J_per_mol)" legend_algebraic[155] = "v_fCa00 in component LCC (fmol_per_sec)" legend_algebraic[109] = "Af_fCa01 in component LCC (J_per_mol)" legend_algebraic[156] = "Ar_fCa01 in component LCC (J_per_mol)" legend_algebraic[159] = "v_fCa01 in component LCC (fmol_per_sec)" legend_algebraic[110] = "Af_fCa02 in component LCC (J_per_mol)" legend_algebraic[160] = "Ar_fCa02 in component LCC (J_per_mol)" legend_algebraic[163] = "v_fCa02 in component LCC (fmol_per_sec)" legend_algebraic[111] = "Af_fCa10 in component LCC (J_per_mol)" legend_algebraic[164] = "Ar_fCa10 in component LCC (J_per_mol)" legend_algebraic[167] = "v_fCa10 in component LCC (fmol_per_sec)" legend_algebraic[112] = "Af_fCa11 in component LCC (J_per_mol)" legend_algebraic[168] = "Ar_fCa11 in component LCC (J_per_mol)" legend_algebraic[171] = "v_fCa11 in component LCC (fmol_per_sec)" legend_algebraic[113] = "Af_fCa12 in component LCC (J_per_mol)" legend_algebraic[172] = "Ar_fCa12 in component LCC (J_per_mol)" legend_algebraic[174] = "v_fCa12 in component LCC (fmol_per_sec)" legend_algebraic[80] = "mu_S000_LCC in component LCC (J_per_mol)" legend_algebraic[157] = "v_S000_LCC in component LCC (fmol_per_sec)" legend_algebraic[81] = "mu_S010_LCC in component LCC (J_per_mol)" legend_algebraic[161] = "v_S010_LCC in component LCC (fmol_per_sec)" legend_algebraic[82] = "mu_S020_LCC in component LCC (J_per_mol)" legend_algebraic[165] = "v_S020_LCC in component LCC (fmol_per_sec)" legend_algebraic[83] = "mu_S100_LCC in component LCC (J_per_mol)" legend_algebraic[169] = "v_S100_LCC in component LCC (fmol_per_sec)" legend_algebraic[84] = "mu_S110_LCC in component LCC (J_per_mol)" legend_algebraic[173] = "v_S110_LCC in component LCC (fmol_per_sec)" legend_algebraic[86] = "mu_S120_LCC in component LCC (J_per_mol)" legend_algebraic[176] = "v_S120_LCC in component LCC (fmol_per_sec)" legend_algebraic[12] = "mu_S001_LCC in component LCC (J_per_mol)" legend_algebraic[158] = "v_S001_LCC in component LCC (fmol_per_sec)" legend_algebraic[13] = "mu_S011_LCC in component LCC (J_per_mol)" legend_algebraic[162] = "v_S011_LCC in component LCC (fmol_per_sec)" legend_algebraic[14] = "mu_S021_LCC in component LCC (J_per_mol)" legend_algebraic[166] = "v_S021_LCC in component LCC (fmol_per_sec)" legend_algebraic[15] = "mu_S101_LCC in component LCC (J_per_mol)" legend_algebraic[170] = "v_S101_LCC in component LCC (fmol_per_sec)" legend_algebraic[85] = "mu_S111_LCC in component LCC (J_per_mol)" legend_algebraic[178] = "v_S111_LCC in component LCC (fmol_per_sec)" legend_algebraic[87] = "mu_S121_LCC in component LCC (J_per_mol)" legend_algebraic[179] = "v_S121_LCC in component LCC (fmol_per_sec)" legend_algebraic[90] = "Af_d000 in component LCC (J_per_mol)" legend_algebraic[118] = "Ar_d000 in component LCC (J_per_mol)" legend_algebraic[119] = "v_d000 in component LCC (fmol_per_sec)" legend_algebraic[92] = "Af_d010 in component LCC (J_per_mol)" legend_algebraic[122] = "Ar_d010 in component LCC (J_per_mol)" legend_algebraic[123] = "v_d010 in component LCC (fmol_per_sec)" legend_algebraic[94] = "Af_d020 in component LCC (J_per_mol)" legend_algebraic[126] = "Ar_d020 in component LCC (J_per_mol)" legend_algebraic[127] = "v_d020 in component LCC (fmol_per_sec)" legend_algebraic[96] = "Af_f1_000 in component LCC (J_per_mol)" legend_algebraic[130] = "Ar_f1_000 in component LCC (J_per_mol)" legend_algebraic[131] = "v_f1_000 in component LCC (fmol_per_sec)" legend_algebraic[100] = "Af_f2_000 in component LCC (J_per_mol)" legend_algebraic[138] = "Ar_f2_000 in component LCC (J_per_mol)" legend_algebraic[139] = "v_f2_000 in component LCC (fmol_per_sec)" legend_algebraic[104] = "Af_f3_010 in component LCC (J_per_mol)" legend_algebraic[146] = "Ar_f3_010 in component LCC (J_per_mol)" legend_algebraic[147] = "v_f3_010 in component LCC (fmol_per_sec)" legend_algebraic[98] = "Af_f1_100 in component LCC (J_per_mol)" legend_algebraic[134] = "Ar_f1_100 in component LCC (J_per_mol)" legend_algebraic[135] = "v_f1_100 in component LCC (fmol_per_sec)" legend_algebraic[102] = "Af_f2_100 in component LCC (J_per_mol)" legend_algebraic[142] = "Ar_f2_100 in component LCC (J_per_mol)" legend_algebraic[143] = "v_f2_100 in component LCC (fmol_per_sec)" legend_algebraic[106] = "Af_f3_110 in component LCC (J_per_mol)" legend_algebraic[150] = "Ar_f3_110 in component LCC (J_per_mol)" legend_algebraic[151] = "v_f3_110 in component LCC (fmol_per_sec)" legend_algebraic[91] = "Af_d001 in component LCC (J_per_mol)" legend_algebraic[120] = "Ar_d001 in component LCC (J_per_mol)" legend_algebraic[121] = "v_d001 in component LCC (fmol_per_sec)" legend_algebraic[93] = "Af_d011 in component LCC (J_per_mol)" legend_algebraic[124] = "Ar_d011 in component LCC (J_per_mol)" legend_algebraic[125] = "v_d011 in component LCC (fmol_per_sec)" legend_algebraic[95] = "Af_d021 in component LCC (J_per_mol)" legend_algebraic[128] = "Ar_d021 in component LCC (J_per_mol)" legend_algebraic[129] = "v_d021 in component LCC (fmol_per_sec)" legend_algebraic[97] = "Af_f1_001 in component LCC (J_per_mol)" legend_algebraic[132] = "Ar_f1_001 in component LCC (J_per_mol)" legend_algebraic[133] = "v_f1_001 in component LCC (fmol_per_sec)" legend_algebraic[101] = "Af_f2_001 in component LCC (J_per_mol)" legend_algebraic[140] = "Ar_f2_001 in component LCC (J_per_mol)" legend_algebraic[141] = "v_f2_001 in component LCC (fmol_per_sec)" legend_algebraic[105] = "Af_f3_011 in component LCC (J_per_mol)" legend_algebraic[148] = "Ar_f3_011 in component LCC (J_per_mol)" legend_algebraic[149] = "v_f3_011 in component LCC (fmol_per_sec)" legend_algebraic[99] = "Af_f1_101 in component LCC (J_per_mol)" legend_algebraic[136] = "Ar_f1_101 in component LCC (J_per_mol)" legend_algebraic[137] = "v_f1_101 in component LCC (fmol_per_sec)" legend_algebraic[103] = "Af_f2_101 in component LCC (J_per_mol)" legend_algebraic[144] = "Ar_f2_101 in component LCC (J_per_mol)" legend_algebraic[145] = "v_f2_101 in component LCC (fmol_per_sec)" legend_algebraic[107] = "Af_f3_111 in component LCC (J_per_mol)" legend_algebraic[152] = "Ar_f3_111 in component LCC (J_per_mol)" legend_algebraic[153] = "v_f3_111 in component LCC (fmol_per_sec)" legend_algebraic[16] = "Ca_tot in component LCC (fmol)" legend_algebraic[184] = "mu_Ca_o in component NCX (J_per_mol)" legend_algebraic[185] = "mu_Ca_i in component NCX (J_per_mol)" legend_algebraic[186] = "mu_Na_o in component NCX (J_per_mol)" legend_algebraic[187] = "mu_Na_i in component NCX (J_per_mol)" legend_algebraic[188] = "mu_P1_NCX in component NCX (J_per_mol)" legend_algebraic[189] = "mu_P2_NCX in component NCX (J_per_mol)" legend_algebraic[190] = "mu_P3_NCX in component NCX (J_per_mol)" legend_algebraic[191] = "mu_P4_NCX in component NCX (J_per_mol)" legend_algebraic[192] = "mu_P5_NCX in component NCX (J_per_mol)" legend_algebraic[193] = "mu_P6_NCX in component NCX (J_per_mol)" legend_algebraic[218] = "v_P1_NCX in component NCX (fmol_per_sec)" legend_algebraic[206] = "v_P2_NCX in component NCX (fmol_per_sec)" legend_algebraic[209] = "v_P3_NCX in component NCX (fmol_per_sec)" legend_algebraic[212] = "v_P4_NCX in component NCX (fmol_per_sec)" legend_algebraic[215] = "v_P5_NCX in component NCX (fmol_per_sec)" legend_algebraic[219] = "v_P6_NCX in component NCX (fmol_per_sec)" legend_algebraic[194] = "Af_r1 in component NCX (J_per_mol)" legend_algebraic[200] = "Ar_r1 in component NCX (J_per_mol)" legend_algebraic[201] = "v_r1 in component NCX (fmol_per_sec)" legend_algebraic[195] = "Af_r2 in component NCX (J_per_mol)" legend_algebraic[202] = "Ar_r2 in component NCX (J_per_mol)" legend_algebraic[204] = "v_r2 in component NCX (fmol_per_sec)" legend_algebraic[196] = "Af_r3 in component NCX (J_per_mol)" legend_algebraic[205] = "Ar_r3 in component NCX (J_per_mol)" legend_algebraic[207] = "v_r3 in component NCX (fmol_per_sec)" legend_algebraic[197] = "Af_r4 in component NCX (J_per_mol)" legend_algebraic[208] = "Ar_r4 in component NCX (J_per_mol)" legend_algebraic[210] = "v_r4 in component NCX (fmol_per_sec)" legend_algebraic[198] = "Af_r5 in component NCX (J_per_mol)" legend_algebraic[211] = "Ar_r5 in component NCX (J_per_mol)" legend_algebraic[213] = "v_r5 in component NCX (fmol_per_sec)" legend_algebraic[199] = "Af_r6 in component NCX (J_per_mol)" legend_algebraic[216] = "Ar_r6 in component NCX (J_per_mol)" legend_algebraic[217] = "v_r6 in component NCX (fmol_per_sec)" legend_algebraic[183] = "V_mem in component NCX (volt)" legend_algebraic[223] = "mu_Ca_i in component pCa (J_per_mol)" legend_algebraic[224] = "mu_Ca_o in component pCa (J_per_mol)" legend_algebraic[225] = "mu_pCa in component pCa (J_per_mol)" legend_algebraic[226] = "mu_pCa_Ca in component pCa (J_per_mol)" legend_algebraic[0] = "mu_PLB in component PLB (J_per_mol)" legend_algebraic[25] = "mu_PKACI in component PLB (J_per_mol)" legend_algebraic[29] = "mu_PLB_PKACI in component PLB (J_per_mol)" legend_algebraic[34] = "mu_PP1 in component PLB (J_per_mol)" legend_algebraic[39] = "mu_PLBp_PP1 in component PLB (J_per_mol)" legend_algebraic[49] = "mu_PLBp in component PLB (J_per_mol)" legend_algebraic[40] = "mu_Ip in component PLB (J_per_mol)" legend_algebraic[50] = "mu_Ip_PP1 in component PLB (J_per_mol)" legend_constants[163] = "nCa_1 in component RyR (dimensionless)" legend_constants[164] = "nCa_2 in component RyR (dimensionless)" legend_algebraic[229] = "mu_Ca_SR in component RyR (J_per_mol)" legend_algebraic[231] = "mu_Ca_D in component RyR (J_per_mol)" legend_algebraic[17] = "q_Ca_gate_complexes in component RyR (fmol)" legend_algebraic[233] = "mu_C_RyR in component RyR (J_per_mol)" legend_algebraic[234] = "mu_CI_RyR in component RyR (J_per_mol)" legend_algebraic[235] = "mu_I_RyR in component RyR (J_per_mol)" legend_algebraic[236] = "mu_O_RyR in component RyR (J_per_mol)" legend_algebraic[238] = "v_OC in component RyR (fmol_per_sec)" legend_algebraic[240] = "v_CCI in component RyR (fmol_per_sec)" legend_algebraic[242] = "v_CII in component RyR (fmol_per_sec)" legend_algebraic[244] = "v_IO in component RyR (fmol_per_sec)" legend_constants[165] = "n_Ca_i in component SERCA (dimensionless)" legend_constants[166] = "n_Ca_SR in component SERCA (dimensionless)" legend_constants[167] = "n_H in component SERCA (dimensionless)" legend_algebraic[18] = "q_SERCA_Ca_complexes in component SERCA (fmol)" legend_algebraic[19] = "c_Ca_i in component SERCA (mM)" legend_algebraic[20] = "c_Ca_SR in component SERCA (mM)" legend_algebraic[21] = "c_H in component SERCA (mM)" legend_algebraic[22] = "c_MgADP in component SERCA (mM)" legend_algebraic[23] = "c_MgATP in component SERCA (mM)" legend_algebraic[24] = "c_Pi in component SERCA (mM)" legend_constants[168] = "vol_i in component SERCA (pL)" legend_constants[174] = "vol_sr in component SERCA (pL)" legend_constants[176] = "vol_isr in component SERCA (pL)" legend_algebraic[239] = "mu_Ca_i in component SERCA (J_per_mol)" legend_algebraic[241] = "mu_Ca_SR in component SERCA (J_per_mol)" legend_algebraic[243] = "mu_H in component SERCA (J_per_mol)" legend_algebraic[245] = "mu_MgADP in component SERCA (J_per_mol)" legend_algebraic[247] = "mu_MgATP in component SERCA (J_per_mol)" legend_algebraic[249] = "mu_P1 in component SERCA (J_per_mol)" legend_algebraic[284] = "v_P1 in component SERCA (fmol_per_sec)" legend_algebraic[248] = "mu_Pi in component SERCA (J_per_mol)" legend_algebraic[250] = "mu_P2 in component SERCA (J_per_mol)" legend_algebraic[291] = "v_P2 in component SERCA (fmol_per_sec)" legend_algebraic[251] = "mu_P2a in component SERCA (J_per_mol)" legend_algebraic[288] = "v_P2a in component SERCA (fmol_per_sec)" legend_algebraic[252] = "mu_P4 in component SERCA (J_per_mol)" legend_algebraic[293] = "v_P4 in component SERCA (fmol_per_sec)" legend_algebraic[253] = "mu_P5 in component SERCA (J_per_mol)" legend_algebraic[294] = "v_P5 in component SERCA (fmol_per_sec)" legend_algebraic[255] = "mu_P6 in component SERCA (J_per_mol)" legend_algebraic[297] = "v_P6 in component SERCA (fmol_per_sec)" legend_algebraic[257] = "mu_P8 in component SERCA (J_per_mol)" legend_algebraic[299] = "v_P8 in component SERCA (fmol_per_sec)" legend_algebraic[258] = "mu_P9 in component SERCA (J_per_mol)" legend_algebraic[301] = "v_P9 in component SERCA (fmol_per_sec)" legend_algebraic[254] = "mu_P10 in component SERCA (J_per_mol)" legend_algebraic[302] = "v_P10 in component SERCA (fmol_per_sec)" legend_algebraic[259] = "mu_PLBp_baseline in component SERCA (J_per_mol)" legend_constants[173] = "mu_PLBp in component SERCA (J_per_mol)" legend_algebraic[260] = "mu_PLBp_eff in component SERCA (J_per_mol)" legend_algebraic[263] = "Af_R1_2 in component SERCA (J_per_mol)" legend_algebraic[264] = "Ar_R1_2 in component SERCA (J_per_mol)" legend_algebraic[280] = "v_SERCA_R1_2 in component SERCA (fmol_per_sec)" legend_algebraic[265] = "Af_R5_6 in component SERCA (J_per_mol)" legend_algebraic[266] = "Ar_R5_6 in component SERCA (J_per_mol)" legend_algebraic[282] = "v_SERCA_R5_6 in component SERCA (fmol_per_sec)" legend_algebraic[267] = "Af_R2_2a in component SERCA (J_per_mol)" legend_algebraic[268] = "Ar_R2_2a in component SERCA (J_per_mol)" legend_algebraic[285] = "v_SERCA_R2_2a in component SERCA (fmol_per_sec)" legend_algebraic[269] = "Af_R2_4 in component SERCA (J_per_mol)" legend_algebraic[270] = "Ar_R2_4 in component SERCA (J_per_mol)" legend_algebraic[287] = "v_SERCA_R2_4 in component SERCA (fmol_per_sec)" legend_algebraic[271] = "Af_R4_5 in component SERCA (J_per_mol)" legend_algebraic[272] = "Ar_R4_5 in component SERCA (J_per_mol)" legend_algebraic[290] = "v_SERCA_R4_5 in component SERCA (fmol_per_sec)" legend_algebraic[273] = "Af_R6_8 in component SERCA (J_per_mol)" legend_algebraic[274] = "Ar_R6_8 in component SERCA (J_per_mol)" legend_algebraic[292] = "v_SERCA_R6_8 in component SERCA (fmol_per_sec)" legend_algebraic[275] = "Af_R8_9 in component SERCA (J_per_mol)" legend_algebraic[276] = "Ar_R8_9 in component SERCA (J_per_mol)" legend_algebraic[296] = "v_SERCA_R8_9 in component SERCA (fmol_per_sec)" legend_algebraic[277] = "Af_R9_10 in component SERCA (J_per_mol)" legend_algebraic[278] = "Ar_R9_10 in component SERCA (J_per_mol)" legend_algebraic[298] = "v_SERCA_R9_10 in component SERCA (fmol_per_sec)" legend_algebraic[261] = "Af_R10_1 in component SERCA (J_per_mol)" legend_algebraic[262] = "Ar_R10_1 in component SERCA (J_per_mol)" legend_algebraic[279] = "v_SERCA_R10_1 in component SERCA (fmol_per_sec)" legend_algebraic[256] = "v_SERCA_R10_1_noPLBp in component SERCA (fmol_per_sec)" legend_rates[20] = "d/dt q_Ca_SR in component environment (fmol)" legend_rates[4] = "d/dt q_Ca_i in component environment (fmol)" legend_rates[2] = "d/dt q_TRPN in component environment (fmol)" legend_rates[3] = "d/dt q_Ca_TRPN in component environment (fmol)" legend_rates[7] = "d/dt q_B_XB in component environment (fmol)" legend_rates[8] = "d/dt q_U_XB in component environment (fmol)" legend_rates[9] = "d/dt q_W_XB in component environment (fmol)" legend_rates[10] = "d/dt q_S_XB in component environment (fmol)" legend_rates[1] = "d/dt q_Ca_D in component environment (fmol)" legend_rates[11] = "d/dt q_Ca_o in component environment (fmol)" legend_rates[5] = "d/dt q_Na_i in component environment (fmol)" legend_rates[6] = "d/dt q_Na_o in component environment (fmol)" legend_rates[21] = "d/dt q_pCa in component environment (fmol)" legend_rates[22] = "d/dt q_pCa_Ca in component environment (fmol)" legend_rates[12] = "d/dt q_PLB in component environment (fmol)" legend_rates[13] = "d/dt q_PKACI in component environment (fmol)" legend_rates[14] = "d/dt q_PLB_PKACI in component environment (fmol)" legend_rates[15] = "d/dt q_PP1 in component environment (fmol)" legend_rates[16] = "d/dt q_PLBp_PP1 in component environment (fmol)" legend_rates[17] = "d/dt q_PLBp in component environment (fmol)" legend_rates[18] = "d/dt q_Ip in component environment (fmol)" legend_rates[19] = "d/dt q_Ip_PP1 in component environment (fmol)" legend_rates[23] = "d/dt q_H in component environment (fmol)" legend_rates[25] = "d/dt q_MgADP in component environment (fmol)" legend_rates[26] = "d/dt q_MgATP in component environment (fmol)" legend_rates[24] = "d/dt q_Pi in component environment (fmol)" legend_rates[0] = "d/dt q_membrane in component environment (fC)" legend_rates[59] = "d/dt G_w in component crossbridge_TRPN (metre)" legend_rates[58] = "d/dt SL in component crossbridge_TRPN (metre)" legend_rates[60] = "d/dt G_s in component crossbridge_TRPN (metre)" legend_rates[61] = "d/dt Cdd in component crossbridge_TRPN (metre)" legend_rates[27] = "d/dt q_000_LCC in component environment (fmol)" legend_rates[28] = "d/dt q_010_LCC in component environment (fmol)" legend_rates[29] = "d/dt q_020_LCC in component environment (fmol)" legend_rates[30] = "d/dt q_100_LCC in component environment (fmol)" legend_rates[31] = "d/dt q_110_LCC in component environment (fmol)" legend_rates[32] = "d/dt q_120_LCC in component environment (fmol)" legend_rates[33] = "d/dt q_001_LCC in component environment (fmol)" legend_rates[34] = "d/dt q_011_LCC in component environment (fmol)" legend_rates[35] = "d/dt q_021_LCC in component environment (fmol)" legend_rates[36] = "d/dt q_101_LCC in component environment (fmol)" legend_rates[37] = "d/dt q_111_LCC in component environment (fmol)" legend_rates[38] = "d/dt q_121_LCC in component environment (fmol)" legend_rates[39] = "d/dt q_P1_NCX in component environment (fmol)" legend_rates[40] = "d/dt q_P2_NCX in component environment (fmol)" legend_rates[41] = "d/dt q_P3_NCX in component environment (fmol)" legend_rates[42] = "d/dt q_P4_NCX in component environment (fmol)" legend_rates[43] = "d/dt q_P5_NCX in component environment (fmol)" legend_rates[44] = "d/dt q_P6_NCX in component environment (fmol)" legend_rates[48] = "d/dt q_O_RyR in component environment (fmol)" legend_rates[45] = "d/dt q_C_RyR in component environment (fmol)" legend_rates[46] = "d/dt q_CI_RyR in component environment (fmol)" legend_rates[47] = "d/dt q_I_RyR in component environment (fmol)" legend_rates[49] = "d/dt q_P1_SERCA in component environment (fmol)" legend_rates[50] = "d/dt q_P2_SERCA in component environment (fmol)" legend_rates[51] = "d/dt q_P2a_SERCA in component environment (fmol)" legend_rates[52] = "d/dt q_P4_SERCA in component environment (fmol)" legend_rates[53] = "d/dt q_P5_SERCA in component environment (fmol)" legend_rates[57] = "d/dt q_P10_SERCA in component environment (fmol)" legend_rates[54] = "d/dt q_P6_SERCA in component environment (fmol)" legend_rates[55] = "d/dt q_P8_SERCA in component environment (fmol)" legend_rates[56] = "d/dt q_P9_SERCA in component environment (fmol)" return (legend_states, legend_algebraic, legend_voi, legend_constants) def initConsts(): constants = [0.0] * sizeConstants; states = [0.0] * sizeStates; 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] return (states, constants) def computeRates(voi, states, constants): rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic algebraic[2] = custom_piecewise([less(voi , 0.0600000), 1.00000*constants[137] , greater(voi , 0.500000) & less(voi , 1.10000), -0.100000*constants[137] , True, 0.00000*constants[137]]) algebraic[26] = algebraic[2] rates[58] = algebraic[26] rootfind_0(voi, constants, rates, states, algebraic) rates[61] = algebraic[44] algebraic[30] = -constants[169]*algebraic[26] algebraic[35] = constants[171]*states[59] algebraic[41] = -algebraic[30]-algebraic[35] algebraic[51] = algebraic[41] rates[59] = algebraic[51] algebraic[31] = -constants[175]*algebraic[26] algebraic[36] = constants[172]*states[60] algebraic[42] = -algebraic[31]-algebraic[36] algebraic[52] = algebraic[42] rates[60] = algebraic[52] algebraic[34] = constants[141]*constants[142]*log(constants[95]*states[15]) algebraic[40] = constants[141]*constants[142]*log(constants[98]*states[18]) algebraic[50] = constants[141]*constants[142]*log(constants[99]*states[19]) algebraic[56] = constants[45]*(exp((algebraic[34]+algebraic[40])/(constants[141]*constants[142]))-exp(algebraic[50]/(constants[141]*constants[142]))) rates[18] = -algebraic[56] rates[19] = algebraic[56] algebraic[0] = constants[141]*constants[142]*log(constants[92]*states[12]) algebraic[25] = constants[141]*constants[142]*log(constants[93]*states[13]) algebraic[29] = constants[141]*constants[142]*log(constants[94]*states[14]) algebraic[55] = constants[41]*(exp((algebraic[0]+algebraic[25])/(constants[141]*constants[142]))-exp(algebraic[29]/(constants[141]*constants[142]))) algebraic[49] = constants[141]*constants[142]*log(constants[97]*states[17]) algebraic[60] = constants[42]*(exp(algebraic[29]/(constants[141]*constants[142]))-exp((algebraic[49]+algebraic[25])/(constants[141]*constants[142]))) rates[13] = algebraic[60]-algebraic[55] rates[14] = algebraic[55]-algebraic[60] rates[15] = algebraic[55]-algebraic[60] algebraic[39] = constants[141]*constants[142]*log(constants[96]*states[16]) algebraic[62] = constants[43]*(exp((algebraic[49]+algebraic[34])/(constants[141]*constants[142]))-exp(algebraic[39]/(constants[141]*constants[142]))) rates[17] = algebraic[60]-algebraic[62] algebraic[58] = constants[141]*constants[142]*log(constants[62]*states[2]) algebraic[61] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[63] = constants[141]*constants[142]*log(constants[63]*states[3]) algebraic[66] = constants[1]*(exp((constants[131]*algebraic[61]+algebraic[58])/(constants[141]*constants[142]))-exp(algebraic[63]/(constants[141]*constants[142]))) rates[2] = -algebraic[66] rates[3] = algebraic[66] algebraic[65] = constants[44]*(exp(algebraic[39]/(constants[141]*constants[142]))-exp((algebraic[0]+algebraic[34])/(constants[141]*constants[142]))) rates[12] = algebraic[65]-algebraic[55] rates[16] = algebraic[62]-algebraic[65] algebraic[64] = constants[141]*constants[142]*log(constants[64]*states[7]) algebraic[67] = constants[141]*constants[142]*log(constants[65]*states[8]) algebraic[73] = constants[2]*(exp((algebraic[64]+constants[130]*algebraic[63])/(constants[141]*constants[142]))-exp((algebraic[67]+constants[130]*algebraic[63])/(constants[141]*constants[142]))) rates[7] = -algebraic[73] algebraic[69] = constants[141]*constants[142]*log(constants[66]*states[9]) algebraic[75] = constants[3]*(exp(algebraic[67]/(constants[141]*constants[142]))-exp(algebraic[69]/(constants[141]*constants[142]))) algebraic[71] = constants[141]*constants[142]*log(constants[67]*states[10]) algebraic[77] = constants[4]*(exp(algebraic[69]/(constants[141]*constants[142]))-exp(algebraic[71]/(constants[141]*constants[142]))) rates[9] = algebraic[75]-algebraic[77] algebraic[79] = constants[5]*(exp(algebraic[71]/(constants[141]*constants[142]))-exp(algebraic[67]/(constants[141]*constants[142]))) rates[8] = (algebraic[73]-algebraic[75])+algebraic[79] rates[10] = algebraic[77]-algebraic[79] algebraic[80] = constants[141]*constants[142]*log(constants[70]*states[27]) algebraic[108] = algebraic[80] algebraic[78] = constants[141]*constants[142]*log(constants[68]*states[1]) algebraic[154] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[117]*algebraic[78] algebraic[155] = constants[27]*(exp(algebraic[108]/(constants[141]*constants[142]))-exp(algebraic[154]/(constants[141]*constants[142]))) algebraic[74] = states[0]/constants[133] algebraic[90] = constants[141]*constants[142]*log(constants[70]*states[27])+constants[118]*constants[143]*algebraic[74] algebraic[118] = constants[141]*constants[142]*log(constants[73]*states[30])+constants[122]*constants[143]*algebraic[74] algebraic[119] = constants[9]*(exp(algebraic[90]/(constants[141]*constants[142]))-exp(algebraic[118]/(constants[141]*constants[142]))) algebraic[96] = constants[141]*constants[142]*log(constants[70]*states[27])+constants[119]*constants[143]*algebraic[74] algebraic[130] = constants[141]*constants[142]*log(constants[71]*states[28])+constants[123]*constants[143]*algebraic[74] algebraic[131] = constants[15]*(exp(algebraic[96]/(constants[141]*constants[142]))-exp(algebraic[130]/(constants[141]*constants[142]))) algebraic[100] = constants[141]*constants[142]*log(constants[70]*states[27])+constants[120]*constants[143]*algebraic[74] algebraic[138] = constants[141]*constants[142]*log(constants[72]*states[29])+constants[124]*constants[143]*algebraic[74] algebraic[139] = constants[19]*(exp(algebraic[100]/(constants[141]*constants[142]))-exp(algebraic[138]/(constants[141]*constants[142]))) algebraic[157] = ((-algebraic[119]-algebraic[131])-algebraic[139])-algebraic[155] rates[27] = algebraic[157] algebraic[91] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[118]*constants[143]*algebraic[74] algebraic[120] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[122]*constants[143]*algebraic[74] algebraic[121] = constants[12]*(exp(algebraic[91]/(constants[141]*constants[142]))-exp(algebraic[120]/(constants[141]*constants[142]))) algebraic[97] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[119]*constants[143]*algebraic[74] algebraic[132] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[123]*constants[143]*algebraic[74] algebraic[133] = constants[17]*(exp(algebraic[97]/(constants[141]*constants[142]))-exp(algebraic[132]/(constants[141]*constants[142]))) algebraic[101] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[120]*constants[143]*algebraic[74] algebraic[140] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[124]*constants[143]*algebraic[74] algebraic[141] = constants[21]*(exp(algebraic[101]/(constants[141]*constants[142]))-exp(algebraic[140]/(constants[141]*constants[142]))) algebraic[158] = ((-algebraic[121]-algebraic[133])-algebraic[141])--algebraic[155] rates[33] = algebraic[158] algebraic[81] = constants[141]*constants[142]*log(constants[71]*states[28]) algebraic[109] = algebraic[81] algebraic[156] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[117]*algebraic[78] algebraic[159] = constants[29]*(exp(algebraic[109]/(constants[141]*constants[142]))-exp(algebraic[156]/(constants[141]*constants[142]))) algebraic[92] = constants[141]*constants[142]*log(constants[71]*states[28])+constants[118]*constants[143]*algebraic[74] algebraic[122] = constants[141]*constants[142]*log(constants[74]*states[31])+constants[122]*constants[143]*algebraic[74] algebraic[123] = constants[10]*(exp(algebraic[92]/(constants[141]*constants[142]))-exp(algebraic[122]/(constants[141]*constants[142]))) algebraic[104] = constants[141]*constants[142]*log(constants[71]*states[28])+constants[121]*constants[143]*algebraic[74] algebraic[146] = constants[141]*constants[142]*log(constants[72]*states[29])+constants[125]*constants[143]*algebraic[74] algebraic[147] = constants[23]*(exp(algebraic[104]/(constants[141]*constants[142]))-exp(algebraic[146]/(constants[141]*constants[142]))) algebraic[161] = ((-algebraic[123]+algebraic[131])-algebraic[159])-algebraic[147] rates[28] = algebraic[161] algebraic[93] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[118]*constants[143]*algebraic[74] algebraic[124] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[122]*constants[143]*algebraic[74] algebraic[125] = constants[13]*(exp(algebraic[93]/(constants[141]*constants[142]))-exp(algebraic[124]/(constants[141]*constants[142]))) algebraic[105] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[121]*constants[143]*algebraic[74] algebraic[148] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[125]*constants[143]*algebraic[74] algebraic[149] = constants[25]*(exp(algebraic[105]/(constants[141]*constants[142]))-exp(algebraic[148]/(constants[141]*constants[142]))) algebraic[162] = ((-algebraic[125]+algebraic[133])--algebraic[159])-algebraic[149] rates[34] = algebraic[162] algebraic[82] = constants[141]*constants[142]*log(constants[72]*states[29]) algebraic[110] = algebraic[82] algebraic[160] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[117]*algebraic[78] algebraic[163] = constants[31]*(exp(algebraic[110]/(constants[141]*constants[142]))-exp(algebraic[160]/(constants[141]*constants[142]))) algebraic[94] = constants[141]*constants[142]*log(constants[72]*states[29])+constants[118]*constants[143]*algebraic[74] algebraic[126] = constants[141]*constants[142]*log(constants[75]*states[32])+constants[122]*constants[143]*algebraic[74] algebraic[127] = constants[11]*(exp(algebraic[94]/(constants[141]*constants[142]))-exp(algebraic[126]/(constants[141]*constants[142]))) algebraic[165] = ((-algebraic[127]+algebraic[139])-algebraic[163])+algebraic[147] rates[29] = algebraic[165] algebraic[95] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[118]*constants[143]*algebraic[74] algebraic[128] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[122]*constants[143]*algebraic[74] algebraic[129] = constants[14]*(exp(algebraic[95]/(constants[141]*constants[142]))-exp(algebraic[128]/(constants[141]*constants[142]))) algebraic[166] = ((-algebraic[129]+algebraic[141])--algebraic[163])+algebraic[149] rates[35] = algebraic[166] algebraic[83] = constants[141]*constants[142]*log(constants[73]*states[30]) algebraic[111] = algebraic[83] algebraic[164] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[117]*algebraic[78] algebraic[167] = constants[28]*(exp(algebraic[111]/(constants[141]*constants[142]))-exp(algebraic[164]/(constants[141]*constants[142]))) algebraic[98] = constants[141]*constants[142]*log(constants[73]*states[30])+constants[119]*constants[143]*algebraic[74] algebraic[134] = constants[141]*constants[142]*log(constants[74]*states[31])+constants[123]*constants[143]*algebraic[74] algebraic[135] = constants[16]*(exp(algebraic[98]/(constants[141]*constants[142]))-exp(algebraic[134]/(constants[141]*constants[142]))) algebraic[102] = constants[141]*constants[142]*log(constants[73]*states[30])+constants[120]*constants[143]*algebraic[74] algebraic[142] = constants[141]*constants[142]*log(constants[75]*states[32])+constants[124]*constants[143]*algebraic[74] algebraic[143] = constants[20]*(exp(algebraic[102]/(constants[141]*constants[142]))-exp(algebraic[142]/(constants[141]*constants[142]))) algebraic[169] = ((-algebraic[143]+algebraic[119])-algebraic[135])-algebraic[167] rates[30] = algebraic[169] algebraic[99] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[119]*constants[143]*algebraic[74] algebraic[136] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[123]*constants[143]*algebraic[74] algebraic[137] = constants[18]*(exp(algebraic[99]/(constants[141]*constants[142]))-exp(algebraic[136]/(constants[141]*constants[142]))) algebraic[103] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[120]*constants[143]*algebraic[74] algebraic[144] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[124]*constants[143]*algebraic[74] algebraic[145] = constants[22]*(exp(algebraic[103]/(constants[141]*constants[142]))-exp(algebraic[144]/(constants[141]*constants[142]))) algebraic[170] = ((-algebraic[145]+algebraic[121])-algebraic[137])--algebraic[167] rates[36] = algebraic[170] algebraic[84] = constants[141]*constants[142]*log(constants[74]*states[31]) algebraic[112] = algebraic[84] algebraic[168] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[117]*algebraic[78] algebraic[171] = constants[30]*(exp(algebraic[112]/(constants[141]*constants[142]))-exp(algebraic[168]/(constants[141]*constants[142]))) algebraic[106] = constants[141]*constants[142]*log(constants[74]*states[31])+constants[121]*constants[143]*algebraic[74] algebraic[150] = constants[141]*constants[142]*log(constants[75]*states[32])+constants[125]*constants[143]*algebraic[74] algebraic[151] = constants[24]*(exp(algebraic[106]/(constants[141]*constants[142]))-exp(algebraic[150]/(constants[141]*constants[142]))) algebraic[173] = ((algebraic[135]+algebraic[123])-algebraic[171])-algebraic[151] rates[31] = algebraic[173] algebraic[86] = constants[141]*constants[142]*log(constants[75]*states[32]) algebraic[113] = algebraic[86] algebraic[172] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[117]*algebraic[78] algebraic[174] = constants[32]*(exp(algebraic[113]/(constants[141]*constants[142]))-exp(algebraic[172]/(constants[141]*constants[142]))) algebraic[176] = ((algebraic[143]+algebraic[127])-algebraic[174])+algebraic[151] rates[32] = algebraic[176] algebraic[85] = constants[141]*constants[142]*log(constants[80]*states[37]) algebraic[88] = algebraic[78]+constants[147]*constants[143]*algebraic[74]+algebraic[85] algebraic[76] = constants[141]*constants[142]*log(constants[69]*states[11]) algebraic[116] = algebraic[76]+algebraic[85] algebraic[114] = constants[147]*constants[143]*algebraic[74] algebraic[175] = custom_piecewise([equal(algebraic[114] , 0.00000), constants[7]*(exp(algebraic[88]/(constants[141]*constants[142]))-exp(algebraic[116]/(constants[141]*constants[142]))) , True, (((constants[7]*algebraic[114])/(constants[141]*constants[142]))/(exp(algebraic[114]/(constants[141]*constants[142]))-1.00000))*(exp(algebraic[88]/(constants[141]*constants[142]))-exp(algebraic[116]/(constants[141]*constants[142])))]) algebraic[107] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[121]*constants[143]*algebraic[74] algebraic[152] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[125]*constants[143]*algebraic[74] algebraic[153] = constants[26]*(exp(algebraic[107]/(constants[141]*constants[142]))-exp(algebraic[152]/(constants[141]*constants[142]))) algebraic[178] = ((algebraic[137]+algebraic[125])-((-algebraic[171]+algebraic[175])-algebraic[175]))-algebraic[153] rates[37] = algebraic[178] algebraic[87] = constants[141]*constants[142]*log(constants[81]*states[38]) algebraic[89] = algebraic[78]+constants[147]*constants[143]*algebraic[74]+algebraic[87] algebraic[117] = algebraic[76]+algebraic[87] algebraic[115] = constants[147]*constants[143]*algebraic[74] algebraic[177] = custom_piecewise([equal(algebraic[115] , 0.00000), constants[8]*(exp(algebraic[89]/(constants[141]*constants[142]))-exp(algebraic[117]/(constants[141]*constants[142]))) , True, (((constants[8]*algebraic[115])/(constants[141]*constants[142]))/(exp(algebraic[115]/(constants[141]*constants[142]))-1.00000))*(exp(algebraic[89]/(constants[141]*constants[142]))-exp(algebraic[117]/(constants[141]*constants[142])))]) algebraic[179] = ((algebraic[145]+algebraic[129])-((-algebraic[174]+algebraic[177])-algebraic[177]))+algebraic[153] rates[38] = algebraic[179] algebraic[188] = constants[141]*constants[142]*log(constants[82]*states[39]) algebraic[194] = algebraic[188] algebraic[187] = constants[141]*constants[142]*log(constants[88]*states[5]) algebraic[189] = constants[141]*constants[142]*log(constants[83]*states[40]) algebraic[200] = constants[127]*algebraic[187]+algebraic[189] algebraic[201] = constants[33]*(exp(algebraic[194]/(constants[141]*constants[142]))-exp(algebraic[200]/(constants[141]*constants[142]))) algebraic[203] = constants[127]*algebraic[201] rates[5] = algebraic[203] algebraic[185] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[195] = algebraic[189]+algebraic[185] algebraic[190] = constants[141]*constants[142]*log(constants[84]*states[41]) algebraic[202] = algebraic[190] algebraic[204] = constants[34]*(exp(algebraic[195]/(constants[141]*constants[142]))-exp(algebraic[202]/(constants[141]*constants[142]))) algebraic[206] = algebraic[201]-algebraic[204] rates[40] = algebraic[206] algebraic[196] = algebraic[190] algebraic[191] = constants[141]*constants[142]*log(constants[85]*states[42]) algebraic[205] = algebraic[191] algebraic[207] = constants[35]*(exp(algebraic[196]/(constants[141]*constants[142]))-exp(algebraic[205]/(constants[141]*constants[142]))) algebraic[209] = algebraic[204]-algebraic[207] rates[41] = algebraic[209] algebraic[197] = algebraic[191] algebraic[184] = constants[141]*constants[142]*log(constants[69]*states[11]) algebraic[192] = constants[141]*constants[142]*log(constants[86]*states[43]) algebraic[208] = algebraic[192]+algebraic[184] algebraic[210] = constants[36]*(exp(algebraic[197]/(constants[141]*constants[142]))-exp(algebraic[208]/(constants[141]*constants[142]))) algebraic[212] = -algebraic[210]+algebraic[207] rates[42] = algebraic[212] algebraic[186] = constants[141]*constants[142]*log(constants[89]*states[6]) algebraic[198] = algebraic[192]+constants[126]*algebraic[186] algebraic[193] = constants[141]*constants[142]*log(constants[87]*states[44]) algebraic[211] = algebraic[193] algebraic[213] = constants[37]*(exp(algebraic[198]/(constants[141]*constants[142]))-exp(algebraic[211]/(constants[141]*constants[142]))) algebraic[214] = -constants[126]*algebraic[213] rates[6] = algebraic[214] algebraic[215] = -algebraic[213]+algebraic[210] rates[43] = algebraic[215] algebraic[183] = states[0]/constants[133] algebraic[199] = constants[141]*constants[142]*log(constants[87]*states[44])+constants[128]*constants[143]*algebraic[183] algebraic[216] = constants[141]*constants[142]*log(constants[82]*states[39])+constants[129]*constants[143]*algebraic[183] algebraic[217] = constants[38]*(exp(algebraic[199]/(constants[141]*constants[142]))-exp(algebraic[216]/(constants[141]*constants[142]))) algebraic[218] = algebraic[217]-algebraic[201] rates[39] = algebraic[218] algebraic[219] = -algebraic[217]+algebraic[213] rates[44] = algebraic[219] algebraic[181] = algebraic[177]+algebraic[175] algebraic[221] = algebraic[210] algebraic[224] = constants[141]*constants[142]*log(constants[69]*states[11]) algebraic[225] = constants[141]*constants[142]*log(constants[90]*states[21]) algebraic[226] = constants[141]*constants[142]*log(constants[91]*states[22]) algebraic[228] = constants[40]*(exp(algebraic[226]/(constants[141]*constants[142]))-exp((algebraic[224]+algebraic[225])/(constants[141]*constants[142]))) rates[11] = algebraic[181]+algebraic[221]+algebraic[228] algebraic[223] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[227] = constants[39]*(exp((algebraic[223]+algebraic[225])/(constants[141]*constants[142]))-exp(algebraic[226]/(constants[141]*constants[142]))) rates[21] = -algebraic[227]+algebraic[228] rates[22] = algebraic[227]-algebraic[228] algebraic[180] = constants[143]*((((((((((((((((((((((((((((((((((((constants[124]*algebraic[141]-constants[120]*algebraic[141])+constants[125]*algebraic[149])-constants[121]*algebraic[149])+constants[123]*algebraic[137])-constants[119]*algebraic[137])+constants[124]*algebraic[145])-constants[120]*algebraic[145])+constants[125]*algebraic[153])-constants[121]*algebraic[153])+constants[123]*algebraic[133])-constants[119]*algebraic[133])+constants[122]*algebraic[129])-constants[118]*algebraic[129])+constants[122]*algebraic[125])-constants[118]*algebraic[125])+constants[122]*algebraic[121])-constants[118]*algebraic[121])+constants[125]*algebraic[151])-constants[121]*algebraic[151])+constants[124]*algebraic[143])-constants[120]*algebraic[143])+constants[123]*algebraic[135])-constants[119]*algebraic[135])+constants[125]*algebraic[147])-constants[121]*algebraic[147])+constants[124]*algebraic[139])-constants[120]*algebraic[139])+constants[123]*algebraic[131])-constants[119]*algebraic[131])+constants[122]*algebraic[127])-constants[118]*algebraic[127])+constants[122]*algebraic[123])-constants[118]*algebraic[123])+constants[122]*algebraic[119])-constants[118]*algebraic[119])-(constants[147]*algebraic[175]+constants[147]*algebraic[177])) algebraic[220] = constants[143]*(-constants[128]*algebraic[217]+constants[129]*algebraic[217]) algebraic[230] = -constants[147]*constants[143]*algebraic[228] algebraic[11] = voi-floor(voi/constants[135])*constants[135] algebraic[28] = custom_piecewise([greater_equal(algebraic[11] , 0.300000) & less_equal(algebraic[11] , 0.300000+constants[136]), (0.00000*0.0300000*constants[133])/constants[136] , True, 0.00000]) algebraic[232] = algebraic[180]+algebraic[220]+algebraic[230]+algebraic[28] rates[0] = algebraic[232] algebraic[231] = constants[141]*constants[142]*log(constants[68]*states[1]) algebraic[233] = constants[141]*constants[142]*log(constants[100]*states[45]) algebraic[236] = constants[141]*constants[142]*log(constants[103]*states[48]) algebraic[238] = constants[47]*(exp(algebraic[236]/(constants[141]*constants[142]))-exp((algebraic[233]+constants[164]*algebraic[231])/(constants[141]*constants[142]))) algebraic[234] = constants[141]*constants[142]*log(constants[101]*states[46]) algebraic[240] = constants[48]*(exp((algebraic[233]+constants[163]*algebraic[231])/(constants[141]*constants[142]))-exp(algebraic[234]/(constants[141]*constants[142]))) rates[45] = algebraic[238]-algebraic[240] algebraic[235] = constants[141]*constants[142]*log(constants[102]*states[47]) algebraic[242] = constants[49]*(exp((algebraic[234]+constants[164]*algebraic[231])/(constants[141]*constants[142]))-exp(algebraic[235]/(constants[141]*constants[142]))) rates[46] = algebraic[240]-algebraic[242] algebraic[244] = constants[50]*(exp(algebraic[235]/(constants[141]*constants[142]))-exp((algebraic[236]+constants[163]*algebraic[231])/(constants[141]*constants[142]))) rates[48] = algebraic[244]-algebraic[238] rates[47] = algebraic[242]-algebraic[244] algebraic[68] = constants[141]*constants[142]*log(constants[68]*states[1]) algebraic[70] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[72] = constants[6]*(exp(algebraic[68]/(constants[141]*constants[142]))-exp(algebraic[70]/(constants[141]*constants[142]))) algebraic[182] = (((constants[117])*algebraic[167]+constants[117]*algebraic[171]+constants[117]*algebraic[163]+constants[117]*algebraic[159]+constants[117]*algebraic[155]+constants[117]*algebraic[174])-algebraic[175])-algebraic[177] algebraic[229] = constants[141]*constants[142]*log(constants[60]*states[20]) algebraic[237] = constants[46]*exp(algebraic[236]/(constants[141]*constants[142]))*(exp(algebraic[229]/(constants[141]*constants[142]))-exp(algebraic[231]/(constants[141]*constants[142]))) algebraic[246] = ((constants[164]*algebraic[238]-constants[163]*algebraic[240])-constants[164]*algebraic[242])+constants[163]*algebraic[244] rates[1] = -algebraic[72]+algebraic[182]+algebraic[237]+algebraic[246] algebraic[254] = constants[141]*constants[142]*log(constants[112]*states[57]) algebraic[259] = constants[141]*constants[142]*log(constants[97]*states[17]) algebraic[260] = algebraic[259]-constants[173] algebraic[261] = algebraic[254]+algebraic[260] algebraic[249] = constants[141]*constants[142]*log(constants[104]*states[49]) algebraic[248] = constants[141]*constants[142]*log(constants[116]*states[24]) algebraic[262] = algebraic[249]+algebraic[248]+algebraic[260] algebraic[279] = constants[59]*(exp(algebraic[261]/(constants[141]*constants[142]))-exp(algebraic[262]/(constants[141]*constants[142]))) algebraic[281] = algebraic[279] rates[24] = algebraic[281] algebraic[247] = constants[141]*constants[142]*log(constants[114]*states[26]) algebraic[263] = algebraic[249]+algebraic[247]+algebraic[260] algebraic[250] = constants[141]*constants[142]*log(constants[105]*states[50]) algebraic[264] = algebraic[250]+algebraic[260] algebraic[280] = constants[51]*(exp(algebraic[263]/(constants[141]*constants[142]))-exp(algebraic[264]/(constants[141]*constants[142]))) algebraic[283] = -algebraic[280] rates[26] = algebraic[283] algebraic[284] = algebraic[279]-algebraic[280] rates[49] = algebraic[284] algebraic[253] = constants[141]*constants[142]*log(constants[108]*states[53]) algebraic[265] = algebraic[253]+algebraic[260] algebraic[245] = constants[141]*constants[142]*log(constants[115]*states[25]) algebraic[255] = constants[141]*constants[142]*log(constants[109]*states[54]) algebraic[266] = algebraic[245]+algebraic[255]+algebraic[260] algebraic[282] = constants[55]*(exp(algebraic[265]/(constants[141]*constants[142]))-exp(algebraic[266]/(constants[141]*constants[142]))) algebraic[286] = algebraic[282] rates[25] = algebraic[286] algebraic[243] = constants[141]*constants[142]*log(constants[113]*states[23]) algebraic[267] = algebraic[250]+algebraic[243]+algebraic[260] algebraic[251] = constants[141]*constants[142]*log(constants[106]*states[51]) algebraic[268] = algebraic[251]+algebraic[260] algebraic[285] = constants[53]*(exp(algebraic[267]/(constants[141]*constants[142]))-exp(algebraic[268]/(constants[141]*constants[142]))) algebraic[288] = algebraic[285] rates[51] = algebraic[288] algebraic[38] = (constants[141]*constants[142])*log(constants[60]*states[20]) algebraic[48] = (constants[141]*constants[142])*log(constants[61]*states[4]) algebraic[54] = constants[0]*(exp(algebraic[38]/(constants[141]*constants[142]))-exp(algebraic[48]/(constants[141]*constants[142]))) algebraic[222] = -algebraic[204] algebraic[239] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[269] = algebraic[250]+constants[165]*algebraic[239]+algebraic[260] algebraic[252] = constants[141]*constants[142]*log(constants[107]*states[52]) algebraic[270] = algebraic[252]+algebraic[260] algebraic[287] = constants[52]*(exp(algebraic[269]/(constants[141]*constants[142]))-exp(algebraic[270]/(constants[141]*constants[142]))) algebraic[289] = -constants[165]*algebraic[287] algebraic[33] = custom_piecewise([greater_equal(algebraic[11] , 0.300000) & less_equal(algebraic[11] , 0.300000+constants[136]), 0.00000*8000.00 , True, 0.00000]) rates[4] = ((-algebraic[66]*constants[131]+algebraic[289]+algebraic[54]+algebraic[72]+algebraic[222])-algebraic[227])+algebraic[33] algebraic[291] = algebraic[280]-algebraic[287] rates[50] = algebraic[291] algebraic[271] = algebraic[252]+algebraic[260] algebraic[272] = algebraic[253]+constants[167]*algebraic[243]+algebraic[260] algebraic[290] = constants[54]*(exp(algebraic[271]/(constants[141]*constants[142]))-exp(algebraic[272]/(constants[141]*constants[142]))) algebraic[293] = algebraic[287]-algebraic[290] rates[52] = algebraic[293] algebraic[294] = algebraic[290]-algebraic[282] rates[53] = algebraic[294] algebraic[273] = algebraic[255]+algebraic[260] algebraic[241] = constants[141]*constants[142]*log(constants[60]*states[20]) algebraic[257] = constants[141]*constants[142]*log(constants[110]*states[55]) algebraic[274] = algebraic[257]+constants[166]*algebraic[241]+algebraic[260] algebraic[292] = constants[56]*(exp(algebraic[273]/(constants[141]*constants[142]))-exp(algebraic[274]/(constants[141]*constants[142]))) algebraic[295] = constants[166]*algebraic[292] rates[20] = (-algebraic[54]-algebraic[237])+algebraic[295] algebraic[297] = algebraic[282]-algebraic[292] rates[54] = algebraic[297] algebraic[275] = algebraic[257]+constants[167]*algebraic[243]+algebraic[260] algebraic[258] = constants[141]*constants[142]*log(constants[111]*states[56]) algebraic[276] = algebraic[258]+algebraic[260] algebraic[296] = constants[57]*(exp(algebraic[275]/(constants[141]*constants[142]))-exp(algebraic[276]/(constants[141]*constants[142]))) algebraic[299] = algebraic[292]-algebraic[296] rates[55] = algebraic[299] algebraic[277] = algebraic[258]+algebraic[260] algebraic[278] = algebraic[243]+algebraic[254]+algebraic[260] algebraic[298] = constants[58]*(exp(algebraic[277]/(constants[141]*constants[142]))-exp(algebraic[278]/(constants[141]*constants[142]))) algebraic[300] = ((constants[167]*algebraic[290]-algebraic[285])-constants[167]*algebraic[296])+algebraic[298] rates[23] = algebraic[300] algebraic[302] = algebraic[298]-algebraic[279] rates[57] = algebraic[302] algebraic[301] = algebraic[296]-algebraic[298] rates[56] = algebraic[301] return(rates) def computeAlgebraic(constants, states, voi): algebraic = array([[0.0] * len(voi)] * sizeAlgebraic) states = array(states) voi = array(voi) algebraic[2] = custom_piecewise([less(voi , 0.0600000), 1.00000*constants[137] , greater(voi , 0.500000) & less(voi , 1.10000), -0.100000*constants[137] , True, 0.00000*constants[137]]) algebraic[26] = algebraic[2] algebraic[30] = -constants[169]*algebraic[26] algebraic[35] = constants[171]*states[59] algebraic[41] = -algebraic[30]-algebraic[35] algebraic[51] = algebraic[41] algebraic[31] = -constants[175]*algebraic[26] algebraic[36] = constants[172]*states[60] algebraic[42] = -algebraic[31]-algebraic[36] algebraic[52] = algebraic[42] algebraic[34] = constants[141]*constants[142]*log(constants[95]*states[15]) algebraic[40] = constants[141]*constants[142]*log(constants[98]*states[18]) algebraic[50] = constants[141]*constants[142]*log(constants[99]*states[19]) algebraic[56] = constants[45]*(exp((algebraic[34]+algebraic[40])/(constants[141]*constants[142]))-exp(algebraic[50]/(constants[141]*constants[142]))) algebraic[0] = constants[141]*constants[142]*log(constants[92]*states[12]) algebraic[25] = constants[141]*constants[142]*log(constants[93]*states[13]) algebraic[29] = constants[141]*constants[142]*log(constants[94]*states[14]) algebraic[55] = constants[41]*(exp((algebraic[0]+algebraic[25])/(constants[141]*constants[142]))-exp(algebraic[29]/(constants[141]*constants[142]))) algebraic[49] = constants[141]*constants[142]*log(constants[97]*states[17]) algebraic[60] = constants[42]*(exp(algebraic[29]/(constants[141]*constants[142]))-exp((algebraic[49]+algebraic[25])/(constants[141]*constants[142]))) algebraic[39] = constants[141]*constants[142]*log(constants[96]*states[16]) algebraic[62] = constants[43]*(exp((algebraic[49]+algebraic[34])/(constants[141]*constants[142]))-exp(algebraic[39]/(constants[141]*constants[142]))) algebraic[58] = constants[141]*constants[142]*log(constants[62]*states[2]) algebraic[61] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[63] = constants[141]*constants[142]*log(constants[63]*states[3]) algebraic[66] = constants[1]*(exp((constants[131]*algebraic[61]+algebraic[58])/(constants[141]*constants[142]))-exp(algebraic[63]/(constants[141]*constants[142]))) algebraic[65] = constants[44]*(exp(algebraic[39]/(constants[141]*constants[142]))-exp((algebraic[0]+algebraic[34])/(constants[141]*constants[142]))) algebraic[64] = constants[141]*constants[142]*log(constants[64]*states[7]) algebraic[67] = constants[141]*constants[142]*log(constants[65]*states[8]) algebraic[73] = constants[2]*(exp((algebraic[64]+constants[130]*algebraic[63])/(constants[141]*constants[142]))-exp((algebraic[67]+constants[130]*algebraic[63])/(constants[141]*constants[142]))) algebraic[69] = constants[141]*constants[142]*log(constants[66]*states[9]) algebraic[75] = constants[3]*(exp(algebraic[67]/(constants[141]*constants[142]))-exp(algebraic[69]/(constants[141]*constants[142]))) algebraic[71] = constants[141]*constants[142]*log(constants[67]*states[10]) algebraic[77] = constants[4]*(exp(algebraic[69]/(constants[141]*constants[142]))-exp(algebraic[71]/(constants[141]*constants[142]))) algebraic[79] = constants[5]*(exp(algebraic[71]/(constants[141]*constants[142]))-exp(algebraic[67]/(constants[141]*constants[142]))) algebraic[80] = constants[141]*constants[142]*log(constants[70]*states[27]) algebraic[108] = algebraic[80] algebraic[78] = constants[141]*constants[142]*log(constants[68]*states[1]) algebraic[154] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[117]*algebraic[78] algebraic[155] = constants[27]*(exp(algebraic[108]/(constants[141]*constants[142]))-exp(algebraic[154]/(constants[141]*constants[142]))) algebraic[74] = states[0]/constants[133] algebraic[90] = constants[141]*constants[142]*log(constants[70]*states[27])+constants[118]*constants[143]*algebraic[74] algebraic[118] = constants[141]*constants[142]*log(constants[73]*states[30])+constants[122]*constants[143]*algebraic[74] algebraic[119] = constants[9]*(exp(algebraic[90]/(constants[141]*constants[142]))-exp(algebraic[118]/(constants[141]*constants[142]))) algebraic[96] = constants[141]*constants[142]*log(constants[70]*states[27])+constants[119]*constants[143]*algebraic[74] algebraic[130] = constants[141]*constants[142]*log(constants[71]*states[28])+constants[123]*constants[143]*algebraic[74] algebraic[131] = constants[15]*(exp(algebraic[96]/(constants[141]*constants[142]))-exp(algebraic[130]/(constants[141]*constants[142]))) algebraic[100] = constants[141]*constants[142]*log(constants[70]*states[27])+constants[120]*constants[143]*algebraic[74] algebraic[138] = constants[141]*constants[142]*log(constants[72]*states[29])+constants[124]*constants[143]*algebraic[74] algebraic[139] = constants[19]*(exp(algebraic[100]/(constants[141]*constants[142]))-exp(algebraic[138]/(constants[141]*constants[142]))) algebraic[157] = ((-algebraic[119]-algebraic[131])-algebraic[139])-algebraic[155] algebraic[91] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[118]*constants[143]*algebraic[74] algebraic[120] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[122]*constants[143]*algebraic[74] algebraic[121] = constants[12]*(exp(algebraic[91]/(constants[141]*constants[142]))-exp(algebraic[120]/(constants[141]*constants[142]))) algebraic[97] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[119]*constants[143]*algebraic[74] algebraic[132] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[123]*constants[143]*algebraic[74] algebraic[133] = constants[17]*(exp(algebraic[97]/(constants[141]*constants[142]))-exp(algebraic[132]/(constants[141]*constants[142]))) algebraic[101] = constants[141]*constants[142]*log(constants[76]*states[33])+constants[120]*constants[143]*algebraic[74] algebraic[140] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[124]*constants[143]*algebraic[74] algebraic[141] = constants[21]*(exp(algebraic[101]/(constants[141]*constants[142]))-exp(algebraic[140]/(constants[141]*constants[142]))) algebraic[158] = ((-algebraic[121]-algebraic[133])-algebraic[141])--algebraic[155] algebraic[81] = constants[141]*constants[142]*log(constants[71]*states[28]) algebraic[109] = algebraic[81] algebraic[156] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[117]*algebraic[78] algebraic[159] = constants[29]*(exp(algebraic[109]/(constants[141]*constants[142]))-exp(algebraic[156]/(constants[141]*constants[142]))) algebraic[92] = constants[141]*constants[142]*log(constants[71]*states[28])+constants[118]*constants[143]*algebraic[74] algebraic[122] = constants[141]*constants[142]*log(constants[74]*states[31])+constants[122]*constants[143]*algebraic[74] algebraic[123] = constants[10]*(exp(algebraic[92]/(constants[141]*constants[142]))-exp(algebraic[122]/(constants[141]*constants[142]))) algebraic[104] = constants[141]*constants[142]*log(constants[71]*states[28])+constants[121]*constants[143]*algebraic[74] algebraic[146] = constants[141]*constants[142]*log(constants[72]*states[29])+constants[125]*constants[143]*algebraic[74] algebraic[147] = constants[23]*(exp(algebraic[104]/(constants[141]*constants[142]))-exp(algebraic[146]/(constants[141]*constants[142]))) algebraic[161] = ((-algebraic[123]+algebraic[131])-algebraic[159])-algebraic[147] algebraic[93] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[118]*constants[143]*algebraic[74] algebraic[124] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[122]*constants[143]*algebraic[74] algebraic[125] = constants[13]*(exp(algebraic[93]/(constants[141]*constants[142]))-exp(algebraic[124]/(constants[141]*constants[142]))) algebraic[105] = constants[141]*constants[142]*log(constants[77]*states[34])+constants[121]*constants[143]*algebraic[74] algebraic[148] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[125]*constants[143]*algebraic[74] algebraic[149] = constants[25]*(exp(algebraic[105]/(constants[141]*constants[142]))-exp(algebraic[148]/(constants[141]*constants[142]))) algebraic[162] = ((-algebraic[125]+algebraic[133])--algebraic[159])-algebraic[149] algebraic[82] = constants[141]*constants[142]*log(constants[72]*states[29]) algebraic[110] = algebraic[82] algebraic[160] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[117]*algebraic[78] algebraic[163] = constants[31]*(exp(algebraic[110]/(constants[141]*constants[142]))-exp(algebraic[160]/(constants[141]*constants[142]))) algebraic[94] = constants[141]*constants[142]*log(constants[72]*states[29])+constants[118]*constants[143]*algebraic[74] algebraic[126] = constants[141]*constants[142]*log(constants[75]*states[32])+constants[122]*constants[143]*algebraic[74] algebraic[127] = constants[11]*(exp(algebraic[94]/(constants[141]*constants[142]))-exp(algebraic[126]/(constants[141]*constants[142]))) algebraic[165] = ((-algebraic[127]+algebraic[139])-algebraic[163])+algebraic[147] algebraic[95] = constants[141]*constants[142]*log(constants[78]*states[35])+constants[118]*constants[143]*algebraic[74] algebraic[128] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[122]*constants[143]*algebraic[74] algebraic[129] = constants[14]*(exp(algebraic[95]/(constants[141]*constants[142]))-exp(algebraic[128]/(constants[141]*constants[142]))) algebraic[166] = ((-algebraic[129]+algebraic[141])--algebraic[163])+algebraic[149] algebraic[83] = constants[141]*constants[142]*log(constants[73]*states[30]) algebraic[111] = algebraic[83] algebraic[164] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[117]*algebraic[78] algebraic[167] = constants[28]*(exp(algebraic[111]/(constants[141]*constants[142]))-exp(algebraic[164]/(constants[141]*constants[142]))) algebraic[98] = constants[141]*constants[142]*log(constants[73]*states[30])+constants[119]*constants[143]*algebraic[74] algebraic[134] = constants[141]*constants[142]*log(constants[74]*states[31])+constants[123]*constants[143]*algebraic[74] algebraic[135] = constants[16]*(exp(algebraic[98]/(constants[141]*constants[142]))-exp(algebraic[134]/(constants[141]*constants[142]))) algebraic[102] = constants[141]*constants[142]*log(constants[73]*states[30])+constants[120]*constants[143]*algebraic[74] algebraic[142] = constants[141]*constants[142]*log(constants[75]*states[32])+constants[124]*constants[143]*algebraic[74] algebraic[143] = constants[20]*(exp(algebraic[102]/(constants[141]*constants[142]))-exp(algebraic[142]/(constants[141]*constants[142]))) algebraic[169] = ((-algebraic[143]+algebraic[119])-algebraic[135])-algebraic[167] algebraic[99] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[119]*constants[143]*algebraic[74] algebraic[136] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[123]*constants[143]*algebraic[74] algebraic[137] = constants[18]*(exp(algebraic[99]/(constants[141]*constants[142]))-exp(algebraic[136]/(constants[141]*constants[142]))) algebraic[103] = constants[141]*constants[142]*log(constants[79]*states[36])+constants[120]*constants[143]*algebraic[74] algebraic[144] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[124]*constants[143]*algebraic[74] algebraic[145] = constants[22]*(exp(algebraic[103]/(constants[141]*constants[142]))-exp(algebraic[144]/(constants[141]*constants[142]))) algebraic[170] = ((-algebraic[145]+algebraic[121])-algebraic[137])--algebraic[167] algebraic[84] = constants[141]*constants[142]*log(constants[74]*states[31]) algebraic[112] = algebraic[84] algebraic[168] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[117]*algebraic[78] algebraic[171] = constants[30]*(exp(algebraic[112]/(constants[141]*constants[142]))-exp(algebraic[168]/(constants[141]*constants[142]))) algebraic[106] = constants[141]*constants[142]*log(constants[74]*states[31])+constants[121]*constants[143]*algebraic[74] algebraic[150] = constants[141]*constants[142]*log(constants[75]*states[32])+constants[125]*constants[143]*algebraic[74] algebraic[151] = constants[24]*(exp(algebraic[106]/(constants[141]*constants[142]))-exp(algebraic[150]/(constants[141]*constants[142]))) algebraic[173] = ((algebraic[135]+algebraic[123])-algebraic[171])-algebraic[151] algebraic[86] = constants[141]*constants[142]*log(constants[75]*states[32]) algebraic[113] = algebraic[86] algebraic[172] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[117]*algebraic[78] algebraic[174] = constants[32]*(exp(algebraic[113]/(constants[141]*constants[142]))-exp(algebraic[172]/(constants[141]*constants[142]))) algebraic[176] = ((algebraic[143]+algebraic[127])-algebraic[174])+algebraic[151] algebraic[85] = constants[141]*constants[142]*log(constants[80]*states[37]) algebraic[88] = algebraic[78]+constants[147]*constants[143]*algebraic[74]+algebraic[85] algebraic[76] = constants[141]*constants[142]*log(constants[69]*states[11]) algebraic[116] = algebraic[76]+algebraic[85] algebraic[114] = constants[147]*constants[143]*algebraic[74] algebraic[175] = custom_piecewise([equal(algebraic[114] , 0.00000), constants[7]*(exp(algebraic[88]/(constants[141]*constants[142]))-exp(algebraic[116]/(constants[141]*constants[142]))) , True, (((constants[7]*algebraic[114])/(constants[141]*constants[142]))/(exp(algebraic[114]/(constants[141]*constants[142]))-1.00000))*(exp(algebraic[88]/(constants[141]*constants[142]))-exp(algebraic[116]/(constants[141]*constants[142])))]) algebraic[107] = constants[141]*constants[142]*log(constants[80]*states[37])+constants[121]*constants[143]*algebraic[74] algebraic[152] = constants[141]*constants[142]*log(constants[81]*states[38])+constants[125]*constants[143]*algebraic[74] algebraic[153] = constants[26]*(exp(algebraic[107]/(constants[141]*constants[142]))-exp(algebraic[152]/(constants[141]*constants[142]))) algebraic[178] = ((algebraic[137]+algebraic[125])-((-algebraic[171]+algebraic[175])-algebraic[175]))-algebraic[153] algebraic[87] = constants[141]*constants[142]*log(constants[81]*states[38]) algebraic[89] = algebraic[78]+constants[147]*constants[143]*algebraic[74]+algebraic[87] algebraic[117] = algebraic[76]+algebraic[87] algebraic[115] = constants[147]*constants[143]*algebraic[74] algebraic[177] = custom_piecewise([equal(algebraic[115] , 0.00000), constants[8]*(exp(algebraic[89]/(constants[141]*constants[142]))-exp(algebraic[117]/(constants[141]*constants[142]))) , True, (((constants[8]*algebraic[115])/(constants[141]*constants[142]))/(exp(algebraic[115]/(constants[141]*constants[142]))-1.00000))*(exp(algebraic[89]/(constants[141]*constants[142]))-exp(algebraic[117]/(constants[141]*constants[142])))]) algebraic[179] = ((algebraic[145]+algebraic[129])-((-algebraic[174]+algebraic[177])-algebraic[177]))+algebraic[153] algebraic[188] = constants[141]*constants[142]*log(constants[82]*states[39]) algebraic[194] = algebraic[188] algebraic[187] = constants[141]*constants[142]*log(constants[88]*states[5]) algebraic[189] = constants[141]*constants[142]*log(constants[83]*states[40]) algebraic[200] = constants[127]*algebraic[187]+algebraic[189] algebraic[201] = constants[33]*(exp(algebraic[194]/(constants[141]*constants[142]))-exp(algebraic[200]/(constants[141]*constants[142]))) algebraic[203] = constants[127]*algebraic[201] algebraic[185] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[195] = algebraic[189]+algebraic[185] algebraic[190] = constants[141]*constants[142]*log(constants[84]*states[41]) algebraic[202] = algebraic[190] algebraic[204] = constants[34]*(exp(algebraic[195]/(constants[141]*constants[142]))-exp(algebraic[202]/(constants[141]*constants[142]))) algebraic[206] = algebraic[201]-algebraic[204] algebraic[196] = algebraic[190] algebraic[191] = constants[141]*constants[142]*log(constants[85]*states[42]) algebraic[205] = algebraic[191] algebraic[207] = constants[35]*(exp(algebraic[196]/(constants[141]*constants[142]))-exp(algebraic[205]/(constants[141]*constants[142]))) algebraic[209] = algebraic[204]-algebraic[207] algebraic[197] = algebraic[191] algebraic[184] = constants[141]*constants[142]*log(constants[69]*states[11]) algebraic[192] = constants[141]*constants[142]*log(constants[86]*states[43]) algebraic[208] = algebraic[192]+algebraic[184] algebraic[210] = constants[36]*(exp(algebraic[197]/(constants[141]*constants[142]))-exp(algebraic[208]/(constants[141]*constants[142]))) algebraic[212] = -algebraic[210]+algebraic[207] algebraic[186] = constants[141]*constants[142]*log(constants[89]*states[6]) algebraic[198] = algebraic[192]+constants[126]*algebraic[186] algebraic[193] = constants[141]*constants[142]*log(constants[87]*states[44]) algebraic[211] = algebraic[193] algebraic[213] = constants[37]*(exp(algebraic[198]/(constants[141]*constants[142]))-exp(algebraic[211]/(constants[141]*constants[142]))) algebraic[214] = -constants[126]*algebraic[213] algebraic[215] = -algebraic[213]+algebraic[210] algebraic[183] = states[0]/constants[133] algebraic[199] = constants[141]*constants[142]*log(constants[87]*states[44])+constants[128]*constants[143]*algebraic[183] algebraic[216] = constants[141]*constants[142]*log(constants[82]*states[39])+constants[129]*constants[143]*algebraic[183] algebraic[217] = constants[38]*(exp(algebraic[199]/(constants[141]*constants[142]))-exp(algebraic[216]/(constants[141]*constants[142]))) algebraic[218] = algebraic[217]-algebraic[201] algebraic[219] = -algebraic[217]+algebraic[213] algebraic[181] = algebraic[177]+algebraic[175] algebraic[221] = algebraic[210] algebraic[224] = constants[141]*constants[142]*log(constants[69]*states[11]) algebraic[225] = constants[141]*constants[142]*log(constants[90]*states[21]) algebraic[226] = constants[141]*constants[142]*log(constants[91]*states[22]) algebraic[228] = constants[40]*(exp(algebraic[226]/(constants[141]*constants[142]))-exp((algebraic[224]+algebraic[225])/(constants[141]*constants[142]))) algebraic[223] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[227] = constants[39]*(exp((algebraic[223]+algebraic[225])/(constants[141]*constants[142]))-exp(algebraic[226]/(constants[141]*constants[142]))) algebraic[180] = constants[143]*((((((((((((((((((((((((((((((((((((constants[124]*algebraic[141]-constants[120]*algebraic[141])+constants[125]*algebraic[149])-constants[121]*algebraic[149])+constants[123]*algebraic[137])-constants[119]*algebraic[137])+constants[124]*algebraic[145])-constants[120]*algebraic[145])+constants[125]*algebraic[153])-constants[121]*algebraic[153])+constants[123]*algebraic[133])-constants[119]*algebraic[133])+constants[122]*algebraic[129])-constants[118]*algebraic[129])+constants[122]*algebraic[125])-constants[118]*algebraic[125])+constants[122]*algebraic[121])-constants[118]*algebraic[121])+constants[125]*algebraic[151])-constants[121]*algebraic[151])+constants[124]*algebraic[143])-constants[120]*algebraic[143])+constants[123]*algebraic[135])-constants[119]*algebraic[135])+constants[125]*algebraic[147])-constants[121]*algebraic[147])+constants[124]*algebraic[139])-constants[120]*algebraic[139])+constants[123]*algebraic[131])-constants[119]*algebraic[131])+constants[122]*algebraic[127])-constants[118]*algebraic[127])+constants[122]*algebraic[123])-constants[118]*algebraic[123])+constants[122]*algebraic[119])-constants[118]*algebraic[119])-(constants[147]*algebraic[175]+constants[147]*algebraic[177])) algebraic[220] = constants[143]*(-constants[128]*algebraic[217]+constants[129]*algebraic[217]) algebraic[230] = -constants[147]*constants[143]*algebraic[228] algebraic[11] = voi-floor(voi/constants[135])*constants[135] algebraic[28] = custom_piecewise([greater_equal(algebraic[11] , 0.300000) & less_equal(algebraic[11] , 0.300000+constants[136]), (0.00000*0.0300000*constants[133])/constants[136] , True, 0.00000]) algebraic[232] = algebraic[180]+algebraic[220]+algebraic[230]+algebraic[28] algebraic[231] = constants[141]*constants[142]*log(constants[68]*states[1]) algebraic[233] = constants[141]*constants[142]*log(constants[100]*states[45]) algebraic[236] = constants[141]*constants[142]*log(constants[103]*states[48]) algebraic[238] = constants[47]*(exp(algebraic[236]/(constants[141]*constants[142]))-exp((algebraic[233]+constants[164]*algebraic[231])/(constants[141]*constants[142]))) algebraic[234] = constants[141]*constants[142]*log(constants[101]*states[46]) algebraic[240] = constants[48]*(exp((algebraic[233]+constants[163]*algebraic[231])/(constants[141]*constants[142]))-exp(algebraic[234]/(constants[141]*constants[142]))) algebraic[235] = constants[141]*constants[142]*log(constants[102]*states[47]) algebraic[242] = constants[49]*(exp((algebraic[234]+constants[164]*algebraic[231])/(constants[141]*constants[142]))-exp(algebraic[235]/(constants[141]*constants[142]))) algebraic[244] = constants[50]*(exp(algebraic[235]/(constants[141]*constants[142]))-exp((algebraic[236]+constants[163]*algebraic[231])/(constants[141]*constants[142]))) algebraic[68] = constants[141]*constants[142]*log(constants[68]*states[1]) algebraic[70] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[72] = constants[6]*(exp(algebraic[68]/(constants[141]*constants[142]))-exp(algebraic[70]/(constants[141]*constants[142]))) algebraic[182] = (((constants[117])*algebraic[167]+constants[117]*algebraic[171]+constants[117]*algebraic[163]+constants[117]*algebraic[159]+constants[117]*algebraic[155]+constants[117]*algebraic[174])-algebraic[175])-algebraic[177] algebraic[229] = constants[141]*constants[142]*log(constants[60]*states[20]) algebraic[237] = constants[46]*exp(algebraic[236]/(constants[141]*constants[142]))*(exp(algebraic[229]/(constants[141]*constants[142]))-exp(algebraic[231]/(constants[141]*constants[142]))) algebraic[246] = ((constants[164]*algebraic[238]-constants[163]*algebraic[240])-constants[164]*algebraic[242])+constants[163]*algebraic[244] algebraic[254] = constants[141]*constants[142]*log(constants[112]*states[57]) algebraic[259] = constants[141]*constants[142]*log(constants[97]*states[17]) algebraic[260] = algebraic[259]-constants[173] algebraic[261] = algebraic[254]+algebraic[260] algebraic[249] = constants[141]*constants[142]*log(constants[104]*states[49]) algebraic[248] = constants[141]*constants[142]*log(constants[116]*states[24]) algebraic[262] = algebraic[249]+algebraic[248]+algebraic[260] algebraic[279] = constants[59]*(exp(algebraic[261]/(constants[141]*constants[142]))-exp(algebraic[262]/(constants[141]*constants[142]))) algebraic[281] = algebraic[279] algebraic[247] = constants[141]*constants[142]*log(constants[114]*states[26]) algebraic[263] = algebraic[249]+algebraic[247]+algebraic[260] algebraic[250] = constants[141]*constants[142]*log(constants[105]*states[50]) algebraic[264] = algebraic[250]+algebraic[260] algebraic[280] = constants[51]*(exp(algebraic[263]/(constants[141]*constants[142]))-exp(algebraic[264]/(constants[141]*constants[142]))) algebraic[283] = -algebraic[280] algebraic[284] = algebraic[279]-algebraic[280] algebraic[253] = constants[141]*constants[142]*log(constants[108]*states[53]) algebraic[265] = algebraic[253]+algebraic[260] algebraic[245] = constants[141]*constants[142]*log(constants[115]*states[25]) algebraic[255] = constants[141]*constants[142]*log(constants[109]*states[54]) algebraic[266] = algebraic[245]+algebraic[255]+algebraic[260] algebraic[282] = constants[55]*(exp(algebraic[265]/(constants[141]*constants[142]))-exp(algebraic[266]/(constants[141]*constants[142]))) algebraic[286] = algebraic[282] algebraic[243] = constants[141]*constants[142]*log(constants[113]*states[23]) algebraic[267] = algebraic[250]+algebraic[243]+algebraic[260] algebraic[251] = constants[141]*constants[142]*log(constants[106]*states[51]) algebraic[268] = algebraic[251]+algebraic[260] algebraic[285] = constants[53]*(exp(algebraic[267]/(constants[141]*constants[142]))-exp(algebraic[268]/(constants[141]*constants[142]))) algebraic[288] = algebraic[285] algebraic[38] = (constants[141]*constants[142])*log(constants[60]*states[20]) algebraic[48] = (constants[141]*constants[142])*log(constants[61]*states[4]) algebraic[54] = constants[0]*(exp(algebraic[38]/(constants[141]*constants[142]))-exp(algebraic[48]/(constants[141]*constants[142]))) algebraic[222] = -algebraic[204] algebraic[239] = constants[141]*constants[142]*log(constants[61]*states[4]) algebraic[269] = algebraic[250]+constants[165]*algebraic[239]+algebraic[260] algebraic[252] = constants[141]*constants[142]*log(constants[107]*states[52]) algebraic[270] = algebraic[252]+algebraic[260] algebraic[287] = constants[52]*(exp(algebraic[269]/(constants[141]*constants[142]))-exp(algebraic[270]/(constants[141]*constants[142]))) algebraic[289] = -constants[165]*algebraic[287] algebraic[33] = custom_piecewise([greater_equal(algebraic[11] , 0.300000) & less_equal(algebraic[11] , 0.300000+constants[136]), 0.00000*8000.00 , True, 0.00000]) algebraic[291] = algebraic[280]-algebraic[287] algebraic[271] = algebraic[252]+algebraic[260] algebraic[272] = algebraic[253]+constants[167]*algebraic[243]+algebraic[260] algebraic[290] = constants[54]*(exp(algebraic[271]/(constants[141]*constants[142]))-exp(algebraic[272]/(constants[141]*constants[142]))) algebraic[293] = algebraic[287]-algebraic[290] algebraic[294] = algebraic[290]-algebraic[282] algebraic[273] = algebraic[255]+algebraic[260] algebraic[241] = constants[141]*constants[142]*log(constants[60]*states[20]) algebraic[257] = constants[141]*constants[142]*log(constants[110]*states[55]) algebraic[274] = algebraic[257]+constants[166]*algebraic[241]+algebraic[260] algebraic[292] = constants[56]*(exp(algebraic[273]/(constants[141]*constants[142]))-exp(algebraic[274]/(constants[141]*constants[142]))) algebraic[295] = constants[166]*algebraic[292] algebraic[297] = algebraic[282]-algebraic[292] algebraic[275] = algebraic[257]+constants[167]*algebraic[243]+algebraic[260] algebraic[258] = constants[141]*constants[142]*log(constants[111]*states[56]) algebraic[276] = algebraic[258]+algebraic[260] algebraic[296] = constants[57]*(exp(algebraic[275]/(constants[141]*constants[142]))-exp(algebraic[276]/(constants[141]*constants[142]))) algebraic[299] = algebraic[292]-algebraic[296] algebraic[277] = algebraic[258]+algebraic[260] algebraic[278] = algebraic[243]+algebraic[254]+algebraic[260] algebraic[298] = constants[58]*(exp(algebraic[277]/(constants[141]*constants[142]))-exp(algebraic[278]/(constants[141]*constants[142]))) algebraic[300] = ((constants[167]*algebraic[290]-algebraic[285])-constants[167]*algebraic[296])+algebraic[298] algebraic[302] = algebraic[298]-algebraic[279] algebraic[301] = algebraic[296]-algebraic[298] algebraic[1] = states[0]/constants[133] algebraic[3] = states[4]/constants[132] algebraic[4] = states[11]/constants[139] algebraic[5] = states[20]/constants[170] algebraic[6] = states[1]/constants[140] algebraic[7] = states[5]/constants[132] algebraic[8] = -log(states[23]/constants[132], 10) algebraic[9] = states[4]+states[11]+states[3]+states[20]+states[1]+states[22]+2.00000*(states[52]+states[53]+states[54])+states[41]+states[42] algebraic[10] = (constants[157]/(constants[148]*constants[149]*constants[150]))*states[10]*(states[60]+constants[148]) algebraic[12] = constants[141]*constants[142]*log(constants[76]*states[33]) algebraic[13] = constants[141]*constants[142]*log(constants[77]*states[34]) algebraic[14] = constants[141]*constants[142]*log(constants[78]*states[35]) algebraic[15] = constants[141]*constants[142]*log(constants[79]*states[36]) algebraic[16] = states[11]+states[1]+states[33]+states[34]+states[35]+states[36]+states[37]+states[38] algebraic[17] = constants[164]*states[48]+constants[163]*states[46]+(constants[163]+constants[164])*states[47] algebraic[18] = constants[165]*states[52]+constants[165]*states[53]+constants[166]*states[54] algebraic[19] = states[4]/constants[168] algebraic[20] = states[20]/constants[174] algebraic[21] = states[23]/constants[176] algebraic[22] = states[25]/constants[168] algebraic[23] = states[26]/constants[168] algebraic[24] = states[24]/constants[168] algebraic[27] = (constants[157]/(constants[148]*constants[149]*constants[150]))*states[9]*states[59] algebraic[32] = algebraic[10]+algebraic[27] algebraic[37] = constants[161]*(exp((constants[162]*states[58])/constants[148])-1.00000) algebraic[53] = ((constants[161]*algebraic[43])/constants[148])*algebraic[44] algebraic[57] = algebraic[37]+algebraic[53] algebraic[59] = algebraic[32]+algebraic[57] algebraic[256] = constants[59]*(exp(algebraic[254]/(constants[141]*constants[142]))-exp((algebraic[249]+algebraic[248])/(constants[141]*constants[142]))) return algebraic initialGuess0 = None def rootfind_0(voi, constants, rates, states, algebraic): """Calculate values of algebraic variables for DAE""" from scipy.optimize import fsolve global initialGuess0 if initialGuess0 is None: initialGuess0 = ones(5)*0.1 if not iterable(voi): soln = fsolve(residualSN_0, initialGuess0, args=(algebraic, voi, constants, rates, states), xtol=1E-6) initialGuess0 = soln algebraic[43] = soln[0] algebraic[44] = soln[1] algebraic[45] = soln[2] algebraic[46] = soln[3] algebraic[47] = soln[4] else: for (i,t) in enumerate(voi): soln = fsolve(residualSN_0, initialGuess0, args=(algebraic[:,i], voi[i], constants, rates[:i], states[:,i]), xtol=1E-6) initialGuess0 = soln algebraic[43][i] = soln[0] algebraic[44][i] = soln[1] algebraic[45][i] = soln[2] algebraic[46][i] = soln[3] algebraic[47][i] = soln[4] def residualSN_0(algebraicCandidate, algebraic, voi, constants, rates, states): resid = array([0.0] * 5) algebraic[43] = algebraicCandidate[0] algebraic[44] = algebraicCandidate[1] algebraic[45] = algebraicCandidate[2] algebraic[46] = algebraicCandidate[3] algebraic[47] = algebraicCandidate[4] resid[0] = (algebraic[43]-(custom_piecewise([greater(algebraic[44] , 0.00000), constants[158] , True, constants[159]]))) resid[1] = (algebraic[44]-algebraic[45]) resid[2] = (algebraic[46]-(constants[160]/algebraic[43])*states[61]) resid[3] = (algebraic[47]-(constants[160]/algebraic[43])*(states[58]-constants[148])) resid[4] = (algebraic[45]-(-algebraic[46]+algebraic[47])) return resid def custom_piecewise(cases): """Compute result of a piecewise function""" return select(cases[0::2],cases[1::2]) def solve_model(): """Solve model with ODE solver""" from scipy.integrate import ode # Initialise constants and state variables (init_states, constants) = initConsts() # Set timespan to solve over voi = linspace(0, 10, 500) # Construct ODE object to solve r = ode(computeRates) r.set_integrator('vode', method='bdf', atol=1e-06, rtol=1e-06, max_step=1) r.set_initial_value(init_states, voi[0]) r.set_f_params(constants) # Solve model states = array([[0.0] * len(voi)] * sizeStates) states[:,0] = init_states for (i,t) in enumerate(voi[1:]): if r.successful(): r.integrate(t) states[:,i+1] = r.y else: break # Compute algebraic variables algebraic = computeAlgebraic(constants, states, voi) return (voi, states, algebraic) def plot_model(voi, states, algebraic): """Plot variables against variable of integration""" import pylab (legend_states, legend_algebraic, legend_voi, legend_constants) = createLegends() pylab.figure(1) pylab.plot(voi,vstack((states,algebraic)).T) pylab.xlabel(legend_voi) pylab.legend(legend_states + legend_algebraic, loc='best') pylab.show() if __name__ == "__main__": (voi, states, algebraic) = solve_model() plot_model(voi, states, algebraic)