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 = 253
sizeStates = 67
sizeConstants = 174
from math import *
from numpy import *
def createLegends():
legend_states = [""] * sizeStates
legend_rates = [""] * sizeStates
legend_algebraic = [""] * sizeAlgebraic
legend_voi = ""
legend_constants = [""] * sizeConstants
legend_voi = "time in component environment (ms)"
legend_constants[143] = "a1 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb)"
legend_constants[148] = "a2 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb)"
legend_constants[0] = "Faraday in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (coulomb_per_millimole)"
legend_constants[1] = "Temp in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (kelvin)"
legend_constants[2] = "Rgas in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (joule_per_mole_kelvin)"
legend_constants[115] = "RT_over_F in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV)"
legend_constants[3] = "Acap in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (cm2)"
legend_constants[123] = "C in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mF)"
legend_constants[4] = "Vmyo in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL)"
legend_constants[5] = "VJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL)"
legend_constants[6] = "VNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL)"
legend_constants[7] = "VSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL)"
legend_states[0] = "Nai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_states[1] = "Ki in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_states[2] = "Cai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_states[3] = "CaSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_states[4] = "CaJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_states[5] = "CaNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_states[6] = "V in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV)"
legend_algebraic[128] = "INa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[225] = "INab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[240] = "INaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_algebraic[236] = "INaK in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_algebraic[203] = "IKv14_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[132] = "IKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[136] = "IKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[223] = "IK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[251] = "ICaK in component COMPUTE_ICa_ICaK (uA_per_uF)"
legend_algebraic[246] = "ICa in component COMPUTE_ICa_ICaK (uA_per_uF)"
legend_algebraic[217] = "Ito1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[11] = "i_Stim in component I_stimulus (uA_per_uF)"
legend_algebraic[141] = "IKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[185] = "IKv14_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_algebraic[227] = "ICab in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_algebraic[229] = "IpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_algebraic[58] = "Jxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms)"
legend_algebraic[42] = "Jup in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms)"
legend_algebraic[88] = "Jtrpn in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM_per_ms)"
legend_algebraic[50] = "Jrel in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms)"
legend_algebraic[59] = "Jtr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms)"
legend_algebraic[98] = "beta_SS in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_algebraic[68] = "beta_JSR in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_algebraic[108] = "beta_i in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_algebraic[252] = "i_tot in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uA_per_uF)"
legend_constants[8] = "Ko in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_constants[9] = "Nao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_constants[10] = "Cao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_constants[11] = "stim_period in component I_stimulus (ms)"
legend_constants[12] = "stim_duration in component I_stimulus (ms)"
legend_constants[13] = "stim_amplitude in component I_stimulus (uA_per_uF)"
legend_constants[14] = "stim_offset in component I_stimulus (ms)"
legend_algebraic[0] = "past in component I_stimulus (ms)"
legend_algebraic[22] = "fb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless)"
legend_constants[15] = "Kfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM)"
legend_constants[16] = "Nfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless)"
legend_algebraic[33] = "rb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless)"
legend_constants[17] = "Krb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM)"
legend_constants[18] = "Nrb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless)"
legend_constants[19] = "KSR in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM)"
legend_constants[20] = "vmaxf in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms)"
legend_constants[21] = "vmaxr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms)"
legend_constants[22] = "v1 in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms)"
legend_states[7] = "O1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_states[8] = "O2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_constants[23] = "tautr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms)"
legend_constants[24] = "tauxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms)"
legend_constants[25] = "LTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[26] = "HTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[27] = "khtrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms)"
legend_constants[28] = "khtrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms)"
legend_constants[29] = "kltrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms)"
legend_constants[30] = "kltrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms)"
legend_constants[31] = "CMDNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[32] = "CSQNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[33] = "EGTAtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[34] = "KmCMDN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[35] = "KmCSQN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_constants[36] = "KmEGTA in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM)"
legend_algebraic[67] = "dLTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms)"
legend_algebraic[78] = "dHTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms)"
legend_states[9] = "LTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_states[10] = "HTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_constants[37] = "kaplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM4_per_ms)"
legend_constants[38] = "kaminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_constants[39] = "kbplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM3_per_ms)"
legend_constants[40] = "kbminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_constants[41] = "kcplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_constants[42] = "kcminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_constants[43] = "ncoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_constants[44] = "mcoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_algebraic[1] = "dC1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_algebraic[12] = "dO2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_algebraic[23] = "dC2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_algebraic[34] = "dO1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms)"
legend_states[11] = "C1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_states[12] = "C2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_constants[45] = "fL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[46] = "gL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[47] = "bL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_constants[48] = "aL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_algebraic[51] = "C0_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[60] = "C1_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[69] = "C2_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[79] = "C3_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[121] = "C1_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[125] = "C2_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[129] = "C3_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[133] = "C4_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[89] = "CCa0_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[99] = "CCa1_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[109] = "CCa2_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[117] = "CCa3_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[137] = "CCa1_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[142] = "CCa2_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[148] = "CCa3_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[156] = "CCa4_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[166] = "C0_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[175] = "C1_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[186] = "C2_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[197] = "C3_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[204] = "C4_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[131] = "CCa0_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[147] = "CCa1_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[150] = "CCa2_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[156] = "CCa3_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[161] = "CCa4_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[2] = "alpha in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[13] = "beta in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[24] = "alpha_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[35] = "beta_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[43] = "gamma in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_constants[116] = "omega in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[176] = "a1_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[187] = "a2_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[188] = "a1_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[198] = "a2_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[199] = "a1_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[205] = "a2_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[206] = "a1_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[211] = "a2_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[212] = "a1_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[218] = "a2_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[177] = "a1_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[189] = "a2_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[190] = "a1_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[200] = "a2_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[201] = "a1_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[207] = "a2_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[208] = "a1_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[213] = "a2_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[214] = "a1_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_algebraic[219] = "a2_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms)"
legend_states[13] = "C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[14] = "C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[15] = "C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[16] = "C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[17] = "C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[18] = "CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[19] = "CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[20] = "CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[21] = "CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[22] = "CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[23] = "Open in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_states[24] = "yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_constants[117] = "a1_Cainf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_algebraic[3] = "yCa_inf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_algebraic[14] = "tau_yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (ms)"
legend_algebraic[4] = "alpha_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[15] = "beta_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[25] = "alpha_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[36] = "beta_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[44] = "C0Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[52] = "C1Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[61] = "C2Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[70] = "C3Kv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[80] = "CI0Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[90] = "CI1Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[100] = "CI2Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[110] = "CI3Kv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[118] = "C1Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[122] = "C2Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[126] = "C3Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[130] = "OKv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[134] = "CI1Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[139] = "CI2Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[144] = "CI3Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[150] = "OIKv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[138] = "C0Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[143] = "C1Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[149] = "C2Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[157] = "C3Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[158] = "OKv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[145] = "CI0Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[151] = "CI1Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[159] = "CI2Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[167] = "CI3Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[168] = "OIKv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_states[25] = "C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[26] = "C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[27] = "C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[28] = "C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[29] = "OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[30] = "CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[31] = "CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[32] = "CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[33] = "CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_states[34] = "OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_algebraic[152] = "a1_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[160] = "a2_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[161] = "a1_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[169] = "a2_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[170] = "a1_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[178] = "a2_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[179] = "a1_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[191] = "a2_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[180] = "a1_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[192] = "a2_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[153] = "a1_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[162] = "a2_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[163] = "a1_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[171] = "a2_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[172] = "a1_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[181] = "a2_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[182] = "a1_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[193] = "a2_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[183] = "a1_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_algebraic[194] = "a2_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_constants[49] = "alphaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_constants[50] = "aaKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV)"
legend_constants[51] = "betaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_constants[52] = "baKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV)"
legend_constants[53] = "alphai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_constants[54] = "aiKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV)"
legend_constants[55] = "betai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms)"
legend_constants[56] = "biKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV)"
legend_constants[57] = "f1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[58] = "f2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[59] = "f3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[60] = "f4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[61] = "b1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[62] = "b2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[63] = "b3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[64] = "b4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_constants[65] = "f1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[66] = "f2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[67] = "f3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[68] = "f4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[69] = "b1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[70] = "b2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[71] = "b3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[72] = "b4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_constants[73] = "alphaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[74] = "aaKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV)"
legend_constants[75] = "betaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[76] = "baKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV)"
legend_constants[77] = "alphai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[78] = "betai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[5] = "alpha_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[16] = "beta_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[118] = "alpha_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[119] = "beta_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[26] = "C0Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[37] = "C1Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[45] = "C2Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[53] = "C3Kv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[27] = "CI0Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[38] = "CI1Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[46] = "CI2Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[54] = "CI3Kv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[62] = "C1Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[71] = "C2Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[81] = "C3Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[91] = "OKv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[63] = "CI1Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[72] = "CI2Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[82] = "CI3Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[92] = "OIKv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[132] = "C0Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[133] = "C1Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[134] = "C2Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[135] = "C3Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[136] = "OKv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[137] = "CI0Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[138] = "CI1Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[139] = "CI2Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[140] = "CI3Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_constants[141] = "OIKv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[73] = "a1_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[83] = "a2_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[35] = "C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[84] = "a1_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[93] = "a2_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[36] = "C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[94] = "a1_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[101] = "a2_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[37] = "C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[102] = "a1_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[111] = "a2_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[38] = "C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[103] = "a1_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[112] = "a2_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[39] = "OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[74] = "a1_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[85] = "a2_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[40] = "CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[86] = "a1_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[95] = "a2_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[41] = "CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[96] = "a1_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[104] = "a2_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[42] = "CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[105] = "a1_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[113] = "a2_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[43] = "CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[106] = "a1_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_algebraic[114] = "a2_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms)"
legend_states[44] = "OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_algebraic[116] = "ENa in component COMPUTE_REVERSAL_POTENTIALS (mV)"
legend_algebraic[120] = "EK in component COMPUTE_REVERSAL_POTENTIALS (mV)"
legend_algebraic[17] = "EKs in component COMPUTE_REVERSAL_POTENTIALS (mV)"
legend_algebraic[124] = "ECa in component COMPUTE_REVERSAL_POTENTIALS (mV)"
legend_constants[142] = "a1 in component COMPUTE_REVERSAL_POTENTIALS (mM)"
legend_algebraic[6] = "a2 in component COMPUTE_REVERSAL_POTENTIALS (mM)"
legend_constants[79] = "GKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF)"
legend_constants[80] = "GKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF)"
legend_constants[81] = "GK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF)"
legend_constants[82] = "GNa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF)"
legend_constants[83] = "GNab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF)"
legend_constants[84] = "KvScale in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless)"
legend_constants[85] = "Kv43Frac in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless)"
legend_constants[120] = "GKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF)"
legend_constants[121] = "PKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (litre_per_farad_second)"
legend_states[45] = "na6 in component INa (dimensionless)"
legend_states[46] = "na7 in component INa (dimensionless)"
legend_states[47] = "OHerg in component IKr (dimensionless)"
legend_states[48] = "O1ks in component IKs (dimensionless)"
legend_states[49] = "O2ks in component IKs (dimensionless)"
legend_constants[122] = "fKo in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless)"
legend_algebraic[147] = "VF_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless)"
legend_algebraic[155] = "VFsq_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (coulomb_per_millimole)"
legend_algebraic[165] = "a1_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM)"
legend_algebraic[174] = "a2 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless)"
legend_algebraic[196] = "a1_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM)"
legend_algebraic[221] = "K1_inf in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless)"
legend_algebraic[210] = "IKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF)"
legend_constants[86] = "kNaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_constants[87] = "KmNa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM)"
legend_constants[88] = "KmCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM)"
legend_constants[89] = "ksat in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_constants[90] = "eta in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_constants[91] = "INaKmax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_constants[92] = "KmNai in component COMPUTE_INaK_INaCa_ICab_IpCa (mM)"
legend_constants[93] = "KmKo in component COMPUTE_INaK_INaCa_ICab_IpCa (mM)"
legend_constants[94] = "IpCamax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF)"
legend_constants[95] = "KmpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM)"
legend_constants[96] = "GCab in component COMPUTE_INaK_INaCa_ICab_IpCa (mS_per_uF)"
legend_algebraic[231] = "VF_over_RT in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_constants[151] = "sigma in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_algebraic[232] = "a1_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_algebraic[233] = "a2_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_algebraic[234] = "fNaK in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_constants[155] = "a1_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_algebraic[235] = "a2_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_algebraic[237] = "a1_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4)"
legend_algebraic[238] = "a2_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4)"
legend_algebraic[239] = "a3_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless)"
legend_constants[160] = "a4_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM)"
legend_constants[163] = "a5_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM3)"
legend_algebraic[248] = "PKprime in component COMPUTE_ICa_ICaK (litre_per_farad_second)"
legend_algebraic[241] = "VF_over_RT in component COMPUTE_ICa_ICaK (dimensionless)"
legend_algebraic[242] = "VFsq_over_RT in component COMPUTE_ICa_ICaK (coulomb_per_millimole)"
legend_algebraic[243] = "a1_Ca in component COMPUTE_ICa_ICaK (mM)"
legend_algebraic[244] = "a2_Ca in component COMPUTE_ICa_ICaK (dimensionless)"
legend_algebraic[249] = "a1_K in component COMPUTE_ICa_ICaK (mM)"
legend_algebraic[250] = "a2_K in component COMPUTE_ICa_ICaK (dimensionless)"
legend_algebraic[245] = "ICamax in component COMPUTE_ICa_ICaK (uA_per_uF)"
legend_algebraic[247] = "Icabar in component COMPUTE_ICa_ICaK (uA_per_uF)"
legend_constants[144] = "PCa in component COMPUTE_ICa_ICaK (litre_per_farad_second)"
legend_constants[124] = "PK in component COMPUTE_ICa_ICaK (litre_per_farad_second)"
legend_constants[97] = "ICahalf in component COMPUTE_ICa_ICaK (uA_per_uF)"
legend_constants[98] = "Pscale in component COMPUTE_ICa_ICaK (dimensionless)"
legend_constants[149] = "Temp_Scale in component INa (dimensionless)"
legend_algebraic[7] = "alpha1 in component INa (per_ms)"
legend_algebraic[18] = "beta1 in component INa (per_ms)"
legend_algebraic[28] = "gamma1 in component INa (per_ms)"
legend_algebraic[39] = "Delta1 in component INa (per_ms)"
legend_algebraic[47] = "On in component INa (per_ms)"
legend_algebraic[55] = "Of in component INa (per_ms)"
legend_algebraic[64] = "GammaGamma in component INa (per_ms)"
legend_algebraic[75] = "DeltaDelta in component INa (per_ms)"
legend_constants[153] = "epsilon in component INa (per_ms)"
legend_constants[154] = "omega_na in component INa (per_ms)"
legend_algebraic[87] = "rho in component INa (per_ms)"
legend_algebraic[97] = "mu in component INa (per_ms)"
legend_constants[152] = "Cn in component INa (per_ms)"
legend_constants[157] = "Cf in component INa (per_ms)"
legend_constants[162] = "parameter_a in component INa (dimensionless)"
legend_algebraic[107] = "k12 in component INa (per_ms)"
legend_algebraic[115] = "k23 in component INa (per_ms)"
legend_algebraic[119] = "k34 in component INa (per_ms)"
legend_algebraic[123] = "k45 in component INa (per_ms)"
legend_algebraic[127] = "k56 in component INa (per_ms)"
legend_constants[158] = "k67 in component INa (per_ms)"
legend_algebraic[131] = "k89 in component INa (per_ms)"
legend_algebraic[135] = "k910 in component INa (per_ms)"
legend_algebraic[140] = "k1011 in component INa (per_ms)"
legend_algebraic[146] = "k1112 in component INa (per_ms)"
legend_algebraic[154] = "k1213 in component INa (per_ms)"
legend_algebraic[164] = "k57 in component INa (per_ms)"
legend_algebraic[173] = "k21 in component INa (per_ms)"
legend_algebraic[184] = "k32 in component INa (per_ms)"
legend_algebraic[195] = "k43 in component INa (per_ms)"
legend_algebraic[202] = "k54 in component INa (per_ms)"
legend_algebraic[209] = "k65 in component INa (per_ms)"
legend_constants[159] = "k76 in component INa (per_ms)"
legend_algebraic[215] = "k98 in component INa (per_ms)"
legend_algebraic[220] = "k109 in component INa (per_ms)"
legend_algebraic[222] = "k1110 in component INa (per_ms)"
legend_algebraic[224] = "k1211 in component INa (per_ms)"
legend_algebraic[226] = "k1312 in component INa (per_ms)"
legend_algebraic[216] = "k75 in component INa (per_ms)"
legend_constants[164] = "k81 in component INa (per_ms)"
legend_constants[166] = "k92 in component INa (per_ms)"
legend_constants[168] = "k103 in component INa (per_ms)"
legend_constants[170] = "k114 in component INa (per_ms)"
legend_constants[172] = "k125 in component INa (per_ms)"
legend_algebraic[228] = "k136 in component INa (per_ms)"
legend_constants[165] = "k18 in component INa (per_ms)"
legend_constants[167] = "k29 in component INa (per_ms)"
legend_constants[169] = "k310 in component INa (per_ms)"
legend_constants[171] = "k411 in component INa (per_ms)"
legend_constants[173] = "k512 in component INa (per_ms)"
legend_algebraic[230] = "k613 in component INa (per_ms)"
legend_states[50] = "na1 in component INa (dimensionless)"
legend_states[51] = "na2 in component INa (dimensionless)"
legend_states[52] = "na3 in component INa (dimensionless)"
legend_states[53] = "na4 in component INa (dimensionless)"
legend_states[54] = "na5 in component INa (dimensionless)"
legend_states[55] = "na8 in component INa (dimensionless)"
legend_states[56] = "na9 in component INa (dimensionless)"
legend_states[57] = "na10 in component INa (dimensionless)"
legend_states[58] = "na11 in component INa (dimensionless)"
legend_states[59] = "na12 in component INa (dimensionless)"
legend_states[60] = "na13 in component INa (dimensionless)"
legend_constants[99] = "TNa in component INa (kelvin)"
legend_constants[125] = "KToverH in component INa (per_ms)"
legend_constants[146] = "FoverRT in component INa (per_mV)"
legend_constants[145] = "RTNa in component INa (joule_per_mole)"
legend_constants[126] = "RTNaF in component INa (mV)"
legend_constants[100] = "T_Const_HERG in component IKr (dimensionless)"
legend_constants[101] = "A0_HERG in component IKr (per_ms)"
legend_constants[102] = "B0_HERG in component IKr (per_mV)"
legend_constants[103] = "A1_HERG in component IKr (per_ms)"
legend_constants[104] = "B1_HERG in component IKr (per_mV)"
legend_constants[105] = "A2_HERG in component IKr (per_ms)"
legend_constants[106] = "B2_HERG in component IKr (per_mV)"
legend_constants[107] = "A3_HERG in component IKr (per_ms)"
legend_constants[108] = "B3_HERG in component IKr (per_mV)"
legend_constants[109] = "A4_HERG in component IKr (per_ms)"
legend_constants[110] = "B4_HERG in component IKr (per_mV)"
legend_constants[111] = "A5_HERG in component IKr (per_ms)"
legend_constants[112] = "B5_HERG in component IKr (per_mV)"
legend_constants[113] = "A6_HERG in component IKr (per_ms)"
legend_constants[114] = "B6_HERG in component IKr (per_mV)"
legend_constants[127] = "C2H_to_C3H in component IKr (per_ms)"
legend_constants[128] = "C3H_to_C2H in component IKr (per_ms)"
legend_algebraic[8] = "C1H_to_C2H in component IKr (per_ms)"
legend_algebraic[19] = "C2H_to_C1H in component IKr (per_ms)"
legend_algebraic[9] = "C3H_to_OH in component IKr (per_ms)"
legend_algebraic[20] = "OH_to_C3H in component IKr (per_ms)"
legend_algebraic[29] = "OH_to_IH in component IKr (per_ms)"
legend_algebraic[40] = "IH_to_OH in component IKr (per_ms)"
legend_algebraic[48] = "C3H_to_IH in component IKr (per_ms)"
legend_algebraic[56] = "IH_to_C3H in component IKr (per_ms)"
legend_states[61] = "C1Herg in component IKr (dimensionless)"
legend_algebraic[30] = "a1_C2 in component IKr (per_ms)"
legend_algebraic[41] = "a2_C2 in component IKr (per_ms)"
legend_states[62] = "C2Herg in component IKr (dimensionless)"
legend_algebraic[65] = "a1_C3 in component IKr (per_ms)"
legend_algebraic[76] = "a2_C3 in component IKr (per_ms)"
legend_states[63] = "C3Herg in component IKr (dimensionless)"
legend_algebraic[49] = "a1_O in component IKr (per_ms)"
legend_algebraic[57] = "a2_O in component IKr (per_ms)"
legend_algebraic[66] = "a1_I in component IKr (per_ms)"
legend_algebraic[77] = "a2_I in component IKr (per_ms)"
legend_states[64] = "IHerg in component IKr (dimensionless)"
legend_constants[129] = "C0ks_C1ks in component IKs (per_ms)"
legend_constants[130] = "C1ks_O1ks in component IKs (per_ms)"
legend_algebraic[10] = "O1ks_O2ks in component IKs (per_ms)"
legend_algebraic[21] = "O1ks_C1ks in component IKs (per_ms)"
legend_algebraic[31] = "O2ks_O1ks in component IKs (per_ms)"
legend_algebraic[32] = "C1ks_C0ks in component IKs (per_ms)"
legend_states[65] = "C0ks in component IKs (dimensionless)"
legend_states[66] = "C1ks in component IKs (dimensionless)"
legend_rates[0] = "d/dt Nai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_rates[1] = "d/dt Ki in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_rates[2] = "d/dt Cai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_rates[3] = "d/dt CaSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_rates[4] = "d/dt CaJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_rates[5] = "d/dt CaNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM)"
legend_rates[6] = "d/dt V in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV)"
legend_rates[9] = "d/dt LTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_rates[10] = "d/dt HTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless)"
legend_rates[11] = "d/dt C1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_rates[8] = "d/dt O2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_rates[12] = "d/dt C2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_rates[7] = "d/dt O1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless)"
legend_rates[13] = "d/dt C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[14] = "d/dt C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[15] = "d/dt C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[16] = "d/dt C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[17] = "d/dt C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[23] = "d/dt Open in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[18] = "d/dt CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[19] = "d/dt CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[20] = "d/dt CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[21] = "d/dt CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[22] = "d/dt CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[24] = "d/dt yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless)"
legend_rates[25] = "d/dt C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[26] = "d/dt C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[27] = "d/dt C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[28] = "d/dt C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[29] = "d/dt OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[30] = "d/dt CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[31] = "d/dt CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[32] = "d/dt CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[33] = "d/dt CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[34] = "d/dt OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless)"
legend_rates[35] = "d/dt C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[36] = "d/dt C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[37] = "d/dt C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[38] = "d/dt C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[39] = "d/dt OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[40] = "d/dt CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[41] = "d/dt CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[42] = "d/dt CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[43] = "d/dt CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[44] = "d/dt OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless)"
legend_rates[50] = "d/dt na1 in component INa (dimensionless)"
legend_rates[51] = "d/dt na2 in component INa (dimensionless)"
legend_rates[52] = "d/dt na3 in component INa (dimensionless)"
legend_rates[53] = "d/dt na4 in component INa (dimensionless)"
legend_rates[54] = "d/dt na5 in component INa (dimensionless)"
legend_rates[45] = "d/dt na6 in component INa (dimensionless)"
legend_rates[46] = "d/dt na7 in component INa (dimensionless)"
legend_rates[55] = "d/dt na8 in component INa (dimensionless)"
legend_rates[56] = "d/dt na9 in component INa (dimensionless)"
legend_rates[57] = "d/dt na10 in component INa (dimensionless)"
legend_rates[58] = "d/dt na11 in component INa (dimensionless)"
legend_rates[59] = "d/dt na12 in component INa (dimensionless)"
legend_rates[60] = "d/dt na13 in component INa (dimensionless)"
legend_rates[61] = "d/dt C1Herg in component IKr (dimensionless)"
legend_rates[62] = "d/dt C2Herg in component IKr (dimensionless)"
legend_rates[63] = "d/dt C3Herg in component IKr (dimensionless)"
legend_rates[47] = "d/dt OHerg in component IKr (dimensionless)"
legend_rates[64] = "d/dt IHerg in component IKr (dimensionless)"
legend_rates[65] = "d/dt C0ks in component IKs (dimensionless)"
legend_rates[66] = "d/dt C1ks in component IKs (dimensionless)"
legend_rates[48] = "d/dt O1ks in component IKs (dimensionless)"
legend_rates[49] = "d/dt O2ks in component IKs (dimensionless)"
return (legend_states, legend_algebraic, legend_voi, legend_constants)
def initConsts():
constants = [0.0] * sizeConstants; states = [0.0] * sizeStates;
constants[0] = 96.5
constants[1] = 310
constants[2] = 8.315
constants[3] = 1.534e-4
constants[4] = 25.84e-6
constants[5] = 0.16e-6
constants[6] = 2.1e-6
constants[7] = 1.2e-9
states[0] = 0.9798304162e1
states[1] = 0.1255589432e3
states[2] = 0.8601192016e-4
states[3] = 0.1420215245e-3
states[4] = 0.2852239446
states[5] = 0.2855294915
states[6] = -0.9065755929e2
constants[8] = 4
constants[9] = 138
constants[10] = 2
constants[11] = 1e3
constants[12] = 3
constants[13] = -15
constants[14] = 0
constants[15] = 0.000168
constants[16] = 1.2
constants[17] = 3.29
constants[18] = 1
constants[19] = 1.2
constants[20] = 0.0748e-3
constants[21] = 0.318e-3
constants[22] = 1.8
states[7] = 0.6601783287e-3
states[8] = 0.397392006e-8
constants[23] = 0.5747
constants[24] = 26.7
constants[25] = 70e-3
constants[26] = 140e-3
constants[27] = 20e0
constants[28] = 0.066e-3
constants[29] = 40e0
constants[30] = 40e-3
constants[31] = 50e-3
constants[32] = 15e0
constants[33] = 0e0
constants[34] = 2.38e-3
constants[35] = 0.8e0
constants[36] = 1.5e-4
states[9] = 0.8046584973e-1
states[10] = 0.9772152481
constants[37] = 0.01215
constants[38] = 0.576
constants[39] = 0.00405
constants[40] = 1.93
constants[41] = 0.1
constants[42] = 0.0008
constants[43] = 4
constants[44] = 3
states[11] = 0.4448038946
states[12] = 0.554535924
constants[45] = 0.3
constants[46] = 4
constants[47] = 2
constants[48] = 2
states[13] = 0.8621934054
states[14] = 0.01141398211
states[15] = 5.666325225e-5
states[16] = 1.25021357e-7
states[17] = 1.034426748e-10
states[18] = 0.1198547081
states[19] = 0.006346794302
states[20] = 1.260326488e-4
states[21] = 1.112315238e-6
states[22] = 3.681315892e-9
states[23] = 7.757981563e-12
states[24] = 0.9997157074
states[25] = 0.9513721351
states[26] = 0.02668288089
states[27] = 2.806380358e-4
states[28] = 1.311837579e-6
states[29] = 2.29955624e-9
states[30] = 0.01513025204
states[31] = 0.005442964601
states[32] = 9.918373359000001e-4
states[33] = 9.514386057000001e-5
states[34] = 2.742677382e-6
constants[49] = 0.543708
constants[50] = 0.028983
constants[51] = 0.080185
constants[52] = 0.0468437
constants[53] = 0.0498424
constants[54] = 0.000373016
constants[55] = 0.000819482
constants[56] = 0.00000005374
constants[57] = 1.8936
constants[58] = 14.224647456
constants[59] = 158.574378389
constants[60] = 142.936645351
constants[61] = 6.77348
constants[62] = 15.6212705152
constants[63] = 28.7532603313
constants[64] = 524.576206679
constants[65] = 0.52465073996
constants[66] = 17.51885408639
constants[67] = 938.58764534556
constants[68] = 54749.19473332601
constants[69] = 1.00947847105
constants[70] = 1.17100540567
constants[71] = 0.63902768758
constants[72] = 2.12035379095
constants[73] = 1.84002414554
constants[74] = 0.00768548031
constants[75] = 0.0108174834
constants[76] = 0.07793378174
constants[77] = 0.00305767916
constants[78] = 0.00000244936
states[35] = 0.5977099765
states[36] = 0.1730990528
states[37] = 0.01881072386
states[38] = 9.16070135e-4
states[39] = 1.975541357e-5
states[40] = 0.03539084346
states[41] = 0.005428824353
states[42] = 0.02287858869
states[43] = 0.03233800003
states[44] = 0.1134082058
constants[79] = 0.0186e0
constants[80] = 0.0035e0
constants[81] = 0.1253051261188079722568586886469e0
constants[82] = 56.32e0
constants[83] = 0.001e0
constants[84] = 0.872
constants[85] = 0.889
states[45] = 1.57442749e-8
states[46] = 2.856607179e-9
states[47] = 9.511789113e-6
states[48] = 2.491710696e-7
states[49] = 1.298547822e-5
constants[86] = 0.44
constants[87] = 87.5
constants[88] = 1.38
constants[89] = 0.2
constants[90] = 0.35
constants[91] = 2.387
constants[92] = 20
constants[93] = 1.5
constants[94] = 0.05
constants[95] = 0.0005
constants[96] = 7.684e-5
constants[97] = -0.265
constants[98] = 7
states[50] = 0.1437575649
states[51] = 0.0417776208
states[52] = 0.004552898364
states[53] = 2.20520743e-4
states[54] = 4.005266484e-6
states[55] = 0.4749950008
states[56] = 0.2707214097
states[57] = 0.05786120057
states[58] = 0.005496292279
states[59] = 1.957874519e-4
states[60] = 4.17679396e-4
constants[99] = 294.16
constants[100] = 5.320000001
constants[101] = 0.017147641733086
constants[102] = 0.03304608038835
constants[103] = 0.03969328381141
constants[104] = -0.0430605416398
constants[105] = 0.02057448605977
constants[106] = 0.02617412715118
constants[107] = 0.00134366604423
constants[108] = -0.02691385498399
constants[109] = 0.10666316491288
constants[110] = 0.00568908859717
constants[111] = 0.00646393910049
constants[112] = -0.04536642959543
constants[113] = 0.00008039374403
constants[114] = 0.00000069808924
states[61] = 0.996697338
states[62] = 4.340879648e-4
states[63] = 7.634099755000001e-5
states[64] = 1.533347007e-6
states[65] = 0.9645606295
states[66] = 0.03542613568
constants[115] = (constants[2]*constants[1])/constants[0]
constants[116] = 0.250000*0.0100000
constants[117] = 0.820000
constants[118] = constants[77]
constants[119] = constants[78]
constants[120] = constants[85]*constants[84]*0.100000
constants[121] = (1.00000-constants[85])*constants[84]*4.29860e-07
constants[122] = power(constants[8]/4.00000, 1.0/2)
constants[123] = constants[3]*0.00100000
constants[124] = constants[98]*4.57400e-07
constants[125] = (1.38100e-23*constants[99])/6.62600e-31
constants[126] = (constants[2]*constants[99])/constants[0]
constants[127] = constants[100]*0.0260836
constants[128] = constants[100]*0.148330
constants[129] = 0.00795601
constants[130] = 0.0396672
constants[131] = constants[116]
constants[132] = constants[119]
constants[133] = constants[65]*constants[119]
constants[134] = constants[66]*constants[119]
constants[135] = constants[67]*constants[119]
constants[136] = constants[68]*constants[119]
constants[137] = constants[118]
constants[138] = constants[118]/constants[69]
constants[139] = constants[118]/constants[70]
constants[140] = constants[118]/constants[71]
constants[141] = constants[118]/constants[72]
constants[142] = constants[8]+0.0183300*constants[9]
constants[143] = constants[123]/(constants[4]*constants[0])
constants[144] = constants[98]*0.000246900
constants[145] = constants[2]*constants[99]
constants[146] = 1.00000/constants[115]
constants[147] = constants[131]/constants[47]
constants[148] = constants[123]/(2.00000*constants[7]*constants[0])
constants[149] = 1.38862
constants[150] = constants[147]/constants[47]
constants[151] = (exp(constants[9]/67.3000)-1.00000)/7.00000
constants[152] = constants[149]*constants[125]*exp(-287913./constants[145]+786.217/constants[2])
constants[153] = constants[149]*constants[125]*exp(-85800.4/constants[145]+70.0780/constants[2])
constants[154] = constants[149]*constants[125]*exp(-121955./constants[145]+225.175/constants[2])
constants[155] = constants[8]/(constants[8]+constants[93])
constants[156] = constants[150]/constants[47]
constants[157] = constants[149]*constants[125]*exp(-59565.2/constants[145]+0.00711000/constants[2])
constants[158] = constants[153]
constants[159] = constants[154]
constants[160] = constants[88]+constants[10]
constants[161] = constants[156]/constants[47]
constants[162] = 1.40043
constants[163] = (power(constants[87], 3.00000)+power(constants[9], 3.00000))/5000.00
constants[164] = constants[157]
constants[165] = constants[152]
constants[166] = constants[164]/constants[162]
constants[167] = constants[165]*constants[162]
constants[168] = constants[166]/constants[162]
constants[169] = constants[167]*constants[162]
constants[170] = constants[168]/constants[162]
constants[171] = constants[169]*constants[162]
constants[172] = constants[170]/constants[162]
constants[173] = constants[171]*constants[162]
return (states, constants)
def computeRates(voi, states, constants):
rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic
rates[23] = constants[45]*states[17]-constants[46]*states[23]
algebraic[1] = -constants[37]*(power(states[3]*1000.00, constants[43]))*states[11]+constants[38]*states[7]
rates[11] = algebraic[1]
algebraic[12] = constants[39]*(power(states[3]*1000.00, constants[44]))*states[7]-constants[40]*states[8]
rates[8] = algebraic[12]
algebraic[3] = (constants[117]/(1.00000+exp((states[6]+28.5000)/7.80000))+1.00000)-constants[117]
algebraic[14] = 1.00000/(0.00336336/(0.500000+exp(states[6]/-5.53900))+0.00779047*exp(states[6]/-49.5104))
rates[24] = (algebraic[3]-states[24])/algebraic[14]
algebraic[8] = constants[100]*constants[101]*exp(constants[102]*states[6])
algebraic[19] = constants[100]*constants[103]*exp(constants[104]*states[6])
rates[61] = algebraic[19]*states[62]-algebraic[8]*states[61]
algebraic[23] = constants[41]*states[7]-constants[42]*states[12]
rates[12] = algebraic[23]
algebraic[32] = 0.216256*exp(-1.88912e-05*states[6])
rates[65] = -constants[129]*states[65]+algebraic[32]*states[66]
algebraic[21] = 0.00700807*exp(-0.149998*states[6])
rates[66] = (constants[129]*states[65]-(algebraic[32]+constants[130])*states[66])+algebraic[21]*states[48]
algebraic[10] = 0.00767254*exp(0.0866295*states[6])
algebraic[31] = 0.00379738*exp(-0.0142567*states[6])
rates[48] = (constants[130]*states[66]-(algebraic[21]+algebraic[10])*states[48])+algebraic[31]*states[49]
rates[49] = algebraic[10]*states[48]-algebraic[31]*states[49]
algebraic[34] = -(algebraic[1]+algebraic[12]+algebraic[23])
rates[7] = algebraic[34]
algebraic[30] = algebraic[8]*states[61]+constants[128]*states[63]
algebraic[41] = (algebraic[19]+constants[127])*states[62]
rates[62] = algebraic[30]-algebraic[41]
algebraic[9] = constants[100]*constants[105]*exp(constants[106]*states[6])
algebraic[40] = constants[100]*constants[111]*exp(constants[112]*states[6])
algebraic[49] = algebraic[9]*states[63]+algebraic[40]*states[64]
algebraic[20] = constants[100]*constants[107]*exp(constants[108]*states[6])
algebraic[29] = constants[100]*constants[109]*exp(constants[110]*states[6])
algebraic[57] = (algebraic[20]+algebraic[29])*states[47]
rates[47] = algebraic[49]-algebraic[57]
algebraic[22] = power(states[2]/constants[15], constants[16])
algebraic[33] = power(states[5]/constants[17], constants[18])
algebraic[42] = (constants[19]*(constants[20]*algebraic[22]-constants[21]*algebraic[33]))/(1.00000+algebraic[22]+algebraic[33])
algebraic[59] = (states[5]-states[4])/constants[23]
rates[5] = (algebraic[42]*constants[4])/constants[6]-(algebraic[59]*constants[5])/constants[6]
algebraic[50] = constants[22]*(states[7]+states[8])*(states[4]-states[3])
algebraic[68] = 1.00000/(1.00000+(constants[32]*constants[35])/(power(states[4]+constants[35], 2.00000)))
rates[4] = algebraic[68]*(algebraic[59]-algebraic[50])
algebraic[67] = constants[29]*states[2]*(1.00000-states[9])-constants[30]*states[9]
rates[9] = algebraic[67]
algebraic[48] = constants[100]*constants[113]*exp(constants[114]*states[6])
algebraic[56] = (algebraic[20]*algebraic[40]*algebraic[48])/(algebraic[9]*algebraic[29])
algebraic[65] = constants[127]*states[62]+algebraic[20]*states[47]+algebraic[56]*states[64]
algebraic[76] = (algebraic[48]+algebraic[9]+constants[128])*states[63]
rates[63] = algebraic[65]-algebraic[76]
algebraic[66] = algebraic[48]*states[63]+algebraic[29]*states[47]
algebraic[77] = (algebraic[56]+algebraic[40])*states[64]
rates[64] = algebraic[66]-algebraic[77]
algebraic[78] = constants[27]*states[2]*(1.00000-states[10])-constants[28]*states[10]
rates[10] = algebraic[78]
algebraic[5] = constants[73]*exp(constants[74]*states[6])
algebraic[26] = 4.00000*algebraic[5]
algebraic[73] = (algebraic[26]+constants[132])*states[35]
algebraic[16] = constants[75]*exp(-constants[76]*states[6])
algebraic[62] = algebraic[16]
algebraic[83] = algebraic[62]*states[36]+constants[137]*states[40]
rates[35] = algebraic[83]-algebraic[73]
algebraic[27] = 4.00000*constants[69]*algebraic[5]
algebraic[74] = (constants[137]+algebraic[27])*states[40]
algebraic[63] = algebraic[16]/constants[65]
algebraic[85] = constants[132]*states[35]+algebraic[63]*states[41]
rates[40] = algebraic[85]-algebraic[74]
algebraic[37] = 3.00000*algebraic[5]
algebraic[84] = (algebraic[37]+algebraic[62]+constants[133])*states[36]
algebraic[71] = 2.00000*algebraic[16]
algebraic[93] = algebraic[71]*states[37]+constants[138]*states[41]+algebraic[26]*states[35]
rates[36] = algebraic[93]-algebraic[84]
algebraic[38] = (3.00000*constants[70]*algebraic[5])/constants[69]
algebraic[86] = (algebraic[38]+constants[138]+algebraic[63])*states[41]
algebraic[72] = (2.00000*constants[65]*algebraic[16])/constants[66]
algebraic[95] = algebraic[72]*states[42]+constants[133]*states[36]+algebraic[27]*states[40]
rates[41] = algebraic[95]-algebraic[86]
algebraic[45] = 2.00000*algebraic[5]
algebraic[94] = (algebraic[45]+algebraic[71]+constants[134])*states[37]
algebraic[81] = 3.00000*algebraic[16]
algebraic[101] = algebraic[81]*states[38]+constants[139]*states[42]+algebraic[37]*states[36]
rates[37] = algebraic[101]-algebraic[94]
algebraic[46] = (2.00000*constants[71]*algebraic[5])/constants[70]
algebraic[96] = (algebraic[46]+constants[139]+algebraic[72])*states[42]
algebraic[82] = (3.00000*constants[66]*algebraic[16])/constants[67]
algebraic[104] = algebraic[82]*states[43]+constants[134]*states[37]+algebraic[38]*states[41]
rates[42] = algebraic[104]-algebraic[96]
algebraic[53] = algebraic[5]
algebraic[102] = (algebraic[53]+algebraic[81]+constants[135])*states[38]
algebraic[91] = 4.00000*algebraic[16]
algebraic[111] = algebraic[91]*states[39]+constants[140]*states[43]+algebraic[45]*states[37]
rates[38] = algebraic[111]-algebraic[102]
algebraic[103] = (algebraic[91]+constants[136])*states[39]
algebraic[112] = algebraic[53]*states[38]+constants[141]*states[44]
rates[39] = algebraic[112]-algebraic[103]
algebraic[54] = (constants[72]*algebraic[5])/constants[71]
algebraic[105] = (algebraic[54]+constants[140]+algebraic[82])*states[43]
algebraic[92] = (4.00000*constants[67]*algebraic[16])/constants[68]
algebraic[113] = algebraic[92]*states[44]+constants[135]*states[38]+algebraic[46]*states[42]
rates[43] = algebraic[113]-algebraic[105]
algebraic[106] = (constants[141]+algebraic[92])*states[44]
algebraic[114] = constants[136]*states[39]+algebraic[54]*states[43]
rates[44] = algebraic[114]-algebraic[106]
algebraic[4] = constants[49]*exp(constants[50]*states[6])
algebraic[44] = 4.00000*algebraic[4]
algebraic[36] = constants[55]*exp(constants[56]*states[6])
algebraic[138] = algebraic[36]
algebraic[152] = (algebraic[44]+algebraic[138])*states[25]
algebraic[15] = constants[51]*exp(-constants[52]*states[6])
algebraic[118] = algebraic[15]
algebraic[25] = constants[53]*exp(-constants[54]*states[6])
algebraic[145] = algebraic[25]
algebraic[160] = algebraic[118]*states[26]+algebraic[145]*states[30]
rates[25] = algebraic[160]-algebraic[152]
algebraic[80] = 4.00000*constants[61]*algebraic[4]
algebraic[153] = (algebraic[145]+algebraic[80])*states[30]
algebraic[134] = algebraic[15]/constants[57]
algebraic[162] = algebraic[138]*states[25]+algebraic[134]*states[31]
rates[30] = algebraic[162]-algebraic[153]
algebraic[52] = 3.00000*algebraic[4]
algebraic[143] = constants[57]*algebraic[36]
algebraic[161] = (algebraic[52]+algebraic[118]+algebraic[143])*states[26]
algebraic[122] = 2.00000*algebraic[15]
algebraic[151] = algebraic[25]/constants[61]
algebraic[169] = algebraic[122]*states[27]+algebraic[151]*states[31]+algebraic[44]*states[25]
rates[26] = algebraic[169]-algebraic[161]
algebraic[90] = (3.00000*constants[62]*algebraic[4])/constants[61]
algebraic[163] = (algebraic[90]+algebraic[151]+algebraic[134])*states[31]
algebraic[139] = (2.00000*constants[57]*algebraic[15])/constants[58]
algebraic[171] = algebraic[139]*states[32]+algebraic[143]*states[26]+algebraic[80]*states[30]
rates[31] = algebraic[171]-algebraic[163]
algebraic[7] = constants[149]*constants[125]*exp(-114007./constants[145]+224.114/constants[2]+(0.286374*states[6])/constants[126])
algebraic[107] = 4.00000*algebraic[7]
algebraic[18] = constants[149]*constants[125]*exp(-272470./constants[145]+708.146/constants[2]+(-2.28528*states[6])/constants[126])
algebraic[173] = algebraic[18]
rates[50] = -(constants[165]+algebraic[107])*states[50]+algebraic[173]*states[51]+constants[164]*states[55]
algebraic[61] = 2.00000*algebraic[4]
algebraic[149] = constants[58]*algebraic[36]
algebraic[170] = (algebraic[61]+algebraic[122]+algebraic[149])*states[27]
algebraic[126] = 3.00000*algebraic[15]
algebraic[159] = algebraic[25]/constants[62]
algebraic[178] = algebraic[126]*states[28]+algebraic[159]*states[32]+algebraic[52]*states[26]
rates[27] = algebraic[178]-algebraic[170]
algebraic[100] = (2.00000*constants[63]*algebraic[4])/constants[62]
algebraic[172] = (algebraic[100]+algebraic[159]+algebraic[139])*states[32]
algebraic[144] = (3.00000*constants[58]*algebraic[15])/constants[59]
algebraic[181] = algebraic[144]*states[33]+algebraic[149]*states[27]+algebraic[90]*states[31]
rates[32] = algebraic[181]-algebraic[172]
algebraic[115] = 3.00000*algebraic[7]
algebraic[184] = 2.00000*algebraic[18]
rates[51] = (algebraic[107]*states[50]-(algebraic[173]+algebraic[115]+constants[167])*states[51])+algebraic[184]*states[52]+constants[166]*states[56]
algebraic[2] = 4.00000*1.20000*0.416000*exp(0.0120000*(states[6]-35.0000))
algebraic[51] = 4.00000*algebraic[2]
algebraic[43] = 0.600000*0.0923300*states[3]
algebraic[166] = algebraic[43]
algebraic[177] = (algebraic[51]+algebraic[166])*states[13]
algebraic[13] = 4.00000*0.450000*0.0490000*exp(-0.0650000*(states[6]-22.0000))
algebraic[121] = algebraic[13]
algebraic[189] = algebraic[121]*states[14]+constants[131]*states[18]
rates[13] = algebraic[189]-algebraic[177]
algebraic[24] = constants[48]*algebraic[2]
algebraic[89] = 4.00000*algebraic[24]
algebraic[176] = (algebraic[89]+constants[131])*states[18]
algebraic[35] = algebraic[13]/constants[47]
algebraic[137] = algebraic[35]
algebraic[187] = algebraic[137]*states[19]+algebraic[166]*states[13]
rates[18] = algebraic[187]-algebraic[176]
algebraic[70] = algebraic[4]
algebraic[157] = constants[59]*algebraic[36]
algebraic[179] = (algebraic[70]+algebraic[126]+algebraic[157])*states[28]
algebraic[130] = 4.00000*algebraic[15]
algebraic[167] = algebraic[25]/constants[63]
algebraic[191] = algebraic[130]*states[29]+algebraic[167]*states[33]+algebraic[61]*states[27]
rates[28] = algebraic[191]-algebraic[179]
algebraic[158] = constants[60]*algebraic[36]
algebraic[180] = (algebraic[130]+algebraic[158])*states[29]
algebraic[168] = algebraic[25]/constants[64]
algebraic[192] = algebraic[70]*states[28]+algebraic[168]*states[34]
rates[29] = algebraic[192]-algebraic[180]
algebraic[110] = (constants[64]*algebraic[4])/constants[63]
algebraic[182] = (algebraic[110]+algebraic[167]+algebraic[144])*states[33]
algebraic[150] = (4.00000*constants[59]*algebraic[15])/constants[60]
algebraic[193] = algebraic[150]*states[34]+algebraic[157]*states[28]+algebraic[100]*states[32]
rates[33] = algebraic[193]-algebraic[182]
algebraic[183] = (algebraic[168]+algebraic[150])*states[34]
algebraic[194] = algebraic[158]*states[29]+algebraic[110]*states[33]
rates[34] = algebraic[194]-algebraic[183]
algebraic[119] = 2.00000*algebraic[7]
algebraic[195] = 3.00000*algebraic[18]
rates[52] = (algebraic[115]*states[51]-(algebraic[184]+algebraic[119]+constants[169])*states[52])+algebraic[195]*states[53]+constants[168]*states[57]
algebraic[60] = 3.00000*algebraic[2]
algebraic[175] = constants[48]*algebraic[166]
algebraic[190] = (algebraic[121]+algebraic[60]+algebraic[175])*states[14]
algebraic[125] = 2.00000*algebraic[13]
algebraic[200] = algebraic[51]*states[13]+algebraic[125]*states[15]+constants[147]*states[19]
rates[14] = algebraic[200]-algebraic[190]
algebraic[99] = 3.00000*algebraic[24]
algebraic[188] = (algebraic[137]+algebraic[99]+constants[147])*states[19]
algebraic[142] = 2.00000*algebraic[35]
algebraic[198] = algebraic[89]*states[18]+algebraic[142]*states[20]+algebraic[175]*states[14]
rates[19] = algebraic[198]-algebraic[188]
algebraic[123] = algebraic[7]
algebraic[202] = 4.00000*algebraic[18]
rates[53] = (algebraic[119]*states[52]-(algebraic[195]+algebraic[123]+constants[171])*states[53])+algebraic[202]*states[54]+constants[170]*states[58]
algebraic[69] = 2.00000*algebraic[2]
algebraic[186] = constants[48]*algebraic[175]
algebraic[201] = (algebraic[125]+algebraic[69]+algebraic[186])*states[15]
algebraic[129] = 3.00000*algebraic[13]
algebraic[207] = algebraic[60]*states[14]+algebraic[129]*states[16]+constants[150]*states[20]
rates[15] = algebraic[207]-algebraic[201]
algebraic[109] = 2.00000*algebraic[24]
algebraic[199] = (algebraic[142]+algebraic[109]+constants[150])*states[20]
algebraic[148] = 3.00000*algebraic[35]
algebraic[205] = algebraic[99]*states[19]+algebraic[148]*states[21]+algebraic[186]*states[15]
rates[20] = algebraic[205]-algebraic[199]
algebraic[79] = algebraic[2]
algebraic[197] = constants[48]*algebraic[186]
algebraic[208] = (algebraic[129]+algebraic[79]+algebraic[197])*states[16]
algebraic[133] = 4.00000*algebraic[13]
algebraic[213] = algebraic[69]*states[15]+algebraic[133]*states[17]+constants[156]*states[21]
rates[16] = algebraic[213]-algebraic[208]
algebraic[117] = algebraic[24]
algebraic[206] = (algebraic[148]+algebraic[117]+constants[156])*states[21]
algebraic[156] = 4.00000*algebraic[35]
algebraic[211] = algebraic[109]*states[20]+algebraic[156]*states[22]+algebraic[197]*states[16]
rates[21] = algebraic[211]-algebraic[206]
algebraic[28] = constants[149]*constants[125]*exp(-196337./constants[145]+529.952/constants[2]+(2.78085*states[6])/constants[126])
algebraic[127] = algebraic[28]
algebraic[87] = constants[149]*constants[125]*exp(-147814./constants[145]+338.915/constants[2]+(2.13600*states[6])/constants[126])
algebraic[164] = algebraic[87]
algebraic[39] = constants[149]*constants[125]*exp(-133690./constants[145]+229.205/constants[2]+(-1.55804*states[6])/constants[126])
algebraic[209] = algebraic[39]
algebraic[97] = constants[149]*constants[125]*exp(-121322./constants[145]+193.265/constants[2]+(-1.74290*states[6])/constants[126])
algebraic[216] = algebraic[97]
rates[54] = (algebraic[123]*states[53]-(algebraic[202]+algebraic[127]+algebraic[164]+constants[173])*states[54])+algebraic[209]*states[45]+algebraic[216]*states[46]+constants[172]*states[59]
rates[46] = (algebraic[164]*states[54]+constants[158]*states[45])-(algebraic[216]+constants[159])*states[46]
algebraic[131] = algebraic[107]*constants[162]
algebraic[215] = algebraic[173]/constants[162]
rates[55] = (constants[165]*states[50]-(constants[164]+algebraic[131])*states[55])+algebraic[215]*states[56]
algebraic[204] = constants[48]*algebraic[197]
algebraic[214] = (algebraic[133]+constants[45]+algebraic[204])*states[17]
algebraic[219] = algebraic[79]*states[16]+constants[46]*states[23]+constants[161]*states[22]
rates[17] = algebraic[219]-algebraic[214]
algebraic[212] = (algebraic[156]+constants[161])*states[22]
algebraic[218] = algebraic[117]*states[21]+algebraic[204]*states[17]
rates[22] = algebraic[218]-algebraic[212]
algebraic[135] = algebraic[115]*constants[162]
algebraic[220] = algebraic[184]/constants[162]
rates[56] = ((constants[167]*states[51]+algebraic[131]*states[55])-(algebraic[215]+constants[166]+algebraic[135])*states[56])+algebraic[220]*states[57]
algebraic[140] = algebraic[119]*constants[162]
algebraic[222] = algebraic[195]/constants[162]
rates[57] = ((constants[169]*states[52]+algebraic[135]*states[56])-(algebraic[140]+constants[168]+algebraic[220])*states[57])+algebraic[222]*states[58]
algebraic[146] = algebraic[123]*constants[162]
algebraic[224] = algebraic[202]/constants[162]
rates[58] = ((constants[171]*states[53]+algebraic[140]*states[57])-(algebraic[222]+constants[170]+algebraic[146])*states[58])+algebraic[224]*states[59]
algebraic[64] = constants[149]*constants[125]*exp(116431./constants[145]+-578.317/constants[2]+(0.764126*states[6])/constants[126])
algebraic[154] = algebraic[64]
algebraic[75] = constants[149]*constants[125]*exp(-55700.7/constants[145]+-130.639/constants[2]+(-3.64982*states[6])/constants[126])
algebraic[226] = algebraic[75]
rates[59] = ((constants[173]*states[54]+algebraic[146]*states[58])-(algebraic[224]+constants[172]+algebraic[154])*states[59])+algebraic[226]*states[60]
algebraic[55] = constants[149]*constants[125]*exp(-97657.8/constants[145]+1.51000/constants[2]+(0.0684862*states[6])/constants[126])
algebraic[228] = algebraic[55]
algebraic[47] = constants[149]*constants[125]*exp(-62123.1/constants[145]+39.2950/constants[2]+(0.288816*states[6])/constants[126])
algebraic[230] = algebraic[47]
rates[45] = (algebraic[127]*states[54]-(algebraic[209]+constants[158]+algebraic[230])*states[45])+constants[159]*states[46]+algebraic[228]*states[60]
rates[60] = (algebraic[230]*states[45]+algebraic[154]*states[59])-(algebraic[226]+algebraic[228])*states[60]
algebraic[116] = constants[115]*log(constants[9]/states[0])
algebraic[128] = constants[82]*(states[45]+states[46])*(states[6]-algebraic[116])
algebraic[225] = constants[83]*(states[6]-algebraic[116])
algebraic[231] = states[6]/constants[115]
algebraic[237] = exp(constants[90]*algebraic[231])*(power(states[0], 3.00000))*constants[10]
algebraic[238] = exp((constants[90]-1.00000)*algebraic[231])*(power(constants[9], 3.00000))*states[2]
algebraic[239] = 1.00000+constants[89]*exp((constants[90]-1.00000)*algebraic[231])
algebraic[240] = (constants[86]*(algebraic[237]-algebraic[238]))/(constants[160]*algebraic[239]*constants[163])
algebraic[232] = 1.00000+0.124500*exp(-0.100000*algebraic[231])
algebraic[233] = 0.0365000*constants[151]*exp(-1.33000*algebraic[231])
algebraic[234] = 1.00000/(algebraic[232]+algebraic[233])
algebraic[235] = 1.00000+power(constants[92]/states[0], 1.50000)
algebraic[236] = (constants[91]*algebraic[234]*constants[155])/algebraic[235]
algebraic[147] = states[6]/constants[115]
algebraic[155] = 1000.00*constants[0]*algebraic[147]
algebraic[174] = exp(algebraic[147])-1.00000
algebraic[196] = states[0]*exp(algebraic[147])-constants[9]
algebraic[203] = (0.0200000*constants[121]*states[39]*algebraic[155]*algebraic[196])/algebraic[174]
rates[0] = -(algebraic[128]+algebraic[225]+3.00000*(algebraic[240]+algebraic[236])+algebraic[203])*constants[143]
algebraic[124] = 0.500000*constants[115]*log(constants[10]/states[2])
algebraic[227] = constants[96]*(states[6]-algebraic[124])
algebraic[229] = (constants[94]*states[2])/(constants[95]+states[2])
algebraic[58] = (states[3]-states[2])/constants[24]
algebraic[88] = constants[25]*algebraic[67]+constants[26]*algebraic[78]
algebraic[108] = 1.00000/(1.00000+(constants[31]*constants[34])/(power(states[2]+constants[34], 2.00000))+(constants[33]*constants[36])/(power(states[2]+constants[36], 2.00000)))
rates[2] = algebraic[108]*(((algebraic[58]-algebraic[42])-algebraic[88])-((algebraic[227]-2.00000*algebraic[240])+algebraic[229])*0.500000*constants[143])
algebraic[241] = states[6]/constants[115]
algebraic[242] = 1000.00*constants[0]*algebraic[241]
algebraic[243] = 0.00100000*exp(2.00000*algebraic[241])-constants[10]*0.341000
algebraic[244] = exp(2.00000*algebraic[241])-1.00000
algebraic[245] = (constants[144]*4.00000*algebraic[242]*algebraic[243])/algebraic[244]
algebraic[246] = algebraic[245]*states[24]*states[23]
algebraic[98] = 1.00000/(1.00000+(constants[31]*constants[34])/(power(states[3]+constants[34], 2.00000))+(constants[33]*constants[36])/(power(states[3]+constants[36], 2.00000)))
rates[3] = algebraic[98]*(((algebraic[50]*constants[5])/constants[7]-(algebraic[58]*constants[4])/constants[7])-algebraic[246]*constants[148])
algebraic[120] = constants[115]*log(constants[8]/states[1])
algebraic[132] = constants[79]*constants[122]*states[47]*(states[6]-algebraic[120])
algebraic[136] = constants[80]*(states[48]+states[49])*(states[6]-algebraic[120])
algebraic[221] = 1.00000/(0.940000+exp((1.26000/constants[115])*(states[6]-algebraic[120])))
algebraic[223] = constants[81]*(power(constants[8]/1.00000, 1.0/2))*algebraic[221]*(states[6]-algebraic[120])
algebraic[247] = custom_piecewise([greater_equal(algebraic[245] , 0.00000), 0.00000 , True, algebraic[245]])
algebraic[248] = constants[124]/(1.00000+algebraic[247]/constants[97])
algebraic[249] = states[1]*exp(algebraic[241])-constants[8]
algebraic[250] = exp(algebraic[241])-1.00000
algebraic[251] = (algebraic[248]*states[23]*states[24]*algebraic[242]*algebraic[249])/algebraic[250]
algebraic[0] = floor(voi/constants[11])*constants[11]
algebraic[11] = custom_piecewise([greater_equal(voi-algebraic[0] , constants[14]) & less_equal(voi-algebraic[0] , constants[14]+constants[12]), constants[13] , True, 0.00000])
algebraic[141] = constants[120]*states[29]*(states[6]-algebraic[120])
algebraic[165] = states[1]*exp(algebraic[147])-constants[8]
algebraic[185] = (constants[121]*states[39]*algebraic[155]*algebraic[165])/algebraic[174]
rates[1] = -(((algebraic[132]+algebraic[136]+algebraic[223]+algebraic[251]+algebraic[11])-2.00000*algebraic[236])+algebraic[141]+algebraic[185])*constants[143]
algebraic[210] = algebraic[185]+algebraic[203]
algebraic[217] = algebraic[141]+algebraic[210]
algebraic[252] = algebraic[128]+algebraic[246]+algebraic[251]+algebraic[132]+algebraic[136]+algebraic[223]+algebraic[240]+algebraic[236]+algebraic[217]+algebraic[229]+algebraic[227]+algebraic[225]+algebraic[11]
rates[6] = -algebraic[252]
return(rates)
def computeAlgebraic(constants, states, voi):
algebraic = array([[0.0] * len(voi)] * sizeAlgebraic)
states = array(states)
voi = array(voi)
algebraic[1] = -constants[37]*(power(states[3]*1000.00, constants[43]))*states[11]+constants[38]*states[7]
algebraic[12] = constants[39]*(power(states[3]*1000.00, constants[44]))*states[7]-constants[40]*states[8]
algebraic[3] = (constants[117]/(1.00000+exp((states[6]+28.5000)/7.80000))+1.00000)-constants[117]
algebraic[14] = 1.00000/(0.00336336/(0.500000+exp(states[6]/-5.53900))+0.00779047*exp(states[6]/-49.5104))
algebraic[8] = constants[100]*constants[101]*exp(constants[102]*states[6])
algebraic[19] = constants[100]*constants[103]*exp(constants[104]*states[6])
algebraic[23] = constants[41]*states[7]-constants[42]*states[12]
algebraic[32] = 0.216256*exp(-1.88912e-05*states[6])
algebraic[21] = 0.00700807*exp(-0.149998*states[6])
algebraic[10] = 0.00767254*exp(0.0866295*states[6])
algebraic[31] = 0.00379738*exp(-0.0142567*states[6])
algebraic[34] = -(algebraic[1]+algebraic[12]+algebraic[23])
algebraic[30] = algebraic[8]*states[61]+constants[128]*states[63]
algebraic[41] = (algebraic[19]+constants[127])*states[62]
algebraic[9] = constants[100]*constants[105]*exp(constants[106]*states[6])
algebraic[40] = constants[100]*constants[111]*exp(constants[112]*states[6])
algebraic[49] = algebraic[9]*states[63]+algebraic[40]*states[64]
algebraic[20] = constants[100]*constants[107]*exp(constants[108]*states[6])
algebraic[29] = constants[100]*constants[109]*exp(constants[110]*states[6])
algebraic[57] = (algebraic[20]+algebraic[29])*states[47]
algebraic[22] = power(states[2]/constants[15], constants[16])
algebraic[33] = power(states[5]/constants[17], constants[18])
algebraic[42] = (constants[19]*(constants[20]*algebraic[22]-constants[21]*algebraic[33]))/(1.00000+algebraic[22]+algebraic[33])
algebraic[59] = (states[5]-states[4])/constants[23]
algebraic[50] = constants[22]*(states[7]+states[8])*(states[4]-states[3])
algebraic[68] = 1.00000/(1.00000+(constants[32]*constants[35])/(power(states[4]+constants[35], 2.00000)))
algebraic[67] = constants[29]*states[2]*(1.00000-states[9])-constants[30]*states[9]
algebraic[48] = constants[100]*constants[113]*exp(constants[114]*states[6])
algebraic[56] = (algebraic[20]*algebraic[40]*algebraic[48])/(algebraic[9]*algebraic[29])
algebraic[65] = constants[127]*states[62]+algebraic[20]*states[47]+algebraic[56]*states[64]
algebraic[76] = (algebraic[48]+algebraic[9]+constants[128])*states[63]
algebraic[66] = algebraic[48]*states[63]+algebraic[29]*states[47]
algebraic[77] = (algebraic[56]+algebraic[40])*states[64]
algebraic[78] = constants[27]*states[2]*(1.00000-states[10])-constants[28]*states[10]
algebraic[5] = constants[73]*exp(constants[74]*states[6])
algebraic[26] = 4.00000*algebraic[5]
algebraic[73] = (algebraic[26]+constants[132])*states[35]
algebraic[16] = constants[75]*exp(-constants[76]*states[6])
algebraic[62] = algebraic[16]
algebraic[83] = algebraic[62]*states[36]+constants[137]*states[40]
algebraic[27] = 4.00000*constants[69]*algebraic[5]
algebraic[74] = (constants[137]+algebraic[27])*states[40]
algebraic[63] = algebraic[16]/constants[65]
algebraic[85] = constants[132]*states[35]+algebraic[63]*states[41]
algebraic[37] = 3.00000*algebraic[5]
algebraic[84] = (algebraic[37]+algebraic[62]+constants[133])*states[36]
algebraic[71] = 2.00000*algebraic[16]
algebraic[93] = algebraic[71]*states[37]+constants[138]*states[41]+algebraic[26]*states[35]
algebraic[38] = (3.00000*constants[70]*algebraic[5])/constants[69]
algebraic[86] = (algebraic[38]+constants[138]+algebraic[63])*states[41]
algebraic[72] = (2.00000*constants[65]*algebraic[16])/constants[66]
algebraic[95] = algebraic[72]*states[42]+constants[133]*states[36]+algebraic[27]*states[40]
algebraic[45] = 2.00000*algebraic[5]
algebraic[94] = (algebraic[45]+algebraic[71]+constants[134])*states[37]
algebraic[81] = 3.00000*algebraic[16]
algebraic[101] = algebraic[81]*states[38]+constants[139]*states[42]+algebraic[37]*states[36]
algebraic[46] = (2.00000*constants[71]*algebraic[5])/constants[70]
algebraic[96] = (algebraic[46]+constants[139]+algebraic[72])*states[42]
algebraic[82] = (3.00000*constants[66]*algebraic[16])/constants[67]
algebraic[104] = algebraic[82]*states[43]+constants[134]*states[37]+algebraic[38]*states[41]
algebraic[53] = algebraic[5]
algebraic[102] = (algebraic[53]+algebraic[81]+constants[135])*states[38]
algebraic[91] = 4.00000*algebraic[16]
algebraic[111] = algebraic[91]*states[39]+constants[140]*states[43]+algebraic[45]*states[37]
algebraic[103] = (algebraic[91]+constants[136])*states[39]
algebraic[112] = algebraic[53]*states[38]+constants[141]*states[44]
algebraic[54] = (constants[72]*algebraic[5])/constants[71]
algebraic[105] = (algebraic[54]+constants[140]+algebraic[82])*states[43]
algebraic[92] = (4.00000*constants[67]*algebraic[16])/constants[68]
algebraic[113] = algebraic[92]*states[44]+constants[135]*states[38]+algebraic[46]*states[42]
algebraic[106] = (constants[141]+algebraic[92])*states[44]
algebraic[114] = constants[136]*states[39]+algebraic[54]*states[43]
algebraic[4] = constants[49]*exp(constants[50]*states[6])
algebraic[44] = 4.00000*algebraic[4]
algebraic[36] = constants[55]*exp(constants[56]*states[6])
algebraic[138] = algebraic[36]
algebraic[152] = (algebraic[44]+algebraic[138])*states[25]
algebraic[15] = constants[51]*exp(-constants[52]*states[6])
algebraic[118] = algebraic[15]
algebraic[25] = constants[53]*exp(-constants[54]*states[6])
algebraic[145] = algebraic[25]
algebraic[160] = algebraic[118]*states[26]+algebraic[145]*states[30]
algebraic[80] = 4.00000*constants[61]*algebraic[4]
algebraic[153] = (algebraic[145]+algebraic[80])*states[30]
algebraic[134] = algebraic[15]/constants[57]
algebraic[162] = algebraic[138]*states[25]+algebraic[134]*states[31]
algebraic[52] = 3.00000*algebraic[4]
algebraic[143] = constants[57]*algebraic[36]
algebraic[161] = (algebraic[52]+algebraic[118]+algebraic[143])*states[26]
algebraic[122] = 2.00000*algebraic[15]
algebraic[151] = algebraic[25]/constants[61]
algebraic[169] = algebraic[122]*states[27]+algebraic[151]*states[31]+algebraic[44]*states[25]
algebraic[90] = (3.00000*constants[62]*algebraic[4])/constants[61]
algebraic[163] = (algebraic[90]+algebraic[151]+algebraic[134])*states[31]
algebraic[139] = (2.00000*constants[57]*algebraic[15])/constants[58]
algebraic[171] = algebraic[139]*states[32]+algebraic[143]*states[26]+algebraic[80]*states[30]
algebraic[7] = constants[149]*constants[125]*exp(-114007./constants[145]+224.114/constants[2]+(0.286374*states[6])/constants[126])
algebraic[107] = 4.00000*algebraic[7]
algebraic[18] = constants[149]*constants[125]*exp(-272470./constants[145]+708.146/constants[2]+(-2.28528*states[6])/constants[126])
algebraic[173] = algebraic[18]
algebraic[61] = 2.00000*algebraic[4]
algebraic[149] = constants[58]*algebraic[36]
algebraic[170] = (algebraic[61]+algebraic[122]+algebraic[149])*states[27]
algebraic[126] = 3.00000*algebraic[15]
algebraic[159] = algebraic[25]/constants[62]
algebraic[178] = algebraic[126]*states[28]+algebraic[159]*states[32]+algebraic[52]*states[26]
algebraic[100] = (2.00000*constants[63]*algebraic[4])/constants[62]
algebraic[172] = (algebraic[100]+algebraic[159]+algebraic[139])*states[32]
algebraic[144] = (3.00000*constants[58]*algebraic[15])/constants[59]
algebraic[181] = algebraic[144]*states[33]+algebraic[149]*states[27]+algebraic[90]*states[31]
algebraic[115] = 3.00000*algebraic[7]
algebraic[184] = 2.00000*algebraic[18]
algebraic[2] = 4.00000*1.20000*0.416000*exp(0.0120000*(states[6]-35.0000))
algebraic[51] = 4.00000*algebraic[2]
algebraic[43] = 0.600000*0.0923300*states[3]
algebraic[166] = algebraic[43]
algebraic[177] = (algebraic[51]+algebraic[166])*states[13]
algebraic[13] = 4.00000*0.450000*0.0490000*exp(-0.0650000*(states[6]-22.0000))
algebraic[121] = algebraic[13]
algebraic[189] = algebraic[121]*states[14]+constants[131]*states[18]
algebraic[24] = constants[48]*algebraic[2]
algebraic[89] = 4.00000*algebraic[24]
algebraic[176] = (algebraic[89]+constants[131])*states[18]
algebraic[35] = algebraic[13]/constants[47]
algebraic[137] = algebraic[35]
algebraic[187] = algebraic[137]*states[19]+algebraic[166]*states[13]
algebraic[70] = algebraic[4]
algebraic[157] = constants[59]*algebraic[36]
algebraic[179] = (algebraic[70]+algebraic[126]+algebraic[157])*states[28]
algebraic[130] = 4.00000*algebraic[15]
algebraic[167] = algebraic[25]/constants[63]
algebraic[191] = algebraic[130]*states[29]+algebraic[167]*states[33]+algebraic[61]*states[27]
algebraic[158] = constants[60]*algebraic[36]
algebraic[180] = (algebraic[130]+algebraic[158])*states[29]
algebraic[168] = algebraic[25]/constants[64]
algebraic[192] = algebraic[70]*states[28]+algebraic[168]*states[34]
algebraic[110] = (constants[64]*algebraic[4])/constants[63]
algebraic[182] = (algebraic[110]+algebraic[167]+algebraic[144])*states[33]
algebraic[150] = (4.00000*constants[59]*algebraic[15])/constants[60]
algebraic[193] = algebraic[150]*states[34]+algebraic[157]*states[28]+algebraic[100]*states[32]
algebraic[183] = (algebraic[168]+algebraic[150])*states[34]
algebraic[194] = algebraic[158]*states[29]+algebraic[110]*states[33]
algebraic[119] = 2.00000*algebraic[7]
algebraic[195] = 3.00000*algebraic[18]
algebraic[60] = 3.00000*algebraic[2]
algebraic[175] = constants[48]*algebraic[166]
algebraic[190] = (algebraic[121]+algebraic[60]+algebraic[175])*states[14]
algebraic[125] = 2.00000*algebraic[13]
algebraic[200] = algebraic[51]*states[13]+algebraic[125]*states[15]+constants[147]*states[19]
algebraic[99] = 3.00000*algebraic[24]
algebraic[188] = (algebraic[137]+algebraic[99]+constants[147])*states[19]
algebraic[142] = 2.00000*algebraic[35]
algebraic[198] = algebraic[89]*states[18]+algebraic[142]*states[20]+algebraic[175]*states[14]
algebraic[123] = algebraic[7]
algebraic[202] = 4.00000*algebraic[18]
algebraic[69] = 2.00000*algebraic[2]
algebraic[186] = constants[48]*algebraic[175]
algebraic[201] = (algebraic[125]+algebraic[69]+algebraic[186])*states[15]
algebraic[129] = 3.00000*algebraic[13]
algebraic[207] = algebraic[60]*states[14]+algebraic[129]*states[16]+constants[150]*states[20]
algebraic[109] = 2.00000*algebraic[24]
algebraic[199] = (algebraic[142]+algebraic[109]+constants[150])*states[20]
algebraic[148] = 3.00000*algebraic[35]
algebraic[205] = algebraic[99]*states[19]+algebraic[148]*states[21]+algebraic[186]*states[15]
algebraic[79] = algebraic[2]
algebraic[197] = constants[48]*algebraic[186]
algebraic[208] = (algebraic[129]+algebraic[79]+algebraic[197])*states[16]
algebraic[133] = 4.00000*algebraic[13]
algebraic[213] = algebraic[69]*states[15]+algebraic[133]*states[17]+constants[156]*states[21]
algebraic[117] = algebraic[24]
algebraic[206] = (algebraic[148]+algebraic[117]+constants[156])*states[21]
algebraic[156] = 4.00000*algebraic[35]
algebraic[211] = algebraic[109]*states[20]+algebraic[156]*states[22]+algebraic[197]*states[16]
algebraic[28] = constants[149]*constants[125]*exp(-196337./constants[145]+529.952/constants[2]+(2.78085*states[6])/constants[126])
algebraic[127] = algebraic[28]
algebraic[87] = constants[149]*constants[125]*exp(-147814./constants[145]+338.915/constants[2]+(2.13600*states[6])/constants[126])
algebraic[164] = algebraic[87]
algebraic[39] = constants[149]*constants[125]*exp(-133690./constants[145]+229.205/constants[2]+(-1.55804*states[6])/constants[126])
algebraic[209] = algebraic[39]
algebraic[97] = constants[149]*constants[125]*exp(-121322./constants[145]+193.265/constants[2]+(-1.74290*states[6])/constants[126])
algebraic[216] = algebraic[97]
algebraic[131] = algebraic[107]*constants[162]
algebraic[215] = algebraic[173]/constants[162]
algebraic[204] = constants[48]*algebraic[197]
algebraic[214] = (algebraic[133]+constants[45]+algebraic[204])*states[17]
algebraic[219] = algebraic[79]*states[16]+constants[46]*states[23]+constants[161]*states[22]
algebraic[212] = (algebraic[156]+constants[161])*states[22]
algebraic[218] = algebraic[117]*states[21]+algebraic[204]*states[17]
algebraic[135] = algebraic[115]*constants[162]
algebraic[220] = algebraic[184]/constants[162]
algebraic[140] = algebraic[119]*constants[162]
algebraic[222] = algebraic[195]/constants[162]
algebraic[146] = algebraic[123]*constants[162]
algebraic[224] = algebraic[202]/constants[162]
algebraic[64] = constants[149]*constants[125]*exp(116431./constants[145]+-578.317/constants[2]+(0.764126*states[6])/constants[126])
algebraic[154] = algebraic[64]
algebraic[75] = constants[149]*constants[125]*exp(-55700.7/constants[145]+-130.639/constants[2]+(-3.64982*states[6])/constants[126])
algebraic[226] = algebraic[75]
algebraic[55] = constants[149]*constants[125]*exp(-97657.8/constants[145]+1.51000/constants[2]+(0.0684862*states[6])/constants[126])
algebraic[228] = algebraic[55]
algebraic[47] = constants[149]*constants[125]*exp(-62123.1/constants[145]+39.2950/constants[2]+(0.288816*states[6])/constants[126])
algebraic[230] = algebraic[47]
algebraic[116] = constants[115]*log(constants[9]/states[0])
algebraic[128] = constants[82]*(states[45]+states[46])*(states[6]-algebraic[116])
algebraic[225] = constants[83]*(states[6]-algebraic[116])
algebraic[231] = states[6]/constants[115]
algebraic[237] = exp(constants[90]*algebraic[231])*(power(states[0], 3.00000))*constants[10]
algebraic[238] = exp((constants[90]-1.00000)*algebraic[231])*(power(constants[9], 3.00000))*states[2]
algebraic[239] = 1.00000+constants[89]*exp((constants[90]-1.00000)*algebraic[231])
algebraic[240] = (constants[86]*(algebraic[237]-algebraic[238]))/(constants[160]*algebraic[239]*constants[163])
algebraic[232] = 1.00000+0.124500*exp(-0.100000*algebraic[231])
algebraic[233] = 0.0365000*constants[151]*exp(-1.33000*algebraic[231])
algebraic[234] = 1.00000/(algebraic[232]+algebraic[233])
algebraic[235] = 1.00000+power(constants[92]/states[0], 1.50000)
algebraic[236] = (constants[91]*algebraic[234]*constants[155])/algebraic[235]
algebraic[147] = states[6]/constants[115]
algebraic[155] = 1000.00*constants[0]*algebraic[147]
algebraic[174] = exp(algebraic[147])-1.00000
algebraic[196] = states[0]*exp(algebraic[147])-constants[9]
algebraic[203] = (0.0200000*constants[121]*states[39]*algebraic[155]*algebraic[196])/algebraic[174]
algebraic[124] = 0.500000*constants[115]*log(constants[10]/states[2])
algebraic[227] = constants[96]*(states[6]-algebraic[124])
algebraic[229] = (constants[94]*states[2])/(constants[95]+states[2])
algebraic[58] = (states[3]-states[2])/constants[24]
algebraic[88] = constants[25]*algebraic[67]+constants[26]*algebraic[78]
algebraic[108] = 1.00000/(1.00000+(constants[31]*constants[34])/(power(states[2]+constants[34], 2.00000))+(constants[33]*constants[36])/(power(states[2]+constants[36], 2.00000)))
algebraic[241] = states[6]/constants[115]
algebraic[242] = 1000.00*constants[0]*algebraic[241]
algebraic[243] = 0.00100000*exp(2.00000*algebraic[241])-constants[10]*0.341000
algebraic[244] = exp(2.00000*algebraic[241])-1.00000
algebraic[245] = (constants[144]*4.00000*algebraic[242]*algebraic[243])/algebraic[244]
algebraic[246] = algebraic[245]*states[24]*states[23]
algebraic[98] = 1.00000/(1.00000+(constants[31]*constants[34])/(power(states[3]+constants[34], 2.00000))+(constants[33]*constants[36])/(power(states[3]+constants[36], 2.00000)))
algebraic[120] = constants[115]*log(constants[8]/states[1])
algebraic[132] = constants[79]*constants[122]*states[47]*(states[6]-algebraic[120])
algebraic[136] = constants[80]*(states[48]+states[49])*(states[6]-algebraic[120])
algebraic[221] = 1.00000/(0.940000+exp((1.26000/constants[115])*(states[6]-algebraic[120])))
algebraic[223] = constants[81]*(power(constants[8]/1.00000, 1.0/2))*algebraic[221]*(states[6]-algebraic[120])
algebraic[247] = custom_piecewise([greater_equal(algebraic[245] , 0.00000), 0.00000 , True, algebraic[245]])
algebraic[248] = constants[124]/(1.00000+algebraic[247]/constants[97])
algebraic[249] = states[1]*exp(algebraic[241])-constants[8]
algebraic[250] = exp(algebraic[241])-1.00000
algebraic[251] = (algebraic[248]*states[23]*states[24]*algebraic[242]*algebraic[249])/algebraic[250]
algebraic[0] = floor(voi/constants[11])*constants[11]
algebraic[11] = custom_piecewise([greater_equal(voi-algebraic[0] , constants[14]) & less_equal(voi-algebraic[0] , constants[14]+constants[12]), constants[13] , True, 0.00000])
algebraic[141] = constants[120]*states[29]*(states[6]-algebraic[120])
algebraic[165] = states[1]*exp(algebraic[147])-constants[8]
algebraic[185] = (constants[121]*states[39]*algebraic[155]*algebraic[165])/algebraic[174]
algebraic[210] = algebraic[185]+algebraic[203]
algebraic[217] = algebraic[141]+algebraic[210]
algebraic[252] = algebraic[128]+algebraic[246]+algebraic[251]+algebraic[132]+algebraic[136]+algebraic[223]+algebraic[240]+algebraic[236]+algebraic[217]+algebraic[229]+algebraic[227]+algebraic[225]+algebraic[11]
algebraic[6] = states[1]+0.0183300*states[0]
algebraic[17] = constants[115]*log(constants[142]/algebraic[6])
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)