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 253 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 174 entries in the constant variable array.
*/
/*
* VOI is time in component environment (ms).
* CONSTANTS[143] is a1 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb).
* CONSTANTS[148] is a2 in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (F_M_per_coulomb).
* CONSTANTS[0] is Faraday in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (coulomb_per_millimole).
* CONSTANTS[1] is Temp in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (kelvin).
* CONSTANTS[2] is Rgas in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (joule_per_mole_kelvin).
* CONSTANTS[115] is RT_over_F in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mV).
* CONSTANTS[3] is Acap in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (cm2).
* CONSTANTS[123] is C in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mF).
* CONSTANTS[4] is Vmyo in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
* CONSTANTS[5] is VJSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
* CONSTANTS[6] is VNSR in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uL).
* CONSTANTS[7] 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[159] is INa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[174] is INab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[186] is INaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
* ALGEBRAIC[182] is INaK in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
* ALGEBRAIC[169] is IKv14_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[160] is IKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[161] is IKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[173] is IK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[197] is ICaK in component COMPUTE_ICa_ICaK (uA_per_uF).
* ALGEBRAIC[192] is ICa in component COMPUTE_ICa_ICaK (uA_per_uF).
* ALGEBRAIC[171] 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[162] is IKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[167] is IKv14_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* ALGEBRAIC[175] is ICab in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
* ALGEBRAIC[176] 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[252] is i_tot in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (uA_per_uF).
* CONSTANTS[8] is Ko in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
* CONSTANTS[9] is Nao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
* CONSTANTS[10] is Cao in component COMPUTE_CONCENTRATION_AND_VOLTAGE_DERIVATIVES (mM).
* CONSTANTS[11] is stim_period in component I_stimulus (ms).
* CONSTANTS[12] is stim_duration in component I_stimulus (ms).
* CONSTANTS[13] is stim_amplitude in component I_stimulus (uA_per_uF).
* CONSTANTS[14] 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[15] is Kfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
* CONSTANTS[16] is Nfb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
* ALGEBRAIC[3] is rb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
* CONSTANTS[17] is Krb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
* CONSTANTS[18] is Nrb in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (dimensionless).
* CONSTANTS[19] is KSR in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM).
* CONSTANTS[20] is vmaxf in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
* CONSTANTS[21] is vmaxr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (per_ms).
* CONSTANTS[22] 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[23] is tautr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms).
* CONSTANTS[24] is tauxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (ms).
* CONSTANTS[25] is LTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[26] is HTRPNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[27] is khtrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms).
* CONSTANTS[28] is khtrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
* CONSTANTS[29] is kltrpn_plus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_mM_per_ms).
* CONSTANTS[30] is kltrpn_minus in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
* CONSTANTS[31] is CMDNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[32] is CSQNtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[33] is EGTAtot in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[34] is KmCMDN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[35] is KmCSQN in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM).
* CONSTANTS[36] 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[37] is kaplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM4_per_ms).
* CONSTANTS[38] is kaminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* CONSTANTS[39] is kbplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_uM3_per_ms).
* CONSTANTS[40] is kbminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* CONSTANTS[41] is kcplus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* CONSTANTS[42] is kcminus in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* CONSTANTS[43] is ncoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
* CONSTANTS[44] is mcoop in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (dimensionless).
* ALGEBRAIC[14] is dC1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* ALGEBRAIC[15] is dO2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* ALGEBRAIC[16] is dC2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
* ALGEBRAIC[17] 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[45] is fL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[46] is gL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[47] is bL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
* CONSTANTS[48] is aL in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
* ALGEBRAIC[23] is C0_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[24] is C1_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[25] is C2_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[26] is C3_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[31] is C1_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[32] is C2_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[33] is C3_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[34] is C4_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[27] is CCa0_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[28] is CCa1_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[29] is CCa2_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[30] is CCa3_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[35] is CCa1_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[36] is CCa2_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[37] is CCa3_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[38] is CCa4_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[39] is C0_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[40] is C1_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[41] is C2_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[42] is C3_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[43] is C4_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[131] is CCa0_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[147] is CCa1_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[150] is CCa2_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[156] is CCa3_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[161] is CCa4_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[18] is alpha in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[19] is beta in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[20] is alpha_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[21] is beta_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[22] is gamma in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* CONSTANTS[116] is omega in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[54] is a1_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[55] is a2_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[56] is a1_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[57] is a2_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[58] is a1_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[59] is a2_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[60] is a1_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[61] is a2_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[62] is a1_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[63] is a2_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[44] is a1_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[45] is a2_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[46] is a1_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[47] is a2_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[48] is a1_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[49] is a2_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[50] is a1_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[51] is a2_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[52] is a1_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
* ALGEBRAIC[53] 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[117] is a1_Cainf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
* ALGEBRAIC[64] is yCa_inf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
* ALGEBRAIC[65] is tau_yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (ms).
* ALGEBRAIC[66] is alpha_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[67] is beta_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[68] is alpha_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[69] is beta_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[70] is C0Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[71] is C1Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[72] is C2Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[73] is C3Kv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[74] is CI0Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[75] is CI1Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[76] is CI2Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[77] is CI3Kv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[78] is C1Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[79] is C2Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[80] is C3Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[81] is OKv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[82] is CI1Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[83] is CI2Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[84] is CI3Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[85] is OIKv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[86] is C0Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[87] is C1Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[88] is C2Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[89] is C3Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[90] is OKv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[91] is CI0Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[92] is CI1Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[93] is CI2Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[94] is CI3Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[95] 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[96] is a1_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[97] is a2_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[98] is a1_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[99] is a2_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[100] is a1_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[101] is a2_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[102] is a1_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[103] is a2_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[104] is a1_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[105] is a2_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[106] is a1_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[107] is a2_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[108] is a1_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[109] is a2_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[110] is a1_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[111] is a2_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[112] is a1_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[113] is a2_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[114] is a1_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* ALGEBRAIC[115] is a2_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* CONSTANTS[49] is alphaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* CONSTANTS[50] is aaKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
* CONSTANTS[51] is betaa0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* CONSTANTS[52] is baKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
* CONSTANTS[53] is alphai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* CONSTANTS[54] is aiKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
* CONSTANTS[55] is betai0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
* CONSTANTS[56] is biKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_mV).
* CONSTANTS[57] is f1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[58] is f2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[59] is f3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[60] is f4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[61] is b1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[62] is b2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[63] is b3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[64] is b4Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (dimensionless).
* CONSTANTS[65] is f1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[66] is f2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[67] is f3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[68] is f4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[69] is b1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[70] is b2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[71] is b3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[72] is b4Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (dimensionless).
* CONSTANTS[73] is alphaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[74] is aaKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV).
* CONSTANTS[75] is betaa0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[76] is baKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_mV).
* CONSTANTS[77] is alphai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[78] is betai0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[116] is alpha_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[117] is beta_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[118] is alpha_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[119] is beta_inact14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[118] is C0Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[119] is C1Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[120] is C2Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[121] is C3Kv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[122] is CI0Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[123] is CI1Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[124] is CI2Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[125] is CI3Kv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[126] is C1Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[127] is C2Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[128] is C3Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[129] is OKv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[130] is CI1Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[131] is CI2Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[132] is CI3Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[133] is OIKv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[132] is C0Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[133] is C1Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[134] is C2Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[135] is C3Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[136] is OKv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[137] is CI0Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[138] is CI1Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[139] is CI2Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[140] is CI3Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* CONSTANTS[141] is OIKv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[134] is a1_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[135] 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[136] is a1_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[137] 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[138] is a1_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[139] 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[140] is a1_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[141] 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[142] is a1_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[143] 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[144] is a1_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[145] 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[146] is a1_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[147] 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[148] is a1_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[149] 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[150] is a1_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[151] 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[152] is a1_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
* ALGEBRAIC[153] 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[154] is ENa in component COMPUTE_REVERSAL_POTENTIALS (mV).
* ALGEBRAIC[155] is EK in component COMPUTE_REVERSAL_POTENTIALS (mV).
* ALGEBRAIC[157] is EKs in component COMPUTE_REVERSAL_POTENTIALS (mV).
* ALGEBRAIC[158] is ECa in component COMPUTE_REVERSAL_POTENTIALS (mV).
* CONSTANTS[142] is a1 in component COMPUTE_REVERSAL_POTENTIALS (mM).
* ALGEBRAIC[156] is a2 in component COMPUTE_REVERSAL_POTENTIALS (mM).
* CONSTANTS[79] is GKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
* CONSTANTS[80] is GKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
* CONSTANTS[81] is GK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
* CONSTANTS[82] is GNa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
* CONSTANTS[83] is GNab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
* CONSTANTS[84] is KvScale in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
* CONSTANTS[85] is Kv43Frac in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
* CONSTANTS[120] is GKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mS_per_uF).
* CONSTANTS[121] 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[122] is fKo in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
* ALGEBRAIC[163] is VF_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
* ALGEBRAIC[164] is VFsq_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (coulomb_per_millimole).
* ALGEBRAIC[165] is a1_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
* ALGEBRAIC[166] is a2 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
* ALGEBRAIC[168] is a1_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
* ALGEBRAIC[172] is K1_inf in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
* ALGEBRAIC[170] is IKv14 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
* CONSTANTS[86] is kNaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
* CONSTANTS[87] is KmNa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
* CONSTANTS[88] is KmCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
* CONSTANTS[89] is ksat in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* CONSTANTS[90] is eta in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* CONSTANTS[91] is INaKmax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
* CONSTANTS[92] is KmNai in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
* CONSTANTS[93] is KmKo in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
* CONSTANTS[94] is IpCamax in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
* CONSTANTS[95] is KmpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
* CONSTANTS[96] is GCab in component COMPUTE_INaK_INaCa_ICab_IpCa (mS_per_uF).
* ALGEBRAIC[177] is VF_over_RT in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* CONSTANTS[151] is sigma in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* ALGEBRAIC[178] is a1_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* ALGEBRAIC[179] is a2_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* ALGEBRAIC[180] is fNaK in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* CONSTANTS[155] is a1_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* ALGEBRAIC[181] is a2_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* ALGEBRAIC[183] is a1_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
* ALGEBRAIC[184] is a2_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
* ALGEBRAIC[185] is a3_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
* CONSTANTS[160] is a4_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM).
* CONSTANTS[163] is a5_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM3).
* ALGEBRAIC[194] is PKprime in component COMPUTE_ICa_ICaK (litre_per_farad_second).
* ALGEBRAIC[187] is VF_over_RT in component COMPUTE_ICa_ICaK (dimensionless).
* ALGEBRAIC[188] is VFsq_over_RT in component COMPUTE_ICa_ICaK (coulomb_per_millimole).
* ALGEBRAIC[189] is a1_Ca in component COMPUTE_ICa_ICaK (mM).
* ALGEBRAIC[190] is a2_Ca in component COMPUTE_ICa_ICaK (dimensionless).
* ALGEBRAIC[195] is a1_K in component COMPUTE_ICa_ICaK (mM).
* ALGEBRAIC[196] is a2_K in component COMPUTE_ICa_ICaK (dimensionless).
* ALGEBRAIC[191] is ICamax in component COMPUTE_ICa_ICaK (uA_per_uF).
* ALGEBRAIC[193] is Icabar in component COMPUTE_ICa_ICaK (uA_per_uF).
* CONSTANTS[144] is PCa in component COMPUTE_ICa_ICaK (litre_per_farad_second).
* CONSTANTS[124] is PK in component COMPUTE_ICa_ICaK (litre_per_farad_second).
* CONSTANTS[97] is ICahalf in component COMPUTE_ICa_ICaK (uA_per_uF).
* CONSTANTS[98] is Pscale in component COMPUTE_ICa_ICaK (dimensionless).
* CONSTANTS[149] is Temp_Scale in component INa (dimensionless).
* ALGEBRAIC[198] is alpha1 in component INa (per_ms).
* ALGEBRAIC[199] is beta1 in component INa (per_ms).
* ALGEBRAIC[200] is gamma1 in component INa (per_ms).
* ALGEBRAIC[201] is Delta1 in component INa (per_ms).
* ALGEBRAIC[202] is On in component INa (per_ms).
* ALGEBRAIC[203] is Of in component INa (per_ms).
* ALGEBRAIC[204] is GammaGamma in component INa (per_ms).
* ALGEBRAIC[205] is DeltaDelta in component INa (per_ms).
* CONSTANTS[153] is epsilon in component INa (per_ms).
* CONSTANTS[154] is omega_na in component INa (per_ms).
* ALGEBRAIC[206] is rho in component INa (per_ms).
* ALGEBRAIC[207] is mu in component INa (per_ms).
* CONSTANTS[152] is Cn in component INa (per_ms).
* CONSTANTS[157] is Cf in component INa (per_ms).
* CONSTANTS[162] is parameter_a in component INa (dimensionless).
* ALGEBRAIC[208] is k12 in component INa (per_ms).
* ALGEBRAIC[209] is k23 in component INa (per_ms).
* ALGEBRAIC[210] is k34 in component INa (per_ms).
* ALGEBRAIC[211] is k45 in component INa (per_ms).
* ALGEBRAIC[212] is k56 in component INa (per_ms).
* CONSTANTS[158] is k67 in component INa (per_ms).
* ALGEBRAIC[213] is k89 in component INa (per_ms).
* ALGEBRAIC[214] is k910 in component INa (per_ms).
* ALGEBRAIC[215] is k1011 in component INa (per_ms).
* ALGEBRAIC[216] is k1112 in component INa (per_ms).
* ALGEBRAIC[217] is k1213 in component INa (per_ms).
* ALGEBRAIC[218] is k57 in component INa (per_ms).
* ALGEBRAIC[219] is k21 in component INa (per_ms).
* ALGEBRAIC[220] is k32 in component INa (per_ms).
* ALGEBRAIC[221] is k43 in component INa (per_ms).
* ALGEBRAIC[222] is k54 in component INa (per_ms).
* ALGEBRAIC[223] is k65 in component INa (per_ms).
* CONSTANTS[159] is k76 in component INa (per_ms).
* ALGEBRAIC[224] is k98 in component INa (per_ms).
* ALGEBRAIC[225] is k109 in component INa (per_ms).
* ALGEBRAIC[226] is k1110 in component INa (per_ms).
* ALGEBRAIC[227] is k1211 in component INa (per_ms).
* ALGEBRAIC[228] is k1312 in component INa (per_ms).
* ALGEBRAIC[229] is k75 in component INa (per_ms).
* CONSTANTS[164] is k81 in component INa (per_ms).
* CONSTANTS[166] is k92 in component INa (per_ms).
* CONSTANTS[168] is k103 in component INa (per_ms).
* CONSTANTS[170] is k114 in component INa (per_ms).
* CONSTANTS[172] is k125 in component INa (per_ms).
* ALGEBRAIC[230] is k136 in component INa (per_ms).
* CONSTANTS[165] is k18 in component INa (per_ms).
* CONSTANTS[167] is k29 in component INa (per_ms).
* CONSTANTS[169] is k310 in component INa (per_ms).
* CONSTANTS[171] is k411 in component INa (per_ms).
* CONSTANTS[173] is k512 in component INa (per_ms).
* ALGEBRAIC[231] 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[99] is TNa in component INa (kelvin).
* CONSTANTS[125] is KToverH in component INa (per_ms).
* CONSTANTS[146] is FoverRT in component INa (per_mV).
* CONSTANTS[145] is RTNa in component INa (joule_per_mole).
* CONSTANTS[126] is RTNaF in component INa (mV).
* CONSTANTS[100] is T_Const_HERG in component IKr (dimensionless).
* CONSTANTS[101] is A0_HERG in component IKr (per_ms).
* CONSTANTS[102] is B0_HERG in component IKr (per_mV).
* CONSTANTS[103] is A1_HERG in component IKr (per_ms).
* CONSTANTS[104] is B1_HERG in component IKr (per_mV).
* CONSTANTS[105] is A2_HERG in component IKr (per_ms).
* CONSTANTS[106] is B2_HERG in component IKr (per_mV).
* CONSTANTS[107] is A3_HERG in component IKr (per_ms).
* CONSTANTS[108] is B3_HERG in component IKr (per_mV).
* CONSTANTS[109] is A4_HERG in component IKr (per_ms).
* CONSTANTS[110] is B4_HERG in component IKr (per_mV).
* CONSTANTS[111] is A5_HERG in component IKr (per_ms).
* CONSTANTS[112] is B5_HERG in component IKr (per_mV).
* CONSTANTS[113] is A6_HERG in component IKr (per_ms).
* CONSTANTS[114] is B6_HERG in component IKr (per_mV).
* CONSTANTS[127] is C2H_to_C3H in component IKr (per_ms).
* CONSTANTS[128] is C3H_to_C2H in component IKr (per_ms).
* ALGEBRAIC[232] is C1H_to_C2H in component IKr (per_ms).
* ALGEBRAIC[233] is C2H_to_C1H in component IKr (per_ms).
* ALGEBRAIC[234] is C3H_to_OH in component IKr (per_ms).
* ALGEBRAIC[235] is OH_to_C3H in component IKr (per_ms).
* ALGEBRAIC[236] is OH_to_IH in component IKr (per_ms).
* ALGEBRAIC[237] is IH_to_OH in component IKr (per_ms).
* ALGEBRAIC[238] is C3H_to_IH in component IKr (per_ms).
* ALGEBRAIC[239] is IH_to_C3H in component IKr (per_ms).
* STATES[61] is C1Herg in component IKr (dimensionless).
* ALGEBRAIC[240] is a1_C2 in component IKr (per_ms).
* ALGEBRAIC[241] is a2_C2 in component IKr (per_ms).
* STATES[62] is C2Herg in component IKr (dimensionless).
* ALGEBRAIC[242] is a1_C3 in component IKr (per_ms).
* ALGEBRAIC[243] is a2_C3 in component IKr (per_ms).
* STATES[63] is C3Herg in component IKr (dimensionless).
* ALGEBRAIC[244] is a1_O in component IKr (per_ms).
* ALGEBRAIC[245] is a2_O in component IKr (per_ms).
* ALGEBRAIC[246] is a1_I in component IKr (per_ms).
* ALGEBRAIC[247] is a2_I in component IKr (per_ms).
* STATES[64] is IHerg in component IKr (dimensionless).
* CONSTANTS[129] is C0ks_C1ks in component IKs (per_ms).
* CONSTANTS[130] is C1ks_O1ks in component IKs (per_ms).
* ALGEBRAIC[248] is O1ks_O2ks in component IKs (per_ms).
* ALGEBRAIC[249] is O1ks_C1ks in component IKs (per_ms).
* ALGEBRAIC[250] is O2ks_O1ks in component IKs (per_ms).
* ALGEBRAIC[251] 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] = 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] = pow((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] = (pow(CONSTANTS[87], 3.00000)+pow(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];
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[159]+ALGEBRAIC[174]+ 3.00000*(ALGEBRAIC[186]+ALGEBRAIC[182])+ALGEBRAIC[169])*CONSTANTS[143];
resid[1] = RATES[1] - - (((ALGEBRAIC[160]+ALGEBRAIC[161]+ALGEBRAIC[173]+ALGEBRAIC[197]+ALGEBRAIC[1]) - 2.00000*ALGEBRAIC[182])+ALGEBRAIC[162]+ALGEBRAIC[167])*CONSTANTS[143];
resid[2] = RATES[2] - ALGEBRAIC[13]*(((ALGEBRAIC[7] - ALGEBRAIC[4]) - ALGEBRAIC[10]) - ((ALGEBRAIC[175] - 2.00000*ALGEBRAIC[186])+ALGEBRAIC[176])*0.500000*CONSTANTS[143]);
resid[3] = RATES[3] - ALGEBRAIC[11]*((( ALGEBRAIC[5]*CONSTANTS[5])/CONSTANTS[7] - ( ALGEBRAIC[7]*CONSTANTS[4])/CONSTANTS[7]) - ALGEBRAIC[192]*CONSTANTS[148]);
resid[4] = RATES[4] - ALGEBRAIC[12]*(ALGEBRAIC[6] - ALGEBRAIC[5]);
resid[5] = RATES[5] - ( ALGEBRAIC[4]*CONSTANTS[4])/CONSTANTS[6] - ( ALGEBRAIC[6]*CONSTANTS[5])/CONSTANTS[6];
resid[6] = RATES[6] - - ALGEBRAIC[252];
resid[7] = RATES[9] - ALGEBRAIC[8];
resid[8] = RATES[10] - ALGEBRAIC[9];
resid[9] = RATES[11] - ALGEBRAIC[14];
resid[10] = RATES[8] - ALGEBRAIC[15];
resid[11] = RATES[12] - ALGEBRAIC[16];
resid[12] = RATES[7] - ALGEBRAIC[17];
resid[13] = RATES[13] - ALGEBRAIC[45] - ALGEBRAIC[44];
resid[14] = RATES[14] - ALGEBRAIC[47] - ALGEBRAIC[46];
resid[15] = RATES[15] - ALGEBRAIC[49] - ALGEBRAIC[48];
resid[16] = RATES[16] - ALGEBRAIC[51] - ALGEBRAIC[50];
resid[17] = RATES[17] - ALGEBRAIC[53] - ALGEBRAIC[52];
resid[18] = RATES[23] - CONSTANTS[45]*STATES[17] - CONSTANTS[46]*STATES[23];
resid[19] = RATES[18] - ALGEBRAIC[55] - ALGEBRAIC[54];
resid[20] = RATES[19] - ALGEBRAIC[57] - ALGEBRAIC[56];
resid[21] = RATES[20] - ALGEBRAIC[59] - ALGEBRAIC[58];
resid[22] = RATES[21] - ALGEBRAIC[61] - ALGEBRAIC[60];
resid[23] = RATES[22] - ALGEBRAIC[63] - ALGEBRAIC[62];
resid[24] = RATES[24] - (ALGEBRAIC[64] - STATES[24])/ALGEBRAIC[65];
resid[25] = RATES[25] - ALGEBRAIC[97] - ALGEBRAIC[96];
resid[26] = RATES[26] - ALGEBRAIC[99] - ALGEBRAIC[98];
resid[27] = RATES[27] - ALGEBRAIC[101] - ALGEBRAIC[100];
resid[28] = RATES[28] - ALGEBRAIC[103] - ALGEBRAIC[102];
resid[29] = RATES[29] - ALGEBRAIC[105] - ALGEBRAIC[104];
resid[30] = RATES[30] - ALGEBRAIC[107] - ALGEBRAIC[106];
resid[31] = RATES[31] - ALGEBRAIC[109] - ALGEBRAIC[108];
resid[32] = RATES[32] - ALGEBRAIC[111] - ALGEBRAIC[110];
resid[33] = RATES[33] - ALGEBRAIC[113] - ALGEBRAIC[112];
resid[34] = RATES[34] - ALGEBRAIC[115] - ALGEBRAIC[114];
resid[35] = RATES[35] - ALGEBRAIC[135] - ALGEBRAIC[134];
resid[36] = RATES[36] - ALGEBRAIC[137] - ALGEBRAIC[136];
resid[37] = RATES[37] - ALGEBRAIC[139] - ALGEBRAIC[138];
resid[38] = RATES[38] - ALGEBRAIC[141] - ALGEBRAIC[140];
resid[39] = RATES[39] - ALGEBRAIC[143] - ALGEBRAIC[142];
resid[40] = RATES[40] - ALGEBRAIC[145] - ALGEBRAIC[144];
resid[41] = RATES[41] - ALGEBRAIC[147] - ALGEBRAIC[146];
resid[42] = RATES[42] - ALGEBRAIC[149] - ALGEBRAIC[148];
resid[43] = RATES[43] - ALGEBRAIC[151] - ALGEBRAIC[150];
resid[44] = RATES[44] - ALGEBRAIC[153] - ALGEBRAIC[152];
resid[45] = RATES[50] - - (CONSTANTS[165]+ALGEBRAIC[208])*STATES[50]+ ALGEBRAIC[219]*STATES[51]+ CONSTANTS[164]*STATES[55];
resid[46] = RATES[51] - ( ALGEBRAIC[208]*STATES[50] - (ALGEBRAIC[219]+ALGEBRAIC[209]+CONSTANTS[167])*STATES[51])+ ALGEBRAIC[220]*STATES[52]+ CONSTANTS[166]*STATES[56];
resid[47] = RATES[52] - ( ALGEBRAIC[209]*STATES[51] - (ALGEBRAIC[220]+ALGEBRAIC[210]+CONSTANTS[169])*STATES[52])+ ALGEBRAIC[221]*STATES[53]+ CONSTANTS[168]*STATES[57];
resid[48] = RATES[53] - ( ALGEBRAIC[210]*STATES[52] - (ALGEBRAIC[221]+ALGEBRAIC[211]+CONSTANTS[171])*STATES[53])+ ALGEBRAIC[222]*STATES[54]+ CONSTANTS[170]*STATES[58];
resid[49] = RATES[54] - ( ALGEBRAIC[211]*STATES[53] - (ALGEBRAIC[222]+ALGEBRAIC[212]+ALGEBRAIC[218]+CONSTANTS[173])*STATES[54])+ ALGEBRAIC[223]*STATES[45]+ ALGEBRAIC[229]*STATES[46]+ CONSTANTS[172]*STATES[59];
resid[50] = RATES[45] - ( ALGEBRAIC[212]*STATES[54] - (ALGEBRAIC[223]+CONSTANTS[158]+ALGEBRAIC[231])*STATES[45])+ CONSTANTS[159]*STATES[46]+ ALGEBRAIC[230]*STATES[60];
resid[51] = RATES[46] - ( ALGEBRAIC[218]*STATES[54]+ CONSTANTS[158]*STATES[45]) - (ALGEBRAIC[229]+CONSTANTS[159])*STATES[46];
resid[52] = RATES[55] - ( CONSTANTS[165]*STATES[50] - (CONSTANTS[164]+ALGEBRAIC[213])*STATES[55])+ ALGEBRAIC[224]*STATES[56];
resid[53] = RATES[56] - (( CONSTANTS[167]*STATES[51]+ ALGEBRAIC[213]*STATES[55]) - (ALGEBRAIC[224]+CONSTANTS[166]+ALGEBRAIC[214])*STATES[56])+ ALGEBRAIC[225]*STATES[57];
resid[54] = RATES[57] - (( CONSTANTS[169]*STATES[52]+ ALGEBRAIC[214]*STATES[56]) - (ALGEBRAIC[215]+CONSTANTS[168]+ALGEBRAIC[225])*STATES[57])+ ALGEBRAIC[226]*STATES[58];
resid[55] = RATES[58] - (( CONSTANTS[171]*STATES[53]+ ALGEBRAIC[215]*STATES[57]) - (ALGEBRAIC[226]+CONSTANTS[170]+ALGEBRAIC[216])*STATES[58])+ ALGEBRAIC[227]*STATES[59];
resid[56] = RATES[59] - (( CONSTANTS[173]*STATES[54]+ ALGEBRAIC[216]*STATES[58]) - (ALGEBRAIC[227]+CONSTANTS[172]+ALGEBRAIC[217])*STATES[59])+ ALGEBRAIC[228]*STATES[60];
resid[57] = RATES[60] - ( ALGEBRAIC[231]*STATES[45]+ ALGEBRAIC[217]*STATES[59]) - (ALGEBRAIC[228]+ALGEBRAIC[230])*STATES[60];
resid[58] = RATES[61] - ALGEBRAIC[233]*STATES[62] - ALGEBRAIC[232]*STATES[61];
resid[59] = RATES[62] - ALGEBRAIC[240] - ALGEBRAIC[241];
resid[60] = RATES[63] - ALGEBRAIC[242] - ALGEBRAIC[243];
resid[61] = RATES[47] - ALGEBRAIC[244] - ALGEBRAIC[245];
resid[62] = RATES[64] - ALGEBRAIC[246] - ALGEBRAIC[247];
resid[63] = RATES[65] - - CONSTANTS[129]*STATES[65]+ ALGEBRAIC[251]*STATES[66];
resid[64] = RATES[66] - ( CONSTANTS[129]*STATES[65] - (ALGEBRAIC[251]+CONSTANTS[130])*STATES[66])+ ALGEBRAIC[249]*STATES[48];
resid[65] = RATES[48] - ( CONSTANTS[130]*STATES[66] - (ALGEBRAIC[249]+ALGEBRAIC[248])*STATES[48])+ ALGEBRAIC[250]*STATES[49];
resid[66] = RATES[49] - ALGEBRAIC[248]*STATES[48] - ALGEBRAIC[250]*STATES[49];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[156] = STATES[1]+ 0.0183300*STATES[0];
ALGEBRAIC[157] = CONSTANTS[115]*log(CONSTANTS[142]/ALGEBRAIC[156]);
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = floor(VOI/CONSTANTS[11])*CONSTANTS[11];
ALGEBRAIC[1] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<=0.00000 ? CONSTANTS[13] : 0.00000);
ALGEBRAIC[2] = pow(STATES[2]/CONSTANTS[15], CONSTANTS[16]);
ALGEBRAIC[3] = pow(STATES[5]/CONSTANTS[17], CONSTANTS[18]);
ALGEBRAIC[4] = ( CONSTANTS[19]*( CONSTANTS[20]*ALGEBRAIC[2] - CONSTANTS[21]*ALGEBRAIC[3]))/(1.00000+ALGEBRAIC[2]+ALGEBRAIC[3]);
ALGEBRAIC[5] = CONSTANTS[22]*(STATES[7]+STATES[8])*(STATES[4] - STATES[3]);
ALGEBRAIC[6] = (STATES[5] - STATES[4])/CONSTANTS[23];
ALGEBRAIC[7] = (STATES[3] - STATES[2])/CONSTANTS[24];
ALGEBRAIC[8] = CONSTANTS[29]*STATES[2]*(1.00000 - STATES[9]) - CONSTANTS[30]*STATES[9];
ALGEBRAIC[9] = CONSTANTS[27]*STATES[2]*(1.00000 - STATES[10]) - CONSTANTS[28]*STATES[10];
ALGEBRAIC[10] = CONSTANTS[25]*ALGEBRAIC[8]+ CONSTANTS[26]*ALGEBRAIC[9];
ALGEBRAIC[11] = 1.00000/(1.00000+( CONSTANTS[31]*CONSTANTS[34])/pow(STATES[3]+CONSTANTS[34], 2.00000)+( CONSTANTS[33]*CONSTANTS[36])/pow(STATES[3]+CONSTANTS[36], 2.00000));
ALGEBRAIC[12] = 1.00000/(1.00000+( CONSTANTS[32]*CONSTANTS[35])/pow(STATES[4]+CONSTANTS[35], 2.00000));
ALGEBRAIC[13] = 1.00000/(1.00000+( CONSTANTS[31]*CONSTANTS[34])/pow(STATES[2]+CONSTANTS[34], 2.00000)+( CONSTANTS[33]*CONSTANTS[36])/pow(STATES[2]+CONSTANTS[36], 2.00000));
ALGEBRAIC[14] = - CONSTANTS[37]*pow( STATES[3]*1000.00, CONSTANTS[43])*STATES[11]+ CONSTANTS[38]*STATES[7];
ALGEBRAIC[15] = CONSTANTS[39]*pow( STATES[3]*1000.00, CONSTANTS[44])*STATES[7] - CONSTANTS[40]*STATES[8];
ALGEBRAIC[16] = CONSTANTS[41]*STATES[7] - CONSTANTS[42]*STATES[12];
ALGEBRAIC[17] = - (ALGEBRAIC[14]+ALGEBRAIC[15]+ALGEBRAIC[16]);
ALGEBRAIC[18] = 4.00000*1.20000*0.416000*exp( 0.0120000*(STATES[6] - 35.0000));
ALGEBRAIC[23] = 4.00000*ALGEBRAIC[18];
ALGEBRAIC[22] = 0.600000*0.0923300*STATES[3];
ALGEBRAIC[39] = ALGEBRAIC[22];
ALGEBRAIC[44] = (ALGEBRAIC[23]+ALGEBRAIC[39])*STATES[13];
ALGEBRAIC[19] = 4.00000*0.450000*0.0490000*exp( - 0.0650000*(STATES[6] - 22.0000));
ALGEBRAIC[31] = ALGEBRAIC[19];
ALGEBRAIC[45] = ALGEBRAIC[31]*STATES[14]+ CONSTANTS[131]*STATES[18];
ALGEBRAIC[24] = 3.00000*ALGEBRAIC[18];
ALGEBRAIC[40] = CONSTANTS[48]*ALGEBRAIC[39];
ALGEBRAIC[46] = (ALGEBRAIC[31]+ALGEBRAIC[24]+ALGEBRAIC[40])*STATES[14];
ALGEBRAIC[32] = 2.00000*ALGEBRAIC[19];
ALGEBRAIC[47] = ALGEBRAIC[23]*STATES[13]+ ALGEBRAIC[32]*STATES[15]+ CONSTANTS[147]*STATES[19];
ALGEBRAIC[25] = 2.00000*ALGEBRAIC[18];
ALGEBRAIC[41] = CONSTANTS[48]*ALGEBRAIC[40];
ALGEBRAIC[48] = (ALGEBRAIC[32]+ALGEBRAIC[25]+ALGEBRAIC[41])*STATES[15];
ALGEBRAIC[33] = 3.00000*ALGEBRAIC[19];
ALGEBRAIC[49] = ALGEBRAIC[24]*STATES[14]+ ALGEBRAIC[33]*STATES[16]+ CONSTANTS[150]*STATES[20];
ALGEBRAIC[26] = ALGEBRAIC[18];
ALGEBRAIC[42] = CONSTANTS[48]*ALGEBRAIC[41];
ALGEBRAIC[50] = (ALGEBRAIC[33]+ALGEBRAIC[26]+ALGEBRAIC[42])*STATES[16];
ALGEBRAIC[34] = 4.00000*ALGEBRAIC[19];
ALGEBRAIC[51] = ALGEBRAIC[25]*STATES[15]+ ALGEBRAIC[34]*STATES[17]+ CONSTANTS[156]*STATES[21];
ALGEBRAIC[43] = CONSTANTS[48]*ALGEBRAIC[42];
ALGEBRAIC[52] = (ALGEBRAIC[34]+CONSTANTS[45]+ALGEBRAIC[43])*STATES[17];
ALGEBRAIC[53] = ALGEBRAIC[26]*STATES[16]+ CONSTANTS[46]*STATES[23]+ CONSTANTS[161]*STATES[22];
ALGEBRAIC[20] = CONSTANTS[48]*ALGEBRAIC[18];
ALGEBRAIC[27] = 4.00000*ALGEBRAIC[20];
ALGEBRAIC[54] = (ALGEBRAIC[27]+CONSTANTS[131])*STATES[18];
ALGEBRAIC[21] = ALGEBRAIC[19]/CONSTANTS[47];
ALGEBRAIC[35] = ALGEBRAIC[21];
ALGEBRAIC[55] = ALGEBRAIC[35]*STATES[19]+ ALGEBRAIC[39]*STATES[13];
ALGEBRAIC[28] = 3.00000*ALGEBRAIC[20];
ALGEBRAIC[56] = (ALGEBRAIC[35]+ALGEBRAIC[28]+CONSTANTS[147])*STATES[19];
ALGEBRAIC[36] = 2.00000*ALGEBRAIC[21];
ALGEBRAIC[57] = ALGEBRAIC[27]*STATES[18]+ ALGEBRAIC[36]*STATES[20]+ ALGEBRAIC[40]*STATES[14];
ALGEBRAIC[29] = 2.00000*ALGEBRAIC[20];
ALGEBRAIC[58] = (ALGEBRAIC[36]+ALGEBRAIC[29]+CONSTANTS[150])*STATES[20];
ALGEBRAIC[37] = 3.00000*ALGEBRAIC[21];
ALGEBRAIC[59] = ALGEBRAIC[28]*STATES[19]+ ALGEBRAIC[37]*STATES[21]+ ALGEBRAIC[41]*STATES[15];
ALGEBRAIC[30] = ALGEBRAIC[20];
ALGEBRAIC[60] = (ALGEBRAIC[37]+ALGEBRAIC[30]+CONSTANTS[156])*STATES[21];
ALGEBRAIC[38] = 4.00000*ALGEBRAIC[21];
ALGEBRAIC[61] = ALGEBRAIC[29]*STATES[20]+ ALGEBRAIC[38]*STATES[22]+ ALGEBRAIC[42]*STATES[16];
ALGEBRAIC[62] = (ALGEBRAIC[38]+CONSTANTS[161])*STATES[22];
ALGEBRAIC[63] = ALGEBRAIC[30]*STATES[21]+ ALGEBRAIC[43]*STATES[17];
ALGEBRAIC[64] = (CONSTANTS[117]/(1.00000+exp((STATES[6]+28.5000)/7.80000))+1.00000) - CONSTANTS[117];
ALGEBRAIC[65] = 1.00000/(0.00653000/(0.500000+exp(STATES[6]/- 7.10000))+ 0.00512000*exp(STATES[6]/- 39.8000));
ALGEBRAIC[66] = CONSTANTS[49]*exp( CONSTANTS[50]*STATES[6]);
ALGEBRAIC[70] = 4.00000*ALGEBRAIC[66];
ALGEBRAIC[69] = CONSTANTS[55]*exp( CONSTANTS[56]*STATES[6]);
ALGEBRAIC[86] = ALGEBRAIC[69];
ALGEBRAIC[96] = (ALGEBRAIC[70]+ALGEBRAIC[86])*STATES[25];
ALGEBRAIC[67] = CONSTANTS[51]*exp( - CONSTANTS[52]*STATES[6]);
ALGEBRAIC[78] = ALGEBRAIC[67];
ALGEBRAIC[68] = CONSTANTS[53]*exp( - CONSTANTS[54]*STATES[6]);
ALGEBRAIC[91] = ALGEBRAIC[68];
ALGEBRAIC[97] = ALGEBRAIC[78]*STATES[26]+ ALGEBRAIC[91]*STATES[30];
ALGEBRAIC[71] = 3.00000*ALGEBRAIC[66];
ALGEBRAIC[87] = CONSTANTS[57]*ALGEBRAIC[69];
ALGEBRAIC[98] = (ALGEBRAIC[71]+ALGEBRAIC[78]+ALGEBRAIC[87])*STATES[26];
ALGEBRAIC[79] = 2.00000*ALGEBRAIC[67];
ALGEBRAIC[92] = ALGEBRAIC[68]/CONSTANTS[61];
ALGEBRAIC[99] = ALGEBRAIC[79]*STATES[27]+ ALGEBRAIC[92]*STATES[31]+ ALGEBRAIC[70]*STATES[25];
ALGEBRAIC[72] = 2.00000*ALGEBRAIC[66];
ALGEBRAIC[88] = CONSTANTS[58]*ALGEBRAIC[69];
ALGEBRAIC[100] = (ALGEBRAIC[72]+ALGEBRAIC[79]+ALGEBRAIC[88])*STATES[27];
ALGEBRAIC[80] = 3.00000*ALGEBRAIC[67];
ALGEBRAIC[93] = ALGEBRAIC[68]/CONSTANTS[62];
ALGEBRAIC[101] = ALGEBRAIC[80]*STATES[28]+ ALGEBRAIC[93]*STATES[32]+ ALGEBRAIC[71]*STATES[26];
ALGEBRAIC[73] = ALGEBRAIC[66];
ALGEBRAIC[89] = CONSTANTS[59]*ALGEBRAIC[69];
ALGEBRAIC[102] = (ALGEBRAIC[73]+ALGEBRAIC[80]+ALGEBRAIC[89])*STATES[28];
ALGEBRAIC[81] = 4.00000*ALGEBRAIC[67];
ALGEBRAIC[94] = ALGEBRAIC[68]/CONSTANTS[63];
ALGEBRAIC[103] = ALGEBRAIC[81]*STATES[29]+ ALGEBRAIC[94]*STATES[33]+ ALGEBRAIC[72]*STATES[27];
ALGEBRAIC[90] = CONSTANTS[60]*ALGEBRAIC[69];
ALGEBRAIC[104] = (ALGEBRAIC[81]+ALGEBRAIC[90])*STATES[29];
ALGEBRAIC[95] = ALGEBRAIC[68]/CONSTANTS[64];
ALGEBRAIC[105] = ALGEBRAIC[73]*STATES[28]+ ALGEBRAIC[95]*STATES[34];
ALGEBRAIC[74] = 4.00000*CONSTANTS[61]*ALGEBRAIC[66];
ALGEBRAIC[106] = (ALGEBRAIC[91]+ALGEBRAIC[74])*STATES[30];
ALGEBRAIC[82] = ALGEBRAIC[67]/CONSTANTS[57];
ALGEBRAIC[107] = ALGEBRAIC[86]*STATES[25]+ ALGEBRAIC[82]*STATES[31];
ALGEBRAIC[75] = ( 3.00000*CONSTANTS[62]*ALGEBRAIC[66])/CONSTANTS[61];
ALGEBRAIC[108] = (ALGEBRAIC[75]+ALGEBRAIC[92]+ALGEBRAIC[82])*STATES[31];
ALGEBRAIC[83] = ( 2.00000*CONSTANTS[57]*ALGEBRAIC[67])/CONSTANTS[58];
ALGEBRAIC[109] = ALGEBRAIC[83]*STATES[32]+ ALGEBRAIC[87]*STATES[26]+ ALGEBRAIC[74]*STATES[30];
ALGEBRAIC[76] = ( 2.00000*CONSTANTS[63]*ALGEBRAIC[66])/CONSTANTS[62];
ALGEBRAIC[110] = (ALGEBRAIC[76]+ALGEBRAIC[93]+ALGEBRAIC[83])*STATES[32];
ALGEBRAIC[84] = ( 3.00000*CONSTANTS[58]*ALGEBRAIC[67])/CONSTANTS[59];
ALGEBRAIC[111] = ALGEBRAIC[84]*STATES[33]+ ALGEBRAIC[88]*STATES[27]+ ALGEBRAIC[75]*STATES[31];
ALGEBRAIC[77] = ( CONSTANTS[64]*ALGEBRAIC[66])/CONSTANTS[63];
ALGEBRAIC[112] = (ALGEBRAIC[77]+ALGEBRAIC[94]+ALGEBRAIC[84])*STATES[33];
ALGEBRAIC[85] = ( 4.00000*CONSTANTS[59]*ALGEBRAIC[67])/CONSTANTS[60];
ALGEBRAIC[113] = ALGEBRAIC[85]*STATES[34]+ ALGEBRAIC[89]*STATES[28]+ ALGEBRAIC[76]*STATES[32];
ALGEBRAIC[114] = (ALGEBRAIC[95]+ALGEBRAIC[85])*STATES[34];
ALGEBRAIC[115] = ALGEBRAIC[90]*STATES[29]+ ALGEBRAIC[77]*STATES[33];
ALGEBRAIC[116] = CONSTANTS[73]*exp( CONSTANTS[74]*STATES[6]);
ALGEBRAIC[118] = 4.00000*ALGEBRAIC[116];
ALGEBRAIC[134] = (ALGEBRAIC[118]+CONSTANTS[132])*STATES[35];
ALGEBRAIC[117] = CONSTANTS[75]*exp( - CONSTANTS[76]*STATES[6]);
ALGEBRAIC[126] = ALGEBRAIC[117];
ALGEBRAIC[135] = ALGEBRAIC[126]*STATES[36]+ CONSTANTS[137]*STATES[40];
ALGEBRAIC[119] = 3.00000*ALGEBRAIC[116];
ALGEBRAIC[136] = (ALGEBRAIC[119]+ALGEBRAIC[126]+CONSTANTS[133])*STATES[36];
ALGEBRAIC[127] = 2.00000*ALGEBRAIC[117];
ALGEBRAIC[137] = ALGEBRAIC[127]*STATES[37]+ CONSTANTS[138]*STATES[41]+ ALGEBRAIC[118]*STATES[35];
ALGEBRAIC[120] = 2.00000*ALGEBRAIC[116];
ALGEBRAIC[138] = (ALGEBRAIC[120]+ALGEBRAIC[127]+CONSTANTS[134])*STATES[37];
ALGEBRAIC[128] = 3.00000*ALGEBRAIC[117];
ALGEBRAIC[139] = ALGEBRAIC[128]*STATES[38]+ CONSTANTS[139]*STATES[42]+ ALGEBRAIC[119]*STATES[36];
ALGEBRAIC[121] = ALGEBRAIC[116];
ALGEBRAIC[140] = (ALGEBRAIC[121]+ALGEBRAIC[128]+CONSTANTS[135])*STATES[38];
ALGEBRAIC[129] = 4.00000*ALGEBRAIC[117];
ALGEBRAIC[141] = ALGEBRAIC[129]*STATES[39]+ CONSTANTS[140]*STATES[43]+ ALGEBRAIC[120]*STATES[37];
ALGEBRAIC[142] = (ALGEBRAIC[129]+CONSTANTS[136])*STATES[39];
ALGEBRAIC[143] = ALGEBRAIC[121]*STATES[38]+ CONSTANTS[141]*STATES[44];
ALGEBRAIC[122] = 4.00000*CONSTANTS[69]*ALGEBRAIC[116];
ALGEBRAIC[144] = (CONSTANTS[137]+ALGEBRAIC[122])*STATES[40];
ALGEBRAIC[130] = ALGEBRAIC[117]/CONSTANTS[65];
ALGEBRAIC[145] = CONSTANTS[132]*STATES[35]+ ALGEBRAIC[130]*STATES[41];
ALGEBRAIC[123] = ( 3.00000*CONSTANTS[70]*ALGEBRAIC[116])/CONSTANTS[69];
ALGEBRAIC[146] = (ALGEBRAIC[123]+CONSTANTS[138]+ALGEBRAIC[130])*STATES[41];
ALGEBRAIC[131] = ( 2.00000*CONSTANTS[65]*ALGEBRAIC[117])/CONSTANTS[66];
ALGEBRAIC[147] = ALGEBRAIC[131]*STATES[42]+ CONSTANTS[133]*STATES[36]+ ALGEBRAIC[122]*STATES[40];
ALGEBRAIC[124] = ( 2.00000*CONSTANTS[71]*ALGEBRAIC[116])/CONSTANTS[70];
ALGEBRAIC[148] = (ALGEBRAIC[124]+CONSTANTS[139]+ALGEBRAIC[131])*STATES[42];
ALGEBRAIC[132] = ( 3.00000*CONSTANTS[66]*ALGEBRAIC[117])/CONSTANTS[67];
ALGEBRAIC[149] = ALGEBRAIC[132]*STATES[43]+ CONSTANTS[134]*STATES[37]+ ALGEBRAIC[123]*STATES[41];
ALGEBRAIC[125] = ( CONSTANTS[72]*ALGEBRAIC[116])/CONSTANTS[71];
ALGEBRAIC[150] = (ALGEBRAIC[125]+CONSTANTS[140]+ALGEBRAIC[132])*STATES[43];
ALGEBRAIC[133] = ( 4.00000*CONSTANTS[67]*ALGEBRAIC[117])/CONSTANTS[68];
ALGEBRAIC[151] = ALGEBRAIC[133]*STATES[44]+ CONSTANTS[135]*STATES[38]+ ALGEBRAIC[124]*STATES[42];
ALGEBRAIC[152] = (CONSTANTS[141]+ALGEBRAIC[133])*STATES[44];
ALGEBRAIC[153] = CONSTANTS[136]*STATES[39]+ ALGEBRAIC[125]*STATES[43];
ALGEBRAIC[154] = CONSTANTS[115]*log(CONSTANTS[9]/STATES[0]);
ALGEBRAIC[159] = CONSTANTS[82]*(STATES[45]+STATES[46])*(STATES[6] - ALGEBRAIC[154]);
ALGEBRAIC[155] = CONSTANTS[115]*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[160] = CONSTANTS[79]*CONSTANTS[122]*STATES[47]*(STATES[6] - ALGEBRAIC[155]);
ALGEBRAIC[161] = CONSTANTS[80]*(STATES[48]+STATES[49])*(STATES[6] - ALGEBRAIC[155]);
ALGEBRAIC[162] = CONSTANTS[120]*STATES[29]*(STATES[6] - ALGEBRAIC[155]);
ALGEBRAIC[163] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[164] = 1000.00*CONSTANTS[0]*ALGEBRAIC[163];
ALGEBRAIC[165] = STATES[1]*exp(ALGEBRAIC[163]) - CONSTANTS[8];
ALGEBRAIC[166] = exp(ALGEBRAIC[163]) - 1.00000;
ALGEBRAIC[167] = ( CONSTANTS[121]*STATES[39]*ALGEBRAIC[164]*ALGEBRAIC[165])/ALGEBRAIC[166];
ALGEBRAIC[168] = STATES[0]*exp(ALGEBRAIC[163]) - CONSTANTS[9];
ALGEBRAIC[169] = ( 0.0200000*CONSTANTS[121]*STATES[39]*ALGEBRAIC[164]*ALGEBRAIC[168])/ALGEBRAIC[166];
ALGEBRAIC[172] = 1.00000/(0.940000+exp( (1.26000/CONSTANTS[115])*(STATES[6] - ALGEBRAIC[155])));
ALGEBRAIC[173] = CONSTANTS[81]* pow((CONSTANTS[8]/1.00000), 1.0 / 2)*ALGEBRAIC[172]*(STATES[6] - ALGEBRAIC[155]);
ALGEBRAIC[174] = CONSTANTS[83]*(STATES[6] - ALGEBRAIC[154]);
ALGEBRAIC[158] = 0.500000*CONSTANTS[115]*log(CONSTANTS[10]/STATES[2]);
ALGEBRAIC[175] = CONSTANTS[96]*(STATES[6] - ALGEBRAIC[158]);
ALGEBRAIC[176] = ( CONSTANTS[94]*STATES[2])/(CONSTANTS[95]+STATES[2]);
ALGEBRAIC[177] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[178] = 1.00000+ 0.124500*exp( - 0.100000*ALGEBRAIC[177]);
ALGEBRAIC[179] = 0.0365000*CONSTANTS[151]*exp( - 1.33000*ALGEBRAIC[177]);
ALGEBRAIC[180] = 1.00000/(ALGEBRAIC[178]+ALGEBRAIC[179]);
ALGEBRAIC[181] = 1.00000+pow(CONSTANTS[92]/STATES[0], 1.50000);
ALGEBRAIC[182] = ( CONSTANTS[91]*ALGEBRAIC[180]*CONSTANTS[155])/ALGEBRAIC[181];
ALGEBRAIC[183] = exp( CONSTANTS[90]*ALGEBRAIC[177])*pow(STATES[0], 3.00000)*CONSTANTS[10];
ALGEBRAIC[184] = exp( (CONSTANTS[90] - 1.00000)*ALGEBRAIC[177])*pow(CONSTANTS[9], 3.00000)*STATES[2];
ALGEBRAIC[185] = 1.00000+ CONSTANTS[89]*exp( (CONSTANTS[90] - 1.00000)*ALGEBRAIC[177]);
ALGEBRAIC[186] = ( CONSTANTS[86]*(ALGEBRAIC[183] - ALGEBRAIC[184]))/( CONSTANTS[160]*ALGEBRAIC[185]*CONSTANTS[163]);
ALGEBRAIC[187] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[188] = 1000.00*CONSTANTS[0]*ALGEBRAIC[187];
ALGEBRAIC[189] = 0.00100000*exp( 2.00000*ALGEBRAIC[187]) - CONSTANTS[10]*0.341000;
ALGEBRAIC[190] = exp( 2.00000*ALGEBRAIC[187]) - 1.00000;
ALGEBRAIC[191] = ( CONSTANTS[144]*4.00000*ALGEBRAIC[188]*ALGEBRAIC[189])/ALGEBRAIC[190];
ALGEBRAIC[192] = ALGEBRAIC[191]*STATES[24]*STATES[23];
ALGEBRAIC[193] = (CONDVAR[2]>=0.00000 ? 0.00000 : ALGEBRAIC[191]);
ALGEBRAIC[194] = CONSTANTS[124]/(1.00000+ALGEBRAIC[193]/CONSTANTS[97]);
ALGEBRAIC[195] = STATES[1]*exp(ALGEBRAIC[187]) - CONSTANTS[8];
ALGEBRAIC[196] = exp(ALGEBRAIC[187]) - 1.00000;
ALGEBRAIC[197] = ( ALGEBRAIC[194]*STATES[23]*STATES[24]*ALGEBRAIC[188]*ALGEBRAIC[195])/ALGEBRAIC[196];
ALGEBRAIC[198] = CONSTANTS[149]*CONSTANTS[125]*exp(- 114007./CONSTANTS[145]+224.114/CONSTANTS[2]+( 0.286374*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[208] = 4.00000*ALGEBRAIC[198];
ALGEBRAIC[209] = 3.00000*ALGEBRAIC[198];
ALGEBRAIC[210] = 2.00000*ALGEBRAIC[198];
ALGEBRAIC[211] = ALGEBRAIC[198];
ALGEBRAIC[200] = CONSTANTS[149]*CONSTANTS[125]*exp(- 196337./CONSTANTS[145]+529.952/CONSTANTS[2]+( 2.78085*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[212] = ALGEBRAIC[200];
ALGEBRAIC[213] = ALGEBRAIC[208]*CONSTANTS[162];
ALGEBRAIC[214] = ALGEBRAIC[209]*CONSTANTS[162];
ALGEBRAIC[215] = ALGEBRAIC[210]*CONSTANTS[162];
ALGEBRAIC[216] = ALGEBRAIC[211]*CONSTANTS[162];
ALGEBRAIC[204] = CONSTANTS[149]*CONSTANTS[125]*exp(116431./CONSTANTS[145]+- 578.317/CONSTANTS[2]+( 0.764126*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[217] = ALGEBRAIC[204];
ALGEBRAIC[206] = CONSTANTS[149]*CONSTANTS[125]*exp(- 147814./CONSTANTS[145]+338.915/CONSTANTS[2]+( 2.13600*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[218] = ALGEBRAIC[206];
ALGEBRAIC[199] = CONSTANTS[149]*CONSTANTS[125]*exp(- 272470./CONSTANTS[145]+708.146/CONSTANTS[2]+( - 2.28528*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[219] = ALGEBRAIC[199];
ALGEBRAIC[220] = 2.00000*ALGEBRAIC[199];
ALGEBRAIC[221] = 3.00000*ALGEBRAIC[199];
ALGEBRAIC[222] = 4.00000*ALGEBRAIC[199];
ALGEBRAIC[201] = CONSTANTS[149]*CONSTANTS[125]*exp(- 133690./CONSTANTS[145]+229.205/CONSTANTS[2]+( - 1.55804*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[223] = ALGEBRAIC[201];
ALGEBRAIC[224] = ALGEBRAIC[219]/CONSTANTS[162];
ALGEBRAIC[225] = ALGEBRAIC[220]/CONSTANTS[162];
ALGEBRAIC[226] = ALGEBRAIC[221]/CONSTANTS[162];
ALGEBRAIC[227] = ALGEBRAIC[222]/CONSTANTS[162];
ALGEBRAIC[205] = CONSTANTS[149]*CONSTANTS[125]*exp(- 55700.7/CONSTANTS[145]+- 130.639/CONSTANTS[2]+( - 3.64982*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[228] = ALGEBRAIC[205];
ALGEBRAIC[207] = CONSTANTS[149]*CONSTANTS[125]*exp(- 121322./CONSTANTS[145]+193.265/CONSTANTS[2]+( - 1.74290*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[229] = ALGEBRAIC[207];
ALGEBRAIC[203] = CONSTANTS[149]*CONSTANTS[125]*exp(- 97657.8/CONSTANTS[145]+1.51000/CONSTANTS[2]+( 0.0684862*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[230] = ALGEBRAIC[203];
ALGEBRAIC[202] = CONSTANTS[149]*CONSTANTS[125]*exp(- 62123.1/CONSTANTS[145]+39.2950/CONSTANTS[2]+( 0.288816*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[231] = ALGEBRAIC[202];
ALGEBRAIC[232] = CONSTANTS[100]*CONSTANTS[101]*exp( CONSTANTS[102]*STATES[6]);
ALGEBRAIC[233] = CONSTANTS[100]*CONSTANTS[103]*exp( CONSTANTS[104]*STATES[6]);
ALGEBRAIC[240] = ALGEBRAIC[232]*STATES[61]+ CONSTANTS[128]*STATES[63];
ALGEBRAIC[241] = (ALGEBRAIC[233]+CONSTANTS[127])*STATES[62];
ALGEBRAIC[235] = CONSTANTS[100]*CONSTANTS[107]*exp( CONSTANTS[108]*STATES[6]);
ALGEBRAIC[234] = CONSTANTS[100]*CONSTANTS[105]*exp( CONSTANTS[106]*STATES[6]);
ALGEBRAIC[236] = CONSTANTS[100]*CONSTANTS[109]*exp( CONSTANTS[110]*STATES[6]);
ALGEBRAIC[237] = CONSTANTS[100]*CONSTANTS[111]*exp( CONSTANTS[112]*STATES[6]);
ALGEBRAIC[238] = CONSTANTS[100]*CONSTANTS[113]*exp( CONSTANTS[114]*STATES[6]);
ALGEBRAIC[239] = ( ALGEBRAIC[235]*ALGEBRAIC[237]*ALGEBRAIC[238])/( ALGEBRAIC[234]*ALGEBRAIC[236]);
ALGEBRAIC[242] = CONSTANTS[127]*STATES[62]+ ALGEBRAIC[235]*STATES[47]+ ALGEBRAIC[239]*STATES[64];
ALGEBRAIC[243] = (ALGEBRAIC[238]+ALGEBRAIC[234]+CONSTANTS[128])*STATES[63];
ALGEBRAIC[244] = ALGEBRAIC[234]*STATES[63]+ ALGEBRAIC[237]*STATES[64];
ALGEBRAIC[245] = (ALGEBRAIC[235]+ALGEBRAIC[236])*STATES[47];
ALGEBRAIC[246] = ALGEBRAIC[238]*STATES[63]+ ALGEBRAIC[236]*STATES[47];
ALGEBRAIC[247] = (ALGEBRAIC[239]+ALGEBRAIC[237])*STATES[64];
ALGEBRAIC[248] = 0.00767254*exp( 0.0866295*STATES[6]);
ALGEBRAIC[249] = 0.00700807*exp( - 0.149998*STATES[6]);
ALGEBRAIC[250] = 0.00379738*exp( - 0.0142567*STATES[6]);
ALGEBRAIC[251] = 0.216256*exp( - 1.88912e-05*STATES[6]);
ALGEBRAIC[170] = ALGEBRAIC[167]+ALGEBRAIC[169];
ALGEBRAIC[171] = ALGEBRAIC[162]+ALGEBRAIC[170];
ALGEBRAIC[252] = ALGEBRAIC[159]+ALGEBRAIC[192]+ALGEBRAIC[197]+ALGEBRAIC[160]+ALGEBRAIC[161]+ALGEBRAIC[173]+ALGEBRAIC[186]+ALGEBRAIC[182]+ALGEBRAIC[171]+ALGEBRAIC[176]+ALGEBRAIC[175]+ALGEBRAIC[174]+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[14];
CONDVAR[1] = (VOI - ALGEBRAIC[0]) - (CONSTANTS[14]+CONSTANTS[12]);
CONDVAR[2] = ALGEBRAIC[191] - 0.00000;
}