Generated Code
The following is python code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
# Size of variable arrays: sizeAlgebraic = 369 sizeStates = 145 sizeConstants = 410 from math import * from numpy import * def createLegends(): legend_states = [""] * sizeStates legend_rates = [""] * sizeStates legend_algebraic = [""] * sizeAlgebraic legend_voi = "" legend_constants = [""] * sizeConstants legend_constants[324] = "AC47_cyt in component ac (uM)" legend_constants[329] = "AC47_eca in component ac (uM)" legend_constants[320] = "AC56_cav in component ac (uM)" legend_constants[328] = "AC56_cyt in component ac (uM)" legend_constants[296] = "AC_tot in component ac (uM)" legend_constants[0] = "ATP in component ac (uM)" legend_constants[1] = "KmATP in component ac (uM)" legend_constants[2] = "KmGiAC56 in component ac (uM)" legend_constants[3] = "KmGsAC47 in component ac (dimensionless)" legend_constants[4] = "KmGsAC56 in component ac (dimensionless)" legend_constants[5] = "KmGsGiAC56 in component ac (dimensionless)" legend_algebraic[0] = "ac_kAC47_cyt_gsa in component ac (dimensionless)" legend_algebraic[1] = "ac_kAC47_eca_gsa in component ac (dimensionless)" legend_algebraic[2] = "ac_kAC56_cav_gsa in component ac (dimensionless)" legend_algebraic[30] = "ac_kAC56_cyt_gsa in component ac (dimensionless)" legend_constants[6] = "afAC47 in component ac (hertz)" legend_constants[7] = "afAC56 in component ac (hertz)" legend_constants[8] = "basalAC47 in component ac (dimensionless)" legend_constants[9] = "basalAC56 in component ac (dimensionless)" legend_constants[294] = "R_b1_tot in component beta (uM)" legend_states[0] = "Gi_bg in component beta_cav (uM)" legend_states[1] = "Gs_aGTP in component beta_cav (uM)" legend_states[2] = "Gs_aGTP in component beta_cyt (uM)" legend_states[3] = "Gs_aGTP in component beta_eca (uM)" legend_algebraic[68] = "dcAMP_AC47_cyt in component ac (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[67] = "dcAMP_AC47_eca in component ac (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[82] = "dcAMP_AC56_cav in component ac (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[106] = "dcAMP_AC56_cyt in component ac (mol_per_m3_per_s_times_1e_minus_3)" legend_constants[227] = "fATP in component ac (dimensionless)" legend_constants[10] = "f_AC47_eca in component ac (dimensionless)" legend_constants[260] = "f_AC56_AC47 in component ac (dimensionless)" legend_constants[11] = "f_AC56_cav in component ac (dimensionless)" legend_algebraic[32] = "gsi in component ac (dimensionless)" legend_constants[12] = "hGsAC47 in component ac (dimensionless)" legend_constants[13] = "hGsAC56 in component ac (dimensionless)" legend_constants[14] = "hGsGiAC56 in component ac (dimensionless)" legend_algebraic[55] = "kAC47_cyt in component ac (hertz)" legend_algebraic[31] = "kAC47_eca in component ac (hertz)" legend_algebraic[56] = "kAC56_cav in component ac (hertz)" legend_algebraic[84] = "kAC56_cyt in component ac (hertz)" legend_constants[15] = "vGsGiAC56 in component ac (dimensionless)" legend_constants[319] = "vr_cav in component cell (dimensionless)" legend_constants[321] = "vr_cyt in component cell (dimensionless)" legend_constants[326] = "vr_eca in component cell (dimensionless)" legend_constants[16] = "ICaL_akap in component akap_sig (uM)" legend_constants[391] = "ICaL_akapf in component akap_sig (uM)" legend_constants[392] = "ICaL_arn in component akap_sig (uM)" legend_constants[393] = "ICaL_arp in component akap_sig (uM)" legend_constants[17] = "ICaL_tot in component akap_sig (uM)" legend_constants[281] = "ICaLf in component akap_sig (uM)" legend_states[4] = "ICaLp in component akap_sig (uM)" legend_constants[18] = "Ka_ical in component akap_sig (uM)" legend_constants[19] = "Ka_ryr in component akap_sig (uM)" legend_constants[20] = "Ki in component akap_sig (uM)" legend_constants[21] = "Kp_ical in component akap_sig (uM)" legend_constants[22] = "Kp_ryr in component akap_sig (uM)" legend_constants[23] = "Kr in component akap_sig (uM)" legend_constants[24] = "Li in component akap_sig (uM)" legend_constants[25] = "Lr in component akap_sig (uM)" legend_constants[26] = "Mi in component akap_sig (uM)" legend_constants[27] = "Mr in component akap_sig (uM)" legend_constants[380] = "PKA_cav in component pka (uM)" legend_constants[390] = "PKAf in component akap_sig (uM)" legend_constants[28] = "PP1_cav in component pp1 (uM)" legend_constants[291] = "PP1f_cav in component akap_sig (uM)" legend_constants[29] = "RyR_akap in component akap_sig (uM)" legend_constants[394] = "RyR_akapf in component akap_sig (uM)" legend_constants[395] = "RyR_arn in component akap_sig (uM)" legend_constants[396] = "RyR_arp in component akap_sig (uM)" legend_constants[30] = "RyR_tot in component akap_sig (uM)" legend_constants[293] = "RyRf in component akap_sig (uM)" legend_states[5] = "RyRp in component akap_sig (uM)" legend_constants[273] = "akap_sig_ICaLf_sum in component akap_sig (dimensionless)" legend_algebraic[3] = "akap_sig_ICaLp_dif in component akap_sig (uM)" legend_constants[387] = "akap_sig_PKAf_arg in component akap_sig (dimensionless)" legend_constants[381] = "akap_sig_PKAf_b in component akap_sig (uM)" legend_constants[382] = "akap_sig_PKAf_c in component akap_sig (mM2_times_1e_minus_6)" legend_constants[383] = "akap_sig_PKAf_d in component akap_sig (mM3_times_1e_minus_9)" legend_constants[388] = "akap_sig_PKAf_mag in component akap_sig (uM)" legend_constants[384] = "akap_sig_PKAf_rr in component akap_sig (mol6_per_m18_times_1e_minus_18)" legend_constants[389] = "akap_sig_PKAf_x in component akap_sig (dimensionless)" legend_constants[385] = "akap_sig_PKAf_yi in component akap_sig (mM3_times_1e_minus_9)" legend_constants[386] = "akap_sig_PKAf_yr in component akap_sig (mM3_times_1e_minus_9)" legend_constants[288] = "akap_sig_PP1f_cav_arg in component akap_sig (dimensionless)" legend_constants[282] = "akap_sig_PP1f_cav_b in component akap_sig (uM)" legend_constants[283] = "akap_sig_PP1f_cav_c in component akap_sig (mM2_times_1e_minus_6)" legend_constants[284] = "akap_sig_PP1f_cav_d in component akap_sig (mM3_times_1e_minus_9)" legend_constants[289] = "akap_sig_PP1f_cav_mag in component akap_sig (uM)" legend_constants[285] = "akap_sig_PP1f_cav_rr in component akap_sig (mol6_per_m18_times_1e_minus_18)" legend_constants[290] = "akap_sig_PP1f_cav_x in component akap_sig (dimensionless)" legend_constants[286] = "akap_sig_PP1f_cav_yi in component akap_sig (mM3_times_1e_minus_9)" legend_constants[287] = "akap_sig_PP1f_cav_yr in component akap_sig (mM3_times_1e_minus_9)" legend_constants[292] = "akap_sig_RyRf_sum in component akap_sig (dimensionless)" legend_algebraic[4] = "akap_sig_RyRp_dif in component akap_sig (uM)" legend_algebraic[10] = "fp_ICaL in component akap_sig (dimensionless)" legend_algebraic[38] = "fp_RyR in component akap_sig (dimensionless)" legend_constants[31] = "ka_ical in component akap_sig (hertz)" legend_constants[32] = "ka_ryr in component akap_sig (hertz)" legend_constants[33] = "kp_ical in component akap_sig (hertz)" legend_constants[34] = "kp_ryr in component akap_sig (hertz)" legend_states[6] = "C in component pka_cav (uM)" legend_voi = "time in component engine (ms)" legend_constants[35] = "Gi_tot in component beta (uM)" legend_constants[295] = "Gs_tot in component beta (uM)" legend_constants[297] = "R_b2_tot in component beta (uM)" legend_constants[36] = "f_Gi_cav in component beta (dimensionless)" legend_constants[215] = "f_Gi_eca in component beta (dimensionless)" legend_constants[37] = "f_Gs_cav in component beta (dimensionless)" legend_constants[216] = "f_Gs_cyt in component beta (dimensionless)" legend_constants[38] = "f_Gs_eca in component beta (dimensionless)" legend_constants[39] = "f_Rb1_cav in component beta (dimensionless)" legend_constants[298] = "f_Rb1_cyt in component beta (dimensionless)" legend_constants[40] = "f_Rb1_eca in component beta (dimensionless)" legend_constants[41] = "f_Rb2_cav in component beta (dimensionless)" legend_constants[299] = "f_Rb2_eca in component beta (dimensionless)" legend_constants[42] = "k_act1_Gi in component beta (hertz)" legend_constants[43] = "k_act1_Gs in component beta (hertz)" legend_constants[44] = "k_act2_Gi in component beta (hertz)" legend_constants[45] = "k_act2_Gs in component beta (hertz)" legend_constants[46] = "k_b1_c in component beta (uM)" legend_constants[47] = "k_b1_h in component beta (uM)" legend_constants[48] = "k_b1_l in component beta (uM)" legend_constants[49] = "k_b2_a in component beta (uM)" legend_constants[50] = "k_b2_c in component beta (uM)" legend_constants[51] = "k_b2_f in component beta (uM)" legend_constants[52] = "k_b2_h in component beta (uM)" legend_constants[53] = "k_b2_l in component beta (uM)" legend_constants[54] = "k_b2_n in component beta (uM)" legend_constants[217] = "k_grk_dp in component beta (hertz)" legend_constants[218] = "k_grk_p in component beta (hertz)" legend_constants[219] = "k_hydr_Gi in component beta (hertz)" legend_constants[55] = "k_hydr_Gs in component beta (hertz)" legend_constants[257] = "k_pka_dp in component beta (hertz)" legend_constants[220] = "k_pka_p in component beta (per_mM_per_ms)" legend_constants[221] = "k_reas_Gi in component beta (per_mM_per_ms)" legend_constants[56] = "k_reas_Gs in component beta (per_mM_per_ms)" legend_constants[57] = "rate_bds in component beta (per_mM_per_ms)" legend_constants[58] = "GRK in component beta_cav (dimensionless)" legend_states[7] = "Gi_aGDP in component beta_cav (uM)" legend_states[8] = "Gi_aGTP in component beta_cav (uM)" legend_algebraic[5] = "Gi_abg in component beta_cav (uM)" legend_algebraic[86] = "Gi_f in component beta_cav (uM)" legend_states[9] = "Gs_aGDP in component beta_cav (uM)" legend_algebraic[6] = "Gs_abg in component beta_cav (uM)" legend_states[10] = "Gs_bg in component beta_cav (uM)" legend_algebraic[156] = "Gs_f in component beta_cav (uM)" legend_algebraic[188] = "LRGs_tot in component beta_cav (uM)" legend_algebraic[164] = "LRb1 in component beta_cav (uM)" legend_algebraic[168] = "LRb1Gs in component beta_cav (uM)" legend_algebraic[180] = "LRb2 in component beta_cav (uM)" legend_algebraic[108] = "LRb2Gi in component beta_cav (uM)" legend_algebraic[184] = "LRb2Gs in component beta_cav (uM)" legend_algebraic[196] = "RGs_tot in component beta_cav (uM)" legend_constants[322] = "R_b1_tot in component beta_cav (uM)" legend_constants[323] = "R_b2_tot in component beta_cav (uM)" legend_algebraic[172] = "Rb1Gs in component beta_cav (uM)" legend_algebraic[160] = "Rb1_f in component beta_cav (uM)" legend_states[11] = "Rb1_grk_tot in component beta_cav (uM)" legend_algebraic[33] = "Rb1_np_tot in component beta_cav (uM)" legend_states[12] = "Rb1_pka_tot in component beta_cav (uM)" legend_algebraic[98] = "Rb2Gi in component beta_cav (uM)" legend_algebraic[192] = "Rb2Gs in component beta_cav (uM)" legend_algebraic[176] = "Rb2_f in component beta_cav (uM)" legend_states[13] = "Rb2_grk_tot in component beta_cav (uM)" legend_algebraic[57] = "Rb2_np_tot in component beta_cav (uM)" legend_algebraic[70] = "Rb2_pka_f in component beta_cav (uM)" legend_states[14] = "Rb2_pka_tot in component beta_cav (uM)" legend_constants[223] = "beta_cav_Gs_f_a in component beta_cav (mol4_per_m12_times_1e_minus_12)" legend_algebraic[131] = "beta_cav_Gs_f_arg in component beta_cav (dimensionless)" legend_algebraic[69] = "beta_cav_Gs_f_b in component beta_cav (uM)" legend_algebraic[85] = "beta_cav_Gs_f_c in component beta_cav (mM2_times_1e_minus_6)" legend_constants[258] = "beta_cav_Gs_f_c11 in component beta_cav (mol5_per_m15_times_1e_minus_15)" legend_constants[271] = "beta_cav_Gs_f_c22 in component beta_cav (mol5_per_m15_times_1e_minus_15)" legend_constants[279] = "beta_cav_Gs_f_c33 in component beta_cav (mol6_per_m18_times_1e_minus_18)" legend_algebraic[97] = "beta_cav_Gs_f_d in component beta_cav (mM3_times_1e_minus_9)" legend_algebraic[147] = "beta_cav_Gs_f_i in component beta_cav (uM)" legend_algebraic[137] = "beta_cav_Gs_f_mag in component beta_cav (uM)" legend_algebraic[151] = "beta_cav_Gs_f_r in component beta_cav (uM)" legend_algebraic[109] = "beta_cav_Gs_f_rr in component beta_cav (mol6_per_m18_times_1e_minus_18)" legend_algebraic[142] = "beta_cav_Gs_f_x in component beta_cav (dimensionless)" legend_algebraic[117] = "beta_cav_Gs_f_yi in component beta_cav (mM3_times_1e_minus_9)" legend_algebraic[125] = "beta_cav_Gs_f_yr in component beta_cav (mM3_times_1e_minus_9)" legend_constants[222] = "beta_cav_Rb2_pka_f_a in component beta_cav (uM)" legend_algebraic[34] = "beta_cav_Rb2_pka_f_b in component beta_cav (mM2_times_1e_minus_6)" legend_algebraic[58] = "beta_cav_Rb2_pka_f_c in component beta_cav (mM3_times_1e_minus_9)" legend_constants[59] = "L in component iso (uM)" legend_constants[60] = "k_GsAct_b2 in component beta_cav (dimensionless)" legend_constants[61] = "GRK in component beta_cyt (dimensionless)" legend_states[15] = "Gs_aGDP in component beta_cyt (uM)" legend_algebraic[7] = "Gs_abg in component beta_cyt (uM)" legend_states[16] = "Gs_bg in component beta_cyt (uM)" legend_algebraic[99] = "Gs_f in component beta_cyt (uM)" legend_algebraic[110] = "LRb1Gs_np in component beta_cyt (uM)" legend_algebraic[118] = "LRb1_np in component beta_cyt (uM)" legend_constants[325] = "R_b1_tot in component beta_cyt (uM)" legend_algebraic[119] = "Rb1Gs_np in component beta_cyt (uM)" legend_states[17] = "Rb1_grk_tot in component beta_cyt (uM)" legend_algebraic[87] = "Rb1_np_f in component beta_cyt (uM)" legend_algebraic[35] = "Rb1_np_tot in component beta_cyt (uM)" legend_states[18] = "Rb1_pka_tot in component beta_cyt (uM)" legend_constants[224] = "beta_cyt_Rb1_np_f_a in component beta_cyt (uM)" legend_algebraic[59] = "beta_cyt_Rb1_np_f_b in component beta_cyt (mM2_times_1e_minus_6)" legend_algebraic[71] = "beta_cyt_Rb1_np_f_c in component beta_cyt (mM3_times_1e_minus_9)" legend_states[19] = "C in component pka_cyt (uM)" legend_constants[62] = "GRK in component beta_eca (dimensionless)" legend_states[20] = "Gi_aGDP in component beta_eca (uM)" legend_states[21] = "Gi_aGTP in component beta_eca (uM)" legend_algebraic[8] = "Gi_abg in component beta_eca (uM)" legend_states[22] = "Gi_bg in component beta_eca (uM)" legend_algebraic[89] = "Gi_f in component beta_eca (uM)" legend_states[23] = "Gs_aGDP in component beta_eca (uM)" legend_algebraic[9] = "Gs_abg in component beta_eca (uM)" legend_states[24] = "Gs_bg in component beta_eca (uM)" legend_algebraic[157] = "Gs_f in component beta_eca (uM)" legend_algebraic[189] = "LRGs_tot in component beta_eca (uM)" legend_algebraic[165] = "LRb1 in component beta_eca (uM)" legend_algebraic[169] = "LRb1Gs in component beta_eca (uM)" legend_algebraic[181] = "LRb2 in component beta_eca (uM)" legend_algebraic[111] = "LRb2Gi in component beta_eca (uM)" legend_algebraic[185] = "LRb2Gs in component beta_eca (uM)" legend_algebraic[197] = "RGs_tot in component beta_eca (uM)" legend_constants[330] = "R_b1_tot in component beta_eca (uM)" legend_constants[331] = "R_b2_tot in component beta_eca (uM)" legend_algebraic[173] = "Rb1Gs in component beta_eca (uM)" legend_algebraic[161] = "Rb1_f in component beta_eca (uM)" legend_states[25] = "Rb1_grk_tot in component beta_eca (uM)" legend_algebraic[36] = "Rb1_np_tot in component beta_eca (uM)" legend_states[26] = "Rb1_pka_tot in component beta_eca (uM)" legend_algebraic[101] = "Rb2Gi in component beta_eca (uM)" legend_algebraic[193] = "Rb2Gs in component beta_eca (uM)" legend_algebraic[177] = "Rb2_f in component beta_eca (uM)" legend_states[27] = "Rb2_grk_tot in component beta_eca (uM)" legend_algebraic[60] = "Rb2_np_tot in component beta_eca (uM)" legend_algebraic[73] = "Rb2_pka_f in component beta_eca (uM)" legend_states[28] = "Rb2_pka_tot in component beta_eca (uM)" legend_constants[226] = "beta_eca_Gs_f_a in component beta_eca (mol4_per_m12_times_1e_minus_12)" legend_algebraic[132] = "beta_eca_Gs_f_arg in component beta_eca (dimensionless)" legend_algebraic[72] = "beta_eca_Gs_f_b in component beta_eca (uM)" legend_algebraic[88] = "beta_eca_Gs_f_c in component beta_eca (mM2_times_1e_minus_6)" legend_constants[259] = "beta_eca_Gs_f_c11 in component beta_eca (mol5_per_m15_times_1e_minus_15)" legend_constants[272] = "beta_eca_Gs_f_c22 in component beta_eca (mol5_per_m15_times_1e_minus_15)" legend_constants[280] = "beta_eca_Gs_f_c33 in component beta_eca (mol6_per_m18_times_1e_minus_18)" legend_algebraic[100] = "beta_eca_Gs_f_d in component beta_eca (mM3_times_1e_minus_9)" legend_algebraic[148] = "beta_eca_Gs_f_i in component beta_eca (uM)" legend_algebraic[138] = "beta_eca_Gs_f_mag in component beta_eca (uM)" legend_algebraic[152] = "beta_eca_Gs_f_r in component beta_eca (uM)" legend_algebraic[112] = "beta_eca_Gs_f_rr in component beta_eca (mol6_per_m18_times_1e_minus_18)" legend_algebraic[143] = "beta_eca_Gs_f_x in component beta_eca (dimensionless)" legend_algebraic[120] = "beta_eca_Gs_f_yi in component beta_eca (mM3_times_1e_minus_9)" legend_algebraic[126] = "beta_eca_Gs_f_yr in component beta_eca (mM3_times_1e_minus_9)" legend_constants[225] = "beta_eca_Rb2_pka_f_a in component beta_eca (uM)" legend_algebraic[37] = "beta_eca_Rb2_pka_f_b in component beta_eca (mM2_times_1e_minus_6)" legend_algebraic[61] = "beta_eca_Rb2_pka_f_c in component beta_eca (mM3_times_1e_minus_9)" legend_constants[63] = "k_GsAct_b2 in component beta_eca (dimensionless)" legend_states[29] = "C in component pka_eca (uM)" legend_constants[302] = "AF in component cell (m2_mol_per_s_per_A_times_1e_minus_4)" legend_algebraic[186] = "Ca in component calcium (mM)" legend_algebraic[102] = "Ca_CaL in component calcium (mM)" legend_algebraic[127] = "Ca_jsr in component calcium (mM)" legend_states[30] = "Ca_nsr in component calcium (mM)" legend_algebraic[149] = "Ca_sr in component calcium (mM)" legend_algebraic[305] = "ICa_tot in component calcium (uA_per_cm2)" legend_algebraic[217] = "ICab in component icab (uA_per_cm2)" legend_algebraic[297] = "INaCa in component inaca (uA_per_cm2)" legend_algebraic[291] = "INaCaSR in component inaca (uA_per_cm2)" legend_algebraic[201] = "Idiff_Ca in component diff (mM_per_ms)" legend_algebraic[209] = "Idiff_sr in component diff (mM_per_ms)" legend_algebraic[304] = "IpCa in component ipca (uA_per_cm2)" legend_algebraic[319] = "Irel in component irel (mM_per_ms)" legend_algebraic[211] = "Itr in component diff (mM_per_ms)" legend_algebraic[328] = "Iup in component iup (mM_per_ms)" legend_constants[64] = "Ka_tni in component calcium (uM)" legend_constants[65] = "Kp_tni in component calcium (uM)" legend_constants[66] = "PP2A in component pp1 (uM)" legend_constants[228] = "bar_sum in component calcium (mM)" legend_constants[67] = "bsl_bar in component calcium (mM)" legend_constants[68] = "bsl_km in component calcium (mM)" legend_constants[69] = "bsr_bar in component calcium (mM)" legend_constants[70] = "bsr_km in component calcium (mM)" legend_algebraic[62] = "calcium_Ca_CaL_b in component calcium (mM)" legend_algebraic[74] = "calcium_Ca_CaL_c in component calcium (mM2)" legend_algebraic[90] = "calcium_Ca_CaL_d in component calcium (mM3)" legend_algebraic[178] = "calcium_Ca_b in component calcium (mM)" legend_algebraic[182] = "calcium_Ca_c in component calcium (mM2)" legend_algebraic[170] = "calcium_Ca_d in component calcium (mM3)" legend_algebraic[113] = "calcium_Ca_jsr_b in component calcium (mM)" legend_algebraic[121] = "calcium_Ca_jsr_c in component calcium (mM2)" legend_constants[312] = "calcium_Ca_nsr_r1 in component calcium (dimensionless)" legend_algebraic[133] = "calcium_Ca_sr_b in component calcium (mM)" legend_algebraic[139] = "calcium_Ca_sr_c in component calcium (mM2)" legend_algebraic[144] = "calcium_Ca_sr_d in component calcium (mM3)" legend_algebraic[153] = "calcium_fhat_val in component calcium (dimensionless)" legend_constants[305] = "calcium_uCa_CaL_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5)" legend_constants[315] = "calcium_uCa_CaL_r2 in component calcium (dimensionless)" legend_constants[308] = "calcium_uCa_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5)" legend_constants[313] = "calcium_uCa_r2 in component calcium (dimensionless)" legend_constants[316] = "calcium_uCa_sr_r1 in component calcium (mol_per_m_per_s_per_A_times_1e5)" legend_constants[317] = "calcium_uCa_sr_r2 in component calcium (dimensionless)" legend_constants[71] = "cbar in component calcium (mM)" legend_constants[72] = "csqn_bar in component calcium (mM)" legend_constants[73] = "csqn_km in component calcium (mM)" legend_states[31] = "f_tni in component calcium (dimensionless)" legend_algebraic[158] = "fhat in component calcium (dimensionless)" legend_algebraic[248] = "ICaL in component ical (uA_per_cm2)" legend_constants[74] = "ka_tni in component calcium (hertz)" legend_constants[75] = "kc in component calcium (mM)" legend_constants[230] = "km_pro in component calcium (mM2)" legend_constants[261] = "km_sum in component calcium (mM)" legend_constants[76] = "kp_tni in component calcium (hertz)" legend_algebraic[166] = "kpro in component calcium (mM2)" legend_algebraic[174] = "ksum in component calcium (mM)" legend_algebraic[162] = "kt in component calcium (mM)" legend_constants[77] = "ktn in component calcium (mM)" legend_constants[229] = "ktp in component calcium (mM)" legend_constants[318] = "r3 in component calcium (dimensionless)" legend_constants[262] = "ss_pro in component calcium (mM2)" legend_constants[274] = "ss_sum in component calcium (mM)" legend_constants[78] = "tbar in component calcium (mM)" legend_states[32] = "uCa in component calcium (mM)" legend_states[33] = "uCa_CaL in component calcium (mM)" legend_states[34] = "uCa_jsr in component calcium (mM)" legend_states[35] = "uCa_sr in component calcium (mM)" legend_constants[304] = "v_CaL in component cell (uL)" legend_constants[307] = "v_cyt in component cell (uL)" legend_constants[310] = "v_jsr in component cell (uL)" legend_constants[311] = "v_nsr in component cell (uL)" legend_constants[314] = "v_sr in component cell (uL)" legend_constants[79] = "CaMK0 in component camk (dimensionless)" legend_constants[80] = "K in component camk (dimensionless)" legend_constants[81] = "Km in component camk (mM)" legend_constants[82] = "PP1_eca in component pp1 (uM)" legend_algebraic[206] = "PP1_tot in component camk (uM)" legend_algebraic[203] = "PP1f_cyt in component pp1 (uM)" legend_algebraic[194] = "active in component camk (dimensionless)" legend_algebraic[190] = "bound in component camk (dimensionless)" legend_algebraic[199] = "c in component camk (dimensionless)" legend_algebraic[200] = "camk_f_ryr_d in component camk (dimensionless)" legend_constants[83] = "camk_trap_alpha in component camk (mS_per_uF)" legend_constants[84] = "camk_trap_beta in component camk (mS_per_uF)" legend_states[36] = "f_ical in component camk (dimensionless)" legend_states[37] = "f_ik1 in component camk (dimensionless)" legend_states[38] = "f_ina in component camk (dimensionless)" legend_states[39] = "f_ito in component camk (dimensionless)" legend_states[40] = "f_plb in component camk (dimensionless)" legend_states[41] = "f_ryr in component camk (dimensionless)" legend_constants[231] = "tau_cal in component camk (ms)" legend_constants[232] = "tau_ik1 in component camk (ms)" legend_constants[233] = "tau_ina in component camk (ms)" legend_constants[234] = "tau_ito in component camk (ms)" legend_constants[85] = "tau_plb in component camk (ms)" legend_constants[86] = "tau_ryr in component camk (ms)" legend_states[42] = "trap in component camk (dimensionless)" legend_states[43] = "cAMP_cav in component camp (uM)" legend_states[44] = "cAMP_cyt in component camp (uM)" legend_states[45] = "cAMP_eca in component camp (uM)" legend_algebraic[95] = "camp_cAMP_cav_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[105] = "camp_cAMP_cav_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[135] = "camp_cAMP_cav_pde in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[115] = "camp_cAMP_cyt_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[123] = "camp_cAMP_cyt_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[146] = "camp_cAMP_cyt_pde in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[83] = "camp_cAMP_eca_j1 in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[96] = "camp_cAMP_eca_j2 in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[124] = "camp_cAMP_eca_pde in component camp (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[114] = "dcAMP_PDE2_cav in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[129] = "dcAMP_PDE2_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[107] = "dcAMP_PDE2_eca in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[122] = "dcAMP_PDE3_cav in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[136] = "dcAMP_PDE3_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[128] = "dcAMP_PDE4_cav in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[141] = "dcAMP_PDE4_cyt in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[116] = "dcAMP_PDE4_eca in component pde (mol_per_m3_per_s_times_1e_minus_3)" legend_constants[235] = "j_cav_cyt in component camp (m3_per_s_times_1e_minus_9)" legend_constants[236] = "j_cav_eca in component camp (m3_per_s_times_1e_minus_9)" legend_constants[237] = "j_eca_cyt in component camp (m3_per_s_times_1e_minus_9)" legend_algebraic[145] = "dcAMP in component pka_cav (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[154] = "dcAMP in component pka_cyt (mol_per_m3_per_s_times_1e_minus_3)" legend_algebraic[134] = "dcAMP in component pka_eca (mol_per_m3_per_s_times_1e_minus_3)" legend_constants[306] = "v_cav in component cell (uL)" legend_constants[309] = "v_eca in component cell (uL)" legend_constants[87] = "F in component phys (C_per_mol)" legend_constants[301] = "capArea in component cell (cm2)" legend_constants[300] = "geoArea in component cell (cm2)" legend_constants[88] = "length in component cell (cm)" legend_constants[89] = "pi in component cell (dimensionless)" legend_constants[90] = "radius in component cell (cm)" legend_constants[303] = "volume in component cell (uL)" legend_algebraic[333] = "CTKCl in component ctkcl (mM_per_ms)" legend_algebraic[354] = "CTNaCl in component ctnacl (mM_per_ms)" legend_states[46] = "Cl in component chloride (mM)" legend_states[47] = "Cl_sr in component chloride (mM)" legend_algebraic[323] = "IClCa in component iclca (uA_per_cm2)" legend_algebraic[330] = "ICl_tot in component chloride (uA_per_cm2)" legend_algebraic[329] = "IClb in component iclb (uA_per_cm2)" legend_algebraic[204] = "Idiff_Cl in component diff (mM_per_ms)" legend_constants[332] = "chloride_Cl_r1 in component chloride (mol_per_m_per_s_per_A_times_1e5)" legend_constants[337] = "chloride_Cl_r2 in component chloride (dimensionless)" legend_constants[341] = "chloride_Cl_sr_r1 in component chloride (mol_per_m_per_s_per_A_times_1e5)" legend_algebraic[327] = "ECl in component nernst (mV)" legend_algebraic[331] = "EK in component nernst (mV)" legend_constants[91] = "KClBar in component ctkcl (mM_per_ms)" legend_algebraic[332] = "ctkcl_CTKCl_z1 in component ctkcl (mV)" legend_constants[92] = "ctkcl_CTKCl_z2 in component ctkcl (mV)" legend_algebraic[352] = "ENa in component nernst (mV)" legend_constants[93] = "NaClBar in component ctnacl (mM_per_ms)" legend_algebraic[353] = "ctnacl_CTNaCl_z1 in component ctnacl (g4_m8_per_s12_per_A4)" legend_constants[346] = "ctnacl_CTNaCl_z2 in component ctnacl (g4_m8_per_s12_per_A4)" legend_algebraic[207] = "Idiff_Na in component diff (mM_per_ms)" legend_states[48] = "Na in component sodium (mM)" legend_states[49] = "Na_sr in component sodium (mM)" legend_constants[94] = "tau in component diff (ms)" legend_constants[95] = "tau_sr in component diff (ms)" legend_constants[96] = "tau_tr in component diff (ms)" legend_constants[97] = "Cao in component extra (mM)" legend_constants[98] = "Clo in component extra (mM)" legend_constants[99] = "Ko in component extra (mM)" legend_constants[100] = "Nao in component extra (mM)" legend_constants[263] = "FRT in component phys (per_mV)" legend_states[50] = "V in component membrane (mV)" legend_algebraic[215] = "efrt in component icab (dimensionless)" legend_constants[101] = "pCab in component icab (cm_per_s)" legend_algebraic[213] = "vfrt in component icab (dimensionless)" legend_algebraic[223] = "f_hat in component ical (dimensionless)" legend_algebraic[235] = "ICaL in component ical_camk (uA_per_cm2)" legend_constants[397] = "ical_f_hat_ratio in component ical (dimensionless)" legend_algebraic[220] = "ical_f_hat_val in component ical (dimensionless)" legend_algebraic[246] = "ICaL in component ical_np (uA_per_cm2)" legend_states[51] = "C in component ical_camk (dimensionless)" legend_states[52] = "CI in component ical_camk (dimensionless)" legend_states[53] = "CIs in component ical_camk (dimensionless)" legend_states[54] = "Cs in component ical_camk (dimensionless)" legend_constants[275] = "FFRT in component phys (s4_A2_per_g_per_m2_per_mol)" legend_algebraic[232] = "IBar in component ical_camk (uA_per_cm2)" legend_states[55] = "O in component ical_camk (dimensionless)" legend_states[56] = "OI in component ical_camk (dimensionless)" legend_states[57] = "OIs in component ical_camk (dimensionless)" legend_states[58] = "Os in component ical_camk (dimensionless)" legend_algebraic[226] = "PCa in component ical_camk (cm_per_s)" legend_algebraic[155] = "ac_inf in component ical_camk (dimensionless)" legend_algebraic[195] = "ac_tau in component ical_np (ms)" legend_algebraic[202] = "alpha in component ical_camk (mS_per_uF)" legend_algebraic[205] = "beta in component ical_camk (mS_per_uF)" legend_algebraic[183] = "delta in component ical_camk (mS_per_uF)" legend_algebraic[242] = "delta1 in component ical_camk (mS_per_uF)" legend_algebraic[159] = "delta_tau in component ical_camk (dimensionless)" legend_algebraic[229] = "ical_camk_IBar_vv in component ical_camk (dimensionless)" legend_algebraic[230] = "ical_camk_delta1_xs_cor in component ical_camk (mS_per_uF)" legend_algebraic[236] = "ical_camk_delta1_y_cor in component ical_camk (mS_per_uF)" legend_algebraic[214] = "in_a in component ical_np (dimensionless)" legend_algebraic[216] = "in_b in component ical_np (dimensionless)" legend_algebraic[167] = "in_hi_inf in component ical_camk (dimensionless)" legend_algebraic[218] = "in_hi_tau in component ical_camk (ms)" legend_algebraic[163] = "in_inf in component ical_camk (dimensionless)" legend_algebraic[171] = "in_lo_inf in component ical_camk (dimensionless)" legend_algebraic[221] = "in_lo_tau in component ical_camk (ms)" legend_algebraic[187] = "inca in component ical_camk (dimensionless)" legend_algebraic[175] = "ss_cal_10 in component ical_camk (dimensionless)" legend_algebraic[179] = "ss_cal_4 in component ical_camk (dimensionless)" legend_constants[102] = "theta in component ical_camk (mS_per_uF)" legend_constants[103] = "theta1 in component ical_camk (mS_per_uF)" legend_algebraic[224] = "x in component ical_camk (mS_per_uF)" legend_algebraic[227] = "xs in component ical_camk (mS_per_uF)" legend_algebraic[233] = "y in component ical_camk (mS_per_uF)" legend_algebraic[239] = "ys in component ical_camk (mS_per_uF)" legend_states[59] = "C in component ical_np (dimensionless)" legend_states[60] = "CI in component ical_np (dimensionless)" legend_states[61] = "CIs in component ical_np (dimensionless)" legend_states[62] = "Cs in component ical_np (dimensionless)" legend_algebraic[244] = "IBar in component ical_np (uA_per_cm2)" legend_states[63] = "O in component ical_np (dimensionless)" legend_states[64] = "OI in component ical_np (dimensionless)" legend_states[65] = "OIs in component ical_np (dimensionless)" legend_states[66] = "Os in component ical_np (dimensionless)" legend_algebraic[238] = "PCa in component ical_np (cm_per_s)" legend_algebraic[191] = "ac_inf in component ical_np (dimensionless)" legend_algebraic[208] = "alpha in component ical_np (mS_per_uF)" legend_algebraic[210] = "beta in component ical_np (mS_per_uF)" legend_algebraic[237] = "delta in component ical_np (mS_per_uF)" legend_algebraic[260] = "delta1 in component ical_np (mS_per_uF)" legend_algebraic[212] = "delta_tau in component ical_np (dimensionless)" legend_algebraic[241] = "ical_np_IBar_vv in component ical_np (dimensionless)" legend_algebraic[249] = "ical_np_delta1_xs_cor in component ical_np (mS_per_uF)" legend_algebraic[254] = "ical_np_delta1_y_cor in component ical_np (mS_per_uF)" legend_algebraic[222] = "in_hi_inf in component ical_np (dimensionless)" legend_algebraic[243] = "in_hi_tau in component ical_np (ms)" legend_algebraic[219] = "in_inf in component ical_np (dimensionless)" legend_algebraic[225] = "in_lo_inf in component ical_np (dimensionless)" legend_algebraic[228] = "in_lo_tau in component ical_np (ms)" legend_algebraic[240] = "inca in component ical_np (dimensionless)" legend_algebraic[231] = "ss_cal_10 in component ical_np (dimensionless)" legend_algebraic[234] = "ss_cal_4 in component ical_np (dimensionless)" legend_constants[104] = "theta in component ical_np (mS_per_uF)" legend_constants[105] = "theta1 in component ical_np (mS_per_uF)" legend_algebraic[245] = "x in component ical_np (mS_per_uF)" legend_algebraic[247] = "xs in component ical_np (mS_per_uF)" legend_algebraic[251] = "y in component ical_np (mS_per_uF)" legend_algebraic[257] = "ys in component ical_np (mS_per_uF)" legend_constants[106] = "Gbar in component iclb (mS_per_cm2)" legend_algebraic[253] = "IClCa_bar in component iclca (uA_per_cm2)" legend_algebraic[317] = "Irel_pure in component irel (mM_per_ms)" legend_algebraic[320] = "KClCa in component iclca (dimensionless)" legend_constants[107] = "PCl in component iclca (cm_per_s)" legend_states[67] = "i2 in component iclca (dimensionless)" legend_algebraic[11] = "iclca_i2_alpha in component iclca (dimensionless)" legend_algebraic[39] = "iclca_i2_beta in component iclca (dimensionless)" legend_constants[108] = "kCaCl in component iclca (mM_per_ms)" legend_constants[109] = "tau in component iclca (ms)" legend_algebraic[250] = "vexp in component iclca (dimensionless)" legend_constants[239] = "Gbar in component ik1 (mS_per_cm2)" legend_algebraic[339] = "IK1 in component ik1 (uA_per_cm2)" legend_algebraic[338] = "IK1_camk in component ik1 (uA_per_cm2)" legend_algebraic[337] = "IK1_np in component ik1 (uA_per_cm2)" legend_algebraic[335] = "ik1_IK1_np_alpha in component ik1 (dimensionless)" legend_algebraic[336] = "ik1_IK1_np_beta in component ik1 (dimensionless)" legend_algebraic[334] = "ik1_IK1_np_vv in component ik1 (mV)" legend_constants[240] = "GKr in component ikr (mS_per_cm2)" legend_algebraic[340] = "IKr in component ikr (uA_per_cm2)" legend_states[68] = "ac in component ikr (dimensionless)" legend_algebraic[12] = "ikr_ac_tau in component ikr (ms)" legend_algebraic[40] = "inf in component ikr (dimensionless)" legend_algebraic[256] = "inx in component ikr (dimensionless)" legend_algebraic[348] = "EKs in component nernst (mV)" legend_algebraic[259] = "G in component iks (mS_per_cm2)" legend_algebraic[351] = "IKs in component iks (uA_per_cm2)" legend_constants[403] = "IKs_arn in component iks_sig (uM)" legend_algebraic[349] = "IKs_np in component iks (uA_per_cm2)" legend_algebraic[350] = "IKs_pka in component iks (uA_per_cm2)" legend_constants[110] = "IKs_tot in component iks_sig (uM)" legend_algebraic[266] = "f_hat in component iks (dimensionless)" legend_algebraic[262] = "fp_iks in component iks_sig (dimensionless)" legend_constants[404] = "iks_f_hat_ratio in component iks (dimensionless)" legend_algebraic[264] = "iks_f_hat_val in component iks (dimensionless)" legend_states[69] = "O1 in component iks_np (dimensionless)" legend_states[70] = "O2 in component iks_np (dimensionless)" legend_states[71] = "O1 in component iks_pka (dimensionless)" legend_states[72] = "O2 in component iks_pka (dimensionless)" legend_states[73] = "C1 in component iks_np (dimensionless)" legend_states[74] = "C10 in component iks_np (dimensionless)" legend_states[75] = "C11 in component iks_np (dimensionless)" legend_states[76] = "C12 in component iks_np (dimensionless)" legend_states[77] = "C13 in component iks_np (dimensionless)" legend_states[78] = "C14 in component iks_np (dimensionless)" legend_states[79] = "C15 in component iks_np (dimensionless)" legend_states[80] = "C2 in component iks_np (dimensionless)" legend_states[81] = "C3 in component iks_np (dimensionless)" legend_states[82] = "C4 in component iks_np (dimensionless)" legend_states[83] = "C5 in component iks_np (dimensionless)" legend_states[84] = "C6 in component iks_np (dimensionless)" legend_states[85] = "C7 in component iks_np (dimensionless)" legend_states[86] = "C8 in component iks_np (dimensionless)" legend_states[87] = "C9 in component iks_np (dimensionless)" legend_algebraic[13] = "a in component iks_np (mS_per_uF)" legend_algebraic[41] = "b in component iks_np (mS_per_uF)" legend_algebraic[63] = "d in component iks_np (mS_per_uF)" legend_algebraic[75] = "e in component iks_np (mS_per_uF)" legend_algebraic[91] = "g in component iks_np (mS_per_uF)" legend_algebraic[92] = "o in component iks_np (mS_per_uF)" legend_algebraic[103] = "p in component iks_np (mS_per_uF)" legend_constants[111] = "t in component iks_np (mS_per_uF)" legend_states[88] = "C1 in component iks_pka (dimensionless)" legend_states[89] = "C10 in component iks_pka (dimensionless)" legend_states[90] = "C11 in component iks_pka (dimensionless)" legend_states[91] = "C12 in component iks_pka (dimensionless)" legend_states[92] = "C13 in component iks_pka (dimensionless)" legend_states[93] = "C14 in component iks_pka (dimensionless)" legend_states[94] = "C15 in component iks_pka (dimensionless)" legend_states[95] = "C2 in component iks_pka (dimensionless)" legend_states[96] = "C3 in component iks_pka (dimensionless)" legend_states[97] = "C4 in component iks_pka (dimensionless)" legend_states[98] = "C5 in component iks_pka (dimensionless)" legend_states[99] = "C6 in component iks_pka (dimensionless)" legend_states[100] = "C7 in component iks_pka (dimensionless)" legend_states[101] = "C8 in component iks_pka (dimensionless)" legend_states[102] = "C9 in component iks_pka (dimensionless)" legend_algebraic[14] = "a in component iks_pka (mS_per_uF)" legend_algebraic[42] = "b in component iks_pka (mS_per_uF)" legend_algebraic[64] = "d in component iks_pka (mS_per_uF)" legend_algebraic[76] = "e in component iks_pka (mS_per_uF)" legend_algebraic[93] = "g in component iks_pka (mS_per_uF)" legend_algebraic[94] = "o in component iks_pka (mS_per_uF)" legend_algebraic[104] = "p in component iks_pka (mS_per_uF)" legend_constants[112] = "t in component iks_pka (mS_per_uF)" legend_constants[405] = "IKs_arp in component iks_sig (uM)" legend_constants[356] = "IKsf in component iks_sig (uM)" legend_states[103] = "IKsp in component iks_sig (uM)" legend_constants[113] = "K in component iks_sig (uM)" legend_constants[114] = "Ka_iks in component iks_sig (uM)" legend_constants[115] = "Kp_iks in component iks_sig (uM)" legend_constants[116] = "L in component iks_sig (uM)" legend_constants[117] = "M in component iks_sig (uM)" legend_constants[399] = "PKA_eca in component pka (uM)" legend_constants[401] = "PKAf in component iks_sig (uM)" legend_constants[360] = "PP1f_eca in component iks_sig (uM)" legend_constants[118] = "Yotiao in component iks_sig (uM)" legend_constants[402] = "Yotiaof in component iks_sig (uM)" legend_constants[353] = "iks_sig_IKsf_sum in component iks_sig (dimensionless)" legend_algebraic[15] = "iks_sig_IKsp_dif in component iks_sig (uM)" legend_constants[400] = "iks_sig_PKAf_sum in component iks_sig (dimensionless)" legend_constants[358] = "iks_sig_PP1f_eca_sum in component iks_sig (dimensionless)" legend_constants[119] = "ka_iks in component iks_sig (hertz)" legend_constants[120] = "kp_iks in component iks_sig (hertz)" legend_algebraic[343] = "IKur in component ikur (uA_per_cm2)" legend_algebraic[341] = "IKur_np in component ikur (uA_per_cm2)" legend_algebraic[342] = "IKur_p in component ikur (uA_per_cm2)" legend_constants[121] = "Ka_ikur in component ikur (uM)" legend_constants[122] = "Kp_ikur in component ikur (uM)" legend_states[104] = "f_ikur in component ikur (dimensionless)" legend_algebraic[270] = "fhat in component ikur (dimensionless)" legend_constants[123] = "gbar_np in component ikur (mS_per_cm2)" legend_algebraic[268] = "ikur_fhat_val in component ikur (dimensionless)" legend_constants[124] = "ka_ikur in component ikur (hertz)" legend_constants[125] = "kp_ikur in component ikur (hertz)" legend_algebraic[359] = "INa in component ina (uA_per_cm2)" legend_algebraic[355] = "INa_both in component ina (uA_per_cm2)" legend_algebraic[356] = "INa_camk in component ina (uA_per_cm2)" legend_algebraic[357] = "INa_np in component ina (uA_per_cm2)" legend_algebraic[358] = "INa_pka in component ina (uA_per_cm2)" legend_constants[126] = "Ka_ina in component ina (uM)" legend_constants[127] = "Kp_ina in component ina (uM)" legend_algebraic[275] = "f_both in component ina (dimensionless)" legend_algebraic[276] = "f_camk_only in component ina (dimensionless)" legend_states[105] = "f_ina in component ina (dimensionless)" legend_algebraic[278] = "f_np in component ina (dimensionless)" legend_algebraic[274] = "f_pka in component ina (dimensionless)" legend_algebraic[277] = "f_pka_only in component ina (dimensionless)" legend_constants[241] = "gNaBar in component ina (mS_per_cm2)" legend_states[106] = "h in component ina_camk (dimensionless)" legend_states[107] = "j in component ina_camk (dimensionless)" legend_states[108] = "m in component ina_camk (dimensionless)" legend_algebraic[273] = "ina_f_pka_val in component ina (dimensionless)" legend_states[109] = "h in component ina_np (dimensionless)" legend_states[110] = "j in component ina_np (dimensionless)" legend_states[111] = "m in component ina_np (dimensionless)" legend_states[112] = "h in component ina_pka (dimensionless)" legend_states[113] = "j in component ina_pka (dimensionless)" legend_states[114] = "m in component ina_pka (dimensionless)" legend_constants[128] = "ka_ina in component ina (hertz)" legend_constants[129] = "kp_ina in component ina (hertz)" legend_constants[130] = "dVIn in component ina_camk (mV)" legend_algebraic[16] = "ina_camk_h_alpha in component ina_camk (mS_per_uF)" legend_algebraic[43] = "ina_camk_h_beta in component ina_camk (mS_per_uF)" legend_algebraic[17] = "ina_camk_j_alpha in component ina_camk (mS_per_uF)" legend_algebraic[44] = "ina_camk_j_beta in component ina_camk (mS_per_uF)" legend_algebraic[18] = "ina_camk_m_alpha in component ina_camk (mS_per_uF)" legend_algebraic[45] = "ina_camk_m_beta in component ina_camk (mS_per_uF)" legend_algebraic[19] = "ina_np_h_alpha in component ina_np (mS_per_uF)" legend_algebraic[46] = "ina_np_h_beta in component ina_np (mS_per_uF)" legend_algebraic[20] = "ina_np_j_alpha in component ina_np (mS_per_uF)" legend_algebraic[47] = "ina_np_j_beta in component ina_np (mS_per_uF)" legend_algebraic[21] = "ina_np_m_alpha in component ina_np (mS_per_uF)" legend_algebraic[48] = "ina_np_m_beta in component ina_np (mS_per_uF)" legend_constants[131] = "dVAc in component ina_pka (mV)" legend_constants[132] = "dVIn in component ina_pka (mV)" legend_algebraic[22] = "ina_pka_h_alpha in component ina_pka (mS_per_uF)" legend_algebraic[49] = "ina_pka_h_beta in component ina_pka (mS_per_uF)" legend_algebraic[23] = "ina_pka_j_alpha in component ina_pka (mS_per_uF)" legend_algebraic[50] = "ina_pka_j_beta in component ina_pka (mS_per_uF)" legend_algebraic[24] = "ina_pka_m_alpha in component ina_pka (mS_per_uF)" legend_algebraic[51] = "ina_pka_m_beta in component ina_pka (mS_per_uF)" legend_algebraic[281] = "INab in component inab (uA_per_cm2)" legend_constants[133] = "P in component inab (cm_per_s)" legend_algebraic[280] = "ePhi in component inab (dimensionless)" legend_algebraic[279] = "inab_INab_phi in component inab (dimensionless)" legend_constants[362] = "KmNai3 in component inaca (mM3)" legend_constants[364] = "KmNao3 in component inaca (mM3)" legend_constants[134] = "Km_Ca in component inaca (mM)" legend_constants[135] = "Km_Cai in component inaca (mM)" legend_constants[136] = "Km_Cao in component inaca (mM)" legend_constants[137] = "Km_Nai in component inaca (mM)" legend_constants[138] = "Km_Nao in component inaca (mM)" legend_algebraic[282] = "Na_i3 in component inaca (mM3)" legend_constants[366] = "Na_o3 in component inaca (mM3)" legend_algebraic[283] = "Na_ss3 in component inaca (mM3)" legend_constants[139] = "eta in component inaca (dimensionless)" legend_algebraic[284] = "exp1 in component inaca (dimensionless)" legend_algebraic[285] = "exp2 in component inaca (dimensionless)" legend_algebraic[286] = "inaca_INaCaSR_denom1 in component inaca (dimensionless)" legend_algebraic[287] = "inaca_INaCaSR_denom2 in component inaca (dimensionless)" legend_algebraic[288] = "inaca_INaCaSR_denom3 in component inaca (mol4_per_m12)" legend_algebraic[289] = "inaca_INaCaSR_denom4 in component inaca (mol4_per_m12)" legend_algebraic[290] = "inaca_INaCaSR_num in component inaca (A_mol4_per_m14_times_1e_minus_2)" legend_algebraic[292] = "inaca_INaCa_denom1 in component inaca (dimensionless)" legend_algebraic[293] = "inaca_INaCa_denom2 in component inaca (dimensionless)" legend_algebraic[294] = "inaca_INaCa_denom3 in component inaca (mol4_per_m12)" legend_algebraic[295] = "inaca_INaCa_denom4 in component inaca (mol4_per_m12)" legend_algebraic[296] = "inaca_INaCa_num in component inaca (A_mol4_per_m14_times_1e_minus_2)" legend_constants[140] = "kSat in component inaca (dimensionless)" legend_constants[141] = "vMax in component inaca (uA_per_cm2)" legend_algebraic[303] = "INaK in component inak (uA_per_cm2)" legend_algebraic[301] = "INaK_np in component inak (uA_per_cm2)" legend_algebraic[302] = "INaK_p in component inak (uA_per_cm2)" legend_constants[142] = "Ka_inak in component inak (uM)" legend_constants[143] = "Kp_inak in component inak (uM)" legend_states[115] = "f_inak in component inak (dimensionless)" legend_algebraic[299] = "fhat in component inak (dimensionless)" legend_constants[144] = "ibar in component inak (uA_per_cm2)" legend_algebraic[298] = "inak_fhat_val in component inak (dimensionless)" legend_constants[145] = "ka_inak in component inak (hertz)" legend_constants[146] = "km_ko in component inak (mM)" legend_constants[147] = "km_np in component inak (mM)" legend_constants[148] = "km_p in component inak (mM)" legend_constants[149] = "kp_inak in component inak (hertz)" legend_algebraic[300] = "phi in component inak (uA_per_cm2)" legend_constants[242] = "pk in component inak (dimensionless)" legend_algebraic[363] = "INaL in component inal (uA_per_cm2)" legend_algebraic[361] = "INaL_camk in component inal (uA_per_cm2)" legend_algebraic[362] = "INaL_np in component inal (uA_per_cm2)" legend_algebraic[360] = "conductance in component inal (mV)" legend_states[116] = "h in component inal (dimensionless)" legend_algebraic[25] = "h_inf in component inal (dimensionless)" legend_algebraic[26] = "inal_m_alpha in component inal (mS_per_uF)" legend_algebraic[52] = "inal_m_beta in component inal (mS_per_uF)" legend_states[117] = "m in component inal (dimensionless)" legend_constants[150] = "tau_h in component inal (ms)" legend_constants[151] = "IpCa_bar in component ipca (uA_per_cm2)" legend_constants[152] = "Km_pCa in component ipca (mM)" legend_algebraic[315] = "Ileak_ryr in component irel (mM_per_ms)" legend_algebraic[307] = "Ileak_ryr_np in component irel (mM_per_ms)" legend_algebraic[309] = "Ileak_ryr_p in component irel (mM_per_ms)" legend_states[118] = "Irel_np in component irel (mM_per_ms)" legend_states[119] = "Irel_p in component irel (mM_per_ms)" legend_constants[153] = "Km_ryr_leak_np in component irel (mM)" legend_constants[154] = "Km_ryr_leak_p in component irel (mM)" legend_algebraic[255] = "alpha_np in component irel (mM_per_ms)" legend_algebraic[261] = "alpha_p in component irel (mM_per_ms)" legend_constants[243] = "beta_0 in component irel (ms)" legend_algebraic[252] = "beta_np in component irel (ms)" legend_algebraic[258] = "beta_p in component irel (ms)" legend_algebraic[313] = "fhat in component irel (dimensionless)" legend_constants[398] = "irel_fhat_ratio in component irel (dimensionless)" legend_algebraic[311] = "irel_fhat_val in component irel (dimensionless)" legend_algebraic[271] = "irel_inf_np in component irel (mM_per_ms)" legend_algebraic[272] = "irel_inf_p in component irel (mM_per_ms)" legend_algebraic[265] = "irel_tau_np in component irel (ms)" legend_algebraic[267] = "irel_tau_p in component irel (ms)" legend_constants[155] = "k_ryr_leak_np in component irel (mS_per_uF)" legend_constants[156] = "k_ryr_leak_p in component irel (mS_per_uF)" legend_algebraic[263] = "x in component irel (dimensionless)" legend_algebraic[269] = "y in component irel (dimensionless)" legend_constants[157] = "Gbar in component ito (mS_per_cm2)" legend_algebraic[347] = "ITo in component ito (uA_per_cm2)" legend_algebraic[345] = "ITo_camk in component ito (uA_per_cm2)" legend_algebraic[346] = "ITo_np in component ito (uA_per_cm2)" legend_algebraic[325] = "R in component ito (dimensionless)" legend_algebraic[27] = "a_inf in component ito (dimensionless)" legend_states[120] = "a_np in component ito (dimensionless)" legend_algebraic[77] = "a_tau in component ito (ms)" legend_algebraic[28] = "alph_if in component ito (dimensionless)" legend_algebraic[54] = "alph_is in component ito (dimensionless)" legend_algebraic[66] = "beta_i in component ito (mS_per_uF)" legend_states[121] = "if_camk in component ito (dimensionless)" legend_states[122] = "if_np in component ito (dimensionless)" legend_states[123] = "is_camk in component ito (dimensionless)" legend_states[124] = "is_np in component ito (dimensionless)" legend_algebraic[53] = "ito_a_np_alpha in component ito (dimensionless)" legend_algebraic[65] = "ito_a_np_beta in component ito (dimensionless)" legend_algebraic[78] = "ito_if_camk_alpha in component ito (mS_per_uF)" legend_algebraic[79] = "ito_if_np_alpha in component ito (mS_per_uF)" legend_algebraic[80] = "ito_is_camk_alpha in component ito (mS_per_uF)" legend_algebraic[81] = "ito_is_np_alpha in component ito (mS_per_uF)" legend_algebraic[344] = "x in component ito (uA_per_cm2)" legend_algebraic[308] = "Imax in component iup (mM_per_ms)" legend_constants[158] = "Ka_plb in component iup (uM)" legend_constants[376] = "Km_both in component iup (mM)" legend_constants[369] = "Km_camk in component iup (mM)" legend_constants[159] = "Km_np in component iup (mM)" legend_constants[374] = "Km_pka in component iup (mM)" legend_algebraic[322] = "Km_up in component iup (mM)" legend_constants[160] = "Kp_plb in component iup (uM)" legend_algebraic[306] = "f_SERCA2a in component iup (dimensionless)" legend_algebraic[314] = "f_both in component iup (dimensionless)" legend_algebraic[316] = "f_camk_only in component iup (dimensionless)" legend_algebraic[321] = "f_np in component iup (dimensionless)" legend_algebraic[312] = "f_pka in component iup (dimensionless)" legend_algebraic[318] = "f_pka_only in component iup (dimensionless)" legend_states[125] = "f_plb in component iup (dimensionless)" legend_algebraic[310] = "iup_f_pka_val in component iup (dimensionless)" legend_constants[161] = "iupmax in component iup (mM_per_ms)" legend_constants[379] = "iupmaxCAMK in component iup (mM_per_ms)" legend_constants[162] = "ka_plb in component iup (hertz)" legend_constants[163] = "kp_plb in component iup (hertz)" legend_algebraic[324] = "leak in component iup (mM_per_ms)" legend_constants[164] = "nsrmax in component iup (mM)" legend_algebraic[326] = "uptake in component iup (mM_per_ms)" legend_algebraic[364] = "IK_tot in component potassium (uA_per_cm2)" legend_algebraic[366] = "INa_tot in component sodium (uA_per_cm2)" legend_algebraic[367] = "i_ion in component membrane (uA_per_cm2)" legend_algebraic[368] = "i_stim in component stimulus (uA_per_cm2)" legend_constants[165] = "PNaK in component nernst (dimensionless)" legend_constants[238] = "RTF in component phys (mV)" legend_states[126] = "K in component potassium (mM)" legend_constants[166] = "KPDEp in component pde (uM)" legend_constants[167] = "KmIbmxPde2 in component pde (dimensionless)" legend_constants[168] = "KmIbmxPde3 in component pde (dimensionless)" legend_constants[169] = "KmIbmxPde4 in component pde (dimensionless)" legend_constants[170] = "KmPDE2 in component pde (uM)" legend_constants[171] = "KmPDE3 in component pde (uM)" legend_constants[172] = "KmPDE4 in component pde (uM)" legend_constants[347] = "PDE2_cav in component pde (uM)" legend_constants[348] = "PDE2_cyt in component pde (uM)" legend_constants[349] = "PDE2_eca in component pde (uM)" legend_constants[173] = "PDE2_tot in component pde (uM)" legend_states[127] = "PDE3_P_cav in component pde (uM)" legend_states[128] = "PDE3_P_cyt in component pde (uM)" legend_constants[377] = "PDE3_cav in component pde (uM)" legend_constants[378] = "PDE3_cyt in component pde (uM)" legend_constants[363] = "PDE3_tot in component pde (uM)" legend_states[129] = "PDE4_P_cav in component pde (uM)" legend_states[130] = "PDE4_P_cyt in component pde (uM)" legend_states[131] = "PDE4_P_eca in component pde (uM)" legend_constants[370] = "PDE4_cav in component pde (uM)" legend_constants[371] = "PDE4_cyt in component pde (uM)" legend_constants[372] = "PDE4_eca in component pde (uM)" legend_constants[365] = "PDE4_tot in component pde (uM)" legend_constants[174] = "delta_k_pde34 in component pde (dimensionless)" legend_constants[175] = "f_pde2_cav in component pde (dimensionless)" legend_constants[333] = "f_pde2_cyt in component pde (dimensionless)" legend_constants[176] = "f_pde2_eca in component pde (dimensionless)" legend_constants[334] = "f_pde2_part in component pde (dimensionless)" legend_constants[367] = "f_pde3_cav in component pde (dimensionless)" legend_constants[373] = "f_pde3_cyt in component pde (dimensionless)" legend_constants[177] = "f_pde4_cav in component pde (dimensionless)" legend_constants[339] = "f_pde4_cyt in component pde (dimensionless)" legend_constants[343] = "f_pde4_eca in component pde (dimensionless)" legend_constants[178] = "f_pde4_part in component pde (dimensionless)" legend_constants[179] = "f_pde_part in component pde (dimensionless)" legend_constants[180] = "ff_pde3_cyt in component pde (dimensionless)" legend_constants[181] = "h_ibmx_pde2 in component pde (dimensionless)" legend_constants[182] = "h_ibmx_pde3 in component pde (dimensionless)" legend_constants[183] = "h_ibmx_pde4 in component pde (dimensionless)" legend_constants[184] = "ibmx in component pde (uM)" legend_constants[342] = "ibmx2 in component pde (uM)" legend_constants[375] = "ibmx3 in component pde (uM)" legend_constants[368] = "ibmx4 in component pde (uM)" legend_constants[338] = "ibmx_h2 in component pde (dimensionless)" legend_constants[350] = "ibmx_h3 in component pde (dimensionless)" legend_constants[354] = "ibmx_h4 in component pde (dimensionless)" legend_constants[185] = "kPDE2 in component pde (hertz)" legend_constants[186] = "kPDE3 in component pde (hertz)" legend_constants[187] = "kPDE4 in component pde (hertz)" legend_constants[244] = "kbPDEp in component pde (hertz)" legend_constants[188] = "kfPDEp in component pde (per_mM_per_ms)" legend_constants[359] = "pde_PDE3_tot_alpha in component pde (dimensionless)" legend_constants[361] = "pde_PDE3_tot_beta in component pde (dimensionless)" legend_constants[189] = "r_pde34_frac in component pde (dimensionless)" legend_constants[357] = "r_pde3_cyt in component pde (dimensionless)" legend_constants[190] = "R in component phys (mJ_per_mol_per_K)" legend_constants[191] = "T in component phys (kelvin)" legend_constants[245] = "K_pki in component pka (uM)" legend_constants[335] = "PKA_cyt in component pka (uM)" legend_constants[192] = "PKA_tot in component pka (uM)" legend_constants[344] = "PKI_cav in component pka (uM)" legend_constants[355] = "PKI_cyt in component pka (uM)" legend_constants[351] = "PKI_eca in component pka (uM)" legend_constants[336] = "PKI_tot in component pka (uM)" legend_constants[264] = "b_pki in component pka (hertz)" legend_constants[193] = "f_cav in component pka (dimensionless)" legend_constants[327] = "f_cyt in component pka (dimensionless)" legend_constants[194] = "f_eca in component pka (dimensionless)" legend_constants[195] = "f_pki in component pka (per_mM_per_ms)" legend_constants[340] = "f_pki_cav in component pka (dimensionless)" legend_constants[352] = "f_pki_cyt in component pka (dimensionless)" legend_constants[345] = "f_pki_eca in component pka (dimensionless)" legend_states[132] = "A2R in component pka_cav (uM)" legend_states[133] = "A2RC in component pka_cav (uM)" legend_states[134] = "ARC in component pka_cav (uM)" legend_constants[196] = "K1 in component pka_cav (uM)" legend_constants[197] = "K2 in component pka_cav (uM)" legend_constants[198] = "K3 in component pka_cav (m3_per_mol_times_1e3)" legend_states[135] = "PKIC in component pka_cav (uM)" legend_algebraic[140] = "RCf in component pka_cav (uM)" legend_constants[246] = "b1 in component pka_cav (hertz)" legend_constants[247] = "b2 in component pka_cav (hertz)" legend_constants[248] = "b3 in component pka_cav (per_mM_per_ms)" legend_constants[199] = "f1 in component pka_cav (per_mM_per_ms)" legend_constants[200] = "f2 in component pka_cav (per_mM_per_ms)" legend_constants[201] = "f3 in component pka_cav (hertz)" legend_states[136] = "A2R in component pka_cyt (uM)" legend_states[137] = "A2RC in component pka_cyt (uM)" legend_states[138] = "ARC in component pka_cyt (uM)" legend_constants[202] = "K1 in component pka_cyt (uM)" legend_constants[203] = "K2 in component pka_cyt (uM)" legend_constants[204] = "K3 in component pka_cyt (m3_per_mol_times_1e3)" legend_states[139] = "PKIC in component pka_cyt (uM)" legend_algebraic[150] = "RCf in component pka_cyt (uM)" legend_constants[265] = "b1 in component pka_cyt (hertz)" legend_constants[266] = "b2 in component pka_cyt (hertz)" legend_constants[267] = "b3 in component pka_cyt (per_mM_per_ms)" legend_constants[249] = "f1 in component pka_cyt (per_mM_per_ms)" legend_constants[250] = "f2 in component pka_cyt (per_mM_per_ms)" legend_constants[251] = "f3 in component pka_cyt (hertz)" legend_states[140] = "A2R in component pka_eca (uM)" legend_states[141] = "A2RC in component pka_eca (uM)" legend_states[142] = "ARC in component pka_eca (uM)" legend_constants[252] = "K1 in component pka_eca (uM)" legend_constants[253] = "K2 in component pka_eca (uM)" legend_constants[254] = "K3 in component pka_eca (m3_per_mol_times_1e3)" legend_states[143] = "PKIC in component pka_eca (uM)" legend_algebraic[130] = "RCf in component pka_eca (uM)" legend_constants[276] = "b1 in component pka_eca (hertz)" legend_constants[277] = "b2 in component pka_eca (hertz)" legend_constants[278] = "b3 in component pka_eca (per_mM_per_ms)" legend_constants[268] = "f1 in component pka_eca (per_mM_per_ms)" legend_constants[269] = "f2 in component pka_eca (per_mM_per_ms)" legend_constants[270] = "f3 in component pka_eca (hertz)" legend_constants[406] = "potassium_K_r1 in component potassium (mol_per_m_per_s_per_A_times_1e5)" legend_constants[205] = "K in component pp1 (uM)" legend_constants[206] = "Kdp in component pp1 (uM)" legend_constants[207] = "Kp in component pp1 (uM)" legend_constants[208] = "PP1_cyt in component pp1 (uM)" legend_algebraic[29] = "di in component pp1 (uM)" legend_constants[209] = "f in component pp1 (dimensionless)" legend_states[144] = "inhib1_p in component pp1 (uM)" legend_constants[255] = "inhib1_tot in component pp1 (uM)" legend_constants[210] = "kdp in component pp1 (hertz)" legend_constants[211] = "kp in component pp1 (hertz)" legend_algebraic[198] = "pp1_PP1f_cyt_sum in component pp1 (uM)" legend_algebraic[365] = "INa_cyt in component sodium (uA_per_cm2)" legend_constants[407] = "sodium_Na_r1 in component sodium (mol_per_m_per_s_per_A_times_1e5)" legend_constants[408] = "sodium_Na_r2 in component sodium (dimensionless)" legend_constants[409] = "sodium_Na_sr_r1 in component sodium (mol_per_m_per_s_per_A_times_1e5)" legend_constants[256] = "amplitude in component stimulus (uA_per_cm2)" legend_constants[212] = "duration in component stimulus (ms)" legend_constants[213] = "offset in component stimulus (ms)" legend_constants[214] = "period in component stimulus (ms)" legend_rates[4] = "d/dt ICaLp in component akap_sig (uM)" legend_rates[5] = "d/dt RyRp in component akap_sig (uM)" legend_rates[7] = "d/dt Gi_aGDP in component beta_cav (uM)" legend_rates[8] = "d/dt Gi_aGTP in component beta_cav (uM)" legend_rates[0] = "d/dt Gi_bg in component beta_cav (uM)" legend_rates[9] = "d/dt Gs_aGDP in component beta_cav (uM)" legend_rates[1] = "d/dt Gs_aGTP in component beta_cav (uM)" legend_rates[10] = "d/dt Gs_bg in component beta_cav (uM)" legend_rates[11] = "d/dt Rb1_grk_tot in component beta_cav (uM)" legend_rates[12] = "d/dt Rb1_pka_tot in component beta_cav (uM)" legend_rates[13] = "d/dt Rb2_grk_tot in component beta_cav (uM)" legend_rates[14] = "d/dt Rb2_pka_tot in component beta_cav (uM)" legend_rates[15] = "d/dt Gs_aGDP in component beta_cyt (uM)" legend_rates[2] = "d/dt Gs_aGTP in component beta_cyt (uM)" legend_rates[16] = "d/dt Gs_bg in component beta_cyt (uM)" legend_rates[17] = "d/dt Rb1_grk_tot in component beta_cyt (uM)" legend_rates[18] = "d/dt Rb1_pka_tot in component beta_cyt (uM)" legend_rates[20] = "d/dt Gi_aGDP in component beta_eca (uM)" legend_rates[21] = "d/dt Gi_aGTP in component beta_eca (uM)" legend_rates[22] = "d/dt Gi_bg in component beta_eca (uM)" legend_rates[23] = "d/dt Gs_aGDP in component beta_eca (uM)" legend_rates[3] = "d/dt Gs_aGTP in component beta_eca (uM)" legend_rates[24] = "d/dt Gs_bg in component beta_eca (uM)" legend_rates[25] = "d/dt Rb1_grk_tot in component beta_eca (uM)" legend_rates[26] = "d/dt Rb1_pka_tot in component beta_eca (uM)" legend_rates[27] = "d/dt Rb2_grk_tot in component beta_eca (uM)" legend_rates[28] = "d/dt Rb2_pka_tot in component beta_eca (uM)" legend_rates[30] = "d/dt Ca_nsr in component calcium (mM)" legend_rates[31] = "d/dt f_tni in component calcium (dimensionless)" legend_rates[32] = "d/dt uCa in component calcium (mM)" legend_rates[33] = "d/dt uCa_CaL in component calcium (mM)" legend_rates[34] = "d/dt uCa_jsr in component calcium (mM)" legend_rates[35] = "d/dt uCa_sr in component calcium (mM)" legend_rates[36] = "d/dt f_ical in component camk (dimensionless)" legend_rates[37] = "d/dt f_ik1 in component camk (dimensionless)" legend_rates[38] = "d/dt f_ina in component camk (dimensionless)" legend_rates[39] = "d/dt f_ito in component camk (dimensionless)" legend_rates[40] = "d/dt f_plb in component camk (dimensionless)" legend_rates[41] = "d/dt f_ryr in component camk (dimensionless)" legend_rates[42] = "d/dt trap in component camk (dimensionless)" legend_rates[43] = "d/dt cAMP_cav in component camp (uM)" legend_rates[44] = "d/dt cAMP_cyt in component camp (uM)" legend_rates[45] = "d/dt cAMP_eca in component camp (uM)" legend_rates[46] = "d/dt Cl in component chloride (mM)" legend_rates[47] = "d/dt Cl_sr in component chloride (mM)" legend_rates[51] = "d/dt C in component ical_camk (dimensionless)" legend_rates[52] = "d/dt CI in component ical_camk (dimensionless)" legend_rates[53] = "d/dt CIs in component ical_camk (dimensionless)" legend_rates[54] = "d/dt Cs in component ical_camk (dimensionless)" legend_rates[55] = "d/dt O in component ical_camk (dimensionless)" legend_rates[56] = "d/dt OI in component ical_camk (dimensionless)" legend_rates[57] = "d/dt OIs in component ical_camk (dimensionless)" legend_rates[58] = "d/dt Os in component ical_camk (dimensionless)" legend_rates[59] = "d/dt C in component ical_np (dimensionless)" legend_rates[60] = "d/dt CI in component ical_np (dimensionless)" legend_rates[61] = "d/dt CIs in component ical_np (dimensionless)" legend_rates[62] = "d/dt Cs in component ical_np (dimensionless)" legend_rates[63] = "d/dt O in component ical_np (dimensionless)" legend_rates[64] = "d/dt OI in component ical_np (dimensionless)" legend_rates[65] = "d/dt OIs in component ical_np (dimensionless)" legend_rates[66] = "d/dt Os in component ical_np (dimensionless)" legend_rates[67] = "d/dt i2 in component iclca (dimensionless)" legend_rates[68] = "d/dt ac in component ikr (dimensionless)" legend_rates[73] = "d/dt C1 in component iks_np (dimensionless)" legend_rates[74] = "d/dt C10 in component iks_np (dimensionless)" legend_rates[75] = "d/dt C11 in component iks_np (dimensionless)" legend_rates[76] = "d/dt C12 in component iks_np (dimensionless)" legend_rates[77] = "d/dt C13 in component iks_np (dimensionless)" legend_rates[78] = "d/dt C14 in component iks_np (dimensionless)" legend_rates[79] = "d/dt C15 in component iks_np (dimensionless)" legend_rates[80] = "d/dt C2 in component iks_np (dimensionless)" legend_rates[81] = "d/dt C3 in component iks_np (dimensionless)" legend_rates[82] = "d/dt C4 in component iks_np (dimensionless)" legend_rates[83] = "d/dt C5 in component iks_np (dimensionless)" legend_rates[84] = "d/dt C6 in component iks_np (dimensionless)" legend_rates[85] = "d/dt C7 in component iks_np (dimensionless)" legend_rates[86] = "d/dt C8 in component iks_np (dimensionless)" legend_rates[87] = "d/dt C9 in component iks_np (dimensionless)" legend_rates[69] = "d/dt O1 in component iks_np (dimensionless)" legend_rates[70] = "d/dt O2 in component iks_np (dimensionless)" legend_rates[88] = "d/dt C1 in component iks_pka (dimensionless)" legend_rates[89] = "d/dt C10 in component iks_pka (dimensionless)" legend_rates[90] = "d/dt C11 in component iks_pka (dimensionless)" legend_rates[91] = "d/dt C12 in component iks_pka (dimensionless)" legend_rates[92] = "d/dt C13 in component iks_pka (dimensionless)" legend_rates[93] = "d/dt C14 in component iks_pka (dimensionless)" legend_rates[94] = "d/dt C15 in component iks_pka (dimensionless)" legend_rates[95] = "d/dt C2 in component iks_pka (dimensionless)" legend_rates[96] = "d/dt C3 in component iks_pka (dimensionless)" legend_rates[97] = "d/dt C4 in component iks_pka (dimensionless)" legend_rates[98] = "d/dt C5 in component iks_pka (dimensionless)" legend_rates[99] = "d/dt C6 in component iks_pka (dimensionless)" legend_rates[100] = "d/dt C7 in component iks_pka (dimensionless)" legend_rates[101] = "d/dt C8 in component iks_pka (dimensionless)" legend_rates[102] = "d/dt C9 in component iks_pka (dimensionless)" legend_rates[71] = "d/dt O1 in component iks_pka (dimensionless)" legend_rates[72] = "d/dt O2 in component iks_pka (dimensionless)" legend_rates[103] = "d/dt IKsp in component iks_sig (uM)" legend_rates[104] = "d/dt f_ikur in component ikur (dimensionless)" legend_rates[105] = "d/dt f_ina in component ina (dimensionless)" legend_rates[106] = "d/dt h in component ina_camk (dimensionless)" legend_rates[107] = "d/dt j in component ina_camk (dimensionless)" legend_rates[108] = "d/dt m in component ina_camk (dimensionless)" legend_rates[109] = "d/dt h in component ina_np (dimensionless)" legend_rates[110] = "d/dt j in component ina_np (dimensionless)" legend_rates[111] = "d/dt m in component ina_np (dimensionless)" legend_rates[112] = "d/dt h in component ina_pka (dimensionless)" legend_rates[113] = "d/dt j in component ina_pka (dimensionless)" legend_rates[114] = "d/dt m in component ina_pka (dimensionless)" legend_rates[115] = "d/dt f_inak in component inak (dimensionless)" legend_rates[116] = "d/dt h in component inal (dimensionless)" legend_rates[117] = "d/dt m in component inal (dimensionless)" legend_rates[118] = "d/dt Irel_np in component irel (mM_per_ms)" legend_rates[119] = "d/dt Irel_p in component irel (mM_per_ms)" legend_rates[120] = "d/dt a_np in component ito (dimensionless)" legend_rates[121] = "d/dt if_camk in component ito (dimensionless)" legend_rates[122] = "d/dt if_np in component ito (dimensionless)" legend_rates[123] = "d/dt is_camk in component ito (dimensionless)" legend_rates[124] = "d/dt is_np in component ito (dimensionless)" legend_rates[125] = "d/dt f_plb in component iup (dimensionless)" legend_rates[50] = "d/dt V in component membrane (mV)" legend_rates[127] = "d/dt PDE3_P_cav in component pde (uM)" legend_rates[128] = "d/dt PDE3_P_cyt in component pde (uM)" legend_rates[129] = "d/dt PDE4_P_cav in component pde (uM)" legend_rates[130] = "d/dt PDE4_P_cyt in component pde (uM)" legend_rates[131] = "d/dt PDE4_P_eca in component pde (uM)" legend_rates[132] = "d/dt A2R in component pka_cav (uM)" legend_rates[133] = "d/dt A2RC in component pka_cav (uM)" legend_rates[134] = "d/dt ARC in component pka_cav (uM)" legend_rates[6] = "d/dt C in component pka_cav (uM)" legend_rates[135] = "d/dt PKIC in component pka_cav (uM)" legend_rates[136] = "d/dt A2R in component pka_cyt (uM)" legend_rates[137] = "d/dt A2RC in component pka_cyt (uM)" legend_rates[138] = "d/dt ARC in component pka_cyt (uM)" legend_rates[19] = "d/dt C in component pka_cyt (uM)" legend_rates[139] = "d/dt PKIC in component pka_cyt (uM)" legend_rates[140] = "d/dt A2R in component pka_eca (uM)" legend_rates[141] = "d/dt A2RC in component pka_eca (uM)" legend_rates[142] = "d/dt ARC in component pka_eca (uM)" legend_rates[29] = "d/dt C in component pka_eca (uM)" legend_rates[143] = "d/dt PKIC in component pka_eca (uM)" legend_rates[126] = "d/dt K in component potassium (mM)" legend_rates[144] = "d/dt inhib1_p in component pp1 (uM)" legend_rates[48] = "d/dt Na in component sodium (mM)" legend_rates[49] = "d/dt Na_sr in component sodium (mM)" return (legend_states, legend_algebraic, legend_voi, legend_constants) def initConsts(): constants = [0.0] * sizeConstants; states = [0.0] * sizeStates; constants[0] = 5000.0 constants[1] = 315.0 constants[2] = 0.0465 constants[3] = 0.031544 constants[4] = 0.0852 constants[5] = 0.4824 constants[6] = 3.3757 constants[7] = 41.32 constants[8] = 0.03135 constants[9] = 0.037696 states[0] = 2.09911481235842013e-3 states[1] = 6.85041638458664965e-3 states[2] = 7.31420577213055985e-4 states[3] = 1.84627603007976003e-2 constants[10] = 0.16479 constants[11] = 0.087459 constants[12] = 1.0043 constants[13] = 1.3574 constants[14] = 0.6623 constants[15] = 0.8569 constants[16] = 0.025 constants[17] = 0.025 states[4] = 6.73713947839316954e-4 constants[18] = 1.27019999999999993e-6 constants[19] = 6.62979999999999944e-5 constants[20] = 0.01 constants[21] = 0.0063064 constants[22] = 0.043003 constants[23] = 0.01 constants[24] = 0.0001 constants[25] = 0.0001 constants[26] = 0.01 constants[27] = 0.01 constants[28] = 0.25 constants[29] = 0.125 constants[30] = 0.125 states[5] = 4.10693810508170991e-3 constants[31] = 5.10090000000000044e-4 constants[32] = 0.0025548 constants[33] = 0.0006903 constants[34] = 0.0038257 states[6] = 3.26565916584702978e-02 constants[35] = 0.5 constants[36] = 0.85 constants[37] = 0.0011071 constants[38] = 0.5664 constants[39] = 0.081161 constants[40] = 0.48744 constants[41] = 0.85 constants[42] = 4.0 constants[43] = 4.9054 constants[44] = 0.05 constants[45] = 0.25945 constants[46] = 2.449 constants[47] = 0.062 constants[48] = 0.567 constants[49] = 1.6655 constants[50] = 1.8463 constants[51] = 0.1 constants[52] = 0.012 constants[53] = 1.053 constants[54] = 1.053 constants[55] = 0.8 constants[56] = 1210.0 constants[57] = 0.35 constants[58] = 1.0 states[7] = 5.02792845976641014e-4 states[8] = 1.59632196638178009e-3 states[9] = 6.07316088556675987e-4 states[10] = 7.45773247314215036e-3 states[11] = 2.49592854373432000e-10 states[12] = 1.49041813757830998e-2 states[13] = 8.91799633266019011e-10 states[14] = 2.75455839709412009e-2 constants[59] = 0.0 constants[60] = 1.0 constants[61] = 1.0 states[15] = 4.19991861054322011e-4 states[16] = 1.15141243826746994e-3 states[17] = 7.07824478944670999e-11 states[18] = 9.44463350378085993e-3 states[19] = 3.62113356111495976e-01 constants[62] = 1.0 states[20] = 3.41341142614041016e-4 states[21] = 3.64315164237569004e-4 states[22] = 7.05656306851923029e-4 states[23] = 6.39038440072506948e-4 states[24] = 1.91017987408719017e-2 states[25] = 1.18055788874765002e-9 states[26] = 2.03016833596287999e-1 states[27] = 1.13428924662652000e-10 states[28] = 1.10248953370551007e-2 constants[63] = 1.0 states[29] = 5.67249910261072965e-01 states[30] = 1.191 constants[64] = 2.71430000000000008e-5 constants[65] = 0.26714 constants[66] = 1.0 constants[67] = 1.124 constants[68] = 0.0087 constants[69] = 0.047 constants[70] = 0.00087 constants[71] = 0.05 constants[72] = 10.0 constants[73] = 0.8 states[31] = 6.73518785672381992e-1 constants[74] = 0.10408 constants[75] = 0.00238 constants[76] = 0.052633 constants[77] = 0.0005 constants[78] = 0.07 states[32] = 0.013394 states[33] = 0.023413 states[34] = 6.8659 states[35] = 0.023413 constants[79] = 0.05 constants[80] = 0.25 constants[81] = 0.0015 constants[82] = 0.1 constants[83] = 0.05 constants[84] = 0.00068 states[36] = 0.0 states[37] = 0.0 states[38] = 0.0 states[39] = 0.0 states[40] = 0.0 states[41] = 0.0 constants[85] = 100000.0 constants[86] = 10000.0 states[42] = 0.0017546 states[43] = 3.47102959606005013e-1 states[44] = 4.74081735738210996e-1 states[45] = 9.62359241535767040e+00 constants[87] = 96487.0 constants[88] = 0.01 constants[89] = 3.14159265358979312e+00 constants[90] = 0.0011 states[46] = 20.273 states[47] = 20.273 constants[91] = 1.77e-5 constants[92] = 87.8251 constants[93] = 2.46108000000000016e-5 states[48] = 6.8909 states[49] = 6.8909 constants[94] = 0.2 constants[95] = 0.02 constants[96] = 75.0 constants[97] = 1.8 constants[98] = 100.0 constants[99] = 5.4 constants[100] = 140.0 states[50] = -87.491 constants[101] = 1.995e-07 states[51] = 1.0 states[52] = 0.0 states[53] = 0.0 states[54] = 0.0 states[55] = 0.0 states[56] = 0.0 states[57] = 0.0 states[58] = 0.0 constants[102] = 1.0 constants[103] = 1e-06 states[59] = 1.0 states[60] = 0.0 states[61] = 0.0 states[62] = 0.0 states[63] = 0.0 states[64] = 0.0 states[65] = 0.0 states[66] = 0.0 constants[104] = 1.0 constants[105] = 1e-06 constants[106] = 0.000225 constants[107] = 9e-07 states[67] = 0.99604 constants[108] = 0.4 constants[109] = 8.0 states[68] = 1.23059999999999995e-08 constants[110] = 0.025 states[69] = 9.37220000000000071e-16 states[70] = 1.65950000000000014e-17 states[71] = 1.12010000000000000e-15 states[72] = 1.66129999999999997e-18 states[73] = 0.91141 states[74] = 5.36959999999999959e-07 states[75] = 2.48610000000000001e-08 states[76] = 2.87759999999999997e-10 states[77] = 1.12169999999999999e-10 states[78] = 2.59670000000000003e-12 states[79] = 8.78740000000000076e-15 states[80] = 0.084395 states[81] = 0.0029306 states[82] = 4.52285000000000024e-05 states[83] = 2.61750000000000024e-07 states[84] = 0.0011424 states[85] = 7.93370000000000029e-05 states[86] = 1.83659999999999992e-06 states[87] = 1.41719999999999995e-08 constants[111] = 0.0027304 states[88] = 0.95624 states[89] = 3.19559999999999984e-07 states[90] = 7.039e-09 states[91] = 3.87629999999999971e-11 states[92] = 5.02769999999999999e-11 states[93] = 5.53740000000000021e-13 states[94] = 2.96639999999999984e-15 states[95] = 0.042127 states[96] = 6.95969999999999965e-04 states[97] = 5.11010000000000042e-06 states[98] = 1.407e-08 states[99] = 9.02690000000000046e-04 states[100] = 2.98259999999999987e-05 states[101] = 3.285e-07 states[102] = 1.206e-09 constants[112] = 0.0046171 states[103] = 7.65988420110534033e-04 constants[113] = 0.01 constants[114] = 9.97940000000000003e-05 constants[115] = 1.11470000000000002e-04 constants[116] = 0.0001 constants[117] = 0.01 constants[118] = 0.025 constants[119] = 0.16305 constants[120] = 1.0542 constants[121] = 0.27623 constants[122] = 0.002331 states[104] = 5.89379755147717982e-02 constants[123] = 0.00384 constants[124] = 0.069537 constants[125] = 0.317 constants[126] = 0.10988 constants[127] = 7.8605 states[105] = 2.39479458960527997e-01 states[106] = 0.83805 states[107] = 0.99281 states[108] = 6.81269999999999988e-04 states[109] = 0.0068172 states[110] = 0.99709 states[111] = 0.90163 states[112] = 0.001236 states[113] = 0.99123 states[114] = 0.79472 constants[128] = 0.01368 constants[129] = 0.052811 constants[130] = 3.25 constants[131] = 3.7 constants[132] = 4.9 constants[133] = 3.2e-09 constants[134] = 0.000125 constants[135] = 0.0036 constants[136] = 1.3 constants[137] = 12.3 constants[138] = 87.5 constants[139] = 0.27 constants[140] = 0.32 constants[141] = 4.5 constants[142] = 0.0011001 constants[143] = 5.7392 states[115] = 1.26345311579565994e-01 constants[144] = 1.4 constants[145] = 0.015265 constants[146] = 1.5 constants[147] = 2.6 constants[148] = 1.846 constants[149] = 0.092455 states[116] = 0.36003 states[117] = 0.0007053 constants[150] = 600.0 constants[151] = 0.0575 constants[152] = 0.0005 states[118] = 3.66750000000000000e-09 states[119] = 7.30739999999999981e-09 constants[153] = 20.0 constants[154] = 1.1 constants[155] = 0.000175 constants[156] = 0.0005 constants[157] = 0.4975 states[120] = 1.76869999999999985e-05 states[121] = 1.0 states[122] = 0.99798 states[123] = 1.0 states[124] = 0.98747 constants[158] = 9.88539999999999992e-04 constants[159] = 0.00092 constants[160] = 0.80737 states[125] = 5.92167467082830967e-01 constants[161] = 0.004375 constants[162] = 0.11348 constants[163] = 0.48302 constants[164] = 15.0 constants[165] = 0.01833 states[126] = 145.62 constants[166] = 0.52218 constants[167] = 21.58 constants[168] = 2.642 constants[169] = 11.89 constants[170] = 50.0 constants[171] = 0.8 constants[172] = 1.4 constants[173] = 0.029268 states[127] = 2.36821659448036986e-02 states[128] = 1.28402905095187994e-02 states[129] = 6.37363047239019025e-03 states[130] = 9.17039986149184062e-03 states[131] = 4.29171113639321980e-05 constants[174] = 3.0 constants[175] = 0.16957 constants[176] = 2.12570000000000006e-04 constants[177] = 0.12481 constants[178] = 0.125 constants[179] = 0.2 constants[180] = 0.35 constants[181] = 1.167 constants[182] = 0.7629 constants[183] = 0.9024 constants[184] = 0.0 constants[185] = 20.0 constants[186] = 2.5 constants[187] = 4.0 constants[188] = 0.0196 constants[189] = 3.71 constants[190] = 8314.0 constants[191] = 310.0 constants[192] = 0.5 constants[193] = 0.0388 constants[194] = 0.1 constants[195] = 50.0 states[132] = 2.25475702283052998e-01 states[133] = 2.76490711096605019e-03 states[134] = 9.04820284659604013e-02 constants[196] = 2.4984 constants[197] = 11.359 constants[198] = 0.3755 states[135] = 1.92819110624504991e-01 constants[199] = 100.0 constants[200] = 100.0 constants[201] = 100.0 states[136] = 4.89063888619455989e-01 states[137] = 6.64997605558790977e-02 states[138] = 6.46928309115710060e-02 constants[202] = 0.1088 constants[203] = 0.4612 constants[204] = 0.3755 states[139] = 1.26950532507959013e-01 states[140] = 8.17161796756963987e-01 states[141] = 1.74057375932567010e-01 states[142] = 2.05444874210056000e-01 states[143] = 2.49911886495889995e-01 constants[205] = 0.001 constants[206] = 1.95259999999999991e-05 constants[207] = 0.001469 constants[208] = 0.2 constants[209] = 0.3 states[144] = 2.82662056977524001e-02 constants[210] = 0.0035731 constants[211] = 0.010145 constants[212] = 0.5 constants[213] = 100.0 constants[214] = 1000.0 constants[215] = 1.00000-constants[36] constants[216] = (1.00000-constants[37])-constants[38] constants[217] = constants[57]*0.000983300 constants[218] = constants[57]*0.00133000 constants[219] = constants[55] constants[220] = constants[57]*0.00650000 constants[221] = constants[56] constants[222] = ((constants[51]+constants[59])*(constants[54]+constants[59]))/constants[54] constants[223] = ((constants[48]*constants[53])*(constants[47]+constants[59]))*(constants[52]+constants[59]) constants[224] = ((constants[47]+constants[59])*(constants[48]+constants[59]))/constants[48] constants[225] = ((constants[51]+constants[59])*(constants[54]+constants[59]))/constants[54] constants[226] = ((constants[48]*constants[53])*(constants[47]+constants[59]))*(constants[52]+constants[59]) constants[227] = constants[0]/(constants[1]+constants[0]) constants[228] = constants[69]+constants[67] constants[229] = 1.50000*constants[77] constants[230] = constants[70]*constants[68] constants[231] = constants[86] constants[232] = constants[85] constants[233] = constants[85] constants[234] = constants[85] constants[235] = 7.50000e-14*1.00000e+06 constants[236] = 5.00000e-15*1.00000e+06 constants[237] = 9.00000e-15*1.00000e+06 constants[238] = (constants[190]*constants[191])/constants[87] constants[239] = 0.500000*(power(constants[99]/5.40000, 1.0/2)) constants[240] = 0.0138542*(power(constants[99]/5.40000, 1.0/2)) constants[241] = (2.15000*8.25000)*1.10000 constants[242] = constants[99]/(constants[99]+constants[146]) constants[243] = 0.666700*4.75000 constants[244] = constants[166]*constants[188] constants[245] = 0.0100000/50.0000 constants[246] = constants[199]*constants[196] constants[247] = constants[200]*constants[197] constants[248] = constants[201]*constants[198] constants[249] = constants[199] constants[250] = constants[200] constants[251] = constants[201] constants[252] = constants[196] constants[253] = constants[197] constants[254] = constants[198] constants[255] = (constants[209]/(1.00000-constants[209]))*constants[205]+constants[209]*constants[208] constants[256] = -80.0000 constants[257] = 0.156290*constants[220] constants[258] = (((constants[46]*constants[47])*constants[53])*(constants[52]+constants[59]))*(constants[48]+constants[59]) constants[259] = (((constants[46]*constants[47])*constants[53])*(constants[52]+constants[59]))*(constants[48]+constants[59]) constants[260] = 1.00000/(1.00000+0.350000) constants[261] = constants[70]+constants[68] constants[262] = (constants[69]*constants[68]+constants[67]*constants[70])+constants[230] constants[263] = 1.00000/constants[238] constants[264] = constants[195]*constants[245] constants[265] = constants[249]*constants[202] constants[266] = constants[250]*constants[203] constants[267] = constants[251]*constants[204] constants[268] = constants[199] constants[269] = constants[200] constants[270] = constants[201] constants[271] = (((constants[50]*constants[52])*constants[48])*(constants[47]+constants[59]))*(constants[53]+constants[59]) constants[272] = (((constants[50]*constants[52])*constants[48])*(constants[47]+constants[59]))*(constants[53]+constants[59]) constants[273] = 1.00000+(constants[16]-constants[17])/constants[24] constants[274] = constants[228]+constants[261] constants[275] = constants[87]*constants[263] constants[276] = constants[268]*constants[252] constants[277] = constants[269]*constants[253] constants[278] = constants[270]*constants[254] constants[279] = ((((constants[46]*constants[50])*constants[47])*constants[52])*(constants[48]+constants[59]))*(constants[53]+constants[59]) constants[280] = ((((constants[46]*constants[50])*constants[47])*constants[52])*(constants[48]+constants[59]))*(constants[53]+constants[59]) constants[281] = (constants[24]/2.00000)*(power(power(constants[273], 2.00000)+(4.00000*constants[17])/constants[24], 1.0/2)-constants[273]) constants[282] = (((constants[16]+constants[29])+constants[20])+constants[23])-constants[28] constants[283] = ((constants[16]*constants[23]+constants[29]*constants[20])+constants[20]*constants[23])-constants[28]*(constants[20]+constants[23]) constants[284] = (constants[28]*constants[20])*constants[23] constants[285] = ((((-constants[284]/27.0000)*(power(constants[282], 3.00000))-(((constants[282]*constants[282])*constants[283])*constants[283])/108.000)+((constants[282]*constants[283])*constants[284])/6.00000)+(power(constants[283], 3.00000))/27.0000)+(constants[284]*constants[284])/4.00000 constants[286] = custom_piecewise([less(constants[285] , 0.00000), power(-constants[285], 1.0/2) , True, 0.00000]) constants[287] = (((custom_piecewise([greater(constants[285] , 0.00000), power(constants[285], 1.0/2) , True, 0.00000]))+constants[284]/2.00000)+(constants[282]*constants[283])/6.00000)-(power(constants[282], 3.00000))/27.0000 constants[288] = arctan(constants[286]/constants[287])/3.00000 constants[289] = power(constants[287]*constants[287]+constants[286]*constants[286], 1.00000/6.00000) constants[290] = (constants[283]/3.00000-(constants[282]*constants[282])/9.00000)/(constants[289]*constants[289]) constants[291] = (constants[289]*cos(constants[288]))*(1.00000-constants[290])-constants[282]/3.00000 constants[292] = 1.00000+(constants[29]-constants[30])/constants[25] constants[293] = (constants[25]/2.00000)*(power(power(constants[292], 2.00000)+(4.00000*constants[30])/constants[25], 1.0/2)-constants[292]) constants[294] = 0.850000*0.0250000 constants[295] = 224.000*constants[294] constants[296] = 3.00000*constants[294] constants[297] = 0.150000*0.0250000 constants[298] = (1.00000-constants[39])-constants[40] constants[299] = 1.00000-constants[41] constants[300] = ((2.00000*constants[89])*constants[90])*(constants[90]+constants[88]) constants[301] = 2.00000*constants[300] constants[302] = constants[301]/constants[87] constants[303] = (((1000.00*constants[89])*constants[90])*constants[90])*constants[88] constants[304] = constants[303]*0.00200000 constants[305] = constants[302]/(2.00000*constants[304]) constants[306] = 0.0200000*constants[303] constants[307] = constants[303]*0.678000 constants[308] = constants[302]/(2.00000*constants[307]) constants[309] = 0.0400000*constants[303] constants[310] = constants[303]*0.00480000 constants[311] = constants[303]*0.0552000 constants[312] = constants[310]/constants[311] constants[313] = constants[311]/constants[307] constants[314] = constants[303]*0.0200000 constants[315] = constants[314]/constants[304] constants[316] = -constants[302]/constants[314] constants[317] = -constants[310]/constants[314] constants[318] = constants[314]/constants[307] constants[319] = constants[303]/constants[306] constants[320] = ((constants[11]*constants[260])*constants[296])*constants[319] constants[321] = constants[303]/constants[307] constants[322] = (constants[39]*constants[294])*constants[319] constants[323] = (constants[41]*constants[297])*constants[319] constants[324] = (((1.00000-constants[10])*(1.00000-constants[260]))*constants[296])*constants[321] constants[325] = (constants[298]*constants[294])*constants[321] constants[326] = constants[303]/constants[309] constants[327] = (1.00000-constants[193])-constants[194] constants[328] = (((1.00000-constants[11])*constants[260])*constants[296])*constants[321] constants[329] = ((constants[10]*(1.00000-constants[260]))*constants[296])*constants[326] constants[330] = (constants[40]*constants[294])*constants[326] constants[331] = (constants[299]*constants[297])*constants[326] constants[332] = constants[302]/constants[307] constants[333] = (1.00000-constants[175])-constants[176] constants[334] = constants[175]+constants[176] constants[335] = (constants[327]*constants[192])*constants[321] constants[336] = 0.200000*constants[192] constants[337] = constants[314]/constants[307] constants[338] = power(constants[184]*1.00000, constants[181]) constants[339] = 1.00000-constants[178] constants[340] = constants[193] constants[341] = constants[302]/constants[314] constants[342] = (1.00000-constants[338]/(constants[167]+constants[338]))*constants[173] constants[343] = constants[178]-constants[177] constants[344] = (constants[340]*constants[336])*constants[319] constants[345] = constants[194] constants[346] = power(87.8251, 4.00000) constants[347] = (constants[342]*constants[175])*constants[319] constants[348] = (constants[342]*constants[333])*constants[321] constants[349] = (constants[342]*constants[176])*constants[326] constants[350] = power(constants[184]*1.00000, constants[182]) constants[351] = (constants[345]*constants[336])*constants[326] constants[352] = (1.00000-constants[340])-constants[345] constants[353] = 1.00000+(constants[118]-constants[110])/constants[116] constants[354] = power(constants[184]*1.00000, constants[183]) constants[355] = (constants[352]*constants[336])*constants[321] constants[356] = (constants[116]/2.00000)*(power(power(constants[353], 2.00000)+(4.00000*constants[110])/constants[116], 1.0/2)-constants[353]) constants[357] = constants[180]/(1.00000-constants[180]) constants[358] = 1.00000+(constants[118]-constants[82])/constants[113] constants[359] = constants[357]*(constants[178]*(((1.00000+constants[189])-constants[189]*constants[334])-constants[179])+constants[334]*(constants[179]-1.00000))+(constants[189]*constants[178])*(constants[179]-constants[334]) constants[360] = (constants[113]/2.00000)*(power(power(constants[358], 2.00000)+(4.00000*constants[82])/constants[113], 1.0/2)-constants[358]) constants[361] = constants[178]*((1.00000+constants[189])+constants[179]*(constants[357]-constants[189]))-constants[179]*(1.00000+constants[357]) constants[362] = power(constants[137], 3.00000) constants[363] = (constants[359]/constants[361])*constants[173] constants[364] = power(constants[138], 3.00000) constants[365] = ((constants[179]-constants[334])*constants[173]+constants[179]*constants[363])/((1.00000+constants[189])*constants[178]-constants[179]) constants[366] = power(constants[100], 3.00000) constants[367] = ((constants[189]*constants[178])*constants[365])/constants[363] constants[368] = (1.00000-constants[354]/(constants[169]+constants[354]))*constants[365] constants[369] = constants[159]-0.000170000 constants[370] = (constants[368]*constants[177])*constants[319] constants[371] = (constants[368]*constants[339])*constants[321] constants[372] = (constants[368]*constants[343])*constants[326] constants[373] = 1.00000-constants[367] constants[374] = constants[159]*(1.00000-0.460000) constants[375] = (1.00000-constants[350]/(constants[168]+constants[350]))*constants[363] constants[376] = constants[374] constants[377] = (constants[375]*constants[367])*constants[319] constants[378] = (constants[375]*constants[373])*constants[321] constants[379] = 3.25000*constants[161] constants[380] = (constants[193]*constants[192])*constants[319] constants[381] = (((constants[16]+constants[29])+constants[26])+constants[27])-constants[380] constants[382] = ((constants[16]*constants[27]+constants[29]*constants[26])+constants[26]*constants[27])-constants[380]*(constants[26]+constants[27]) constants[383] = (constants[380]*constants[26])*constants[27] constants[384] = ((((-constants[383]/27.0000)*(power(constants[381], 3.00000))-(((constants[381]*constants[381])*constants[382])*constants[382])/108.000)+((constants[381]*constants[382])*constants[383])/6.00000)+(power(constants[382], 3.00000))/27.0000)+(constants[383]*constants[383])/4.00000 constants[385] = custom_piecewise([less(constants[384] , 0.00000), power(-constants[384], 1.0/2) , True, 0.00000]) constants[386] = (((custom_piecewise([greater(constants[384] , 0.00000), power(constants[384], 1.0/2) , True, 0.00000]))+constants[383]/2.00000)+(constants[381]*constants[382])/6.00000)-(power(constants[381], 3.00000))/27.0000 constants[387] = arctan(constants[385]/constants[386])/3.00000 constants[388] = power(constants[386]*constants[386]+constants[385]*constants[385], 1.00000/6.00000) constants[389] = (constants[382]/3.00000-(constants[381]*constants[381])/9.00000)/(constants[388]*constants[388]) constants[390] = (constants[388]*cos(constants[387]))*(1.00000-constants[389])-constants[381]/3.00000 constants[391] = ((constants[16]-constants[17])+constants[281])/((constants[291]/constants[20]+1.00000)*(constants[390]/constants[26]+1.00000)) constants[392] = ((constants[281]*constants[391])*constants[390])/(constants[24]*constants[26]) constants[393] = (constants[392]*constants[291])/constants[20] constants[394] = ((constants[29]-constants[30])+constants[293])/((constants[291]/constants[23]+1.00000)*(constants[390]/constants[27]+1.00000)) constants[395] = ((constants[293]*constants[394])*constants[390])/(constants[25]*constants[27]) constants[396] = (constants[395]*constants[291])/constants[23] constants[397] = 0.0269000+constants[392]/constants[17] constants[398] = 0.0329000+constants[395]/constants[30] constants[399] = (constants[194]*constants[192])*constants[326] constants[400] = 1.00000+(constants[118]-constants[399])/constants[117] constants[401] = (constants[117]/2.00000)*(power(power(constants[400], 2.00000)+(4.00000*constants[399])/constants[117], 1.0/2)-constants[400]) constants[402] = ((constants[118]-constants[110])+constants[356])/((1.00000+constants[360]/constants[113])*(1.00000+constants[401]/constants[117])) constants[403] = ((constants[356]*constants[402])*constants[401])/(constants[116]*constants[117]) constants[404] = 0.0306000+constants[403]/constants[110] constants[405] = (constants[403]*constants[360])/constants[113] constants[406] = -constants[302]/constants[307] constants[407] = -constants[302]/constants[307] constants[408] = constants[314]/constants[307] constants[409] = (3.00000*constants[302])/constants[314] return (states, constants) def computeRates(voi, states, constants): rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic rates[7] = 0.00100000*(constants[219]*states[8]-(constants[221]*states[0])*states[7]) rates[9] = 0.00100000*(constants[55]*states[1]-(constants[56]*states[10])*states[9]) rates[15] = 0.00100000*(constants[55]*states[2]-(constants[56]*states[16])*states[15]) rates[20] = 0.00100000*(constants[219]*states[21]-(constants[221]*states[22])*states[20]) rates[23] = 0.00100000*(constants[55]*states[3]-(constants[56]*states[24])*states[23]) rates[31] = 0.00100000*(((constants[74]*states[19])*(1.00000-states[31]))/(constants[64]+(1.00000-states[31])*1.00000)-((constants[76]*constants[66])*states[31])/(constants[65]+states[31]*1.00000)) rates[104] = 0.00100000*(((constants[124]*states[29])*(1.00000-states[104]))/(constants[121]+(1.00000-states[104])*1.00000)-((constants[125]*constants[82])*states[104])/(constants[122]+states[104]*1.00000)) rates[105] = 0.00100000*(((constants[128]*states[6])*(1.00000-states[105]))/(constants[126]+(1.00000-states[105])*1.00000)-((constants[129]*constants[28])*states[105])/(constants[127]+states[105]*1.00000)) rates[115] = 0.00100000*(((constants[145]*states[6])*(1.00000-states[115]))/(constants[142]+(1.00000-states[115])*1.00000)-((constants[149]*constants[28])*states[115])/(constants[143]+states[115]*1.00000)) rates[127] = 0.00100000*((constants[188]*states[6])*(constants[377]-states[127])-constants[244]*states[127]) rates[128] = 0.00100000*((constants[188]*states[19])*(constants[378]-states[128])-constants[244]*states[128]) rates[129] = 0.00100000*((constants[188]*states[6])*(constants[370]-states[129])-constants[244]*states[129]) rates[130] = 0.00100000*((constants[188]*states[19])*(constants[371]-states[130])-constants[244]*states[130]) rates[131] = 0.00100000*((constants[188]*states[29])*(constants[372]-states[131])-constants[244]*states[131]) rates[132] = 0.00100000*(constants[201]*states[133]-(constants[248]*states[132])*states[6]) rates[133] = 0.00100000*(((constants[200]*states[134])*states[43]-(constants[247]+constants[201])*states[133])+(constants[248]*states[132])*states[6]) rates[6] = 0.00100000*(((constants[201]*states[133]-(constants[248]*states[132])*states[6])+constants[264]*states[135])-(constants[195]*(constants[344]-states[135]))*states[6]) rates[135] = 0.00100000*((constants[195]*(constants[344]-states[135]))*states[6]-constants[264]*states[135]) rates[136] = 0.00100000*(constants[251]*states[137]-(constants[267]*states[136])*states[19]) rates[137] = 0.00100000*(((constants[250]*states[138])*states[44]-(constants[266]+constants[251])*states[137])+(constants[267]*states[136])*states[19]) rates[19] = 0.00100000*(((constants[251]*states[137]-(constants[267]*states[136])*states[19])+constants[264]*states[139])-(constants[195]*(constants[355]-states[139]))*states[19]) rates[139] = 0.00100000*((constants[195]*(constants[355]-states[139]))*states[19]-constants[264]*states[139]) rates[140] = 0.00100000*(constants[270]*states[141]-(constants[278]*states[140])*states[29]) rates[141] = 0.00100000*(((constants[269]*states[142])*states[45]-(constants[277]+constants[270])*states[141])+(constants[278]*states[140])*states[29]) rates[29] = 0.00100000*(((constants[270]*states[141]-(constants[278]*states[140])*states[29])+constants[264]*states[143])-(constants[195]*(constants[351]-states[143]))*states[29]) rates[143] = 0.00100000*((constants[195]*(constants[351]-states[143]))*states[29]-constants[264]*states[143]) algebraic[3] = constants[393]-states[4] rates[4] = 0.00100000*(((constants[31]*states[6])*algebraic[3])/(constants[18]+algebraic[3])-((constants[33]*constants[28])*states[4])/(constants[21]+states[4])) algebraic[4] = constants[396]-states[5] rates[5] = 0.00100000*(((constants[32]*states[6])*algebraic[4])/(constants[19]+algebraic[4])-((constants[34]*constants[28])*states[5])/(constants[22]+states[5])) algebraic[15] = constants[405]-states[103] rates[103] = 0.00100000*(((constants[119]*states[29])*algebraic[15])/(constants[114]+algebraic[15])-((constants[120]*constants[82])*states[103])/(constants[115]+states[103])) algebraic[25] = 1.00000/(1.00000+exp((states[50]+91.0000)/6.10000)) rates[116] = (algebraic[25]-states[116])/constants[150] algebraic[29] = constants[255]-states[144] rates[144] = 0.00100000*(((constants[211]*states[19])*algebraic[29])/(constants[207]+algebraic[29])-((constants[210]*constants[66])*states[144])/(constants[206]+states[144])) algebraic[33] = (constants[322]-states[12])-states[11] rates[12] = 0.00100000*((constants[220]*states[6])*algebraic[33]-constants[257]*states[12]) algebraic[35] = (constants[325]-states[18])-states[17] rates[18] = 0.00100000*((constants[220]*states[19])*algebraic[35]-constants[257]*states[18]) algebraic[36] = (constants[330]-states[26])-states[25] rates[26] = 0.00100000*((constants[220]*states[29])*algebraic[36]-constants[257]*states[26]) algebraic[11] = 0.0250000/(1.00000+exp((states[50]+58.0000)/5.00000)) algebraic[39] = 0.200000/(1.00000+exp((states[50]+19.0000)/-9.00000)) rates[67] = (algebraic[11]/(algebraic[11]+algebraic[39])-states[67])/constants[109] algebraic[12] = 1.00000/((0.000600000*(states[50]-1.73840))/(1.00000-exp(-0.136000*(states[50]-1.73840)))-(0.000300000*(states[50]+38.3608))/(1.00000-exp(0.152200*(states[50]+38.3608)))) algebraic[40] = 1.00000/(1.00000+exp((states[50]+10.0850)/-4.25000)) rates[68] = (algebraic[40]-states[68])/algebraic[12] algebraic[13] = 0.00739900/(1.00000+exp((constants[263]*(states[50]-0.0311960))/-0.800190)) algebraic[41] = 0.00569920/(1.00000+exp((constants[263]*(states[50]-0.0415200))/1.34890)) rates[73] = algebraic[41]*states[80]-states[73]*(4.00000*algebraic[13]) algebraic[14] = 0.00994150/(1.00000+exp((constants[263]*(states[50]-0.0448090))/-0.581720)) algebraic[42] = 0.00332010/(1.00000+exp((constants[263]*(states[50]-0.0942170))/0.953640)) rates[88] = algebraic[42]*states[95]-states[88]*(4.00000*algebraic[14]) algebraic[16] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), 0.00000 , True, 0.135000*exp(((87.0000+states[50])+constants[130])/-6.80000)]) algebraic[43] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), 1.00000/(0.130000*(1.00000+exp(((states[50]+constants[130])+27.4034)/-11.1000))) , True, 3.56000*exp(0.0790000*((states[50]+constants[130])+7.00000))+310000.*exp(0.350000*((states[50]+constants[130])+7.00000))]) rates[106] = algebraic[16]*(1.00000-states[106])-algebraic[43]*states[106] algebraic[17] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), 0.00000 , True, ((-127140.*exp(0.244400*(states[50]+constants[130]))+-6.94800e-05*exp(-0.0439100*(states[50]+constants[130])))*((states[50]+constants[130])+37.7800))/(1.00000+exp(0.311000*((states[50]+constants[130])+79.2300)))]) algebraic[44] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), (0.300000*exp(-2.53500e-07*(states[50]+constants[130])))/(1.00000+exp(-0.100000*((states[50]+constants[130])+32.0000))) , True, (0.121200*exp(-0.0105200*(states[50]+constants[130])))/(1.00000+exp(-0.137800*((states[50]+constants[130])+40.1400)))]) rates[107] = algebraic[17]*(1.00000-states[107])-algebraic[44]*states[107] algebraic[18] = (0.320000*(states[50]+58.4729))/(1.00000-exp(-0.100000*(states[50]+58.4729))) algebraic[45] = 0.0800000*exp((13.7299-states[50])/11.0000) rates[108] = algebraic[18]*(1.00000-states[108])-algebraic[45]*states[108] algebraic[19] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, 0.135000*exp((87.0000+states[50])/-6.80000)]) algebraic[46] = custom_piecewise([greater_equal(states[50] , -40.0000), 1.00000/(0.130000*(1.00000+exp((states[50]+27.4034)/-11.1000))) , True, 3.56000*exp(0.0790000*(states[50]+7.00000))+310000.*exp(0.350000*(states[50]+7.00000))]) rates[109] = algebraic[19]*(1.00000-states[109])-algebraic[46]*states[109] algebraic[20] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, ((-127140.*exp(0.244400*states[50])-6.94800e-05*exp(-0.0439100*states[50]))*(states[50]+37.7800))/(1.00000+exp(0.311000*(states[50]+79.2300)))]) algebraic[47] = custom_piecewise([greater_equal(states[50] , -40.0000), (0.300000*exp(-2.53500e-07*states[50]))/(1.00000+exp(-0.100000*(states[50]+32.0000))) , True, (0.121200*exp(-0.0105200*states[50]))/(1.00000+exp(-0.137800*(states[50]+40.1400)))]) rates[110] = algebraic[20]*(1.00000-states[110])-algebraic[47]*states[110] algebraic[21] = (0.320000*(states[50]+58.4729))/(1.00000-exp(-0.100000*(states[50]+58.4729))) algebraic[48] = 0.0800000*exp((13.7299-states[50])/11.0000) rates[111] = algebraic[21]*(1.00000-states[111])-algebraic[48]*states[111] algebraic[22] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, 0.135000*exp(((87.0000+states[50])+constants[132])/-6.80000)]) algebraic[49] = custom_piecewise([greater_equal(states[50] , -40.0000), (1.00000/0.130000)/(1.00000+exp(((states[50]+constants[132])+27.4034)/-11.1000)) , True, 3.56000*exp(0.0790000*((states[50]+constants[132])+7.00000))+310000.*exp(0.350000*((states[50]+constants[132])+7.00000))]) rates[112] = algebraic[22]*(1.00000-states[112])-algebraic[49]*states[112] algebraic[23] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, ((-127140.*exp(0.244400*(states[50]+constants[132]))-6.94800e-05*exp(-0.0439100*(states[50]+constants[132])))*((states[50]+constants[132])+37.7800))/(1.00000+exp(0.311000*((states[50]+constants[132])+79.2300)))]) algebraic[50] = custom_piecewise([greater_equal(states[50] , -40.0000), (0.300000*exp(-2.53500e-07*(states[50]+constants[132])))/(1.00000+exp(-0.100000*((states[50]+constants[132])+32.0000))) , True, (0.121200*exp(-0.0105200*(states[50]+constants[132])))/(1.00000+exp(-0.137800*((states[50]+constants[132])+40.1400)))]) rates[113] = algebraic[23]*(1.00000-states[113])-algebraic[50]*states[113] algebraic[24] = (0.320000*((states[50]+constants[131])+58.4729))/(1.00000-exp(-0.100000*((states[50]+constants[131])+58.4729))) algebraic[51] = 0.0800000*exp(((states[50]+constants[131])-13.7299)/-11.0000) rates[114] = algebraic[24]*(1.00000-states[114])-algebraic[51]*states[114] algebraic[26] = (0.320000*(states[50]+47.1300))/(1.00000-exp(-0.100000*(states[50]+47.1300))) algebraic[52] = 0.0800000*exp(states[50]/-11.0000) rates[117] = algebraic[26]*(1.00000-states[117])-algebraic[52]*states[117] algebraic[57] = (constants[323]-states[14])-states[13] rates[14] = 0.00100000*((constants[220]*states[6])*algebraic[57]-constants[257]*states[14]) algebraic[60] = (constants[331]-states[28])-states[27] rates[28] = 0.00100000*((constants[220]*states[29])*algebraic[60]-constants[257]*states[28]) algebraic[27] = 1.00000/(1.00000+exp((states[50]+9.43700)/-7.13300)) algebraic[53] = 1.00000/(1.00000+exp((states[50]-18.4099)/-29.3814)) algebraic[65] = 1.00000/(1.00000+exp((states[50]+100.000)/29.3814)) algebraic[77] = 1.00000/(algebraic[53]/1.20890+3.50000*algebraic[65]) rates[120] = (algebraic[27]-states[120])/algebraic[77] algebraic[66] = ((1.00000/(1.00000+exp((states[50]+19.0000)/-9.00000)))/0.500000)/9.79530 algebraic[28] = 1.00000/(1.00000+exp((states[50]+58.0000)/5.00000)) algebraic[78] = 0.0479600*algebraic[28] rates[121] = algebraic[78]*(1.00000-states[121])-algebraic[66]*states[121] algebraic[79] = 0.0214400*algebraic[28] rates[122] = algebraic[79]*(1.00000-states[122])-algebraic[66]*states[122] algebraic[54] = (1.00000/(1.00000+exp((states[50]+60.0000)/5.00000)))/250.000 algebraic[80] = 2.46000*algebraic[54] rates[123] = algebraic[80]*(1.00000-states[123])-algebraic[66]*states[123] algebraic[81] = 0.560340*algebraic[54] rates[124] = algebraic[81]*(1.00000-states[124])-algebraic[66]*states[124] algebraic[63] = 0.0906540*exp((-0.111570*states[50])*constants[263]) algebraic[91] = 0.388390/(1.00000+exp((constants[263]*(states[50]+0.150190))/-0.606930)) rates[74] = (algebraic[41]*states[75]+algebraic[91]*states[85])-states[74]*(2.00000*algebraic[13]+2.00000*algebraic[63]) rates[75] = ((((2.00000*algebraic[13])*states[74]+(2.00000*algebraic[41])*states[76])+(2.00000*algebraic[91])*states[86])+(3.00000*algebraic[63])*states[77])-states[75]*(((algebraic[13]+algebraic[41])+algebraic[91])+2.00000*algebraic[63]) rates[76] = ((algebraic[13]*states[75]+(3.00000*algebraic[91])*states[87])+(3.00000*algebraic[63])*states[78])-states[76]*((2.00000*algebraic[41]+2.00000*algebraic[91])+2.00000*algebraic[63]) rates[77] = (algebraic[41]*states[78]+algebraic[91]*states[75])-states[77]*(algebraic[13]+3.00000*algebraic[63]) rates[78] = ((algebraic[13]*states[77]+(2.00000*algebraic[91])*states[76])+(4.00000*algebraic[63])*states[79])-states[78]*((algebraic[41]+algebraic[91])+3.00000*algebraic[63]) algebraic[75] = 0.00311240+(0.0283300-0.00311240)/(1.00000+exp((constants[263]*(states[50]+0.0516600))/1.55220)) rates[79] = (algebraic[91]*states[78]-states[79]*(4.00000*algebraic[63]+constants[111]))+algebraic[75]*states[69] rates[80] = (((4.00000*algebraic[13])*states[73]+(2.00000*algebraic[41])*states[81])+algebraic[63]*states[84])-states[80]*((3.00000*algebraic[13]+algebraic[41])+algebraic[91]) rates[81] = (((3.00000*algebraic[13])*states[80]+(3.00000*algebraic[41])*states[82])+algebraic[63]*states[85])-states[81]*((2.00000*algebraic[13]+2.00000*algebraic[41])+2.00000*algebraic[91]) rates[82] = (((2.00000*algebraic[13])*states[81]+(4.00000*algebraic[41])*states[83])+algebraic[63]*states[86])-states[82]*((algebraic[13]+3.00000*algebraic[41])+3.00000*algebraic[91]) rates[83] = (algebraic[13]*states[82]+algebraic[63]*states[87])-states[83]*(4.00000*algebraic[41]+4.00000*algebraic[91]) rates[84] = (algebraic[41]*states[85]+algebraic[91]*states[80])-states[84]*(3.00000*algebraic[13]+algebraic[63]) rates[85] = ((((3.00000*algebraic[13])*states[84]+(2.00000*algebraic[41])*states[86])+(2.00000*algebraic[91])*states[81])+(2.00000*algebraic[63])*states[74])-states[85]*(((2.00000*algebraic[13]+algebraic[41])+algebraic[91])+algebraic[63]) rates[86] = ((((2.00000*algebraic[13])*states[85]+(3.00000*algebraic[41])*states[87])+(3.00000*algebraic[91])*states[82])+(2.00000*algebraic[63])*states[75])-states[86]*(((algebraic[13]+2.00000*algebraic[41])+2.00000*algebraic[91])+algebraic[63]) rates[87] = ((algebraic[13]*states[86]+(4.00000*algebraic[91])*states[83])+(2.00000*algebraic[63])*states[76])-states[87]*((3.00000*algebraic[41]+3.00000*algebraic[91])+algebraic[63]) algebraic[64] = 0.0657000*exp((-0.118990*states[50])*constants[263]) algebraic[93] = 0.563560/(1.00000+exp((constants[263]*(states[50]+0.179860))/-0.583810)) rates[89] = (algebraic[42]*states[90]+algebraic[93]*states[100])-states[89]*(2.00000*algebraic[14]+2.00000*algebraic[64]) rates[90] = ((((2.00000*algebraic[14])*states[89]+(2.00000*algebraic[42])*states[91])+(2.00000*algebraic[93])*states[101])+(3.00000*algebraic[64])*states[92])-states[90]*(((algebraic[14]+algebraic[42])+algebraic[93])+2.00000*algebraic[64]) rates[91] = ((algebraic[14]*states[90]+(3.00000*algebraic[93])*states[102])+(3.00000*algebraic[64])*states[93])-states[91]*((2.00000*algebraic[42]+2.00000*algebraic[93])+2.00000*algebraic[64]) rates[92] = (algebraic[42]*states[93]+algebraic[93]*states[90])-states[92]*(algebraic[14]+3.00000*algebraic[64]) rates[93] = ((algebraic[14]*states[92]+(2.00000*algebraic[93])*states[91])+(4.00000*algebraic[64])*states[94])-states[93]*((algebraic[42]+algebraic[93])+3.00000*algebraic[64]) algebraic[76] = 0.000385250+(0.0124060-0.000385250)/(1.00000+exp((constants[263]*(states[50]+0.0641180))/0.779920)) rates[94] = (algebraic[93]*states[93]-states[94]*(4.00000*algebraic[64]+constants[112]))+algebraic[76]*states[71] rates[95] = (((4.00000*algebraic[14])*states[88]+(2.00000*algebraic[42])*states[96])+algebraic[64]*states[99])-states[95]*((3.00000*algebraic[14]+algebraic[42])+algebraic[93]) rates[96] = (((3.00000*algebraic[14])*states[95]+(3.00000*algebraic[42])*states[97])+algebraic[64]*states[100])-states[96]*((2.00000*algebraic[14]+2.00000*algebraic[42])+2.00000*algebraic[93]) rates[97] = (((2.00000*algebraic[14])*states[96]+(4.00000*algebraic[42])*states[98])+algebraic[64]*states[101])-states[97]*((algebraic[14]+3.00000*algebraic[42])+3.00000*algebraic[93]) rates[98] = (algebraic[14]*states[97]+algebraic[64]*states[102])-states[98]*(4.00000*algebraic[42]+4.00000*algebraic[93]) rates[99] = (algebraic[42]*states[100]+algebraic[93]*states[95])-states[99]*(3.00000*algebraic[14]+algebraic[64]) rates[100] = ((((3.00000*algebraic[14])*states[99]+(2.00000*algebraic[42])*states[101])+(2.00000*algebraic[93])*states[96])+(2.00000*algebraic[64])*states[89])-states[100]*(((2.00000*algebraic[14]+algebraic[42])+algebraic[93])+algebraic[64]) rates[101] = ((((2.00000*algebraic[14])*states[100]+(3.00000*algebraic[42])*states[102])+(3.00000*algebraic[93])*states[97])+(2.00000*algebraic[64])*states[90])-states[101]*(((algebraic[14]+2.00000*algebraic[42])+2.00000*algebraic[93])+algebraic[64]) rates[102] = ((algebraic[14]*states[101]+(4.00000*algebraic[93])*states[98])+(2.00000*algebraic[64])*states[91])-states[102]*((3.00000*algebraic[42]+3.00000*algebraic[93])+algebraic[64]) algebraic[92] = 0.000441980*exp((-1.20220*states[50])*constants[263]) algebraic[103] = 0.000401730*exp((0.000208730*states[50])*constants[263]) rates[69] = (-(algebraic[75]+algebraic[103])*states[69]+algebraic[92]*states[70])+constants[111]*states[79] rates[70] = algebraic[103]*states[69]-algebraic[92]*states[70] algebraic[94] = 0.000237300*exp((-1.97420*states[50])*constants[263]) algebraic[104] = 0.000226520*exp((0.000246900*states[50])*constants[263]) rates[71] = (-(algebraic[76]+algebraic[104])*states[71]+algebraic[94]*states[72])+constants[112]*states[94] rates[72] = algebraic[104]*states[71]-algebraic[94]*states[72] algebraic[5] = ((constants[36]*constants[35])*constants[319]-states[8])-states[7] algebraic[34] = (algebraic[5]*(constants[59]+constants[51])-states[14]*(constants[51]+constants[59]))+(constants[49]*constants[51])*(1.00000+constants[59]/constants[54]) algebraic[58] = (-states[14]*constants[49])*constants[51] algebraic[70] = (-algebraic[34]+power(algebraic[34]*algebraic[34]-(4.00000*constants[222])*algebraic[58], 1.0/2))/(2.00000*constants[222]) algebraic[86] = algebraic[5]/(1.00000+(algebraic[70]/constants[49])*(1.00000+constants[59]/constants[51])) algebraic[98] = (algebraic[70]*algebraic[86])/constants[49] algebraic[108] = (algebraic[98]*constants[59])/constants[51] rates[8] = 0.00100000*((constants[44]*algebraic[98]+constants[42]*algebraic[108])-constants[219]*states[8]) rates[0] = 0.00100000*((constants[44]*algebraic[98]+constants[42]*algebraic[108])-(constants[221]*states[0])*states[7]) algebraic[8] = ((constants[215]*constants[35])*constants[326]-states[21])-states[20] algebraic[37] = (algebraic[8]*(constants[59]+constants[51])-states[28]*(constants[51]+constants[59]))+(constants[49]*constants[51])*(1.00000+constants[59]/constants[54]) algebraic[61] = (-states[28]*constants[49])*constants[51] algebraic[73] = (-algebraic[37]+power(algebraic[37]*algebraic[37]-(4.00000*constants[225])*algebraic[61], 1.0/2))/(2.00000*constants[225]) algebraic[89] = algebraic[8]/(1.00000+(algebraic[73]/constants[49])*(1.00000+constants[59]/constants[51])) algebraic[101] = (algebraic[73]*algebraic[89])/constants[49] algebraic[111] = (constants[59]/constants[51])*algebraic[101] rates[21] = 0.00100000*((constants[44]*algebraic[101]+constants[42]*algebraic[111])-constants[219]*states[21]) rates[22] = 0.00100000*((constants[44]*algebraic[101]+constants[42]*algebraic[111])-(constants[221]*states[22])*states[20]) algebraic[7] = ((constants[216]*constants[295])*constants[321]-states[2])-states[15] algebraic[59] = (algebraic[7]*(constants[47]+constants[59])-algebraic[35]*(constants[47]+constants[59]))+(constants[46]*constants[47])*(1.00000+constants[59]/constants[48]) algebraic[71] = (-algebraic[35]*constants[47])*constants[46] algebraic[87] = (-algebraic[59]+power(algebraic[59]*algebraic[59]-(4.00000*constants[224])*algebraic[71], 1.0/2))/(2.00000*constants[224]) algebraic[99] = algebraic[7]/(1.00000+(algebraic[87]/constants[46])*(1.00000+constants[59]/constants[47])) algebraic[110] = ((constants[59]*algebraic[87])*algebraic[99])/(constants[46]*constants[47]) algebraic[119] = (algebraic[99]*algebraic[87])/constants[46] rates[2] = 0.00100000*((constants[45]*algebraic[119]+constants[43]*algebraic[110])-constants[55]*states[2]) rates[16] = 0.00100000*((constants[45]*algebraic[119]+constants[43]*algebraic[110])-(constants[56]*states[16])*states[15]) algebraic[118] = (constants[59]*algebraic[87])/constants[48] rates[17] = 0.00100000*((constants[218]*constants[61])*(algebraic[118]+algebraic[110])-constants[217]*states[17]) algebraic[130] = ((constants[399]-states[142])-states[141])-states[140] rates[142] = 0.00100000*((((constants[268]*algebraic[130])*states[45]-constants[276]*states[142])-(constants[269]*states[142])*states[45])+constants[277]*states[141]) algebraic[1] = power(states[3]*1.00000, constants[12]) algebraic[31] = constants[6]*(constants[8]+algebraic[1]/(constants[3]+algebraic[1])) algebraic[67] = (algebraic[31]*constants[329])*constants[227] algebraic[83] = (constants[236]*(states[43]-states[45]))/constants[309] algebraic[96] = (constants[237]*(states[45]-states[44]))/constants[309] algebraic[107] = (constants[349]*constants[185])/(1.00000+constants[170]/states[45]) algebraic[116] = ((constants[372]+(constants[174]-1.00000)*states[131])*constants[187])/(1.00000+constants[172]/states[45]) algebraic[124] = algebraic[107]+algebraic[116] algebraic[134] = (((-constants[268]*algebraic[130])*states[45]+constants[276]*states[142])-(constants[269]*states[142])*states[45])+constants[277]*states[141] rates[45] = 0.00100000*((((algebraic[134]+algebraic[67])-algebraic[124])+algebraic[83])-algebraic[96]) algebraic[140] = ((constants[380]-states[134])-states[133])-states[132] rates[134] = 0.00100000*((((constants[199]*algebraic[140])*states[43]-constants[246]*states[134])-(constants[200]*states[134])*states[43])+constants[247]*states[133]) algebraic[2] = power(states[1]*1.00000, constants[13]) algebraic[32] = power(states[1]*1.00000, constants[14]) algebraic[56] = (constants[7]*(constants[9]+algebraic[2]/(constants[4]+algebraic[2])))*(1.00000-((1.00000-(constants[15]*algebraic[32])/(constants[5]+algebraic[32]))*states[0])/(constants[2]+states[0])) algebraic[82] = (algebraic[56]*constants[320])*constants[227] algebraic[95] = (constants[236]*(states[43]-states[45]))/constants[306] algebraic[105] = (constants[235]*(states[43]-states[44]))/constants[306] algebraic[114] = (constants[347]*constants[185])/(1.00000+constants[170]/states[43]) algebraic[122] = ((constants[377]+(constants[174]-1.00000)*states[127])*constants[186])/(1.00000+constants[171]/states[43]) algebraic[128] = ((constants[370]+(constants[174]-1.00000)*states[129])*constants[187])/(1.00000+constants[172]/states[43]) algebraic[135] = (algebraic[114]+algebraic[122])+algebraic[128] algebraic[145] = (((-constants[199]*algebraic[140])*states[43]+constants[246]*states[134])-(constants[200]*states[134])*states[43])+constants[247]*states[133] rates[43] = 0.00100000*((((algebraic[145]+algebraic[82])-algebraic[135])-algebraic[95])-algebraic[105]) algebraic[150] = ((constants[335]-states[138])-states[137])-states[136] rates[138] = 0.00100000*((((constants[249]*algebraic[150])*states[44]-constants[265]*states[138])-(constants[250]*states[138])*states[44])+constants[266]*states[137]) algebraic[0] = power(states[2]*1.00000, constants[12]) algebraic[55] = constants[6]*(constants[8]+algebraic[0]/(constants[3]+algebraic[0])) algebraic[68] = (algebraic[55]*constants[324])*constants[227] algebraic[30] = power(states[2]*1.00000, constants[13]) algebraic[84] = constants[7]*(constants[9]+algebraic[30]/(constants[4]+algebraic[30])) algebraic[106] = (algebraic[84]*constants[328])*constants[227] algebraic[115] = (constants[235]*(states[43]-states[44]))/constants[307] algebraic[123] = (constants[237]*(states[45]-states[44]))/constants[307] algebraic[129] = (constants[348]*constants[185])/(1.00000+constants[170]/states[44]) algebraic[136] = ((constants[378]+(constants[174]-1.00000)*states[128])*constants[186])/(1.00000+constants[171]/states[44]) algebraic[141] = ((constants[371]+(constants[174]-1.00000)*states[130])*constants[187])/(1.00000+constants[172]/states[44]) algebraic[146] = (algebraic[129]+algebraic[136])+algebraic[141] algebraic[154] = (((-constants[249]*algebraic[150])*states[44]+constants[265]*states[138])-(constants[250]*states[138])*states[44])+constants[266]*states[137] rates[44] = 0.00100000*(((((algebraic[154]+algebraic[68])+algebraic[106])-algebraic[146])+algebraic[115])+algebraic[123]) algebraic[6] = ((constants[37]*constants[295])*constants[319]-states[1])-states[9] algebraic[69] = (((constants[258]+constants[271])/constants[223]+algebraic[33])+algebraic[57])-algebraic[6] algebraic[85] = ((constants[271]*(algebraic[33]-algebraic[6])+constants[258]*(algebraic[57]-algebraic[6]))+constants[279])/constants[223] algebraic[97] = (algebraic[6]*constants[279])/constants[223] algebraic[109] = ((((-algebraic[97]/27.0000)*(power(algebraic[69], 3.00000))-(((algebraic[69]*algebraic[69])*algebraic[85])*algebraic[85])/108.000)+((algebraic[69]*algebraic[85])*algebraic[97])/6.00000)+(power(algebraic[85], 3.00000))/27.0000)+(algebraic[97]*algebraic[97])/4.00000 algebraic[117] = custom_piecewise([less(algebraic[109] , 0.00000), power(-algebraic[109], 1.0/2) , True, 0.00000]) algebraic[125] = (((custom_piecewise([greater(algebraic[109] , 0.00000), power(algebraic[109], 1.0/2) , True, 0.00000]))+algebraic[97]/2.00000)+(algebraic[69]*algebraic[85])/6.00000)-(power(algebraic[69], 3.00000))/27.0000 algebraic[131] = arctan(algebraic[117]/algebraic[125])/3.00000 algebraic[137] = power(algebraic[125]*algebraic[125]+algebraic[117]*algebraic[117], 1.00000/6.00000) algebraic[142] = (algebraic[85]/3.00000-(algebraic[69]*algebraic[69])/9.00000)/(algebraic[137]*algebraic[137]) algebraic[147] = (algebraic[137]*sin(algebraic[131]))*(1.00000+algebraic[142]) algebraic[151] = (algebraic[137]*cos(algebraic[131]))*(1.00000-algebraic[142])-algebraic[69]/3.00000 algebraic[156] = power(algebraic[151]*algebraic[151]+algebraic[147]*algebraic[147], 1.0/2) algebraic[160] = algebraic[33]/((1.00000+constants[59]/constants[48])+(algebraic[156]*(constants[47]+constants[59]))/(constants[46]*constants[47])) algebraic[164] = (constants[59]*algebraic[160])/constants[48] algebraic[168] = ((constants[59]*algebraic[160])*algebraic[156])/(constants[46]*constants[47]) rates[11] = 0.00100000*((constants[218]*constants[58])*(algebraic[164]+algebraic[168])-constants[217]*states[11]) algebraic[9] = ((constants[38]*constants[295])*constants[326]-states[3])-states[23] algebraic[72] = (((constants[259]+constants[272])/constants[226]+algebraic[36])+algebraic[60])-algebraic[9] algebraic[88] = ((constants[272]*(algebraic[36]-algebraic[9])+constants[259]*(algebraic[60]-algebraic[9]))+constants[280])/constants[226] algebraic[100] = (algebraic[9]*constants[280])/constants[226] algebraic[112] = ((((-algebraic[100]/27.0000)*(power(algebraic[72], 3.00000))-(((algebraic[72]*algebraic[72])*algebraic[88])*algebraic[88])/108.000)+((algebraic[72]*algebraic[88])*algebraic[100])/6.00000)+(power(algebraic[88], 3.00000))/27.0000)+(algebraic[100]*algebraic[100])/4.00000 algebraic[120] = custom_piecewise([less(algebraic[112] , 0.00000), power(-algebraic[112], 1.0/2) , True, 0.00000]) algebraic[126] = (((custom_piecewise([greater(algebraic[112] , 0.00000), power(algebraic[112], 1.0/2) , True, 0.00000]))+algebraic[100]/2.00000)+(algebraic[72]*algebraic[88])/6.00000)-(power(algebraic[72], 3.00000))/27.0000 algebraic[132] = arctan(algebraic[120]/algebraic[126])/3.00000 algebraic[138] = power(algebraic[126]*algebraic[126]+algebraic[120]*algebraic[120], 1.00000/6.00000) algebraic[143] = (algebraic[88]/3.00000-(algebraic[72]*algebraic[72])/9.00000)/(algebraic[138]*algebraic[138]) algebraic[148] = (algebraic[138]*sin(algebraic[132]))*(1.00000+algebraic[143]) algebraic[152] = (algebraic[138]*cos(algebraic[132]))*(1.00000-algebraic[143])-algebraic[72]/3.00000 algebraic[157] = power(algebraic[152]*algebraic[152]+algebraic[148]*algebraic[148], 1.0/2) algebraic[161] = algebraic[36]/((1.00000+constants[59]/constants[48])+(algebraic[157]*(constants[47]+constants[59]))/(constants[46]*constants[47])) algebraic[165] = (constants[59]*algebraic[161])/constants[48] algebraic[169] = ((constants[59]*algebraic[161])*algebraic[157])/(constants[46]*constants[47]) rates[25] = 0.00100000*((constants[218]*constants[62])*(algebraic[165]+algebraic[169])-constants[217]*states[25]) algebraic[176] = algebraic[57]/((1.00000+constants[59]/constants[53])+(algebraic[156]*(constants[52]+constants[59]))/(constants[50]*constants[52])) algebraic[180] = (constants[59]*algebraic[176])/constants[53] algebraic[184] = ((constants[59]*algebraic[176])*algebraic[156])/(constants[50]*constants[52]) rates[13] = 0.00100000*((constants[218]*constants[58])*(algebraic[180]+algebraic[184])-constants[217]*states[13]) algebraic[177] = algebraic[60]/((1.00000+constants[59]/constants[53])+(algebraic[157]*(constants[52]+constants[59]))/(constants[50]*constants[52])) algebraic[181] = (constants[59]*algebraic[177])/constants[53] algebraic[185] = ((constants[59]*algebraic[177])*algebraic[157])/(constants[50]*constants[52]) rates[27] = 0.00100000*((constants[218]*constants[62])*(algebraic[181]+algebraic[185])-constants[217]*states[27]) algebraic[188] = algebraic[168]+constants[60]*algebraic[184] algebraic[172] = (algebraic[160]*algebraic[156])/constants[46] algebraic[192] = (algebraic[176]*algebraic[156])/constants[50] algebraic[196] = algebraic[172]+constants[60]*algebraic[192] rates[1] = 0.00100000*((constants[45]*algebraic[196]+constants[43]*algebraic[188])-constants[55]*states[1]) rates[10] = 0.00100000*((constants[45]*algebraic[196]+constants[43]*algebraic[188])-(constants[56]*states[10])*states[9]) algebraic[189] = algebraic[169]+constants[63]*algebraic[185] algebraic[173] = (algebraic[161]*algebraic[157])/constants[46] algebraic[193] = (algebraic[177]*algebraic[157])/constants[50] algebraic[197] = algebraic[173]+constants[63]*algebraic[193] rates[3] = 0.00100000*((constants[45]*algebraic[197]+constants[43]*algebraic[189])-constants[55]*states[3]) rates[24] = 0.00100000*((constants[45]*algebraic[197]+constants[43]*algebraic[189])-(constants[56]*states[24])*states[23]) algebraic[133] = constants[274]-states[35] algebraic[139] = constants[262]-states[35]*constants[261] algebraic[144] = -constants[230]*states[35] algebraic[149] = -algebraic[133]/3.00000+((2.00000/3.00000)*(power(algebraic[133]*algebraic[133]-3.00000*algebraic[139], 1.0/2)))*cos(arccos((((9.00000*algebraic[133])*algebraic[139]-((2.00000*algebraic[133])*algebraic[133])*algebraic[133])-27.0000*algebraic[144])/(2.00000*(power(algebraic[133]*algebraic[133]-3.00000*algebraic[139], 1.50000))))/3.00000) algebraic[190] = (constants[79]*(1.00000-states[42]))/(1.00000+constants[81]/algebraic[149]) algebraic[194] = algebraic[190]+states[42] algebraic[199] = algebraic[194]/(algebraic[194]+constants[80]) rates[36] = (algebraic[199]-states[36])/constants[231] rates[37] = (algebraic[199]-states[37])/constants[232] rates[38] = (algebraic[199]-states[38])/constants[233] rates[39] = (algebraic[199]-states[39])/constants[234] rates[40] = (algebraic[199]-states[40])/constants[85] algebraic[200] = 1.00000/(1.00000+power(constants[80]/algebraic[194], 2.00000)) rates[41] = (algebraic[200]-states[41])/constants[86] algebraic[198] = (constants[205]-constants[208])+states[144] algebraic[203] = 0.500000*(power(power(algebraic[198], 2.00000)+(4.00000*constants[205])*constants[208], 1.0/2)-algebraic[198]) rates[125] = 0.00100000*(((constants[162]*states[19])*(1.00000-states[125]))/(constants[158]+(1.00000-states[125])*1.00000)-((constants[163]*algebraic[203])*states[125])/(constants[160]+states[125]*1.00000)) algebraic[206] = (constants[28]/constants[319]+constants[82]/constants[326])+algebraic[203]/constants[321] rates[42] = (constants[83]*algebraic[190])*algebraic[194]-(constants[84]*states[42])*(0.100000+(0.900000*algebraic[206])/0.137100) algebraic[155] = 1.00000/((1.00000+exp((4.79800+states[50])/-7.56990))*(1.00000+exp((25.0000+states[50])/-5.00000))) algebraic[195] = 0.590000+(0.800000*exp(0.0520000*(states[50]+13.0000)))/(1.00000+exp(0.132000*(states[50]+13.0000))) algebraic[202] = algebraic[155]/algebraic[195] algebraic[205] = (1.00000-algebraic[155])/algebraic[195] algebraic[179] = 1.00000+power(0.00200000/algebraic[149], 4.00000) algebraic[183] = 6.00000/algebraic[179] algebraic[163] = 1.00000/(1.00000+exp((29.9790+states[50])/3.17750)) algebraic[171] = (0.100000+algebraic[163])/1.10000 algebraic[214] = 1.00000/((70.0000*(1.00000-0.500000*states[36]))*(1.00000+exp((states[50]+49.1000)/10.3490))) algebraic[216] = 1.00000/(1.00000+exp((states[50]+0.213000)/-10.8070)) algebraic[221] = 1.00000/(algebraic[214]+algebraic[216]/38.4940) algebraic[224] = algebraic[171]/algebraic[221] algebraic[233] = (1.00000-algebraic[171])/algebraic[221] rates[51] = ((-((algebraic[202]+algebraic[183])+algebraic[233])*states[51]+algebraic[205]*states[55])+constants[102]*states[54])+algebraic[224]*states[52] rates[55] = ((-((algebraic[205]+algebraic[183])+algebraic[233])*states[55]+algebraic[202]*states[51])+constants[102]*states[58])+algebraic[224]*states[56] algebraic[167] = (0.000100000+algebraic[163])/1.00010 algebraic[159] = 0.100000*states[36] algebraic[175] = 1.00000+power(0.0100000/algebraic[149], 10.0000) algebraic[187] = (32.5000-(18.0000-algebraic[159])/algebraic[179])-10.0000/algebraic[175] algebraic[218] = 1.00000/(algebraic[214]+algebraic[216]/algebraic[187]) algebraic[227] = algebraic[167]/algebraic[218] algebraic[239] = (1.00000-algebraic[167])/algebraic[218] rates[54] = ((-((algebraic[202]+constants[102])+algebraic[239])*states[54]+algebraic[183]*states[51])+algebraic[205]*states[58])+algebraic[227]*states[53] rates[58] = ((-((algebraic[205]+constants[102])+algebraic[239])*states[58]+algebraic[183]*states[55])+algebraic[202]*states[54])+algebraic[227]*states[57] algebraic[230] = custom_piecewise([less(fabs(algebraic[227]) , 1.00000e-12), 1.00000e-12 , True, algebraic[227]]) algebraic[236] = custom_piecewise([less(fabs(algebraic[233]) , 1.00000e-12), 1.00000e-12 , True, algebraic[233]]) algebraic[242] = (constants[103]*((algebraic[224]*algebraic[239])*algebraic[183]))/((algebraic[236]*algebraic[230])*constants[102]) rates[52] = ((-((algebraic[202]+algebraic[242])+algebraic[224])*states[52]+algebraic[233]*states[51])+constants[103]*states[53])+algebraic[205]*states[56] rates[53] = ((-((algebraic[202]+constants[103])+algebraic[227])*states[53]+algebraic[239]*states[54])+algebraic[242]*states[52])+algebraic[205]*states[57] rates[56] = ((-((algebraic[205]+algebraic[242])+algebraic[224])*states[56]+algebraic[233]*states[55])+constants[103]*states[57])+algebraic[202]*states[52] rates[57] = ((-((algebraic[205]+constants[103])+algebraic[227])*states[57]+algebraic[239]*states[58])+algebraic[242]*states[56])+algebraic[202]*states[53] algebraic[62] = constants[274]-states[33] algebraic[74] = constants[262]-states[33]*constants[261] algebraic[90] = -constants[230]*states[33] algebraic[102] = -algebraic[62]/3.00000+((2.00000/3.00000)*(power(algebraic[62]*algebraic[62]-3.00000*algebraic[74], 1.0/2)))*cos(arccos((((9.00000*algebraic[62])*algebraic[74]-((2.00000*algebraic[62])*algebraic[62])*algebraic[62])-27.0000*algebraic[90])/(2.00000*(power(algebraic[62]*algebraic[62]-3.00000*algebraic[74], 1.50000))))/3.00000) algebraic[209] = (algebraic[149]-algebraic[102])/constants[95] algebraic[10] = (states[4]+constants[392])/constants[17] algebraic[220] = (algebraic[10]-constants[397])/(0.927300-constants[397]) algebraic[223] = custom_piecewise([less(algebraic[220] , 0.00000), 0.00000 , True, algebraic[220]]) algebraic[226] = 0.000257900*(1.00000+0.100000*states[36]) algebraic[229] = exp((2.00000*states[50])*constants[263]) algebraic[232] = ((((algebraic[226]*4.00000)*states[50])*constants[275])*(algebraic[149]*algebraic[229]-0.341000*constants[97]))/(algebraic[229]-1.00000) algebraic[235] = algebraic[232]*(states[55]+states[58]) algebraic[238] = 0.000155200*(1.00000+0.400000*states[36]) algebraic[241] = exp((2.00000*states[50])*constants[263]) algebraic[244] = ((((algebraic[238]*4.00000)*states[50])*constants[275])*(algebraic[149]*algebraic[241]-0.341000*constants[97]))/(algebraic[241]-1.00000) algebraic[246] = algebraic[244]*(states[63]+states[66]) algebraic[248] = (1.00000-algebraic[223])*algebraic[246]+algebraic[223]*algebraic[235] rates[33] = -constants[305]*algebraic[248]+constants[315]*algebraic[209] algebraic[191] = 1.00000/((1.00000+exp((13.5600-states[50])/9.45000))*(1.00000+exp((25.0000+states[50])/-5.00000))) algebraic[208] = algebraic[191]/algebraic[195] algebraic[210] = (1.00000-algebraic[191])/algebraic[195] algebraic[234] = 1.00000+power(0.00110000/algebraic[149], 4.00000) algebraic[237] = 14.9186/algebraic[234] algebraic[219] = 1.00000/(1.00000+exp((17.5000+states[50])/3.00000)) algebraic[225] = (0.247400+algebraic[219])/1.24740 algebraic[228] = 1.00000/(algebraic[214]+algebraic[216]/26.5530) algebraic[245] = algebraic[225]/algebraic[228] algebraic[251] = (1.00000-algebraic[225])/algebraic[228] rates[59] = ((-((algebraic[208]+algebraic[237])+algebraic[251])*states[59]+algebraic[210]*states[63])+constants[104]*states[62])+algebraic[245]*states[60] rates[63] = ((-((algebraic[210]+algebraic[237])+algebraic[251])*states[63]+algebraic[208]*states[59])+constants[104]*states[66])+algebraic[245]*states[64] algebraic[222] = (0.00100000+algebraic[219])/1.00100 algebraic[212] = 5.00000*states[36] algebraic[231] = 1.00000+power(0.0120000/algebraic[149], 10.0000) algebraic[240] = (13.8250-(6.38360-algebraic[212])/algebraic[234])-3.36960/algebraic[231] algebraic[243] = 1.00000/(algebraic[214]+algebraic[216]/algebraic[240]) algebraic[247] = algebraic[222]/algebraic[243] algebraic[257] = (1.00000-algebraic[222])/algebraic[243] rates[62] = ((-((algebraic[208]+constants[104])+algebraic[257])*states[62]+algebraic[237]*states[59])+algebraic[210]*states[66])+algebraic[247]*states[61] rates[66] = ((-((algebraic[210]+constants[104])+algebraic[257])*states[66]+algebraic[237]*states[63])+algebraic[208]*states[62])+algebraic[247]*states[65] algebraic[249] = custom_piecewise([less(fabs(algebraic[247]) , 1.00000e-12), 1.00000e-12 , True, algebraic[247]]) algebraic[254] = custom_piecewise([less(fabs(algebraic[251]) , 1.00000e-12), 1.00000e-12 , True, algebraic[251]]) algebraic[260] = (constants[105]*((algebraic[245]*algebraic[257])*algebraic[237]))/((algebraic[254]*algebraic[249])*constants[104]) rates[60] = ((-((algebraic[208]+algebraic[260])+algebraic[245])*states[60]+algebraic[251]*states[59])+constants[105]*states[61])+algebraic[210]*states[64] rates[61] = ((-((algebraic[208]+constants[105])+algebraic[247])*states[61]+algebraic[257]*states[62])+algebraic[260]*states[60])+algebraic[210]*states[65] rates[64] = ((-((algebraic[210]+algebraic[260])+algebraic[245])*states[64]+algebraic[251]*states[63])+constants[105]*states[65])+algebraic[208]*states[60] rates[65] = ((-((algebraic[210]+constants[105])+algebraic[247])*states[65]+algebraic[257]*states[66])+algebraic[260]*states[64])+algebraic[208]*states[61] algebraic[252] = constants[243]*(1.00000+2.00000*states[41]) algebraic[255] = 0.112500*algebraic[252] algebraic[113] = (constants[72]+constants[73])-states[34] algebraic[121] = states[34]*constants[73] algebraic[127] = (power(algebraic[113]*algebraic[113]+4.00000*algebraic[121], 1.0/2)-algebraic[113])/2.00000 algebraic[269] = (algebraic[248]*1.00000)/(1.00000+power(1.00000/algebraic[127], 8.00000)) algebraic[271] = algebraic[255]*algebraic[269] algebraic[263] = 1.00000+0.0123000/algebraic[127] algebraic[265] = algebraic[252]/algebraic[263] rates[118] = -(algebraic[271]+states[118])/algebraic[265] algebraic[258] = constants[243]*(1.00000+0.00000*states[41]) algebraic[261] = 0.112500*algebraic[258] algebraic[272] = (1.99250*algebraic[261])*algebraic[269] algebraic[267] = (0.535700*algebraic[258])/algebraic[263] rates[119] = -(algebraic[272]+states[119])/algebraic[267] algebraic[286] = 1.00000+power(constants[134]/algebraic[149], 2.00000) algebraic[285] = exp(((constants[139]-1.00000)*states[50])*constants[263]) algebraic[287] = 1.00000+constants[140]*algebraic[285] algebraic[283] = power(states[49], 3.00000) algebraic[288] = (constants[136]*algebraic[283]+constants[364]*algebraic[149])+(constants[362]*constants[97])*(1.00000+algebraic[149]/constants[135]) algebraic[289] = ((constants[135]*constants[366])*(1.00000+algebraic[283]/constants[362])+algebraic[283]*constants[97])+constants[366]*algebraic[149] algebraic[284] = exp((constants[139]*states[50])*constants[263]) algebraic[290] = (0.200000*constants[141])*((algebraic[283]*constants[97])*algebraic[284]-(constants[366]*algebraic[149])*algebraic[285]) algebraic[291] = algebraic[290]/((algebraic[286]*algebraic[287])*(algebraic[288]+algebraic[289])) algebraic[207] = (states[49]-states[48])/constants[94] rates[49] = -(constants[409]*algebraic[291]+algebraic[207]) algebraic[38] = (states[5]+constants[395])/constants[30] algebraic[311] = (algebraic[38]-constants[398])/(0.958600-constants[398]) algebraic[313] = custom_piecewise([less(algebraic[311] , 0.00000), 0.00000 , True, algebraic[311]]) algebraic[317] = (1.00000-algebraic[313])*states[118]+algebraic[313]*states[119] algebraic[307] = (constants[155]*exp(algebraic[127]/constants[153]))*(algebraic[127]-algebraic[149]) algebraic[309] = (constants[156]*exp(algebraic[127]/constants[154]))*(algebraic[127]-algebraic[149]) algebraic[315] = (1.00000-algebraic[313])*algebraic[307]+algebraic[313]*algebraic[309] algebraic[319] = algebraic[315]+algebraic[317] algebraic[211] = (states[30]-algebraic[127])/constants[96] rates[34] = algebraic[211]-algebraic[319] algebraic[153] = (states[31]-0.673519)/(0.999180-0.673519) algebraic[158] = custom_piecewise([less(algebraic[153] , 0.00000), 0.00000 , True, algebraic[153]]) algebraic[162] = (1.00000-algebraic[158])*constants[77]+algebraic[158]*constants[229] algebraic[174] = algebraic[162]+constants[75] algebraic[178] = ((algebraic[174]-states[32])+constants[71])+constants[78] algebraic[166] = algebraic[162]*constants[75] algebraic[182] = ((algebraic[166]-states[32]*algebraic[174])+constants[78]*constants[75])+constants[71]*algebraic[162] algebraic[170] = -algebraic[166]*states[32] algebraic[186] = -algebraic[178]/3.00000+((2.00000/3.00000)*(power(algebraic[178]*algebraic[178]-3.00000*algebraic[182], 1.0/2)))*cos(arccos((((9.00000*algebraic[178])*algebraic[182]-((2.00000*algebraic[178])*algebraic[178])*algebraic[178])-27.0000*algebraic[170])/(2.00000*(power(algebraic[178]*algebraic[178]-3.00000*algebraic[182], 1.50000))))/3.00000) algebraic[201] = (algebraic[149]-algebraic[186])/constants[94] rates[35] = -(((constants[316]*algebraic[291]+constants[317]*algebraic[319])+algebraic[201])+algebraic[209]) algebraic[250] = exp(states[50]*constants[263]) algebraic[253] = (((constants[107]*states[50])*constants[275])*(states[46]-constants[98]*algebraic[250]))/(1.00000-algebraic[250]) algebraic[320] = 1.00000-1.00000/(1.00000+power(algebraic[317]/constants[108], 2.00000)) algebraic[323] = (algebraic[253]*algebraic[320])*states[67] algebraic[204] = (states[47]-states[46])/constants[94] rates[47] = constants[341]*algebraic[323]-algebraic[204] algebraic[306] = 1.00000/(1.00000+power(0.0300000/algebraic[194], 2.00000)) algebraic[308] = (1.00000-algebraic[306])*constants[161]+algebraic[306]*constants[379] algebraic[324] = (algebraic[308]*states[30])/constants[164] algebraic[310] = (states[125]-0.659100)/(0.994500-0.659100) algebraic[312] = custom_piecewise([less(algebraic[310] , 0.00000), 0.00000 , True, algebraic[310]]) algebraic[314] = algebraic[312]*states[40] algebraic[316] = states[40]-algebraic[314] algebraic[318] = algebraic[312]-algebraic[314] algebraic[321] = ((1.00000-algebraic[318])-algebraic[316])-algebraic[314] algebraic[322] = ((algebraic[321]*constants[159]+algebraic[318]*constants[374])+algebraic[316]*constants[369])+algebraic[314]*constants[376] algebraic[326] = (algebraic[308]*algebraic[186])/(algebraic[186]+algebraic[322]) algebraic[328] = algebraic[326]-algebraic[324] rates[30] = algebraic[328]-constants[312]*algebraic[211] algebraic[213] = (2.00000*states[50])*constants[263] algebraic[215] = exp(algebraic[213]) algebraic[217] = ((((constants[101]*2.00000)*constants[87])*algebraic[213])*(algebraic[186]*algebraic[215]-0.341000*constants[97]))/(algebraic[215]-1.00000) algebraic[292] = 1.00000+power(constants[134]/algebraic[186], 2.00000) algebraic[293] = 1.00000+constants[140]*algebraic[285] algebraic[282] = power(states[48], 3.00000) algebraic[294] = (constants[136]*algebraic[282]+constants[364]*algebraic[186])+(constants[362]*constants[97])*(1.00000+algebraic[186]/constants[135]) algebraic[295] = ((constants[135]*constants[366])*(1.00000+algebraic[282]/constants[362])+algebraic[282]*constants[97])+constants[366]*algebraic[186] algebraic[296] = (0.800000*constants[141])*((algebraic[282]*constants[97])*algebraic[284]-(constants[366]*algebraic[186])*algebraic[285]) algebraic[297] = algebraic[296]/((algebraic[292]*algebraic[293])*(algebraic[294]+algebraic[295])) algebraic[304] = (constants[151]*algebraic[186])/(constants[152]+algebraic[186]) rates[32] = (-constants[308]*((algebraic[217]+algebraic[304])-2.00000*algebraic[297])-constants[313]*algebraic[328])+constants[318]*algebraic[201] algebraic[327] = -constants[238]*log(constants[98]/states[46]) algebraic[331] = constants[238]*log(constants[99]/states[126]) algebraic[332] = algebraic[331]-algebraic[327] algebraic[333] = (constants[91]*algebraic[332])/(algebraic[332]+constants[92]) algebraic[352] = constants[238]*log(constants[100]/states[48]) algebraic[353] = power(algebraic[352]-algebraic[327], 4.00000) algebraic[354] = (constants[93]*algebraic[353])/(algebraic[353]+constants[346]) algebraic[329] = constants[106]*(states[50]-algebraic[327]) rates[46] = ((constants[332]*algebraic[329]+algebraic[354])+algebraic[333])+constants[337]*algebraic[204] algebraic[355] = (((constants[241]*(power(states[114], 3.00000)))*states[106])*states[107])*(states[50]-algebraic[352]) algebraic[356] = (((constants[241]*(power(states[108], 3.00000)))*states[106])*states[107])*(states[50]-algebraic[352]) algebraic[357] = (((constants[241]*(power(states[111], 3.00000)))*states[109])*states[110])*(states[50]-algebraic[352]) algebraic[358] = ((((constants[241]*(power(states[114], 3.00000)))*states[112])*states[113])*(states[50]-algebraic[352]))*1.25000 algebraic[273] = (states[105]-0.239480)/(0.950143-0.239480) algebraic[274] = custom_piecewise([less(algebraic[273] , 0.00000), 0.00000 , True, algebraic[273]]) algebraic[275] = algebraic[274]*states[38] algebraic[276] = states[38]-algebraic[275] algebraic[277] = algebraic[274]-algebraic[275] algebraic[278] = ((1.00000-algebraic[277])-algebraic[276])-algebraic[275] algebraic[359] = ((algebraic[278]*algebraic[357]+algebraic[277]*algebraic[358])+algebraic[276]*algebraic[356])+algebraic[275]*algebraic[355] algebraic[279] = states[50]*constants[263] algebraic[280] = exp(algebraic[279]) algebraic[281] = (((constants[133]*constants[87])*algebraic[279])*(states[48]*algebraic[280]-constants[100]))/(algebraic[280]-1.00000) algebraic[300] = (constants[144]*constants[242])/(1.00000+exp(-(states[50]+92.0000)*constants[263])) algebraic[301] = algebraic[300]*(power(states[48]/(states[48]+constants[147]), 3.00000)) algebraic[302] = algebraic[300]*(power(states[48]/(states[48]+constants[148]), 3.00000)) algebraic[298] = (states[115]-0.126345)/(0.998014-0.126345) algebraic[299] = custom_piecewise([less(algebraic[298] , 0.00000), 0.00000 , True, algebraic[298]]) algebraic[303] = (1.00000-algebraic[299])*algebraic[301]+algebraic[299]*algebraic[302] algebraic[360] = ((power(states[117], 3.00000))*states[116])*(states[50]-algebraic[352]) algebraic[361] = 0.0160000*algebraic[360] algebraic[362] = 0.00650000*algebraic[360] algebraic[363] = (1.00000-states[38])*algebraic[362]+states[38]*algebraic[361] algebraic[365] = (((algebraic[359]+algebraic[281])+algebraic[363])+algebraic[303]*3.00000)+algebraic[297]*3.00000 rates[48] = (constants[407]*algebraic[365]+constants[408]*algebraic[207])+algebraic[354] algebraic[305] = (((algebraic[248]+algebraic[217])+algebraic[304])-algebraic[297]*2.00000)-algebraic[291]*2.00000 algebraic[330] = algebraic[329]+algebraic[323] algebraic[334] = states[50]-algebraic[331] algebraic[335] = 1.02000/(1.00000+exp(0.238500*(algebraic[334]-59.2150))) algebraic[336] = (0.491240*exp(0.0803200*(algebraic[334]+5.47600))+exp(0.0617500*(algebraic[334]-594.310)))/(1.00000+exp(-0.514300*(algebraic[334]+4.75300))) algebraic[337] = (constants[239]*(algebraic[335]/(algebraic[335]+algebraic[336])))*algebraic[334] algebraic[338] = algebraic[337]*1.20000 algebraic[339] = (1.00000-states[37])*algebraic[337]+states[37]*algebraic[338] algebraic[256] = 1.00000/(1.00000+exp((states[50]+10.0000)/15.4000)) algebraic[340] = ((constants[240]*states[68])*algebraic[256])*(states[50]-algebraic[331]) algebraic[348] = constants[238]*log((constants[99]+constants[165]*constants[100])/(states[126]+constants[165]*states[48])) algebraic[259] = 0.195610*(1.00000+0.600000/(1.00000+power(3.80000e-05/algebraic[186], 1.40000))) algebraic[349] = (algebraic[259]*(states[69]+states[70]))*(states[50]-algebraic[348]) algebraic[350] = (algebraic[259]*(states[71]+states[72]))*(states[50]-algebraic[348]) algebraic[262] = (states[103]+constants[403])/constants[110] algebraic[264] = (algebraic[262]-constants[404])/(0.785000-constants[404]) algebraic[266] = custom_piecewise([less(algebraic[264] , 0.00000), 0.00000 , True, algebraic[264]]) algebraic[351] = algebraic[266]*algebraic[350]+(1.00000-algebraic[266])*algebraic[349] algebraic[341] = (constants[123]*(states[50]-algebraic[331]))/(1.00000+exp((15.0000-states[50])/17.0000)) algebraic[342] = ((constants[123]*(states[50]-algebraic[331]))/(1.00000+exp((36.0000-states[50])/17.0000)))*3.62000 algebraic[268] = (states[104]-0.0589380)/(0.393747-0.0589380) algebraic[270] = custom_piecewise([less(algebraic[268] , 0.00000), 0.00000 , True, algebraic[268]]) algebraic[343] = (1.00000-algebraic[270])*algebraic[341]+algebraic[270]*algebraic[342] algebraic[325] = exp(states[50]/550.000) algebraic[344] = ((constants[157]*(power(states[120], 3.00000)))*algebraic[325])*(states[50]-algebraic[331]) algebraic[345] = algebraic[344]*(0.735600*states[121]+0.264400*states[123]) algebraic[346] = algebraic[344]*(0.735600*states[122]+0.264400*states[124]) algebraic[347] = (1.00000-states[39])*algebraic[346]+states[39]*algebraic[345] algebraic[364] = ((((algebraic[339]+algebraic[340])+algebraic[351])+algebraic[343])+algebraic[347])-2.00000*algebraic[303] algebraic[366] = algebraic[365]+algebraic[291]*3.00000 algebraic[367] = ((algebraic[366]+algebraic[364])+algebraic[305])+algebraic[330] algebraic[368] = (custom_piecewise([less((voi-constants[213])-constants[214]*floor((voi-constants[213])/constants[214]) , constants[212]), 1.00000 , True, 0.00000]))*constants[256] rates[50] = -(algebraic[367]+algebraic[368])/1.00000 rates[126] = constants[406]*(algebraic[364]+algebraic[368])+algebraic[333] return(rates) def computeAlgebraic(constants, states, voi): algebraic = array([[0.0] * len(voi)] * sizeAlgebraic) states = array(states) voi = array(voi) algebraic[3] = constants[393]-states[4] algebraic[4] = constants[396]-states[5] algebraic[15] = constants[405]-states[103] algebraic[25] = 1.00000/(1.00000+exp((states[50]+91.0000)/6.10000)) algebraic[29] = constants[255]-states[144] algebraic[33] = (constants[322]-states[12])-states[11] algebraic[35] = (constants[325]-states[18])-states[17] algebraic[36] = (constants[330]-states[26])-states[25] algebraic[11] = 0.0250000/(1.00000+exp((states[50]+58.0000)/5.00000)) algebraic[39] = 0.200000/(1.00000+exp((states[50]+19.0000)/-9.00000)) algebraic[12] = 1.00000/((0.000600000*(states[50]-1.73840))/(1.00000-exp(-0.136000*(states[50]-1.73840)))-(0.000300000*(states[50]+38.3608))/(1.00000-exp(0.152200*(states[50]+38.3608)))) algebraic[40] = 1.00000/(1.00000+exp((states[50]+10.0850)/-4.25000)) algebraic[13] = 0.00739900/(1.00000+exp((constants[263]*(states[50]-0.0311960))/-0.800190)) algebraic[41] = 0.00569920/(1.00000+exp((constants[263]*(states[50]-0.0415200))/1.34890)) algebraic[14] = 0.00994150/(1.00000+exp((constants[263]*(states[50]-0.0448090))/-0.581720)) algebraic[42] = 0.00332010/(1.00000+exp((constants[263]*(states[50]-0.0942170))/0.953640)) algebraic[16] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), 0.00000 , True, 0.135000*exp(((87.0000+states[50])+constants[130])/-6.80000)]) algebraic[43] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), 1.00000/(0.130000*(1.00000+exp(((states[50]+constants[130])+27.4034)/-11.1000))) , True, 3.56000*exp(0.0790000*((states[50]+constants[130])+7.00000))+310000.*exp(0.350000*((states[50]+constants[130])+7.00000))]) algebraic[17] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), 0.00000 , True, ((-127140.*exp(0.244400*(states[50]+constants[130]))+-6.94800e-05*exp(-0.0439100*(states[50]+constants[130])))*((states[50]+constants[130])+37.7800))/(1.00000+exp(0.311000*((states[50]+constants[130])+79.2300)))]) algebraic[44] = custom_piecewise([greater_equal(states[50]+constants[130] , -40.0000), (0.300000*exp(-2.53500e-07*(states[50]+constants[130])))/(1.00000+exp(-0.100000*((states[50]+constants[130])+32.0000))) , True, (0.121200*exp(-0.0105200*(states[50]+constants[130])))/(1.00000+exp(-0.137800*((states[50]+constants[130])+40.1400)))]) algebraic[18] = (0.320000*(states[50]+58.4729))/(1.00000-exp(-0.100000*(states[50]+58.4729))) algebraic[45] = 0.0800000*exp((13.7299-states[50])/11.0000) algebraic[19] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, 0.135000*exp((87.0000+states[50])/-6.80000)]) algebraic[46] = custom_piecewise([greater_equal(states[50] , -40.0000), 1.00000/(0.130000*(1.00000+exp((states[50]+27.4034)/-11.1000))) , True, 3.56000*exp(0.0790000*(states[50]+7.00000))+310000.*exp(0.350000*(states[50]+7.00000))]) algebraic[20] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, ((-127140.*exp(0.244400*states[50])-6.94800e-05*exp(-0.0439100*states[50]))*(states[50]+37.7800))/(1.00000+exp(0.311000*(states[50]+79.2300)))]) algebraic[47] = custom_piecewise([greater_equal(states[50] , -40.0000), (0.300000*exp(-2.53500e-07*states[50]))/(1.00000+exp(-0.100000*(states[50]+32.0000))) , True, (0.121200*exp(-0.0105200*states[50]))/(1.00000+exp(-0.137800*(states[50]+40.1400)))]) algebraic[21] = (0.320000*(states[50]+58.4729))/(1.00000-exp(-0.100000*(states[50]+58.4729))) algebraic[48] = 0.0800000*exp((13.7299-states[50])/11.0000) algebraic[22] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, 0.135000*exp(((87.0000+states[50])+constants[132])/-6.80000)]) algebraic[49] = custom_piecewise([greater_equal(states[50] , -40.0000), (1.00000/0.130000)/(1.00000+exp(((states[50]+constants[132])+27.4034)/-11.1000)) , True, 3.56000*exp(0.0790000*((states[50]+constants[132])+7.00000))+310000.*exp(0.350000*((states[50]+constants[132])+7.00000))]) algebraic[23] = custom_piecewise([greater_equal(states[50] , -40.0000), 0.00000 , True, ((-127140.*exp(0.244400*(states[50]+constants[132]))-6.94800e-05*exp(-0.0439100*(states[50]+constants[132])))*((states[50]+constants[132])+37.7800))/(1.00000+exp(0.311000*((states[50]+constants[132])+79.2300)))]) algebraic[50] = custom_piecewise([greater_equal(states[50] , -40.0000), (0.300000*exp(-2.53500e-07*(states[50]+constants[132])))/(1.00000+exp(-0.100000*((states[50]+constants[132])+32.0000))) , True, (0.121200*exp(-0.0105200*(states[50]+constants[132])))/(1.00000+exp(-0.137800*((states[50]+constants[132])+40.1400)))]) algebraic[24] = (0.320000*((states[50]+constants[131])+58.4729))/(1.00000-exp(-0.100000*((states[50]+constants[131])+58.4729))) algebraic[51] = 0.0800000*exp(((states[50]+constants[131])-13.7299)/-11.0000) algebraic[26] = (0.320000*(states[50]+47.1300))/(1.00000-exp(-0.100000*(states[50]+47.1300))) algebraic[52] = 0.0800000*exp(states[50]/-11.0000) algebraic[57] = (constants[323]-states[14])-states[13] algebraic[60] = (constants[331]-states[28])-states[27] algebraic[27] = 1.00000/(1.00000+exp((states[50]+9.43700)/-7.13300)) algebraic[53] = 1.00000/(1.00000+exp((states[50]-18.4099)/-29.3814)) algebraic[65] = 1.00000/(1.00000+exp((states[50]+100.000)/29.3814)) algebraic[77] = 1.00000/(algebraic[53]/1.20890+3.50000*algebraic[65]) algebraic[66] = ((1.00000/(1.00000+exp((states[50]+19.0000)/-9.00000)))/0.500000)/9.79530 algebraic[28] = 1.00000/(1.00000+exp((states[50]+58.0000)/5.00000)) algebraic[78] = 0.0479600*algebraic[28] algebraic[79] = 0.0214400*algebraic[28] algebraic[54] = (1.00000/(1.00000+exp((states[50]+60.0000)/5.00000)))/250.000 algebraic[80] = 2.46000*algebraic[54] algebraic[81] = 0.560340*algebraic[54] algebraic[63] = 0.0906540*exp((-0.111570*states[50])*constants[263]) algebraic[91] = 0.388390/(1.00000+exp((constants[263]*(states[50]+0.150190))/-0.606930)) algebraic[75] = 0.00311240+(0.0283300-0.00311240)/(1.00000+exp((constants[263]*(states[50]+0.0516600))/1.55220)) algebraic[64] = 0.0657000*exp((-0.118990*states[50])*constants[263]) algebraic[93] = 0.563560/(1.00000+exp((constants[263]*(states[50]+0.179860))/-0.583810)) algebraic[76] = 0.000385250+(0.0124060-0.000385250)/(1.00000+exp((constants[263]*(states[50]+0.0641180))/0.779920)) algebraic[92] = 0.000441980*exp((-1.20220*states[50])*constants[263]) algebraic[103] = 0.000401730*exp((0.000208730*states[50])*constants[263]) algebraic[94] = 0.000237300*exp((-1.97420*states[50])*constants[263]) algebraic[104] = 0.000226520*exp((0.000246900*states[50])*constants[263]) algebraic[5] = ((constants[36]*constants[35])*constants[319]-states[8])-states[7] algebraic[34] = (algebraic[5]*(constants[59]+constants[51])-states[14]*(constants[51]+constants[59]))+(constants[49]*constants[51])*(1.00000+constants[59]/constants[54]) algebraic[58] = (-states[14]*constants[49])*constants[51] algebraic[70] = (-algebraic[34]+power(algebraic[34]*algebraic[34]-(4.00000*constants[222])*algebraic[58], 1.0/2))/(2.00000*constants[222]) algebraic[86] = algebraic[5]/(1.00000+(algebraic[70]/constants[49])*(1.00000+constants[59]/constants[51])) algebraic[98] = (algebraic[70]*algebraic[86])/constants[49] algebraic[108] = (algebraic[98]*constants[59])/constants[51] algebraic[8] = ((constants[215]*constants[35])*constants[326]-states[21])-states[20] algebraic[37] = (algebraic[8]*(constants[59]+constants[51])-states[28]*(constants[51]+constants[59]))+(constants[49]*constants[51])*(1.00000+constants[59]/constants[54]) algebraic[61] = (-states[28]*constants[49])*constants[51] algebraic[73] = (-algebraic[37]+power(algebraic[37]*algebraic[37]-(4.00000*constants[225])*algebraic[61], 1.0/2))/(2.00000*constants[225]) algebraic[89] = algebraic[8]/(1.00000+(algebraic[73]/constants[49])*(1.00000+constants[59]/constants[51])) algebraic[101] = (algebraic[73]*algebraic[89])/constants[49] algebraic[111] = (constants[59]/constants[51])*algebraic[101] algebraic[7] = ((constants[216]*constants[295])*constants[321]-states[2])-states[15] algebraic[59] = (algebraic[7]*(constants[47]+constants[59])-algebraic[35]*(constants[47]+constants[59]))+(constants[46]*constants[47])*(1.00000+constants[59]/constants[48]) algebraic[71] = (-algebraic[35]*constants[47])*constants[46] algebraic[87] = (-algebraic[59]+power(algebraic[59]*algebraic[59]-(4.00000*constants[224])*algebraic[71], 1.0/2))/(2.00000*constants[224]) algebraic[99] = algebraic[7]/(1.00000+(algebraic[87]/constants[46])*(1.00000+constants[59]/constants[47])) algebraic[110] = ((constants[59]*algebraic[87])*algebraic[99])/(constants[46]*constants[47]) algebraic[119] = (algebraic[99]*algebraic[87])/constants[46] algebraic[118] = (constants[59]*algebraic[87])/constants[48] algebraic[130] = ((constants[399]-states[142])-states[141])-states[140] algebraic[1] = power(states[3]*1.00000, constants[12]) algebraic[31] = constants[6]*(constants[8]+algebraic[1]/(constants[3]+algebraic[1])) algebraic[67] = (algebraic[31]*constants[329])*constants[227] algebraic[83] = (constants[236]*(states[43]-states[45]))/constants[309] algebraic[96] = (constants[237]*(states[45]-states[44]))/constants[309] algebraic[107] = (constants[349]*constants[185])/(1.00000+constants[170]/states[45]) algebraic[116] = ((constants[372]+(constants[174]-1.00000)*states[131])*constants[187])/(1.00000+constants[172]/states[45]) algebraic[124] = algebraic[107]+algebraic[116] algebraic[134] = (((-constants[268]*algebraic[130])*states[45]+constants[276]*states[142])-(constants[269]*states[142])*states[45])+constants[277]*states[141] algebraic[140] = ((constants[380]-states[134])-states[133])-states[132] algebraic[2] = power(states[1]*1.00000, constants[13]) algebraic[32] = power(states[1]*1.00000, constants[14]) algebraic[56] = (constants[7]*(constants[9]+algebraic[2]/(constants[4]+algebraic[2])))*(1.00000-((1.00000-(constants[15]*algebraic[32])/(constants[5]+algebraic[32]))*states[0])/(constants[2]+states[0])) algebraic[82] = (algebraic[56]*constants[320])*constants[227] algebraic[95] = (constants[236]*(states[43]-states[45]))/constants[306] algebraic[105] = (constants[235]*(states[43]-states[44]))/constants[306] algebraic[114] = (constants[347]*constants[185])/(1.00000+constants[170]/states[43]) algebraic[122] = ((constants[377]+(constants[174]-1.00000)*states[127])*constants[186])/(1.00000+constants[171]/states[43]) algebraic[128] = ((constants[370]+(constants[174]-1.00000)*states[129])*constants[187])/(1.00000+constants[172]/states[43]) algebraic[135] = (algebraic[114]+algebraic[122])+algebraic[128] algebraic[145] = (((-constants[199]*algebraic[140])*states[43]+constants[246]*states[134])-(constants[200]*states[134])*states[43])+constants[247]*states[133] algebraic[150] = ((constants[335]-states[138])-states[137])-states[136] algebraic[0] = power(states[2]*1.00000, constants[12]) algebraic[55] = constants[6]*(constants[8]+algebraic[0]/(constants[3]+algebraic[0])) algebraic[68] = (algebraic[55]*constants[324])*constants[227] algebraic[30] = power(states[2]*1.00000, constants[13]) algebraic[84] = constants[7]*(constants[9]+algebraic[30]/(constants[4]+algebraic[30])) algebraic[106] = (algebraic[84]*constants[328])*constants[227] algebraic[115] = (constants[235]*(states[43]-states[44]))/constants[307] algebraic[123] = (constants[237]*(states[45]-states[44]))/constants[307] algebraic[129] = (constants[348]*constants[185])/(1.00000+constants[170]/states[44]) algebraic[136] = ((constants[378]+(constants[174]-1.00000)*states[128])*constants[186])/(1.00000+constants[171]/states[44]) algebraic[141] = ((constants[371]+(constants[174]-1.00000)*states[130])*constants[187])/(1.00000+constants[172]/states[44]) algebraic[146] = (algebraic[129]+algebraic[136])+algebraic[141] algebraic[154] = (((-constants[249]*algebraic[150])*states[44]+constants[265]*states[138])-(constants[250]*states[138])*states[44])+constants[266]*states[137] algebraic[6] = ((constants[37]*constants[295])*constants[319]-states[1])-states[9] algebraic[69] = (((constants[258]+constants[271])/constants[223]+algebraic[33])+algebraic[57])-algebraic[6] algebraic[85] = ((constants[271]*(algebraic[33]-algebraic[6])+constants[258]*(algebraic[57]-algebraic[6]))+constants[279])/constants[223] algebraic[97] = (algebraic[6]*constants[279])/constants[223] algebraic[109] = ((((-algebraic[97]/27.0000)*(power(algebraic[69], 3.00000))-(((algebraic[69]*algebraic[69])*algebraic[85])*algebraic[85])/108.000)+((algebraic[69]*algebraic[85])*algebraic[97])/6.00000)+(power(algebraic[85], 3.00000))/27.0000)+(algebraic[97]*algebraic[97])/4.00000 algebraic[117] = custom_piecewise([less(algebraic[109] , 0.00000), power(-algebraic[109], 1.0/2) , True, 0.00000]) algebraic[125] = (((custom_piecewise([greater(algebraic[109] , 0.00000), power(algebraic[109], 1.0/2) , True, 0.00000]))+algebraic[97]/2.00000)+(algebraic[69]*algebraic[85])/6.00000)-(power(algebraic[69], 3.00000))/27.0000 algebraic[131] = arctan(algebraic[117]/algebraic[125])/3.00000 algebraic[137] = power(algebraic[125]*algebraic[125]+algebraic[117]*algebraic[117], 1.00000/6.00000) algebraic[142] = (algebraic[85]/3.00000-(algebraic[69]*algebraic[69])/9.00000)/(algebraic[137]*algebraic[137]) algebraic[147] = (algebraic[137]*sin(algebraic[131]))*(1.00000+algebraic[142]) algebraic[151] = (algebraic[137]*cos(algebraic[131]))*(1.00000-algebraic[142])-algebraic[69]/3.00000 algebraic[156] = power(algebraic[151]*algebraic[151]+algebraic[147]*algebraic[147], 1.0/2) algebraic[160] = algebraic[33]/((1.00000+constants[59]/constants[48])+(algebraic[156]*(constants[47]+constants[59]))/(constants[46]*constants[47])) algebraic[164] = (constants[59]*algebraic[160])/constants[48] algebraic[168] = ((constants[59]*algebraic[160])*algebraic[156])/(constants[46]*constants[47]) algebraic[9] = ((constants[38]*constants[295])*constants[326]-states[3])-states[23] algebraic[72] = (((constants[259]+constants[272])/constants[226]+algebraic[36])+algebraic[60])-algebraic[9] algebraic[88] = ((constants[272]*(algebraic[36]-algebraic[9])+constants[259]*(algebraic[60]-algebraic[9]))+constants[280])/constants[226] algebraic[100] = (algebraic[9]*constants[280])/constants[226] algebraic[112] = ((((-algebraic[100]/27.0000)*(power(algebraic[72], 3.00000))-(((algebraic[72]*algebraic[72])*algebraic[88])*algebraic[88])/108.000)+((algebraic[72]*algebraic[88])*algebraic[100])/6.00000)+(power(algebraic[88], 3.00000))/27.0000)+(algebraic[100]*algebraic[100])/4.00000 algebraic[120] = custom_piecewise([less(algebraic[112] , 0.00000), power(-algebraic[112], 1.0/2) , True, 0.00000]) algebraic[126] = (((custom_piecewise([greater(algebraic[112] , 0.00000), power(algebraic[112], 1.0/2) , True, 0.00000]))+algebraic[100]/2.00000)+(algebraic[72]*algebraic[88])/6.00000)-(power(algebraic[72], 3.00000))/27.0000 algebraic[132] = arctan(algebraic[120]/algebraic[126])/3.00000 algebraic[138] = power(algebraic[126]*algebraic[126]+algebraic[120]*algebraic[120], 1.00000/6.00000) algebraic[143] = (algebraic[88]/3.00000-(algebraic[72]*algebraic[72])/9.00000)/(algebraic[138]*algebraic[138]) algebraic[148] = (algebraic[138]*sin(algebraic[132]))*(1.00000+algebraic[143]) algebraic[152] = (algebraic[138]*cos(algebraic[132]))*(1.00000-algebraic[143])-algebraic[72]/3.00000 algebraic[157] = power(algebraic[152]*algebraic[152]+algebraic[148]*algebraic[148], 1.0/2) algebraic[161] = algebraic[36]/((1.00000+constants[59]/constants[48])+(algebraic[157]*(constants[47]+constants[59]))/(constants[46]*constants[47])) algebraic[165] = (constants[59]*algebraic[161])/constants[48] algebraic[169] = ((constants[59]*algebraic[161])*algebraic[157])/(constants[46]*constants[47]) algebraic[176] = algebraic[57]/((1.00000+constants[59]/constants[53])+(algebraic[156]*(constants[52]+constants[59]))/(constants[50]*constants[52])) algebraic[180] = (constants[59]*algebraic[176])/constants[53] algebraic[184] = ((constants[59]*algebraic[176])*algebraic[156])/(constants[50]*constants[52]) algebraic[177] = algebraic[60]/((1.00000+constants[59]/constants[53])+(algebraic[157]*(constants[52]+constants[59]))/(constants[50]*constants[52])) algebraic[181] = (constants[59]*algebraic[177])/constants[53] algebraic[185] = ((constants[59]*algebraic[177])*algebraic[157])/(constants[50]*constants[52]) algebraic[188] = algebraic[168]+constants[60]*algebraic[184] algebraic[172] = (algebraic[160]*algebraic[156])/constants[46] algebraic[192] = (algebraic[176]*algebraic[156])/constants[50] algebraic[196] = algebraic[172]+constants[60]*algebraic[192] algebraic[189] = algebraic[169]+constants[63]*algebraic[185] algebraic[173] = (algebraic[161]*algebraic[157])/constants[46] algebraic[193] = (algebraic[177]*algebraic[157])/constants[50] algebraic[197] = algebraic[173]+constants[63]*algebraic[193] algebraic[133] = constants[274]-states[35] algebraic[139] = constants[262]-states[35]*constants[261] algebraic[144] = -constants[230]*states[35] algebraic[149] = -algebraic[133]/3.00000+((2.00000/3.00000)*(power(algebraic[133]*algebraic[133]-3.00000*algebraic[139], 1.0/2)))*cos(arccos((((9.00000*algebraic[133])*algebraic[139]-((2.00000*algebraic[133])*algebraic[133])*algebraic[133])-27.0000*algebraic[144])/(2.00000*(power(algebraic[133]*algebraic[133]-3.00000*algebraic[139], 1.50000))))/3.00000) algebraic[190] = (constants[79]*(1.00000-states[42]))/(1.00000+constants[81]/algebraic[149]) algebraic[194] = algebraic[190]+states[42] algebraic[199] = algebraic[194]/(algebraic[194]+constants[80]) algebraic[200] = 1.00000/(1.00000+power(constants[80]/algebraic[194], 2.00000)) algebraic[198] = (constants[205]-constants[208])+states[144] algebraic[203] = 0.500000*(power(power(algebraic[198], 2.00000)+(4.00000*constants[205])*constants[208], 1.0/2)-algebraic[198]) algebraic[206] = (constants[28]/constants[319]+constants[82]/constants[326])+algebraic[203]/constants[321] algebraic[155] = 1.00000/((1.00000+exp((4.79800+states[50])/-7.56990))*(1.00000+exp((25.0000+states[50])/-5.00000))) algebraic[195] = 0.590000+(0.800000*exp(0.0520000*(states[50]+13.0000)))/(1.00000+exp(0.132000*(states[50]+13.0000))) algebraic[202] = algebraic[155]/algebraic[195] algebraic[205] = (1.00000-algebraic[155])/algebraic[195] algebraic[179] = 1.00000+power(0.00200000/algebraic[149], 4.00000) algebraic[183] = 6.00000/algebraic[179] algebraic[163] = 1.00000/(1.00000+exp((29.9790+states[50])/3.17750)) algebraic[171] = (0.100000+algebraic[163])/1.10000 algebraic[214] = 1.00000/((70.0000*(1.00000-0.500000*states[36]))*(1.00000+exp((states[50]+49.1000)/10.3490))) algebraic[216] = 1.00000/(1.00000+exp((states[50]+0.213000)/-10.8070)) algebraic[221] = 1.00000/(algebraic[214]+algebraic[216]/38.4940) algebraic[224] = algebraic[171]/algebraic[221] algebraic[233] = (1.00000-algebraic[171])/algebraic[221] algebraic[167] = (0.000100000+algebraic[163])/1.00010 algebraic[159] = 0.100000*states[36] algebraic[175] = 1.00000+power(0.0100000/algebraic[149], 10.0000) algebraic[187] = (32.5000-(18.0000-algebraic[159])/algebraic[179])-10.0000/algebraic[175] algebraic[218] = 1.00000/(algebraic[214]+algebraic[216]/algebraic[187]) algebraic[227] = algebraic[167]/algebraic[218] algebraic[239] = (1.00000-algebraic[167])/algebraic[218] algebraic[230] = custom_piecewise([less(fabs(algebraic[227]) , 1.00000e-12), 1.00000e-12 , True, algebraic[227]]) algebraic[236] = custom_piecewise([less(fabs(algebraic[233]) , 1.00000e-12), 1.00000e-12 , True, algebraic[233]]) algebraic[242] = (constants[103]*((algebraic[224]*algebraic[239])*algebraic[183]))/((algebraic[236]*algebraic[230])*constants[102]) algebraic[62] = constants[274]-states[33] algebraic[74] = constants[262]-states[33]*constants[261] algebraic[90] = -constants[230]*states[33] algebraic[102] = -algebraic[62]/3.00000+((2.00000/3.00000)*(power(algebraic[62]*algebraic[62]-3.00000*algebraic[74], 1.0/2)))*cos(arccos((((9.00000*algebraic[62])*algebraic[74]-((2.00000*algebraic[62])*algebraic[62])*algebraic[62])-27.0000*algebraic[90])/(2.00000*(power(algebraic[62]*algebraic[62]-3.00000*algebraic[74], 1.50000))))/3.00000) algebraic[209] = (algebraic[149]-algebraic[102])/constants[95] algebraic[10] = (states[4]+constants[392])/constants[17] algebraic[220] = (algebraic[10]-constants[397])/(0.927300-constants[397]) algebraic[223] = custom_piecewise([less(algebraic[220] , 0.00000), 0.00000 , True, algebraic[220]]) algebraic[226] = 0.000257900*(1.00000+0.100000*states[36]) algebraic[229] = exp((2.00000*states[50])*constants[263]) algebraic[232] = ((((algebraic[226]*4.00000)*states[50])*constants[275])*(algebraic[149]*algebraic[229]-0.341000*constants[97]))/(algebraic[229]-1.00000) algebraic[235] = algebraic[232]*(states[55]+states[58]) algebraic[238] = 0.000155200*(1.00000+0.400000*states[36]) algebraic[241] = exp((2.00000*states[50])*constants[263]) algebraic[244] = ((((algebraic[238]*4.00000)*states[50])*constants[275])*(algebraic[149]*algebraic[241]-0.341000*constants[97]))/(algebraic[241]-1.00000) algebraic[246] = algebraic[244]*(states[63]+states[66]) algebraic[248] = (1.00000-algebraic[223])*algebraic[246]+algebraic[223]*algebraic[235] algebraic[191] = 1.00000/((1.00000+exp((13.5600-states[50])/9.45000))*(1.00000+exp((25.0000+states[50])/-5.00000))) algebraic[208] = algebraic[191]/algebraic[195] algebraic[210] = (1.00000-algebraic[191])/algebraic[195] algebraic[234] = 1.00000+power(0.00110000/algebraic[149], 4.00000) algebraic[237] = 14.9186/algebraic[234] algebraic[219] = 1.00000/(1.00000+exp((17.5000+states[50])/3.00000)) algebraic[225] = (0.247400+algebraic[219])/1.24740 algebraic[228] = 1.00000/(algebraic[214]+algebraic[216]/26.5530) algebraic[245] = algebraic[225]/algebraic[228] algebraic[251] = (1.00000-algebraic[225])/algebraic[228] algebraic[222] = (0.00100000+algebraic[219])/1.00100 algebraic[212] = 5.00000*states[36] algebraic[231] = 1.00000+power(0.0120000/algebraic[149], 10.0000) algebraic[240] = (13.8250-(6.38360-algebraic[212])/algebraic[234])-3.36960/algebraic[231] algebraic[243] = 1.00000/(algebraic[214]+algebraic[216]/algebraic[240]) algebraic[247] = algebraic[222]/algebraic[243] algebraic[257] = (1.00000-algebraic[222])/algebraic[243] algebraic[249] = custom_piecewise([less(fabs(algebraic[247]) , 1.00000e-12), 1.00000e-12 , True, algebraic[247]]) algebraic[254] = custom_piecewise([less(fabs(algebraic[251]) , 1.00000e-12), 1.00000e-12 , True, algebraic[251]]) algebraic[260] = (constants[105]*((algebraic[245]*algebraic[257])*algebraic[237]))/((algebraic[254]*algebraic[249])*constants[104]) algebraic[252] = constants[243]*(1.00000+2.00000*states[41]) algebraic[255] = 0.112500*algebraic[252] algebraic[113] = (constants[72]+constants[73])-states[34] algebraic[121] = states[34]*constants[73] algebraic[127] = (power(algebraic[113]*algebraic[113]+4.00000*algebraic[121], 1.0/2)-algebraic[113])/2.00000 algebraic[269] = (algebraic[248]*1.00000)/(1.00000+power(1.00000/algebraic[127], 8.00000)) algebraic[271] = algebraic[255]*algebraic[269] algebraic[263] = 1.00000+0.0123000/algebraic[127] algebraic[265] = algebraic[252]/algebraic[263] algebraic[258] = constants[243]*(1.00000+0.00000*states[41]) algebraic[261] = 0.112500*algebraic[258] algebraic[272] = (1.99250*algebraic[261])*algebraic[269] algebraic[267] = (0.535700*algebraic[258])/algebraic[263] algebraic[286] = 1.00000+power(constants[134]/algebraic[149], 2.00000) algebraic[285] = exp(((constants[139]-1.00000)*states[50])*constants[263]) algebraic[287] = 1.00000+constants[140]*algebraic[285] algebraic[283] = power(states[49], 3.00000) algebraic[288] = (constants[136]*algebraic[283]+constants[364]*algebraic[149])+(constants[362]*constants[97])*(1.00000+algebraic[149]/constants[135]) algebraic[289] = ((constants[135]*constants[366])*(1.00000+algebraic[283]/constants[362])+algebraic[283]*constants[97])+constants[366]*algebraic[149] algebraic[284] = exp((constants[139]*states[50])*constants[263]) algebraic[290] = (0.200000*constants[141])*((algebraic[283]*constants[97])*algebraic[284]-(constants[366]*algebraic[149])*algebraic[285]) algebraic[291] = algebraic[290]/((algebraic[286]*algebraic[287])*(algebraic[288]+algebraic[289])) algebraic[207] = (states[49]-states[48])/constants[94] algebraic[38] = (states[5]+constants[395])/constants[30] algebraic[311] = (algebraic[38]-constants[398])/(0.958600-constants[398]) algebraic[313] = custom_piecewise([less(algebraic[311] , 0.00000), 0.00000 , True, algebraic[311]]) algebraic[317] = (1.00000-algebraic[313])*states[118]+algebraic[313]*states[119] algebraic[307] = (constants[155]*exp(algebraic[127]/constants[153]))*(algebraic[127]-algebraic[149]) algebraic[309] = (constants[156]*exp(algebraic[127]/constants[154]))*(algebraic[127]-algebraic[149]) algebraic[315] = (1.00000-algebraic[313])*algebraic[307]+algebraic[313]*algebraic[309] algebraic[319] = algebraic[315]+algebraic[317] algebraic[211] = (states[30]-algebraic[127])/constants[96] algebraic[153] = (states[31]-0.673519)/(0.999180-0.673519) algebraic[158] = custom_piecewise([less(algebraic[153] , 0.00000), 0.00000 , True, algebraic[153]]) algebraic[162] = (1.00000-algebraic[158])*constants[77]+algebraic[158]*constants[229] algebraic[174] = algebraic[162]+constants[75] algebraic[178] = ((algebraic[174]-states[32])+constants[71])+constants[78] algebraic[166] = algebraic[162]*constants[75] algebraic[182] = ((algebraic[166]-states[32]*algebraic[174])+constants[78]*constants[75])+constants[71]*algebraic[162] algebraic[170] = -algebraic[166]*states[32] algebraic[186] = -algebraic[178]/3.00000+((2.00000/3.00000)*(power(algebraic[178]*algebraic[178]-3.00000*algebraic[182], 1.0/2)))*cos(arccos((((9.00000*algebraic[178])*algebraic[182]-((2.00000*algebraic[178])*algebraic[178])*algebraic[178])-27.0000*algebraic[170])/(2.00000*(power(algebraic[178]*algebraic[178]-3.00000*algebraic[182], 1.50000))))/3.00000) algebraic[201] = (algebraic[149]-algebraic[186])/constants[94] algebraic[250] = exp(states[50]*constants[263]) algebraic[253] = (((constants[107]*states[50])*constants[275])*(states[46]-constants[98]*algebraic[250]))/(1.00000-algebraic[250]) algebraic[320] = 1.00000-1.00000/(1.00000+power(algebraic[317]/constants[108], 2.00000)) algebraic[323] = (algebraic[253]*algebraic[320])*states[67] algebraic[204] = (states[47]-states[46])/constants[94] algebraic[306] = 1.00000/(1.00000+power(0.0300000/algebraic[194], 2.00000)) algebraic[308] = (1.00000-algebraic[306])*constants[161]+algebraic[306]*constants[379] algebraic[324] = (algebraic[308]*states[30])/constants[164] algebraic[310] = (states[125]-0.659100)/(0.994500-0.659100) algebraic[312] = custom_piecewise([less(algebraic[310] , 0.00000), 0.00000 , True, algebraic[310]]) algebraic[314] = algebraic[312]*states[40] algebraic[316] = states[40]-algebraic[314] algebraic[318] = algebraic[312]-algebraic[314] algebraic[321] = ((1.00000-algebraic[318])-algebraic[316])-algebraic[314] algebraic[322] = ((algebraic[321]*constants[159]+algebraic[318]*constants[374])+algebraic[316]*constants[369])+algebraic[314]*constants[376] algebraic[326] = (algebraic[308]*algebraic[186])/(algebraic[186]+algebraic[322]) algebraic[328] = algebraic[326]-algebraic[324] algebraic[213] = (2.00000*states[50])*constants[263] algebraic[215] = exp(algebraic[213]) algebraic[217] = ((((constants[101]*2.00000)*constants[87])*algebraic[213])*(algebraic[186]*algebraic[215]-0.341000*constants[97]))/(algebraic[215]-1.00000) algebraic[292] = 1.00000+power(constants[134]/algebraic[186], 2.00000) algebraic[293] = 1.00000+constants[140]*algebraic[285] algebraic[282] = power(states[48], 3.00000) algebraic[294] = (constants[136]*algebraic[282]+constants[364]*algebraic[186])+(constants[362]*constants[97])*(1.00000+algebraic[186]/constants[135]) algebraic[295] = ((constants[135]*constants[366])*(1.00000+algebraic[282]/constants[362])+algebraic[282]*constants[97])+constants[366]*algebraic[186] algebraic[296] = (0.800000*constants[141])*((algebraic[282]*constants[97])*algebraic[284]-(constants[366]*algebraic[186])*algebraic[285]) algebraic[297] = algebraic[296]/((algebraic[292]*algebraic[293])*(algebraic[294]+algebraic[295])) algebraic[304] = (constants[151]*algebraic[186])/(constants[152]+algebraic[186]) algebraic[327] = -constants[238]*log(constants[98]/states[46]) algebraic[331] = constants[238]*log(constants[99]/states[126]) algebraic[332] = algebraic[331]-algebraic[327] algebraic[333] = (constants[91]*algebraic[332])/(algebraic[332]+constants[92]) algebraic[352] = constants[238]*log(constants[100]/states[48]) algebraic[353] = power(algebraic[352]-algebraic[327], 4.00000) algebraic[354] = (constants[93]*algebraic[353])/(algebraic[353]+constants[346]) algebraic[329] = constants[106]*(states[50]-algebraic[327]) algebraic[355] = (((constants[241]*(power(states[114], 3.00000)))*states[106])*states[107])*(states[50]-algebraic[352]) algebraic[356] = (((constants[241]*(power(states[108], 3.00000)))*states[106])*states[107])*(states[50]-algebraic[352]) algebraic[357] = (((constants[241]*(power(states[111], 3.00000)))*states[109])*states[110])*(states[50]-algebraic[352]) algebraic[358] = ((((constants[241]*(power(states[114], 3.00000)))*states[112])*states[113])*(states[50]-algebraic[352]))*1.25000 algebraic[273] = (states[105]-0.239480)/(0.950143-0.239480) algebraic[274] = custom_piecewise([less(algebraic[273] , 0.00000), 0.00000 , True, algebraic[273]]) algebraic[275] = algebraic[274]*states[38] algebraic[276] = states[38]-algebraic[275] algebraic[277] = algebraic[274]-algebraic[275] algebraic[278] = ((1.00000-algebraic[277])-algebraic[276])-algebraic[275] algebraic[359] = ((algebraic[278]*algebraic[357]+algebraic[277]*algebraic[358])+algebraic[276]*algebraic[356])+algebraic[275]*algebraic[355] algebraic[279] = states[50]*constants[263] algebraic[280] = exp(algebraic[279]) algebraic[281] = (((constants[133]*constants[87])*algebraic[279])*(states[48]*algebraic[280]-constants[100]))/(algebraic[280]-1.00000) algebraic[300] = (constants[144]*constants[242])/(1.00000+exp(-(states[50]+92.0000)*constants[263])) algebraic[301] = algebraic[300]*(power(states[48]/(states[48]+constants[147]), 3.00000)) algebraic[302] = algebraic[300]*(power(states[48]/(states[48]+constants[148]), 3.00000)) algebraic[298] = (states[115]-0.126345)/(0.998014-0.126345) algebraic[299] = custom_piecewise([less(algebraic[298] , 0.00000), 0.00000 , True, algebraic[298]]) algebraic[303] = (1.00000-algebraic[299])*algebraic[301]+algebraic[299]*algebraic[302] algebraic[360] = ((power(states[117], 3.00000))*states[116])*(states[50]-algebraic[352]) algebraic[361] = 0.0160000*algebraic[360] algebraic[362] = 0.00650000*algebraic[360] algebraic[363] = (1.00000-states[38])*algebraic[362]+states[38]*algebraic[361] algebraic[365] = (((algebraic[359]+algebraic[281])+algebraic[363])+algebraic[303]*3.00000)+algebraic[297]*3.00000 algebraic[305] = (((algebraic[248]+algebraic[217])+algebraic[304])-algebraic[297]*2.00000)-algebraic[291]*2.00000 algebraic[330] = algebraic[329]+algebraic[323] algebraic[334] = states[50]-algebraic[331] algebraic[335] = 1.02000/(1.00000+exp(0.238500*(algebraic[334]-59.2150))) algebraic[336] = (0.491240*exp(0.0803200*(algebraic[334]+5.47600))+exp(0.0617500*(algebraic[334]-594.310)))/(1.00000+exp(-0.514300*(algebraic[334]+4.75300))) algebraic[337] = (constants[239]*(algebraic[335]/(algebraic[335]+algebraic[336])))*algebraic[334] algebraic[338] = algebraic[337]*1.20000 algebraic[339] = (1.00000-states[37])*algebraic[337]+states[37]*algebraic[338] algebraic[256] = 1.00000/(1.00000+exp((states[50]+10.0000)/15.4000)) algebraic[340] = ((constants[240]*states[68])*algebraic[256])*(states[50]-algebraic[331]) algebraic[348] = constants[238]*log((constants[99]+constants[165]*constants[100])/(states[126]+constants[165]*states[48])) algebraic[259] = 0.195610*(1.00000+0.600000/(1.00000+power(3.80000e-05/algebraic[186], 1.40000))) algebraic[349] = (algebraic[259]*(states[69]+states[70]))*(states[50]-algebraic[348]) algebraic[350] = (algebraic[259]*(states[71]+states[72]))*(states[50]-algebraic[348]) algebraic[262] = (states[103]+constants[403])/constants[110] algebraic[264] = (algebraic[262]-constants[404])/(0.785000-constants[404]) algebraic[266] = custom_piecewise([less(algebraic[264] , 0.00000), 0.00000 , True, algebraic[264]]) algebraic[351] = algebraic[266]*algebraic[350]+(1.00000-algebraic[266])*algebraic[349] algebraic[341] = (constants[123]*(states[50]-algebraic[331]))/(1.00000+exp((15.0000-states[50])/17.0000)) algebraic[342] = ((constants[123]*(states[50]-algebraic[331]))/(1.00000+exp((36.0000-states[50])/17.0000)))*3.62000 algebraic[268] = (states[104]-0.0589380)/(0.393747-0.0589380) algebraic[270] = custom_piecewise([less(algebraic[268] , 0.00000), 0.00000 , True, algebraic[268]]) algebraic[343] = (1.00000-algebraic[270])*algebraic[341]+algebraic[270]*algebraic[342] algebraic[325] = exp(states[50]/550.000) algebraic[344] = ((constants[157]*(power(states[120], 3.00000)))*algebraic[325])*(states[50]-algebraic[331]) algebraic[345] = algebraic[344]*(0.735600*states[121]+0.264400*states[123]) algebraic[346] = algebraic[344]*(0.735600*states[122]+0.264400*states[124]) algebraic[347] = (1.00000-states[39])*algebraic[346]+states[39]*algebraic[345] algebraic[364] = ((((algebraic[339]+algebraic[340])+algebraic[351])+algebraic[343])+algebraic[347])-2.00000*algebraic[303] algebraic[366] = algebraic[365]+algebraic[291]*3.00000 algebraic[367] = ((algebraic[366]+algebraic[364])+algebraic[305])+algebraic[330] algebraic[368] = (custom_piecewise([less((voi-constants[213])-constants[214]*floor((voi-constants[213])/constants[214]) , constants[212]), 1.00000 , True, 0.00000]))*constants[256] return algebraic 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)