Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/* There are a total of 256 entries in the algebraic variable array. There are a total of 67 entries in each of the rate and state variable arrays. There are a total of 186 entries in the constant variable array. */ /* * VOI is time in component environment (ms). * CONSTANTS[0] is iso in component environment (dimensionless). * CONSTANTS[1] is CSQN2 in component environment (dimensionless). * CONSTANTS[2] is RyR2 in component environment (dimensionless). * CONSTANTS[155] is a1 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb). * CONSTANTS[160] is a2 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb). * CONSTANTS[3] is Faraday in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (coulomb_per_millimole). * CONSTANTS[4] is Temp in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (kelvin). * CONSTANTS[5] is Rgas in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (joule_per_mole_kelvin). * CONSTANTS[124] is RT_over_F in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV). * CONSTANTS[6] is Acap in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (cm2). * CONSTANTS[135] is C in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mF). * CONSTANTS[7] is Vmyo in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL). * CONSTANTS[8] is VJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL). * CONSTANTS[9] is VNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL). * CONSTANTS[10] is VSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL). * STATES[0] is Nai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * STATES[1] is Ki in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * STATES[2] is Cai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * STATES[3] is CaSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * STATES[4] is CaJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * STATES[5] is CaNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * STATES[6] is V in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV). * ALGEBRAIC[162] is INa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[177] is INab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[189] is INaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * ALGEBRAIC[185] is INaK in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * ALGEBRAIC[172] is IKv14_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[163] is IKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[164] is IKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[176] is IK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[200] is ICaK in component COMPUTE_ICa_ICaK (uA_per_uF). * ALGEBRAIC[195] is ICa in component COMPUTE_ICa_ICaK (uA_per_uF). * ALGEBRAIC[174] is Ito1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[1] is i_Stim in component I_stimulus (uA_per_uF). * ALGEBRAIC[165] is IKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[170] is IKv14_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * ALGEBRAIC[178] is ICab in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * ALGEBRAIC[179] is IpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * ALGEBRAIC[7] is Jxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms). * ALGEBRAIC[4] is Jup in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms). * ALGEBRAIC[10] is Jtrpn in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM_per_ms). * ALGEBRAIC[5] is Jrel in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms). * ALGEBRAIC[6] is Jtr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms). * ALGEBRAIC[11] is beta_SS in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * ALGEBRAIC[12] is beta_JSR in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * ALGEBRAIC[13] is beta_i in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * ALGEBRAIC[255] is i_tot in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uA_per_uF). * CONSTANTS[11] is Ko in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * CONSTANTS[12] is Nao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * CONSTANTS[13] is Cao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * CONSTANTS[14] is stim_period in component I_stimulus (ms). * CONSTANTS[15] is stim_duration in component I_stimulus (ms). * CONSTANTS[16] is stim_amplitude in component I_stimulus (uA_per_uF). * CONSTANTS[17] is stim_offset in component I_stimulus (ms). * ALGEBRAIC[0] is past in component I_stimulus (ms). * ALGEBRAIC[2] is fb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless). * CONSTANTS[18] is Kfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM). * CONSTANTS[19] is Nfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless). * ALGEBRAIC[3] is rb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless). * CONSTANTS[20] is Krb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM). * CONSTANTS[21] is Nrb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless). * CONSTANTS[22] is KSR in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM). * CONSTANTS[23] is vmaxf in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms). * CONSTANTS[24] is vmaxr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms). * CONSTANTS[25] is v1 in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms). * STATES[7] is O1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * STATES[8] is O2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[26] is tautr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms). * CONSTANTS[27] is tauxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms). * CONSTANTS[28] is LTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[29] is HTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[30] is khtrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms). * CONSTANTS[31] is khtrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms). * CONSTANTS[32] is kltrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms). * CONSTANTS[33] is kltrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms). * CONSTANTS[34] is CMDNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[125] is CSQNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[35] is EGTAtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[36] is KmCMDN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[37] is KmCSQN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * CONSTANTS[38] is KmEGTA in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM). * ALGEBRAIC[8] is dLTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms). * ALGEBRAIC[9] is dHTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms). * STATES[9] is LTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * STATES[10] is HTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * CONSTANTS[39] is kaplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM4_per_ms). * CONSTANTS[40] is kaminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * CONSTANTS[41] is kbplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM3_per_ms). * CONSTANTS[42] is kbminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * CONSTANTS[43] is kcplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * CONSTANTS[44] is kcminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * CONSTANTS[45] is ncoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[46] is mcoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[47] is kryr in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * ALGEBRAIC[14] is klumenC1O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * ALGEBRAIC[15] is klumenO1C2 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * ALGEBRAIC[16] is klumenC2O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[48] is HmaxC1O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[126] is HmaxO1C2 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[49] is HmaxC2O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[50] is Hmin in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[51] is H50C1O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (mM). * CONSTANTS[52] is H50O1C2 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (mM). * CONSTANTS[127] is H50C2O1 in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (mM). * CONSTANTS[53] is HN in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * ALGEBRAIC[17] is dC1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * ALGEBRAIC[18] is dO2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * ALGEBRAIC[19] is dC2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * ALGEBRAIC[20] is dO1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms). * STATES[11] is C1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * STATES[12] is C2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * CONSTANTS[54] is fL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[55] is gL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[56] is bL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * CONSTANTS[57] is aL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * ALGEBRAIC[26] is C0_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[27] is C1_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[28] is C2_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[29] is C3_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[34] is C1_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[35] is C2_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[36] is C3_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[37] is C4_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[30] is CCa0_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[31] is CCa1_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[32] is CCa2_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[33] is CCa3_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[38] is CCa1_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[39] is CCa2_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[40] is CCa3_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[41] is CCa4_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[42] is C0_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[43] is C1_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[44] is C2_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[45] is C3_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[46] is C4_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[143] is CCa0_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[159] is CCa1_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[162] is CCa2_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[168] is CCa3_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[173] is CCa4_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[21] is alpha in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[22] is beta in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[23] is alpha_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[24] is beta_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[25] is gamma in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * CONSTANTS[128] is omega in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[57] is a1_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[58] is a2_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[59] is a1_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[60] is a2_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[61] is a1_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[62] is a2_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[63] is a1_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[64] is a2_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[65] is a1_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[66] is a2_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[47] is a1_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[48] is a2_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[49] is a1_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[50] is a2_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[51] is a1_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[52] is a2_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[53] is a1_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[54] is a2_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[55] is a1_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * ALGEBRAIC[56] is a2_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms). * STATES[13] is C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[14] is C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[15] is C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[16] is C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[17] is C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[18] is CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[19] is CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[20] is CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[21] is CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[22] is CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[23] is Open in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * STATES[24] is yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * CONSTANTS[129] is a1_Cainf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * ALGEBRAIC[67] is yCa_inf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * ALGEBRAIC[68] is tau_yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (ms). * ALGEBRAIC[69] is alpha_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[70] is beta_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[71] is alpha_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[72] is beta_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[73] is C0Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[74] is C1Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[75] is C2Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[76] is C3Kv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[77] is CI0Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[78] is CI1Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[79] is CI2Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[80] is CI3Kv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[81] is C1Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[82] is C2Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[83] is C3Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[84] is OKv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[85] is CI1Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[86] is CI2Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[87] is CI3Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[88] is OIKv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[89] is C0Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[90] is C1Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[91] is C2Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[92] is C3Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[93] is OKv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[94] is CI0Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[95] is CI1Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[96] is CI2Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[97] is CI3Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[98] is OIKv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * STATES[25] is C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[26] is C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[27] is C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[28] is C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[29] is OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[30] is CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[31] is CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[32] is CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[33] is CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * STATES[34] is OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * ALGEBRAIC[99] is a1_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[100] is a2_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[101] is a1_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[102] is a2_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[103] is a1_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[104] is a2_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[105] is a1_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[106] is a2_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[107] is a1_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[108] is a2_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[109] is a1_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[110] is a2_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[111] is a1_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[112] is a2_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[113] is a1_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[114] is a2_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[115] is a1_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[116] is a2_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[117] is a1_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * ALGEBRAIC[118] is a2_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * CONSTANTS[58] is alphaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * CONSTANTS[59] is aaKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV). * CONSTANTS[60] is betaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * CONSTANTS[61] is baKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV). * CONSTANTS[62] is alphai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * CONSTANTS[63] is aiKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV). * CONSTANTS[64] is betai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms). * CONSTANTS[65] is biKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV). * CONSTANTS[66] is f1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[67] is f2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[68] is f3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[69] is f4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[70] is b1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[71] is b2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[72] is b3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[73] is b4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * CONSTANTS[74] is f1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[75] is f2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[76] is f3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[77] is f4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[78] is b1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[79] is b2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[80] is b3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[81] is b4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * CONSTANTS[82] is alphaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[83] is aaKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV). * CONSTANTS[84] is betaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[85] is baKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV). * CONSTANTS[86] is alphai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[87] is betai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[119] is alpha_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[120] is beta_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[130] is alpha_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[131] is beta_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[121] is C0Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[122] is C1Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[123] is C2Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[124] is C3Kv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[125] is CI0Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[126] is CI1Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[127] is CI2Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[128] is CI3Kv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[129] is C1Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[130] is C2Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[131] is C3Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[132] is OKv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[133] is CI1Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[134] is CI2Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[135] is CI3Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[136] is OIKv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[144] is C0Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[145] is C1Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[146] is C2Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[147] is C3Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[148] is OKv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[149] is CI0Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[150] is CI1Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[151] is CI2Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[152] is CI3Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * CONSTANTS[153] is OIKv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[137] is a1_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[138] is a2_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[35] is C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[139] is a1_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[140] is a2_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[36] is C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[141] is a1_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[142] is a2_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[37] is C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[143] is a1_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[144] is a2_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[38] is C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[145] is a1_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[146] is a2_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[39] is OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[147] is a1_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[148] is a2_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[40] is CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[149] is a1_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[150] is a2_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[41] is CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[151] is a1_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[152] is a2_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[42] is CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[153] is a1_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[154] is a2_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[43] is CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[155] is a1_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * ALGEBRAIC[156] is a2_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms). * STATES[44] is OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * ALGEBRAIC[157] is ENa in component COMPUTE_REVERSAL_POTENTIALS (mV). * ALGEBRAIC[158] is EK in component COMPUTE_REVERSAL_POTENTIALS (mV). * ALGEBRAIC[160] is EKs in component COMPUTE_REVERSAL_POTENTIALS (mV). * ALGEBRAIC[161] is ECa in component COMPUTE_REVERSAL_POTENTIALS (mV). * CONSTANTS[154] is a1 in component COMPUTE_REVERSAL_POTENTIALS (mM). * ALGEBRAIC[159] is a2 in component COMPUTE_REVERSAL_POTENTIALS (mM). * CONSTANTS[88] is GKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF). * CONSTANTS[89] is GKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF). * CONSTANTS[90] is GK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF). * CONSTANTS[91] is GNa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF). * CONSTANTS[92] is GNab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF). * CONSTANTS[93] is KvScale in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless). * CONSTANTS[94] is Kv43Frac in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless). * CONSTANTS[132] is GKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF). * CONSTANTS[133] is PKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (litre_per_farad_second). * STATES[45] is na6 in component INa (dimensionless). * STATES[46] is na7 in component INa (dimensionless). * STATES[47] is OHerg in component IKr (dimensionless). * STATES[48] is O1ks in component IKs (dimensionless). * STATES[49] is O2ks in component IKs (dimensionless). * CONSTANTS[134] is fKo in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless). * ALGEBRAIC[166] is VF_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless). * ALGEBRAIC[167] is VFsq_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (coulomb_per_millimole). * ALGEBRAIC[168] is a1_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM). * ALGEBRAIC[169] is a2 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless). * ALGEBRAIC[171] is a1_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM). * ALGEBRAIC[175] is K1_inf in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless). * ALGEBRAIC[173] is IKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF). * CONSTANTS[95] is kNaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * CONSTANTS[96] is KmNa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM). * CONSTANTS[97] is KmCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM). * CONSTANTS[98] is ksat in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * CONSTANTS[99] is eta in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * CONSTANTS[100] is INaKmax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * CONSTANTS[101] is KmNai in component COMPUTE_INaK_INaCa_ICab_IpCa (mM). * CONSTANTS[102] is KmKo in component COMPUTE_INaK_INaCa_ICab_IpCa (mM). * CONSTANTS[103] is IpCamax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF). * CONSTANTS[104] is KmpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM). * CONSTANTS[105] is GCab in component COMPUTE_INaK_INaCa_ICab_IpCa (mS_per_uF). * ALGEBRAIC[180] is VF_over_RT in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * CONSTANTS[163] is sigma in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * ALGEBRAIC[181] is a1_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * ALGEBRAIC[182] is a2_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * ALGEBRAIC[183] is fNaK in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * CONSTANTS[167] is a1_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * ALGEBRAIC[184] is a2_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * ALGEBRAIC[186] is a1_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4). * ALGEBRAIC[187] is a2_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4). * ALGEBRAIC[188] is a3_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless). * CONSTANTS[172] is a4_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM). * CONSTANTS[175] is a5_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM3). * ALGEBRAIC[197] is PKprime in component COMPUTE_ICa_ICaK (litre_per_farad_second). * ALGEBRAIC[190] is VF_over_RT in component COMPUTE_ICa_ICaK (dimensionless). * ALGEBRAIC[191] is VFsq_over_RT in component COMPUTE_ICa_ICaK (coulomb_per_millimole). * ALGEBRAIC[192] is a1_Ca in component COMPUTE_ICa_ICaK (mM). * ALGEBRAIC[193] is a2_Ca in component COMPUTE_ICa_ICaK (dimensionless). * ALGEBRAIC[198] is a1_K in component COMPUTE_ICa_ICaK (mM). * ALGEBRAIC[199] is a2_K in component COMPUTE_ICa_ICaK (dimensionless). * ALGEBRAIC[194] is ICamax in component COMPUTE_ICa_ICaK (uA_per_uF). * ALGEBRAIC[196] is Icabar in component COMPUTE_ICa_ICaK (uA_per_uF). * CONSTANTS[156] is PCa in component COMPUTE_ICa_ICaK (litre_per_farad_second). * CONSTANTS[136] is PK in component COMPUTE_ICa_ICaK (litre_per_farad_second). * CONSTANTS[106] is ICahalf in component COMPUTE_ICa_ICaK (uA_per_uF). * CONSTANTS[107] is Pscale in component COMPUTE_ICa_ICaK (dimensionless). * CONSTANTS[161] is Temp_Scale in component INa (dimensionless). * ALGEBRAIC[201] is alpha1 in component INa (per_ms). * ALGEBRAIC[202] is beta1 in component INa (per_ms). * ALGEBRAIC[203] is gamma1 in component INa (per_ms). * ALGEBRAIC[204] is Delta1 in component INa (per_ms). * ALGEBRAIC[205] is On in component INa (per_ms). * ALGEBRAIC[206] is Of in component INa (per_ms). * ALGEBRAIC[207] is GammaGamma in component INa (per_ms). * ALGEBRAIC[208] is DeltaDelta in component INa (per_ms). * CONSTANTS[165] is epsilon in component INa (per_ms). * CONSTANTS[166] is omega_na in component INa (per_ms). * ALGEBRAIC[209] is rho in component INa (per_ms). * ALGEBRAIC[210] is mu in component INa (per_ms). * CONSTANTS[164] is Cn in component INa (per_ms). * CONSTANTS[169] is Cf in component INa (per_ms). * CONSTANTS[174] is parameter_a in component INa (dimensionless). * ALGEBRAIC[211] is k12 in component INa (per_ms). * ALGEBRAIC[212] is k23 in component INa (per_ms). * ALGEBRAIC[213] is k34 in component INa (per_ms). * ALGEBRAIC[214] is k45 in component INa (per_ms). * ALGEBRAIC[215] is k56 in component INa (per_ms). * CONSTANTS[170] is k67 in component INa (per_ms). * ALGEBRAIC[216] is k89 in component INa (per_ms). * ALGEBRAIC[217] is k910 in component INa (per_ms). * ALGEBRAIC[218] is k1011 in component INa (per_ms). * ALGEBRAIC[219] is k1112 in component INa (per_ms). * ALGEBRAIC[220] is k1213 in component INa (per_ms). * ALGEBRAIC[221] is k57 in component INa (per_ms). * ALGEBRAIC[222] is k21 in component INa (per_ms). * ALGEBRAIC[223] is k32 in component INa (per_ms). * ALGEBRAIC[224] is k43 in component INa (per_ms). * ALGEBRAIC[225] is k54 in component INa (per_ms). * ALGEBRAIC[226] is k65 in component INa (per_ms). * CONSTANTS[171] is k76 in component INa (per_ms). * ALGEBRAIC[227] is k98 in component INa (per_ms). * ALGEBRAIC[228] is k109 in component INa (per_ms). * ALGEBRAIC[229] is k1110 in component INa (per_ms). * ALGEBRAIC[230] is k1211 in component INa (per_ms). * ALGEBRAIC[231] is k1312 in component INa (per_ms). * ALGEBRAIC[232] is k75 in component INa (per_ms). * CONSTANTS[176] is k81 in component INa (per_ms). * CONSTANTS[178] is k92 in component INa (per_ms). * CONSTANTS[180] is k103 in component INa (per_ms). * CONSTANTS[182] is k114 in component INa (per_ms). * CONSTANTS[184] is k125 in component INa (per_ms). * ALGEBRAIC[233] is k136 in component INa (per_ms). * CONSTANTS[177] is k18 in component INa (per_ms). * CONSTANTS[179] is k29 in component INa (per_ms). * CONSTANTS[181] is k310 in component INa (per_ms). * CONSTANTS[183] is k411 in component INa (per_ms). * CONSTANTS[185] is k512 in component INa (per_ms). * ALGEBRAIC[234] is k613 in component INa (per_ms). * STATES[50] is na1 in component INa (dimensionless). * STATES[51] is na2 in component INa (dimensionless). * STATES[52] is na3 in component INa (dimensionless). * STATES[53] is na4 in component INa (dimensionless). * STATES[54] is na5 in component INa (dimensionless). * STATES[55] is na8 in component INa (dimensionless). * STATES[56] is na9 in component INa (dimensionless). * STATES[57] is na10 in component INa (dimensionless). * STATES[58] is na11 in component INa (dimensionless). * STATES[59] is na12 in component INa (dimensionless). * STATES[60] is na13 in component INa (dimensionless). * CONSTANTS[108] is TNa in component INa (kelvin). * CONSTANTS[137] is KToverH in component INa (per_ms). * CONSTANTS[158] is FoverRT in component INa (per_mV). * CONSTANTS[157] is RTNa in component INa (joule_per_mole). * CONSTANTS[138] is RTNaF in component INa (mV). * CONSTANTS[109] is T_Const_HERG in component IKr (dimensionless). * CONSTANTS[110] is A0_HERG in component IKr (per_ms). * CONSTANTS[111] is B0_HERG in component IKr (per_mV). * CONSTANTS[112] is A1_HERG in component IKr (per_ms). * CONSTANTS[113] is B1_HERG in component IKr (per_mV). * CONSTANTS[114] is A2_HERG in component IKr (per_ms). * CONSTANTS[115] is B2_HERG in component IKr (per_mV). * CONSTANTS[116] is A3_HERG in component IKr (per_ms). * CONSTANTS[117] is B3_HERG in component IKr (per_mV). * CONSTANTS[118] is A4_HERG in component IKr (per_ms). * CONSTANTS[119] is B4_HERG in component IKr (per_mV). * CONSTANTS[120] is A5_HERG in component IKr (per_ms). * CONSTANTS[121] is B5_HERG in component IKr (per_mV). * CONSTANTS[122] is A6_HERG in component IKr (per_ms). * CONSTANTS[123] is B6_HERG in component IKr (per_mV). * CONSTANTS[139] is C2H_to_C3H in component IKr (per_ms). * CONSTANTS[140] is C3H_to_C2H in component IKr (per_ms). * ALGEBRAIC[235] is C1H_to_C2H in component IKr (per_ms). * ALGEBRAIC[236] is C2H_to_C1H in component IKr (per_ms). * ALGEBRAIC[237] is C3H_to_OH in component IKr (per_ms). * ALGEBRAIC[238] is OH_to_C3H in component IKr (per_ms). * ALGEBRAIC[239] is OH_to_IH in component IKr (per_ms). * ALGEBRAIC[240] is IH_to_OH in component IKr (per_ms). * ALGEBRAIC[241] is C3H_to_IH in component IKr (per_ms). * ALGEBRAIC[242] is IH_to_C3H in component IKr (per_ms). * STATES[61] is C1Herg in component IKr (dimensionless). * ALGEBRAIC[243] is a1_C2 in component IKr (per_ms). * ALGEBRAIC[244] is a2_C2 in component IKr (per_ms). * STATES[62] is C2Herg in component IKr (dimensionless). * ALGEBRAIC[245] is a1_C3 in component IKr (per_ms). * ALGEBRAIC[246] is a2_C3 in component IKr (per_ms). * STATES[63] is C3Herg in component IKr (dimensionless). * ALGEBRAIC[247] is a1_O in component IKr (per_ms). * ALGEBRAIC[248] is a2_O in component IKr (per_ms). * ALGEBRAIC[249] is a1_I in component IKr (per_ms). * ALGEBRAIC[250] is a2_I in component IKr (per_ms). * STATES[64] is IHerg in component IKr (dimensionless). * CONSTANTS[141] is C0ks_C1ks in component IKs (per_ms). * CONSTANTS[142] is C1ks_O1ks in component IKs (per_ms). * ALGEBRAIC[251] is O1ks_O2ks in component IKs (per_ms). * ALGEBRAIC[252] is O1ks_C1ks in component IKs (per_ms). * ALGEBRAIC[253] is O2ks_O1ks in component IKs (per_ms). * ALGEBRAIC[254] is C1ks_C0ks in component IKs (per_ms). * STATES[65] is C0ks in component IKs (dimensionless). * STATES[66] is C1ks in component IKs (dimensionless). * RATES[0] is d/dt Nai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * RATES[1] is d/dt Ki in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * RATES[2] is d/dt Cai in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * RATES[3] is d/dt CaSS in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * RATES[4] is d/dt CaJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * RATES[5] is d/dt CaNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM). * RATES[6] is d/dt V in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV). * RATES[9] is d/dt LTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * RATES[10] is d/dt HTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless). * RATES[11] is d/dt C1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * RATES[8] is d/dt O2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * RATES[12] is d/dt C2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * RATES[7] is d/dt O1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless). * RATES[13] is d/dt C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[14] is d/dt C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[15] is d/dt C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[16] is d/dt C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[17] is d/dt C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[23] is d/dt Open in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[18] is d/dt CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[19] is d/dt CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[20] is d/dt CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[21] is d/dt CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[22] is d/dt CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[24] is d/dt yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless). * RATES[25] is d/dt C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[26] is d/dt C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[27] is d/dt C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[28] is d/dt C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[29] is d/dt OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[30] is d/dt CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[31] is d/dt CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[32] is d/dt CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[33] is d/dt CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[34] is d/dt OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless). * RATES[35] is d/dt C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[36] is d/dt C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[37] is d/dt C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[38] is d/dt C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[39] is d/dt OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[40] is d/dt CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[41] is d/dt CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[42] is d/dt CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[43] is d/dt CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[44] is d/dt OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless). * RATES[50] is d/dt na1 in component INa (dimensionless). * RATES[51] is d/dt na2 in component INa (dimensionless). * RATES[52] is d/dt na3 in component INa (dimensionless). * RATES[53] is d/dt na4 in component INa (dimensionless). * RATES[54] is d/dt na5 in component INa (dimensionless). * RATES[45] is d/dt na6 in component INa (dimensionless). * RATES[46] is d/dt na7 in component INa (dimensionless). * RATES[55] is d/dt na8 in component INa (dimensionless). * RATES[56] is d/dt na9 in component INa (dimensionless). * RATES[57] is d/dt na10 in component INa (dimensionless). * RATES[58] is d/dt na11 in component INa (dimensionless). * RATES[59] is d/dt na12 in component INa (dimensionless). * RATES[60] is d/dt na13 in component INa (dimensionless). * RATES[61] is d/dt C1Herg in component IKr (dimensionless). * RATES[62] is d/dt C2Herg in component IKr (dimensionless). * RATES[63] is d/dt C3Herg in component IKr (dimensionless). * RATES[47] is d/dt OHerg in component IKr (dimensionless). * RATES[64] is d/dt IHerg in component IKr (dimensionless). * RATES[65] is d/dt C0ks in component IKs (dimensionless). * RATES[66] is d/dt C1ks in component IKs (dimensionless). * RATES[48] is d/dt O1ks in component IKs (dimensionless). * RATES[49] is d/dt O2ks in component IKs (dimensionless). * There are a total of 3 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 0; CONSTANTS[1] = 0; CONSTANTS[2] = 0; CONSTANTS[3] = 96.5; CONSTANTS[4] = 310; CONSTANTS[5] = 8.315; CONSTANTS[6] = 0.0001534; CONSTANTS[7] = 2.584e-5; CONSTANTS[8] = 1.6e-7; CONSTANTS[9] = 2.1e-6; CONSTANTS[10] = 1.2e-9; STATES[0] = 9.85573275838928; STATES[1] = 125.427082712469; STATES[2] = 0.000363968672182656; STATES[3] = 0.000506604278037024; STATES[4] = 0.421936980515042; STATES[5] = 0.423551621440241; STATES[6] = -86.7261544519706; CONSTANTS[11] = 4; CONSTANTS[12] = 138; CONSTANTS[13] = 2; CONSTANTS[14] = 500; CONSTANTS[15] = 3; CONSTANTS[16] = -15; CONSTANTS[17] = 100; CONSTANTS[18] = 0.000168; CONSTANTS[19] = 1.2; CONSTANTS[20] = 3.29; CONSTANTS[21] = 1; CONSTANTS[22] = 1.2; CONSTANTS[23] = 7.48e-5; CONSTANTS[24] = 0.000318; CONSTANTS[25] = 1.8; STATES[7] = 0.00113684728532807; STATES[8] = 3.11350788541838e-7; CONSTANTS[26] = 0.5747; CONSTANTS[27] = 26.7; CONSTANTS[28] = 0.07; CONSTANTS[29] = 0.14; CONSTANTS[30] = 20; CONSTANTS[31] = 6.6e-5; CONSTANTS[32] = 40; CONSTANTS[33] = 0.04; CONSTANTS[34] = 0.05; CONSTANTS[35] = 0; CONSTANTS[36] = 0.00238; CONSTANTS[37] = 0.8; CONSTANTS[38] = 0.00015; STATES[9] = 0.280466039150394; STATES[10] = 0.99347761599363; CONSTANTS[39] = 0.01215; CONSTANTS[40] = 0.576; CONSTANTS[41] = 0.00405; CONSTANTS[42] = 1.93; CONSTANTS[43] = 0.1; CONSTANTS[44] = 0.0008; CONSTANTS[45] = 4; CONSTANTS[46] = 3; CONSTANTS[47] = 1; CONSTANTS[48] = 0.5; CONSTANTS[49] = 0.5; CONSTANTS[50] = 5; CONSTANTS[51] = 1; CONSTANTS[52] = 1; CONSTANTS[53] = 2.5; STATES[11] = 0.132070890861418; STATES[12] = 0.866791951404883; CONSTANTS[54] = 0.3; CONSTANTS[55] = 4; CONSTANTS[56] = 2; CONSTANTS[57] = 2; STATES[13] = 0.465679150104636; STATES[14] = 0.00834457719966281; STATES[15] = 5.60736209083906e-5; STATES[16] = 1.6747092904465e-7; STATES[17] = 1.87571666668874e-10; STATES[18] = 0.489846779190386; STATES[19] = 0.035111008610982; STATES[20] = 0.000943745917866092; STATES[21] = 1.12741202215634e-5; STATES[22] = 5.05056944609524e-8; STATES[23] = 1.40806027419488e-11; STATES[24] = 0.995434385054729; STATES[25] = 0.908189132330738; STATES[26] = 0.0343385704915328; STATES[27] = 0.000487654173162347; STATES[28] = 3.11550715247964e-6; STATES[29] = 7.42911977991342e-9; STATES[30] = 0.0349937004781012; STATES[31] = 0.0171163265867255; STATES[32] = 0.00428471710061031; STATES[33] = 0.000564724236640674; STATES[34] = 2.19603439704397e-5; CONSTANTS[58] = 0.543708; CONSTANTS[59] = 0.028983; CONSTANTS[60] = 0.080185; CONSTANTS[61] = 0.0468437; CONSTANTS[62] = 0.0498424; CONSTANTS[63] = 0.000373016; CONSTANTS[64] = 0.000819482; CONSTANTS[65] = 5.374e-8; CONSTANTS[66] = 1.8936; CONSTANTS[67] = 14.224647456; CONSTANTS[68] = 158.574378389; CONSTANTS[69] = 142.936645351; CONSTANTS[70] = 6.77348; CONSTANTS[71] = 15.6212705152; CONSTANTS[72] = 28.7532603313; CONSTANTS[73] = 524.576206679; CONSTANTS[74] = 0.52465073996; CONSTANTS[75] = 17.51885408639; CONSTANTS[76] = 938.58764534556; CONSTANTS[77] = 54749.194733326; CONSTANTS[78] = 1.00947847105; CONSTANTS[79] = 1.17100540567; CONSTANTS[80] = 0.63902768758; CONSTANTS[81] = 2.12035379095; CONSTANTS[82] = 1.84002414554; CONSTANTS[83] = 0.00768548031; CONSTANTS[84] = 0.0108174834; CONSTANTS[85] = 0.07793378174; CONSTANTS[86] = 0.00305767916; CONSTANTS[87] = 2.44936e-6; STATES[35] = 0.149374350989705; STATES[36] = 0.0606794865684932; STATES[37] = 0.00930314185504921; STATES[38] = 0.000676403999474111; STATES[39] = 3.85187206387239e-5; STATES[40] = 0.0442722560882536; STATES[41] = 0.00952432663172288; STATES[42] = 0.0567396669678271; STATES[43] = 0.113122845136053; STATES[44] = 0.556269044084734; CONSTANTS[88] = 0.0186; CONSTANTS[89] = 0.0035; CONSTANTS[90] = 0.125305126118808; CONSTANTS[91] = 56.32; CONSTANTS[92] = 0.001; CONSTANTS[93] = 0.872; CONSTANTS[94] = 0.889; STATES[45] = 1.02118700961583e-7; STATES[46] = 1.93499158844817e-8; STATES[47] = 0.00120284688677794; STATES[48] = 5.65460174551007e-7; STATES[49] = 0.0258818770122187; CONSTANTS[95] = 0.44; CONSTANTS[96] = 87.5; CONSTANTS[97] = 1.38; CONSTANTS[98] = 0.2; CONSTANTS[99] = 0.35; CONSTANTS[100] = 2.387; CONSTANTS[101] = 20; CONSTANTS[102] = 1.5; CONSTANTS[103] = 0.05; CONSTANTS[104] = 0.0005; CONSTANTS[105] = 7.684e-5; CONSTANTS[106] = -0.265; CONSTANTS[107] = 7; STATES[50] = 0.111284526171411; STATES[51] = 0.0481019786429977; STATES[52] = 0.00779692701457915; STATES[53] = 0.000561699600929369; STATES[54] = 1.51746424723121e-5; STATES[55] = 0.368582741846592; STATES[56] = 0.312463212648791; STATES[57] = 0.0993398770493615; STATES[58] = 0.0140431688972267; STATES[59] = 0.000750073829883749; STATES[60] = 0.0370604970714329; CONSTANTS[108] = 294.16; CONSTANTS[109] = 5.320000001; CONSTANTS[110] = 0.017147641733086; CONSTANTS[111] = 0.03304608038835; CONSTANTS[112] = 0.03969328381141; CONSTANTS[113] = -0.0430605416398; CONSTANTS[114] = 0.02057448605977; CONSTANTS[115] = 0.02617412715118; CONSTANTS[116] = 0.00134366604423; CONSTANTS[117] = -0.02691385498399; CONSTANTS[118] = 0.10666316491288; CONSTANTS[119] = 0.00568908859717; CONSTANTS[120] = 0.00646393910049; CONSTANTS[121] = -0.04536642959543; CONSTANTS[122] = 8.039374403e-5; CONSTANTS[123] = 6.9808924e-7; STATES[61] = 0.994948338598163; STATES[62] = 0.000595653663190548; STATES[63] = 0.000228183228829573; STATES[64] = 0.000243789721526602; STATES[65] = 0.938064990549233; STATES[66] = 0.0360525668093578; CONSTANTS[124] = ( CONSTANTS[5]*CONSTANTS[4])/CONSTANTS[3]; CONSTANTS[125] = (CONSTANTS[2]==1.00000 ? 15.0000 : 0.500000*15.0000); CONSTANTS[126] = (CONSTANTS[1]==1.00000 ? 0.200000 : 0.500000); CONSTANTS[127] = (CONSTANTS[1]==1.00000 ? 0.500000 : 1.00000); CONSTANTS[128] = 0.250000*0.0100000; CONSTANTS[129] = 0.820000; CONSTANTS[130] = CONSTANTS[86]; CONSTANTS[131] = CONSTANTS[87]; CONSTANTS[132] = CONSTANTS[94]*CONSTANTS[93]*0.100000; CONSTANTS[133] = (1.00000 - CONSTANTS[94])*CONSTANTS[93]*4.29860e-07; CONSTANTS[134] = pow((CONSTANTS[11]/4.00000), 1.0 / 2); CONSTANTS[135] = CONSTANTS[6]*0.00100000; CONSTANTS[136] = CONSTANTS[107]*4.57400e-07; CONSTANTS[137] = ( 1.38100e-23*CONSTANTS[108])/6.62600e-31; CONSTANTS[138] = ( CONSTANTS[5]*CONSTANTS[108])/CONSTANTS[3]; CONSTANTS[139] = CONSTANTS[109]*0.0260836; CONSTANTS[140] = CONSTANTS[109]*0.148330; CONSTANTS[141] = 0.00795601; CONSTANTS[142] = 0.0396672; CONSTANTS[143] = CONSTANTS[128]; CONSTANTS[144] = CONSTANTS[131]; CONSTANTS[145] = CONSTANTS[74]*CONSTANTS[131]; CONSTANTS[146] = CONSTANTS[75]*CONSTANTS[131]; CONSTANTS[147] = CONSTANTS[76]*CONSTANTS[131]; CONSTANTS[148] = CONSTANTS[77]*CONSTANTS[131]; CONSTANTS[149] = CONSTANTS[130]; CONSTANTS[150] = CONSTANTS[130]/CONSTANTS[78]; CONSTANTS[151] = CONSTANTS[130]/CONSTANTS[79]; CONSTANTS[152] = CONSTANTS[130]/CONSTANTS[80]; CONSTANTS[153] = CONSTANTS[130]/CONSTANTS[81]; CONSTANTS[154] = CONSTANTS[11]+ 0.0183300*CONSTANTS[12]; CONSTANTS[155] = CONSTANTS[135]/( CONSTANTS[7]*CONSTANTS[3]); CONSTANTS[156] = (CONSTANTS[0]==0.00000 ? CONSTANTS[107]*0.000246900 : 1.50000*CONSTANTS[107]*0.000246900); CONSTANTS[157] = CONSTANTS[5]*CONSTANTS[108]; CONSTANTS[158] = 1.00000/CONSTANTS[124]; CONSTANTS[159] = CONSTANTS[143]/CONSTANTS[56]; CONSTANTS[160] = CONSTANTS[135]/( 2.00000*CONSTANTS[10]*CONSTANTS[3]); CONSTANTS[161] = 1.38862; CONSTANTS[162] = CONSTANTS[159]/CONSTANTS[56]; CONSTANTS[163] = (exp(CONSTANTS[12]/67.3000) - 1.00000)/7.00000; CONSTANTS[164] = CONSTANTS[161]*CONSTANTS[137]*exp(- 287913./CONSTANTS[157]+786.217/CONSTANTS[5]); CONSTANTS[165] = CONSTANTS[161]*CONSTANTS[137]*exp(- 85800.4/CONSTANTS[157]+70.0780/CONSTANTS[5]); CONSTANTS[166] = CONSTANTS[161]*CONSTANTS[137]*exp(- 121955./CONSTANTS[157]+225.175/CONSTANTS[5]); CONSTANTS[167] = CONSTANTS[11]/(CONSTANTS[11]+CONSTANTS[102]); CONSTANTS[168] = CONSTANTS[162]/CONSTANTS[56]; CONSTANTS[169] = CONSTANTS[161]*CONSTANTS[137]*exp(- 59565.2/CONSTANTS[157]+0.00711000/CONSTANTS[5]); CONSTANTS[170] = CONSTANTS[165]; CONSTANTS[171] = CONSTANTS[166]; CONSTANTS[172] = CONSTANTS[97]+CONSTANTS[13]; CONSTANTS[173] = CONSTANTS[168]/CONSTANTS[56]; CONSTANTS[174] = 1.40043; CONSTANTS[175] = (pow(CONSTANTS[96], 3.00000)+pow(CONSTANTS[12], 3.00000))/5000.00; CONSTANTS[176] = CONSTANTS[169]; CONSTANTS[177] = CONSTANTS[164]; CONSTANTS[178] = CONSTANTS[176]/CONSTANTS[174]; CONSTANTS[179] = CONSTANTS[177]*CONSTANTS[174]; CONSTANTS[180] = CONSTANTS[178]/CONSTANTS[174]; CONSTANTS[181] = CONSTANTS[179]*CONSTANTS[174]; CONSTANTS[182] = CONSTANTS[180]/CONSTANTS[174]; CONSTANTS[183] = CONSTANTS[181]*CONSTANTS[174]; CONSTANTS[184] = CONSTANTS[182]/CONSTANTS[174]; CONSTANTS[185] = CONSTANTS[183]*CONSTANTS[174]; RATES[0] = 0.1001; RATES[1] = 0.1001; RATES[2] = 0.1001; RATES[3] = 0.1001; RATES[4] = 0.1001; RATES[5] = 0.1001; RATES[6] = 0.1001; RATES[9] = 0.1001; RATES[10] = 0.1001; RATES[11] = 0.1001; RATES[8] = 0.1001; RATES[12] = 0.1001; RATES[7] = 0.1001; RATES[13] = 0.1001; RATES[14] = 0.1001; RATES[15] = 0.1001; RATES[16] = 0.1001; RATES[17] = 0.1001; RATES[23] = 0.1001; RATES[18] = 0.1001; RATES[19] = 0.1001; RATES[20] = 0.1001; RATES[21] = 0.1001; RATES[22] = 0.1001; RATES[24] = 0.1001; RATES[25] = 0.1001; RATES[26] = 0.1001; RATES[27] = 0.1001; RATES[28] = 0.1001; RATES[29] = 0.1001; RATES[30] = 0.1001; RATES[31] = 0.1001; RATES[32] = 0.1001; RATES[33] = 0.1001; RATES[34] = 0.1001; RATES[35] = 0.1001; RATES[36] = 0.1001; RATES[37] = 0.1001; RATES[38] = 0.1001; RATES[39] = 0.1001; RATES[40] = 0.1001; RATES[41] = 0.1001; RATES[42] = 0.1001; RATES[43] = 0.1001; RATES[44] = 0.1001; RATES[50] = 0.1001; RATES[51] = 0.1001; RATES[52] = 0.1001; RATES[53] = 0.1001; RATES[54] = 0.1001; RATES[45] = 0.1001; RATES[46] = 0.1001; RATES[55] = 0.1001; RATES[56] = 0.1001; RATES[57] = 0.1001; RATES[58] = 0.1001; RATES[59] = 0.1001; RATES[60] = 0.1001; RATES[61] = 0.1001; RATES[62] = 0.1001; RATES[63] = 0.1001; RATES[47] = 0.1001; RATES[64] = 0.1001; RATES[65] = 0.1001; RATES[66] = 0.1001; RATES[48] = 0.1001; RATES[49] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[0] - - (ALGEBRAIC[162]+ALGEBRAIC[177]+ 3.00000*(ALGEBRAIC[189]+ALGEBRAIC[185])+ALGEBRAIC[172])*CONSTANTS[155]; resid[1] = RATES[1] - - (((ALGEBRAIC[163]+ALGEBRAIC[164]+ALGEBRAIC[176]+ALGEBRAIC[200]+ALGEBRAIC[1]) - 2.00000*ALGEBRAIC[185])+ALGEBRAIC[165]+ALGEBRAIC[170])*CONSTANTS[155]; resid[2] = RATES[2] - ALGEBRAIC[13]*(((ALGEBRAIC[7] - ALGEBRAIC[4]) - ALGEBRAIC[10]) - ((ALGEBRAIC[178] - 2.00000*ALGEBRAIC[189])+ALGEBRAIC[179])*0.500000*CONSTANTS[155]); resid[3] = RATES[3] - ALGEBRAIC[11]*((( ALGEBRAIC[5]*CONSTANTS[8])/CONSTANTS[10] - ( ALGEBRAIC[7]*CONSTANTS[7])/CONSTANTS[10]) - ALGEBRAIC[195]*CONSTANTS[160]); resid[4] = RATES[4] - ALGEBRAIC[12]*(ALGEBRAIC[6] - ALGEBRAIC[5]); resid[5] = RATES[5] - ( ALGEBRAIC[4]*CONSTANTS[7])/CONSTANTS[9] - ( ALGEBRAIC[6]*CONSTANTS[8])/CONSTANTS[9]; resid[6] = RATES[6] - - ALGEBRAIC[255]; resid[7] = RATES[9] - ALGEBRAIC[8]; resid[8] = RATES[10] - ALGEBRAIC[9]; resid[9] = RATES[11] - ALGEBRAIC[17]; resid[10] = RATES[8] - ALGEBRAIC[18]; resid[11] = RATES[12] - ALGEBRAIC[19]; resid[12] = RATES[7] - ALGEBRAIC[20]; resid[13] = RATES[13] - ALGEBRAIC[48] - ALGEBRAIC[47]; resid[14] = RATES[14] - ALGEBRAIC[50] - ALGEBRAIC[49]; resid[15] = RATES[15] - ALGEBRAIC[52] - ALGEBRAIC[51]; resid[16] = RATES[16] - ALGEBRAIC[54] - ALGEBRAIC[53]; resid[17] = RATES[17] - ALGEBRAIC[56] - ALGEBRAIC[55]; resid[18] = RATES[23] - CONSTANTS[54]*STATES[17] - CONSTANTS[55]*STATES[23]; resid[19] = RATES[18] - ALGEBRAIC[58] - ALGEBRAIC[57]; resid[20] = RATES[19] - ALGEBRAIC[60] - ALGEBRAIC[59]; resid[21] = RATES[20] - ALGEBRAIC[62] - ALGEBRAIC[61]; resid[22] = RATES[21] - ALGEBRAIC[64] - ALGEBRAIC[63]; resid[23] = RATES[22] - ALGEBRAIC[66] - ALGEBRAIC[65]; resid[24] = RATES[24] - (ALGEBRAIC[67] - STATES[24])/ALGEBRAIC[68]; resid[25] = RATES[25] - ALGEBRAIC[100] - ALGEBRAIC[99]; resid[26] = RATES[26] - ALGEBRAIC[102] - ALGEBRAIC[101]; resid[27] = RATES[27] - ALGEBRAIC[104] - ALGEBRAIC[103]; resid[28] = RATES[28] - ALGEBRAIC[106] - ALGEBRAIC[105]; resid[29] = RATES[29] - ALGEBRAIC[108] - ALGEBRAIC[107]; resid[30] = RATES[30] - ALGEBRAIC[110] - ALGEBRAIC[109]; resid[31] = RATES[31] - ALGEBRAIC[112] - ALGEBRAIC[111]; resid[32] = RATES[32] - ALGEBRAIC[114] - ALGEBRAIC[113]; resid[33] = RATES[33] - ALGEBRAIC[116] - ALGEBRAIC[115]; resid[34] = RATES[34] - ALGEBRAIC[118] - ALGEBRAIC[117]; resid[35] = RATES[35] - ALGEBRAIC[138] - ALGEBRAIC[137]; resid[36] = RATES[36] - ALGEBRAIC[140] - ALGEBRAIC[139]; resid[37] = RATES[37] - ALGEBRAIC[142] - ALGEBRAIC[141]; resid[38] = RATES[38] - ALGEBRAIC[144] - ALGEBRAIC[143]; resid[39] = RATES[39] - ALGEBRAIC[146] - ALGEBRAIC[145]; resid[40] = RATES[40] - ALGEBRAIC[148] - ALGEBRAIC[147]; resid[41] = RATES[41] - ALGEBRAIC[150] - ALGEBRAIC[149]; resid[42] = RATES[42] - ALGEBRAIC[152] - ALGEBRAIC[151]; resid[43] = RATES[43] - ALGEBRAIC[154] - ALGEBRAIC[153]; resid[44] = RATES[44] - ALGEBRAIC[156] - ALGEBRAIC[155]; resid[45] = RATES[50] - - (CONSTANTS[177]+ALGEBRAIC[211])*STATES[50]+ ALGEBRAIC[222]*STATES[51]+ CONSTANTS[176]*STATES[55]; resid[46] = RATES[51] - ( ALGEBRAIC[211]*STATES[50] - (ALGEBRAIC[222]+ALGEBRAIC[212]+CONSTANTS[179])*STATES[51])+ ALGEBRAIC[223]*STATES[52]+ CONSTANTS[178]*STATES[56]; resid[47] = RATES[52] - ( ALGEBRAIC[212]*STATES[51] - (ALGEBRAIC[223]+ALGEBRAIC[213]+CONSTANTS[181])*STATES[52])+ ALGEBRAIC[224]*STATES[53]+ CONSTANTS[180]*STATES[57]; resid[48] = RATES[53] - ( ALGEBRAIC[213]*STATES[52] - (ALGEBRAIC[224]+ALGEBRAIC[214]+CONSTANTS[183])*STATES[53])+ ALGEBRAIC[225]*STATES[54]+ CONSTANTS[182]*STATES[58]; resid[49] = RATES[54] - ( ALGEBRAIC[214]*STATES[53] - (ALGEBRAIC[225]+ALGEBRAIC[215]+ALGEBRAIC[221]+CONSTANTS[185])*STATES[54])+ ALGEBRAIC[226]*STATES[45]+ ALGEBRAIC[232]*STATES[46]+ CONSTANTS[184]*STATES[59]; resid[50] = RATES[45] - ( ALGEBRAIC[215]*STATES[54] - (ALGEBRAIC[226]+CONSTANTS[170]+ALGEBRAIC[234])*STATES[45])+ CONSTANTS[171]*STATES[46]+ ALGEBRAIC[233]*STATES[60]; resid[51] = RATES[46] - ( ALGEBRAIC[221]*STATES[54]+ CONSTANTS[170]*STATES[45]) - (ALGEBRAIC[232]+CONSTANTS[171])*STATES[46]; resid[52] = RATES[55] - ( CONSTANTS[177]*STATES[50] - (CONSTANTS[176]+ALGEBRAIC[216])*STATES[55])+ ALGEBRAIC[227]*STATES[56]; resid[53] = RATES[56] - (( CONSTANTS[179]*STATES[51]+ ALGEBRAIC[216]*STATES[55]) - (ALGEBRAIC[227]+CONSTANTS[178]+ALGEBRAIC[217])*STATES[56])+ ALGEBRAIC[228]*STATES[57]; resid[54] = RATES[57] - (( CONSTANTS[181]*STATES[52]+ ALGEBRAIC[217]*STATES[56]) - (ALGEBRAIC[218]+CONSTANTS[180]+ALGEBRAIC[228])*STATES[57])+ ALGEBRAIC[229]*STATES[58]; resid[55] = RATES[58] - (( CONSTANTS[183]*STATES[53]+ ALGEBRAIC[218]*STATES[57]) - (ALGEBRAIC[229]+CONSTANTS[182]+ALGEBRAIC[219])*STATES[58])+ ALGEBRAIC[230]*STATES[59]; resid[56] = RATES[59] - (( CONSTANTS[185]*STATES[54]+ ALGEBRAIC[219]*STATES[58]) - (ALGEBRAIC[230]+CONSTANTS[184]+ALGEBRAIC[220])*STATES[59])+ ALGEBRAIC[231]*STATES[60]; resid[57] = RATES[60] - ( ALGEBRAIC[234]*STATES[45]+ ALGEBRAIC[220]*STATES[59]) - (ALGEBRAIC[231]+ALGEBRAIC[233])*STATES[60]; resid[58] = RATES[61] - ALGEBRAIC[236]*STATES[62] - ALGEBRAIC[235]*STATES[61]; resid[59] = RATES[62] - ALGEBRAIC[243] - ALGEBRAIC[244]; resid[60] = RATES[63] - ALGEBRAIC[245] - ALGEBRAIC[246]; resid[61] = RATES[47] - ALGEBRAIC[247] - ALGEBRAIC[248]; resid[62] = RATES[64] - ALGEBRAIC[249] - ALGEBRAIC[250]; resid[63] = RATES[65] - - CONSTANTS[141]*STATES[65]+ ALGEBRAIC[254]*STATES[66]; resid[64] = RATES[66] - ( CONSTANTS[141]*STATES[65] - (ALGEBRAIC[254]+CONSTANTS[142])*STATES[66])+ ALGEBRAIC[252]*STATES[48]; resid[65] = RATES[48] - ( CONSTANTS[142]*STATES[66] - (ALGEBRAIC[252]+ALGEBRAIC[251])*STATES[48])+ ALGEBRAIC[253]*STATES[49]; resid[66] = RATES[49] - ALGEBRAIC[251]*STATES[48] - ALGEBRAIC[253]*STATES[49]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[159] = STATES[1]+ 0.0183300*STATES[0]; ALGEBRAIC[160] = CONSTANTS[124]*log(CONSTANTS[154]/ALGEBRAIC[159]); } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = floor(VOI/CONSTANTS[14])*CONSTANTS[14]; ALGEBRAIC[1] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<=0.00000 ? CONSTANTS[16] : 0.00000); ALGEBRAIC[2] = pow(STATES[2]/CONSTANTS[18], CONSTANTS[19]); ALGEBRAIC[3] = pow(STATES[5]/CONSTANTS[20], CONSTANTS[21]); ALGEBRAIC[4] = (CONSTANTS[0]==0.00000 ? ( CONSTANTS[22]*( CONSTANTS[23]*ALGEBRAIC[2] - CONSTANTS[24]*ALGEBRAIC[3]))/(1.00000+ALGEBRAIC[2]+ALGEBRAIC[3]) : ( 1.50000*CONSTANTS[22]*( CONSTANTS[23]*ALGEBRAIC[2] - CONSTANTS[24]*ALGEBRAIC[3]))/(1.00000+ALGEBRAIC[2]+ALGEBRAIC[3])); ALGEBRAIC[5] = CONSTANTS[25]*(STATES[7]+STATES[8])*(STATES[4] - STATES[3]); ALGEBRAIC[6] = (STATES[5] - STATES[4])/CONSTANTS[26]; ALGEBRAIC[7] = (STATES[3] - STATES[2])/CONSTANTS[27]; ALGEBRAIC[8] = CONSTANTS[32]*STATES[2]*(1.00000 - STATES[9]) - CONSTANTS[33]*STATES[9]; ALGEBRAIC[9] = CONSTANTS[30]*STATES[2]*(1.00000 - STATES[10]) - CONSTANTS[31]*STATES[10]; ALGEBRAIC[10] = CONSTANTS[28]*ALGEBRAIC[8]+ CONSTANTS[29]*ALGEBRAIC[9]; ALGEBRAIC[11] = 1.00000/(1.00000+( CONSTANTS[34]*CONSTANTS[36])/pow(STATES[3]+CONSTANTS[36], 2.00000)+( CONSTANTS[35]*CONSTANTS[38])/pow(STATES[3]+CONSTANTS[38], 2.00000)); ALGEBRAIC[12] = 1.00000/(1.00000+( CONSTANTS[125]*CONSTANTS[37])/pow(STATES[4]+CONSTANTS[37], 2.00000)); ALGEBRAIC[13] = 1.00000/(1.00000+( CONSTANTS[34]*CONSTANTS[36])/pow(STATES[2]+CONSTANTS[36], 2.00000)+( CONSTANTS[35]*CONSTANTS[38])/pow(STATES[2]+CONSTANTS[38], 2.00000)); ALGEBRAIC[14] = CONSTANTS[48] - (CONSTANTS[48] - CONSTANTS[50])/(1.00000+pow(CONSTANTS[51]/STATES[4], CONSTANTS[53])); ALGEBRAIC[17] = - CONSTANTS[39]*ALGEBRAIC[14]*pow( STATES[3]*1000.00, CONSTANTS[45]*CONSTANTS[47])*STATES[11]+ CONSTANTS[40]*STATES[7]; ALGEBRAIC[18] = CONSTANTS[41]*pow( STATES[3]*1000.00, CONSTANTS[46]*CONSTANTS[47])*STATES[7] - CONSTANTS[42]*STATES[8]; ALGEBRAIC[15] = CONSTANTS[126] - (CONSTANTS[126] - CONSTANTS[50])/(1.00000+pow(CONSTANTS[52]/STATES[4], CONSTANTS[53])); ALGEBRAIC[16] = CONSTANTS[49] - (CONSTANTS[49] - CONSTANTS[50])/(1.00000+pow(CONSTANTS[127]/STATES[4], CONSTANTS[53])); ALGEBRAIC[19] = (CONSTANTS[43]/ALGEBRAIC[15])*STATES[7] - CONSTANTS[44]*ALGEBRAIC[16]*STATES[12]; ALGEBRAIC[20] = - (ALGEBRAIC[17]+ALGEBRAIC[18]+ALGEBRAIC[19]); ALGEBRAIC[21] = 4.00000*1.20000*0.416000*exp( 0.0120000*(STATES[6] - 35.0000)); ALGEBRAIC[26] = 4.00000*ALGEBRAIC[21]; ALGEBRAIC[25] = 0.600000*0.0923300*STATES[3]; ALGEBRAIC[42] = ALGEBRAIC[25]; ALGEBRAIC[47] = (ALGEBRAIC[26]+ALGEBRAIC[42])*STATES[13]; ALGEBRAIC[22] = 4.00000*0.450000*0.0490000*exp( - 0.0650000*(STATES[6] - 22.0000)); ALGEBRAIC[34] = ALGEBRAIC[22]; ALGEBRAIC[48] = ALGEBRAIC[34]*STATES[14]+ CONSTANTS[143]*STATES[18]; ALGEBRAIC[27] = 3.00000*ALGEBRAIC[21]; ALGEBRAIC[43] = CONSTANTS[57]*ALGEBRAIC[42]; ALGEBRAIC[49] = (ALGEBRAIC[34]+ALGEBRAIC[27]+ALGEBRAIC[43])*STATES[14]; ALGEBRAIC[35] = 2.00000*ALGEBRAIC[22]; ALGEBRAIC[50] = ALGEBRAIC[26]*STATES[13]+ ALGEBRAIC[35]*STATES[15]+ CONSTANTS[159]*STATES[19]; ALGEBRAIC[28] = 2.00000*ALGEBRAIC[21]; ALGEBRAIC[44] = CONSTANTS[57]*ALGEBRAIC[43]; ALGEBRAIC[51] = (ALGEBRAIC[35]+ALGEBRAIC[28]+ALGEBRAIC[44])*STATES[15]; ALGEBRAIC[36] = 3.00000*ALGEBRAIC[22]; ALGEBRAIC[52] = ALGEBRAIC[27]*STATES[14]+ ALGEBRAIC[36]*STATES[16]+ CONSTANTS[162]*STATES[20]; ALGEBRAIC[29] = ALGEBRAIC[21]; ALGEBRAIC[45] = CONSTANTS[57]*ALGEBRAIC[44]; ALGEBRAIC[53] = (ALGEBRAIC[36]+ALGEBRAIC[29]+ALGEBRAIC[45])*STATES[16]; ALGEBRAIC[37] = 4.00000*ALGEBRAIC[22]; ALGEBRAIC[54] = ALGEBRAIC[28]*STATES[15]+ ALGEBRAIC[37]*STATES[17]+ CONSTANTS[168]*STATES[21]; ALGEBRAIC[46] = CONSTANTS[57]*ALGEBRAIC[45]; ALGEBRAIC[55] = (ALGEBRAIC[37]+CONSTANTS[54]+ALGEBRAIC[46])*STATES[17]; ALGEBRAIC[56] = ALGEBRAIC[29]*STATES[16]+ CONSTANTS[55]*STATES[23]+ CONSTANTS[173]*STATES[22]; ALGEBRAIC[23] = CONSTANTS[57]*ALGEBRAIC[21]; ALGEBRAIC[30] = 4.00000*ALGEBRAIC[23]; ALGEBRAIC[57] = (ALGEBRAIC[30]+CONSTANTS[143])*STATES[18]; ALGEBRAIC[24] = ALGEBRAIC[22]/CONSTANTS[56]; ALGEBRAIC[38] = ALGEBRAIC[24]; ALGEBRAIC[58] = ALGEBRAIC[38]*STATES[19]+ ALGEBRAIC[42]*STATES[13]; ALGEBRAIC[31] = 3.00000*ALGEBRAIC[23]; ALGEBRAIC[59] = (ALGEBRAIC[38]+ALGEBRAIC[31]+CONSTANTS[159])*STATES[19]; ALGEBRAIC[39] = 2.00000*ALGEBRAIC[24]; ALGEBRAIC[60] = ALGEBRAIC[30]*STATES[18]+ ALGEBRAIC[39]*STATES[20]+ ALGEBRAIC[43]*STATES[14]; ALGEBRAIC[32] = 2.00000*ALGEBRAIC[23]; ALGEBRAIC[61] = (ALGEBRAIC[39]+ALGEBRAIC[32]+CONSTANTS[162])*STATES[20]; ALGEBRAIC[40] = 3.00000*ALGEBRAIC[24]; ALGEBRAIC[62] = ALGEBRAIC[31]*STATES[19]+ ALGEBRAIC[40]*STATES[21]+ ALGEBRAIC[44]*STATES[15]; ALGEBRAIC[33] = ALGEBRAIC[23]; ALGEBRAIC[63] = (ALGEBRAIC[40]+ALGEBRAIC[33]+CONSTANTS[168])*STATES[21]; ALGEBRAIC[41] = 4.00000*ALGEBRAIC[24]; ALGEBRAIC[64] = ALGEBRAIC[32]*STATES[20]+ ALGEBRAIC[41]*STATES[22]+ ALGEBRAIC[45]*STATES[16]; ALGEBRAIC[65] = (ALGEBRAIC[41]+CONSTANTS[173])*STATES[22]; ALGEBRAIC[66] = ALGEBRAIC[33]*STATES[21]+ ALGEBRAIC[46]*STATES[17]; ALGEBRAIC[67] = (CONSTANTS[129]/(1.00000+exp((STATES[6]+28.5000)/7.80000))+1.00000) - CONSTANTS[129]; ALGEBRAIC[68] = 1.00000/(0.00336336/(0.500000+exp(STATES[6]/- 5.53900))+ 0.00779047*exp(STATES[6]/- 49.5104)); ALGEBRAIC[69] = CONSTANTS[58]*exp( CONSTANTS[59]*STATES[6]); ALGEBRAIC[73] = 4.00000*ALGEBRAIC[69]; ALGEBRAIC[72] = CONSTANTS[64]*exp( CONSTANTS[65]*STATES[6]); ALGEBRAIC[89] = ALGEBRAIC[72]; ALGEBRAIC[99] = (ALGEBRAIC[73]+ALGEBRAIC[89])*STATES[25]; ALGEBRAIC[70] = CONSTANTS[60]*exp( - CONSTANTS[61]*STATES[6]); ALGEBRAIC[81] = ALGEBRAIC[70]; ALGEBRAIC[71] = CONSTANTS[62]*exp( - CONSTANTS[63]*STATES[6]); ALGEBRAIC[94] = ALGEBRAIC[71]; ALGEBRAIC[100] = ALGEBRAIC[81]*STATES[26]+ ALGEBRAIC[94]*STATES[30]; ALGEBRAIC[74] = 3.00000*ALGEBRAIC[69]; ALGEBRAIC[90] = CONSTANTS[66]*ALGEBRAIC[72]; ALGEBRAIC[101] = (ALGEBRAIC[74]+ALGEBRAIC[81]+ALGEBRAIC[90])*STATES[26]; ALGEBRAIC[82] = 2.00000*ALGEBRAIC[70]; ALGEBRAIC[95] = ALGEBRAIC[71]/CONSTANTS[70]; ALGEBRAIC[102] = ALGEBRAIC[82]*STATES[27]+ ALGEBRAIC[95]*STATES[31]+ ALGEBRAIC[73]*STATES[25]; ALGEBRAIC[75] = 2.00000*ALGEBRAIC[69]; ALGEBRAIC[91] = CONSTANTS[67]*ALGEBRAIC[72]; ALGEBRAIC[103] = (ALGEBRAIC[75]+ALGEBRAIC[82]+ALGEBRAIC[91])*STATES[27]; ALGEBRAIC[83] = 3.00000*ALGEBRAIC[70]; ALGEBRAIC[96] = ALGEBRAIC[71]/CONSTANTS[71]; ALGEBRAIC[104] = ALGEBRAIC[83]*STATES[28]+ ALGEBRAIC[96]*STATES[32]+ ALGEBRAIC[74]*STATES[26]; ALGEBRAIC[76] = ALGEBRAIC[69]; ALGEBRAIC[92] = CONSTANTS[68]*ALGEBRAIC[72]; ALGEBRAIC[105] = (ALGEBRAIC[76]+ALGEBRAIC[83]+ALGEBRAIC[92])*STATES[28]; ALGEBRAIC[84] = 4.00000*ALGEBRAIC[70]; ALGEBRAIC[97] = ALGEBRAIC[71]/CONSTANTS[72]; ALGEBRAIC[106] = ALGEBRAIC[84]*STATES[29]+ ALGEBRAIC[97]*STATES[33]+ ALGEBRAIC[75]*STATES[27]; ALGEBRAIC[93] = CONSTANTS[69]*ALGEBRAIC[72]; ALGEBRAIC[107] = (ALGEBRAIC[84]+ALGEBRAIC[93])*STATES[29]; ALGEBRAIC[98] = ALGEBRAIC[71]/CONSTANTS[73]; ALGEBRAIC[108] = ALGEBRAIC[76]*STATES[28]+ ALGEBRAIC[98]*STATES[34]; ALGEBRAIC[77] = 4.00000*CONSTANTS[70]*ALGEBRAIC[69]; ALGEBRAIC[109] = (ALGEBRAIC[94]+ALGEBRAIC[77])*STATES[30]; ALGEBRAIC[85] = ALGEBRAIC[70]/CONSTANTS[66]; ALGEBRAIC[110] = ALGEBRAIC[89]*STATES[25]+ ALGEBRAIC[85]*STATES[31]; ALGEBRAIC[78] = ( 3.00000*CONSTANTS[71]*ALGEBRAIC[69])/CONSTANTS[70]; ALGEBRAIC[111] = (ALGEBRAIC[78]+ALGEBRAIC[95]+ALGEBRAIC[85])*STATES[31]; ALGEBRAIC[86] = ( 2.00000*CONSTANTS[66]*ALGEBRAIC[70])/CONSTANTS[67]; ALGEBRAIC[112] = ALGEBRAIC[86]*STATES[32]+ ALGEBRAIC[90]*STATES[26]+ ALGEBRAIC[77]*STATES[30]; ALGEBRAIC[79] = ( 2.00000*CONSTANTS[72]*ALGEBRAIC[69])/CONSTANTS[71]; ALGEBRAIC[113] = (ALGEBRAIC[79]+ALGEBRAIC[96]+ALGEBRAIC[86])*STATES[32]; ALGEBRAIC[87] = ( 3.00000*CONSTANTS[67]*ALGEBRAIC[70])/CONSTANTS[68]; ALGEBRAIC[114] = ALGEBRAIC[87]*STATES[33]+ ALGEBRAIC[91]*STATES[27]+ ALGEBRAIC[78]*STATES[31]; ALGEBRAIC[80] = ( CONSTANTS[73]*ALGEBRAIC[69])/CONSTANTS[72]; ALGEBRAIC[115] = (ALGEBRAIC[80]+ALGEBRAIC[97]+ALGEBRAIC[87])*STATES[33]; ALGEBRAIC[88] = ( 4.00000*CONSTANTS[68]*ALGEBRAIC[70])/CONSTANTS[69]; ALGEBRAIC[116] = ALGEBRAIC[88]*STATES[34]+ ALGEBRAIC[92]*STATES[28]+ ALGEBRAIC[79]*STATES[32]; ALGEBRAIC[117] = (ALGEBRAIC[98]+ALGEBRAIC[88])*STATES[34]; ALGEBRAIC[118] = ALGEBRAIC[93]*STATES[29]+ ALGEBRAIC[80]*STATES[33]; ALGEBRAIC[119] = CONSTANTS[82]*exp( CONSTANTS[83]*STATES[6]); ALGEBRAIC[121] = 4.00000*ALGEBRAIC[119]; ALGEBRAIC[137] = (ALGEBRAIC[121]+CONSTANTS[144])*STATES[35]; ALGEBRAIC[120] = CONSTANTS[84]*exp( - CONSTANTS[85]*STATES[6]); ALGEBRAIC[129] = ALGEBRAIC[120]; ALGEBRAIC[138] = ALGEBRAIC[129]*STATES[36]+ CONSTANTS[149]*STATES[40]; ALGEBRAIC[122] = 3.00000*ALGEBRAIC[119]; ALGEBRAIC[139] = (ALGEBRAIC[122]+ALGEBRAIC[129]+CONSTANTS[145])*STATES[36]; ALGEBRAIC[130] = 2.00000*ALGEBRAIC[120]; ALGEBRAIC[140] = ALGEBRAIC[130]*STATES[37]+ CONSTANTS[150]*STATES[41]+ ALGEBRAIC[121]*STATES[35]; ALGEBRAIC[123] = 2.00000*ALGEBRAIC[119]; ALGEBRAIC[141] = (ALGEBRAIC[123]+ALGEBRAIC[130]+CONSTANTS[146])*STATES[37]; ALGEBRAIC[131] = 3.00000*ALGEBRAIC[120]; ALGEBRAIC[142] = ALGEBRAIC[131]*STATES[38]+ CONSTANTS[151]*STATES[42]+ ALGEBRAIC[122]*STATES[36]; ALGEBRAIC[124] = ALGEBRAIC[119]; ALGEBRAIC[143] = (ALGEBRAIC[124]+ALGEBRAIC[131]+CONSTANTS[147])*STATES[38]; ALGEBRAIC[132] = 4.00000*ALGEBRAIC[120]; ALGEBRAIC[144] = ALGEBRAIC[132]*STATES[39]+ CONSTANTS[152]*STATES[43]+ ALGEBRAIC[123]*STATES[37]; ALGEBRAIC[145] = (ALGEBRAIC[132]+CONSTANTS[148])*STATES[39]; ALGEBRAIC[146] = ALGEBRAIC[124]*STATES[38]+ CONSTANTS[153]*STATES[44]; ALGEBRAIC[125] = 4.00000*CONSTANTS[78]*ALGEBRAIC[119]; ALGEBRAIC[147] = (CONSTANTS[149]+ALGEBRAIC[125])*STATES[40]; ALGEBRAIC[133] = ALGEBRAIC[120]/CONSTANTS[74]; ALGEBRAIC[148] = CONSTANTS[144]*STATES[35]+ ALGEBRAIC[133]*STATES[41]; ALGEBRAIC[126] = ( 3.00000*CONSTANTS[79]*ALGEBRAIC[119])/CONSTANTS[78]; ALGEBRAIC[149] = (ALGEBRAIC[126]+CONSTANTS[150]+ALGEBRAIC[133])*STATES[41]; ALGEBRAIC[134] = ( 2.00000*CONSTANTS[74]*ALGEBRAIC[120])/CONSTANTS[75]; ALGEBRAIC[150] = ALGEBRAIC[134]*STATES[42]+ CONSTANTS[145]*STATES[36]+ ALGEBRAIC[125]*STATES[40]; ALGEBRAIC[127] = ( 2.00000*CONSTANTS[80]*ALGEBRAIC[119])/CONSTANTS[79]; ALGEBRAIC[151] = (ALGEBRAIC[127]+CONSTANTS[151]+ALGEBRAIC[134])*STATES[42]; ALGEBRAIC[135] = ( 3.00000*CONSTANTS[75]*ALGEBRAIC[120])/CONSTANTS[76]; ALGEBRAIC[152] = ALGEBRAIC[135]*STATES[43]+ CONSTANTS[146]*STATES[37]+ ALGEBRAIC[126]*STATES[41]; ALGEBRAIC[128] = ( CONSTANTS[81]*ALGEBRAIC[119])/CONSTANTS[80]; ALGEBRAIC[153] = (ALGEBRAIC[128]+CONSTANTS[152]+ALGEBRAIC[135])*STATES[43]; ALGEBRAIC[136] = ( 4.00000*CONSTANTS[76]*ALGEBRAIC[120])/CONSTANTS[77]; ALGEBRAIC[154] = ALGEBRAIC[136]*STATES[44]+ CONSTANTS[147]*STATES[38]+ ALGEBRAIC[127]*STATES[42]; ALGEBRAIC[155] = (CONSTANTS[153]+ALGEBRAIC[136])*STATES[44]; ALGEBRAIC[156] = CONSTANTS[148]*STATES[39]+ ALGEBRAIC[128]*STATES[43]; ALGEBRAIC[157] = CONSTANTS[124]*log(CONSTANTS[12]/STATES[0]); ALGEBRAIC[162] = CONSTANTS[91]*(STATES[45]+STATES[46])*(STATES[6] - ALGEBRAIC[157]); ALGEBRAIC[158] = CONSTANTS[124]*log(CONSTANTS[11]/STATES[1]); ALGEBRAIC[163] = CONSTANTS[88]*CONSTANTS[134]*STATES[47]*(STATES[6] - ALGEBRAIC[158]); ALGEBRAIC[164] = CONSTANTS[89]*(STATES[48]+STATES[49])*(STATES[6] - ALGEBRAIC[158]); ALGEBRAIC[165] = CONSTANTS[132]*STATES[29]*(STATES[6] - ALGEBRAIC[158]); ALGEBRAIC[166] = STATES[6]/CONSTANTS[124]; ALGEBRAIC[167] = 1000.00*CONSTANTS[3]*ALGEBRAIC[166]; ALGEBRAIC[168] = STATES[1]*exp(ALGEBRAIC[166]) - CONSTANTS[11]; ALGEBRAIC[169] = exp(ALGEBRAIC[166]) - 1.00000; ALGEBRAIC[170] = ( CONSTANTS[133]*STATES[39]*ALGEBRAIC[167]*ALGEBRAIC[168])/ALGEBRAIC[169]; ALGEBRAIC[171] = STATES[0]*exp(ALGEBRAIC[166]) - CONSTANTS[12]; ALGEBRAIC[172] = ( 0.0200000*CONSTANTS[133]*STATES[39]*ALGEBRAIC[167]*ALGEBRAIC[171])/ALGEBRAIC[169]; ALGEBRAIC[175] = 1.00000/(0.940000+exp( (1.26000/CONSTANTS[124])*(STATES[6] - ALGEBRAIC[158]))); ALGEBRAIC[176] = CONSTANTS[90]* pow((CONSTANTS[11]/1.00000), 1.0 / 2)*ALGEBRAIC[175]*(STATES[6] - ALGEBRAIC[158]); ALGEBRAIC[177] = CONSTANTS[92]*(STATES[6] - ALGEBRAIC[157]); ALGEBRAIC[161] = 0.500000*CONSTANTS[124]*log(CONSTANTS[13]/STATES[2]); ALGEBRAIC[178] = CONSTANTS[105]*(STATES[6] - ALGEBRAIC[161]); ALGEBRAIC[179] = ( CONSTANTS[103]*STATES[2])/(CONSTANTS[104]+STATES[2]); ALGEBRAIC[180] = STATES[6]/CONSTANTS[124]; ALGEBRAIC[181] = 1.00000+ 0.124500*exp( - 0.100000*ALGEBRAIC[180]); ALGEBRAIC[182] = 0.0365000*CONSTANTS[163]*exp( - 1.33000*ALGEBRAIC[180]); ALGEBRAIC[183] = 1.00000/(ALGEBRAIC[181]+ALGEBRAIC[182]); ALGEBRAIC[184] = 1.00000+pow(CONSTANTS[101]/STATES[0], 1.50000); ALGEBRAIC[185] = ( CONSTANTS[100]*ALGEBRAIC[183]*CONSTANTS[167])/ALGEBRAIC[184]; ALGEBRAIC[186] = exp( CONSTANTS[99]*ALGEBRAIC[180])*pow(STATES[0], 3.00000)*CONSTANTS[13]; ALGEBRAIC[187] = exp( (CONSTANTS[99] - 1.00000)*ALGEBRAIC[180])*pow(CONSTANTS[12], 3.00000)*STATES[2]; ALGEBRAIC[188] = 1.00000+ CONSTANTS[98]*exp( (CONSTANTS[99] - 1.00000)*ALGEBRAIC[180]); ALGEBRAIC[189] = ( CONSTANTS[95]*(ALGEBRAIC[186] - ALGEBRAIC[187]))/( CONSTANTS[172]*ALGEBRAIC[188]*CONSTANTS[175]); ALGEBRAIC[190] = STATES[6]/CONSTANTS[124]; ALGEBRAIC[191] = 1000.00*CONSTANTS[3]*ALGEBRAIC[190]; ALGEBRAIC[192] = 0.00100000*exp( 2.00000*ALGEBRAIC[190]) - CONSTANTS[13]*0.341000; ALGEBRAIC[193] = exp( 2.00000*ALGEBRAIC[190]) - 1.00000; ALGEBRAIC[194] = ( CONSTANTS[156]*4.00000*ALGEBRAIC[191]*ALGEBRAIC[192])/ALGEBRAIC[193]; ALGEBRAIC[195] = ALGEBRAIC[194]*STATES[24]*STATES[23]; ALGEBRAIC[196] = (CONDVAR[2]>=0.00000 ? 0.00000 : ALGEBRAIC[194]); ALGEBRAIC[197] = CONSTANTS[136]/(1.00000+ALGEBRAIC[196]/CONSTANTS[106]); ALGEBRAIC[198] = STATES[1]*exp(ALGEBRAIC[190]) - CONSTANTS[11]; ALGEBRAIC[199] = exp(ALGEBRAIC[190]) - 1.00000; ALGEBRAIC[200] = ( ALGEBRAIC[197]*STATES[23]*STATES[24]*ALGEBRAIC[191]*ALGEBRAIC[198])/ALGEBRAIC[199]; ALGEBRAIC[201] = CONSTANTS[161]*CONSTANTS[137]*exp(- 114007./CONSTANTS[157]+224.114/CONSTANTS[5]+( 0.286374*STATES[6])/CONSTANTS[138]); ALGEBRAIC[211] = 4.00000*ALGEBRAIC[201]; ALGEBRAIC[212] = 3.00000*ALGEBRAIC[201]; ALGEBRAIC[213] = 2.00000*ALGEBRAIC[201]; ALGEBRAIC[214] = ALGEBRAIC[201]; ALGEBRAIC[203] = CONSTANTS[161]*CONSTANTS[137]*exp(- 196337./CONSTANTS[157]+529.952/CONSTANTS[5]+( 2.78085*STATES[6])/CONSTANTS[138]); ALGEBRAIC[215] = ALGEBRAIC[203]; ALGEBRAIC[216] = ALGEBRAIC[211]*CONSTANTS[174]; ALGEBRAIC[217] = ALGEBRAIC[212]*CONSTANTS[174]; ALGEBRAIC[218] = ALGEBRAIC[213]*CONSTANTS[174]; ALGEBRAIC[219] = ALGEBRAIC[214]*CONSTANTS[174]; ALGEBRAIC[207] = CONSTANTS[161]*CONSTANTS[137]*exp(116431./CONSTANTS[157]+- 578.317/CONSTANTS[5]+( 0.764126*STATES[6])/CONSTANTS[138]); ALGEBRAIC[220] = ALGEBRAIC[207]; ALGEBRAIC[209] = CONSTANTS[161]*CONSTANTS[137]*exp(- 147814./CONSTANTS[157]+338.915/CONSTANTS[5]+( 2.13600*STATES[6])/CONSTANTS[138]); ALGEBRAIC[221] = ALGEBRAIC[209]; ALGEBRAIC[202] = CONSTANTS[161]*CONSTANTS[137]*exp(- 272470./CONSTANTS[157]+708.146/CONSTANTS[5]+( - 2.28528*STATES[6])/CONSTANTS[138]); ALGEBRAIC[222] = ALGEBRAIC[202]; ALGEBRAIC[223] = 2.00000*ALGEBRAIC[202]; ALGEBRAIC[224] = 3.00000*ALGEBRAIC[202]; ALGEBRAIC[225] = 4.00000*ALGEBRAIC[202]; ALGEBRAIC[204] = CONSTANTS[161]*CONSTANTS[137]*exp(- 133690./CONSTANTS[157]+229.205/CONSTANTS[5]+( - 1.55804*STATES[6])/CONSTANTS[138]); ALGEBRAIC[226] = ALGEBRAIC[204]; ALGEBRAIC[227] = ALGEBRAIC[222]/CONSTANTS[174]; ALGEBRAIC[228] = ALGEBRAIC[223]/CONSTANTS[174]; ALGEBRAIC[229] = ALGEBRAIC[224]/CONSTANTS[174]; ALGEBRAIC[230] = ALGEBRAIC[225]/CONSTANTS[174]; ALGEBRAIC[208] = CONSTANTS[161]*CONSTANTS[137]*exp(- 55700.7/CONSTANTS[157]+- 130.639/CONSTANTS[5]+( - 3.64982*STATES[6])/CONSTANTS[138]); ALGEBRAIC[231] = ALGEBRAIC[208]; ALGEBRAIC[210] = CONSTANTS[161]*CONSTANTS[137]*exp(- 121322./CONSTANTS[157]+193.265/CONSTANTS[5]+( - 1.74290*STATES[6])/CONSTANTS[138]); ALGEBRAIC[232] = ALGEBRAIC[210]; ALGEBRAIC[206] = CONSTANTS[161]*CONSTANTS[137]*exp(- 97657.8/CONSTANTS[157]+1.51000/CONSTANTS[5]+( 0.0684862*STATES[6])/CONSTANTS[138]); ALGEBRAIC[233] = ALGEBRAIC[206]; ALGEBRAIC[205] = CONSTANTS[161]*CONSTANTS[137]*exp(- 62123.1/CONSTANTS[157]+39.2950/CONSTANTS[5]+( 0.288816*STATES[6])/CONSTANTS[138]); ALGEBRAIC[234] = ALGEBRAIC[205]; ALGEBRAIC[235] = CONSTANTS[109]*CONSTANTS[110]*exp( CONSTANTS[111]*STATES[6]); ALGEBRAIC[236] = CONSTANTS[109]*CONSTANTS[112]*exp( CONSTANTS[113]*STATES[6]); ALGEBRAIC[243] = ALGEBRAIC[235]*STATES[61]+ CONSTANTS[140]*STATES[63]; ALGEBRAIC[244] = (ALGEBRAIC[236]+CONSTANTS[139])*STATES[62]; ALGEBRAIC[238] = CONSTANTS[109]*CONSTANTS[116]*exp( CONSTANTS[117]*STATES[6]); ALGEBRAIC[237] = CONSTANTS[109]*CONSTANTS[114]*exp( CONSTANTS[115]*STATES[6]); ALGEBRAIC[239] = CONSTANTS[109]*CONSTANTS[118]*exp( CONSTANTS[119]*STATES[6]); ALGEBRAIC[240] = CONSTANTS[109]*CONSTANTS[120]*exp( CONSTANTS[121]*STATES[6]); ALGEBRAIC[241] = CONSTANTS[109]*CONSTANTS[122]*exp( CONSTANTS[123]*STATES[6]); ALGEBRAIC[242] = ( ALGEBRAIC[238]*ALGEBRAIC[240]*ALGEBRAIC[241])/( ALGEBRAIC[237]*ALGEBRAIC[239]); ALGEBRAIC[245] = CONSTANTS[139]*STATES[62]+ ALGEBRAIC[238]*STATES[47]+ ALGEBRAIC[242]*STATES[64]; ALGEBRAIC[246] = (ALGEBRAIC[241]+ALGEBRAIC[237]+CONSTANTS[140])*STATES[63]; ALGEBRAIC[247] = ALGEBRAIC[237]*STATES[63]+ ALGEBRAIC[240]*STATES[64]; ALGEBRAIC[248] = (ALGEBRAIC[238]+ALGEBRAIC[239])*STATES[47]; ALGEBRAIC[249] = ALGEBRAIC[241]*STATES[63]+ ALGEBRAIC[239]*STATES[47]; ALGEBRAIC[250] = (ALGEBRAIC[242]+ALGEBRAIC[240])*STATES[64]; ALGEBRAIC[251] = 0.00767254*exp( 0.0866295*STATES[6]); ALGEBRAIC[252] = 0.00700807*exp( - 0.149998*STATES[6]); ALGEBRAIC[253] = 0.00379738*exp( - 0.0142567*STATES[6]); ALGEBRAIC[254] = 0.216256*exp( - 1.88912e-05*STATES[6]); ALGEBRAIC[173] = ALGEBRAIC[170]+ALGEBRAIC[172]; ALGEBRAIC[174] = ALGEBRAIC[165]+ALGEBRAIC[173]; ALGEBRAIC[255] = ALGEBRAIC[162]+ALGEBRAIC[195]+ALGEBRAIC[200]+ALGEBRAIC[163]+ALGEBRAIC[164]+ALGEBRAIC[176]+ALGEBRAIC[189]+ALGEBRAIC[185]+ALGEBRAIC[174]+ALGEBRAIC[179]+ALGEBRAIC[178]+ALGEBRAIC[177]+ALGEBRAIC[1]; } void getStateInformation(double* SI) { SI[0] = 1.0; SI[1] = 1.0; SI[2] = 1.0; SI[3] = 1.0; SI[4] = 1.0; SI[5] = 1.0; SI[6] = 1.0; SI[7] = 1.0; SI[8] = 1.0; SI[9] = 1.0; SI[10] = 1.0; SI[11] = 1.0; SI[12] = 1.0; SI[13] = 1.0; SI[14] = 1.0; SI[15] = 1.0; SI[16] = 1.0; SI[17] = 1.0; SI[18] = 1.0; SI[19] = 1.0; SI[20] = 1.0; SI[21] = 1.0; SI[22] = 1.0; SI[23] = 1.0; SI[24] = 1.0; SI[25] = 1.0; SI[26] = 1.0; SI[27] = 1.0; SI[28] = 1.0; SI[29] = 1.0; SI[30] = 1.0; SI[31] = 1.0; SI[32] = 1.0; SI[33] = 1.0; SI[34] = 1.0; SI[35] = 1.0; SI[36] = 1.0; SI[37] = 1.0; SI[38] = 1.0; SI[39] = 1.0; SI[40] = 1.0; SI[41] = 1.0; SI[42] = 1.0; SI[43] = 1.0; SI[44] = 1.0; SI[45] = 1.0; SI[46] = 1.0; SI[47] = 1.0; SI[48] = 1.0; SI[49] = 1.0; SI[50] = 1.0; SI[51] = 1.0; SI[52] = 1.0; SI[53] = 1.0; SI[54] = 1.0; SI[55] = 1.0; SI[56] = 1.0; SI[57] = 1.0; SI[58] = 1.0; SI[59] = 1.0; SI[60] = 1.0; SI[61] = 1.0; SI[62] = 1.0; SI[63] = 1.0; SI[64] = 1.0; SI[65] = 1.0; SI[66] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = (VOI - ALGEBRAIC[0]) - CONSTANTS[17]; CONDVAR[1] = (VOI - ALGEBRAIC[0]) - (CONSTANTS[17]+CONSTANTS[15]); CONDVAR[2] = ALGEBRAIC[194] - 0.00000; }