Generated Code

The following is c 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[128] is INa in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[225] is INab in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[240] is INaCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
 * ALGEBRAIC[236] is INaK in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
 * ALGEBRAIC[203] is IKv14_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[132] is IKr in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[136] is IKs in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[223] is IK1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[251] is ICaK in component COMPUTE_ICa_ICaK (uA_per_uF).
 * ALGEBRAIC[246] is ICa in component COMPUTE_ICa_ICaK (uA_per_uF).
 * ALGEBRAIC[217] is Ito1 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[11] is i_Stim in component I_stimulus (uA_per_uF).
 * ALGEBRAIC[141] is IKv43 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[185] is IKv14_K in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (uA_per_uF).
 * ALGEBRAIC[227] is ICab in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
 * ALGEBRAIC[229] is IpCa in component COMPUTE_INaK_INaCa_ICab_IpCa (uA_per_uF).
 * ALGEBRAIC[58] is Jxfer in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
 * ALGEBRAIC[42] is Jup in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
 * ALGEBRAIC[88] is Jtrpn in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (mM_per_ms).
 * ALGEBRAIC[50] is Jrel in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
 * ALGEBRAIC[59] is Jtr in component COMPUTE_INTRACELLULAR_CALCIUM_FLUXES (mM_per_ms).
 * ALGEBRAIC[98] is beta_SS in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
 * ALGEBRAIC[68] is beta_JSR in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (dimensionless).
 * ALGEBRAIC[108] 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[22] 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[33] 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[67] is dLTRPNCa in component COMPUTE_Jtrpn_and_BUFFER_SCALE_FACTORS (per_ms).
 * ALGEBRAIC[78] 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[1] is dC1_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
 * ALGEBRAIC[12] is dO2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
 * ALGEBRAIC[23] is dC2_RyR in component COMPUTE_DERIVATIVES_OF_RyR_RECEPTOR_STATES (per_ms).
 * ALGEBRAIC[34] 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[51] is C0_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[60] is C1_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[69] is C2_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[79] is C3_to_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[121] is C1_to_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[125] is C2_to_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[129] is C3_to_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[133] is C4_to_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[89] is CCa0_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[99] is CCa1_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[109] is CCa2_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[117] is CCa3_to_CCa4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[137] is CCa1_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[142] is CCa2_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[148] is CCa3_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[156] is CCa4_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[166] is C0_to_CCa0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[175] is C1_to_CCa1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[186] is C2_to_CCa2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[197] is C3_to_CCa3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[204] 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[2] is alpha in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[13] is beta in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[24] is alpha_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[35] is beta_prime in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[43] 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[176] is a1_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[187] is a2_Ca0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[188] is a1_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[198] is a2_Ca1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[199] is a1_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[205] is a2_Ca2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[206] is a1_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[211] is a2_Ca3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[212] is a1_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[218] is a2_Ca4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[177] is a1_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[189] is a2_C0 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[190] is a1_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[200] is a2_C1 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[201] is a1_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[207] is a2_C2 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[208] is a1_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[213] is a2_C3 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[214] is a1_C4 in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[219] 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[3] is yCa_inf in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (dimensionless).
 * ALGEBRAIC[14] is tau_yCa in component COMPUTE_DERIVATIVES_OF_LTYPE_CHANNEL_STATES (ms).
 * ALGEBRAIC[4] is alpha_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[15] is beta_act43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[25] is alpha_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[36] is beta_inact43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[44] is C0Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[52] is C1Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[61] is C2Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[70] is C3Kv43_to_OKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[80] is CI0Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[90] is CI1Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[100] is CI2Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[110] is CI3Kv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[118] is C1Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[122] is C2Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[126] is C3Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[130] is OKv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[134] is CI1Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[139] is CI2Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[144] is CI3Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[150] is OIKv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[138] is C0Kv43_to_CI0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[143] is C1Kv43_to_CI1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[149] is C2Kv43_to_CI2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[157] is C3Kv43_to_CI3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[158] is OKv43_to_OIKv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[145] is CI0Kv43_to_C0Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[151] is CI1Kv43_to_C1Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[159] is CI2Kv43_to_C2Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[167] is CI3Kv43_to_C3Kv43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[168] 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[152] is a1_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[160] is a2_C043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[161] is a1_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[169] is a2_C143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[170] is a1_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[178] is a2_C243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[179] is a1_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[191] is a2_C343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[180] is a1_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[192] is a2_O43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[153] is a1_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[162] is a2_I043 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[163] is a1_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[171] is a2_I143 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[172] is a1_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[181] is a2_I243 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[182] is a1_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[193] is a2_I343 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[183] is a1_OI43 in component COMPUTE_DERIVATIVES_OF_Kv4_3_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[194] 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[5] is alpha_act14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[16] 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[26] is C0Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[37] is C1Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[45] is C2Kv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[53] is C3Kv14_to_OKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[27] is CI0Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[38] is CI1Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[46] is CI2Kv14_to_CI3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[54] is CI3Kv14_to_OIKv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[62] is C1Kv14_to_C0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[71] is C2Kv14_to_C1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[81] is C3Kv14_to_C2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[91] is OKv14_to_C3Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[63] is CI1Kv14_to_CI0Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[72] is CI2Kv14_to_CI1Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[82] is CI3Kv14_to_CI2Kv14 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[92] 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[73] is a1_C0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[83] 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[84] is a1_C1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[93] 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[94] is a1_C2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[101] 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[102] is a1_C3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[111] 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[103] is a1_O in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[112] 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[74] is a1_CI0 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[85] 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[86] is a1_CI1 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[95] 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[96] is a1_CI2 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[104] 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[105] is a1_CI3 in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[113] 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[106] is a1_OI in component COMPUTE_DERIVATIVES_OF_Kv1_4_CHANNEL_STATES (per_ms).
 * ALGEBRAIC[114] 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[116] is ENa in component COMPUTE_REVERSAL_POTENTIALS (mV).
 * ALGEBRAIC[120] is EK in component COMPUTE_REVERSAL_POTENTIALS (mV).
 * ALGEBRAIC[17] is EKs in component COMPUTE_REVERSAL_POTENTIALS (mV).
 * ALGEBRAIC[124] is ECa in component COMPUTE_REVERSAL_POTENTIALS (mV).
 * CONSTANTS[142] is a1 in component COMPUTE_REVERSAL_POTENTIALS (mM).
 * ALGEBRAIC[6] 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[147] is VF_over_RT in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
 * ALGEBRAIC[155] 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[174] is a2 in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
 * ALGEBRAIC[196] is a1_Na in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (mM).
 * ALGEBRAIC[221] is K1_inf in component COMPUTE_INa_IKr_IKs_Ito1_IK1_INab_IKp (dimensionless).
 * ALGEBRAIC[210] 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[231] 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[232] is a1_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
 * ALGEBRAIC[233] is a2_Na in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
 * ALGEBRAIC[234] is fNaK in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
 * CONSTANTS[155] is a1_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
 * ALGEBRAIC[235] is a2_K in component COMPUTE_INaK_INaCa_ICab_IpCa (dimensionless).
 * ALGEBRAIC[237] is a1_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
 * ALGEBRAIC[238] is a2_ncx in component COMPUTE_INaK_INaCa_ICab_IpCa (mM4).
 * ALGEBRAIC[239] 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[248] is PKprime in component COMPUTE_ICa_ICaK (litre_per_farad_second).
 * ALGEBRAIC[241] is VF_over_RT in component COMPUTE_ICa_ICaK (dimensionless).
 * ALGEBRAIC[242] is VFsq_over_RT in component COMPUTE_ICa_ICaK (coulomb_per_millimole).
 * ALGEBRAIC[243] is a1_Ca in component COMPUTE_ICa_ICaK (mM).
 * ALGEBRAIC[244] is a2_Ca in component COMPUTE_ICa_ICaK (dimensionless).
 * ALGEBRAIC[249] is a1_K in component COMPUTE_ICa_ICaK (mM).
 * ALGEBRAIC[250] is a2_K in component COMPUTE_ICa_ICaK (dimensionless).
 * ALGEBRAIC[245] is ICamax in component COMPUTE_ICa_ICaK (uA_per_uF).
 * ALGEBRAIC[247] 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[7] is alpha1 in component INa (per_ms).
 * ALGEBRAIC[18] is beta1 in component INa (per_ms).
 * ALGEBRAIC[28] is gamma1 in component INa (per_ms).
 * ALGEBRAIC[39] is Delta1 in component INa (per_ms).
 * ALGEBRAIC[47] is On in component INa (per_ms).
 * ALGEBRAIC[55] is Of in component INa (per_ms).
 * ALGEBRAIC[64] is GammaGamma in component INa (per_ms).
 * ALGEBRAIC[75] 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[87] is rho in component INa (per_ms).
 * ALGEBRAIC[97] 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[107] is k12 in component INa (per_ms).
 * ALGEBRAIC[115] is k23 in component INa (per_ms).
 * ALGEBRAIC[119] is k34 in component INa (per_ms).
 * ALGEBRAIC[123] is k45 in component INa (per_ms).
 * ALGEBRAIC[127] is k56 in component INa (per_ms).
 * CONSTANTS[158] is k67 in component INa (per_ms).
 * ALGEBRAIC[131] is k89 in component INa (per_ms).
 * ALGEBRAIC[135] is k910 in component INa (per_ms).
 * ALGEBRAIC[140] is k1011 in component INa (per_ms).
 * ALGEBRAIC[146] is k1112 in component INa (per_ms).
 * ALGEBRAIC[154] is k1213 in component INa (per_ms).
 * ALGEBRAIC[164] is k57 in component INa (per_ms).
 * ALGEBRAIC[173] is k21 in component INa (per_ms).
 * ALGEBRAIC[184] is k32 in component INa (per_ms).
 * ALGEBRAIC[195] is k43 in component INa (per_ms).
 * ALGEBRAIC[202] is k54 in component INa (per_ms).
 * ALGEBRAIC[209] is k65 in component INa (per_ms).
 * CONSTANTS[159] is k76 in component INa (per_ms).
 * ALGEBRAIC[215] is k98 in component INa (per_ms).
 * ALGEBRAIC[220] is k109 in component INa (per_ms).
 * ALGEBRAIC[222] is k1110 in component INa (per_ms).
 * ALGEBRAIC[224] is k1211 in component INa (per_ms).
 * ALGEBRAIC[226] is k1312 in component INa (per_ms).
 * ALGEBRAIC[216] 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[228] 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[230] 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[8] is C1H_to_C2H in component IKr (per_ms).
 * ALGEBRAIC[19] is C2H_to_C1H in component IKr (per_ms).
 * ALGEBRAIC[9] is C3H_to_OH in component IKr (per_ms).
 * ALGEBRAIC[20] is OH_to_C3H in component IKr (per_ms).
 * ALGEBRAIC[29] is OH_to_IH in component IKr (per_ms).
 * ALGEBRAIC[40] is IH_to_OH in component IKr (per_ms).
 * ALGEBRAIC[48] is C3H_to_IH in component IKr (per_ms).
 * ALGEBRAIC[56] is IH_to_C3H in component IKr (per_ms).
 * STATES[61] is C1Herg in component IKr (dimensionless).
 * ALGEBRAIC[30] is a1_C2 in component IKr (per_ms).
 * ALGEBRAIC[41] is a2_C2 in component IKr (per_ms).
 * STATES[62] is C2Herg in component IKr (dimensionless).
 * ALGEBRAIC[65] is a1_C3 in component IKr (per_ms).
 * ALGEBRAIC[76] is a2_C3 in component IKr (per_ms).
 * STATES[63] is C3Herg in component IKr (dimensionless).
 * ALGEBRAIC[49] is a1_O in component IKr (per_ms).
 * ALGEBRAIC[57] is a2_O in component IKr (per_ms).
 * ALGEBRAIC[66] is a1_I in component IKr (per_ms).
 * ALGEBRAIC[77] 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[10] is O1ks_O2ks in component IKs (per_ms).
 * ALGEBRAIC[21] is O1ks_C1ks in component IKs (per_ms).
 * ALGEBRAIC[31] is O2ks_O1ks in component IKs (per_ms).
 * ALGEBRAIC[32] 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).
 */
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];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[23] =  CONSTANTS[45]*STATES[17] -  CONSTANTS[46]*STATES[23];
ALGEBRAIC[1] =  - CONSTANTS[37]*pow( STATES[3]*1000.00, CONSTANTS[43])*STATES[11]+ CONSTANTS[38]*STATES[7];
RATES[11] = ALGEBRAIC[1];
ALGEBRAIC[12] =  CONSTANTS[39]*pow( STATES[3]*1000.00, CONSTANTS[44])*STATES[7] -  CONSTANTS[40]*STATES[8];
RATES[8] = ALGEBRAIC[12];
ALGEBRAIC[3] = (CONSTANTS[117]/(1.00000+exp((STATES[6]+28.5000)/7.80000))+1.00000) - CONSTANTS[117];
ALGEBRAIC[14] = 1.00000/(0.00653000/(0.500000+exp(STATES[6]/- 7.10000))+ 0.00512000*exp(STATES[6]/- 39.8000));
RATES[24] = (ALGEBRAIC[3] - STATES[24])/ALGEBRAIC[14];
ALGEBRAIC[8] =  CONSTANTS[100]*CONSTANTS[101]*exp( CONSTANTS[102]*STATES[6]);
ALGEBRAIC[19] =  CONSTANTS[100]*CONSTANTS[103]*exp( CONSTANTS[104]*STATES[6]);
RATES[61] =  ALGEBRAIC[19]*STATES[62] -  ALGEBRAIC[8]*STATES[61];
ALGEBRAIC[23] =  CONSTANTS[41]*STATES[7] -  CONSTANTS[42]*STATES[12];
RATES[12] = ALGEBRAIC[23];
ALGEBRAIC[32] =  0.216256*exp( - 1.88912e-05*STATES[6]);
RATES[65] =  - CONSTANTS[129]*STATES[65]+ ALGEBRAIC[32]*STATES[66];
ALGEBRAIC[21] =  0.00700807*exp( - 0.149998*STATES[6]);
RATES[66] = ( CONSTANTS[129]*STATES[65] -  (ALGEBRAIC[32]+CONSTANTS[130])*STATES[66])+ ALGEBRAIC[21]*STATES[48];
ALGEBRAIC[10] =  0.00767254*exp( 0.0866295*STATES[6]);
ALGEBRAIC[31] =  0.00379738*exp( - 0.0142567*STATES[6]);
RATES[48] = ( CONSTANTS[130]*STATES[66] -  (ALGEBRAIC[21]+ALGEBRAIC[10])*STATES[48])+ ALGEBRAIC[31]*STATES[49];
RATES[49] =  ALGEBRAIC[10]*STATES[48] -  ALGEBRAIC[31]*STATES[49];
ALGEBRAIC[34] = - (ALGEBRAIC[1]+ALGEBRAIC[12]+ALGEBRAIC[23]);
RATES[7] = ALGEBRAIC[34];
ALGEBRAIC[30] =  ALGEBRAIC[8]*STATES[61]+ CONSTANTS[128]*STATES[63];
ALGEBRAIC[41] =  (ALGEBRAIC[19]+CONSTANTS[127])*STATES[62];
RATES[62] = ALGEBRAIC[30] - ALGEBRAIC[41];
ALGEBRAIC[9] =  CONSTANTS[100]*CONSTANTS[105]*exp( CONSTANTS[106]*STATES[6]);
ALGEBRAIC[40] =  CONSTANTS[100]*CONSTANTS[111]*exp( CONSTANTS[112]*STATES[6]);
ALGEBRAIC[49] =  ALGEBRAIC[9]*STATES[63]+ ALGEBRAIC[40]*STATES[64];
ALGEBRAIC[20] =  CONSTANTS[100]*CONSTANTS[107]*exp( CONSTANTS[108]*STATES[6]);
ALGEBRAIC[29] =  CONSTANTS[100]*CONSTANTS[109]*exp( CONSTANTS[110]*STATES[6]);
ALGEBRAIC[57] =  (ALGEBRAIC[20]+ALGEBRAIC[29])*STATES[47];
RATES[47] = ALGEBRAIC[49] - ALGEBRAIC[57];
ALGEBRAIC[22] = pow(STATES[2]/CONSTANTS[15], CONSTANTS[16]);
ALGEBRAIC[33] = pow(STATES[5]/CONSTANTS[17], CONSTANTS[18]);
ALGEBRAIC[42] = ( CONSTANTS[19]*( CONSTANTS[20]*ALGEBRAIC[22] -  CONSTANTS[21]*ALGEBRAIC[33]))/(1.00000+ALGEBRAIC[22]+ALGEBRAIC[33]);
ALGEBRAIC[59] = (STATES[5] - STATES[4])/CONSTANTS[23];
RATES[5] = ( ALGEBRAIC[42]*CONSTANTS[4])/CONSTANTS[6] - ( ALGEBRAIC[59]*CONSTANTS[5])/CONSTANTS[6];
ALGEBRAIC[50] =  CONSTANTS[22]*(STATES[7]+STATES[8])*(STATES[4] - STATES[3]);
ALGEBRAIC[68] = 1.00000/(1.00000+( CONSTANTS[32]*CONSTANTS[35])/pow(STATES[4]+CONSTANTS[35], 2.00000));
RATES[4] =  ALGEBRAIC[68]*(ALGEBRAIC[59] - ALGEBRAIC[50]);
ALGEBRAIC[67] =  CONSTANTS[29]*STATES[2]*(1.00000 - STATES[9]) -  CONSTANTS[30]*STATES[9];
RATES[9] = ALGEBRAIC[67];
ALGEBRAIC[48] =  CONSTANTS[100]*CONSTANTS[113]*exp( CONSTANTS[114]*STATES[6]);
ALGEBRAIC[56] = ( ALGEBRAIC[20]*ALGEBRAIC[40]*ALGEBRAIC[48])/( ALGEBRAIC[9]*ALGEBRAIC[29]);
ALGEBRAIC[65] =  CONSTANTS[127]*STATES[62]+ ALGEBRAIC[20]*STATES[47]+ ALGEBRAIC[56]*STATES[64];
ALGEBRAIC[76] =  (ALGEBRAIC[48]+ALGEBRAIC[9]+CONSTANTS[128])*STATES[63];
RATES[63] = ALGEBRAIC[65] - ALGEBRAIC[76];
ALGEBRAIC[66] =  ALGEBRAIC[48]*STATES[63]+ ALGEBRAIC[29]*STATES[47];
ALGEBRAIC[77] =  (ALGEBRAIC[56]+ALGEBRAIC[40])*STATES[64];
RATES[64] = ALGEBRAIC[66] - ALGEBRAIC[77];
ALGEBRAIC[78] =  CONSTANTS[27]*STATES[2]*(1.00000 - STATES[10]) -  CONSTANTS[28]*STATES[10];
RATES[10] = ALGEBRAIC[78];
ALGEBRAIC[5] =  CONSTANTS[73]*exp( CONSTANTS[74]*STATES[6]);
ALGEBRAIC[26] =  4.00000*ALGEBRAIC[5];
ALGEBRAIC[73] =  (ALGEBRAIC[26]+CONSTANTS[132])*STATES[35];
ALGEBRAIC[16] =  CONSTANTS[75]*exp( - CONSTANTS[76]*STATES[6]);
ALGEBRAIC[62] = ALGEBRAIC[16];
ALGEBRAIC[83] =  ALGEBRAIC[62]*STATES[36]+ CONSTANTS[137]*STATES[40];
RATES[35] = ALGEBRAIC[83] - ALGEBRAIC[73];
ALGEBRAIC[27] =  4.00000*CONSTANTS[69]*ALGEBRAIC[5];
ALGEBRAIC[74] =  (CONSTANTS[137]+ALGEBRAIC[27])*STATES[40];
ALGEBRAIC[63] = ALGEBRAIC[16]/CONSTANTS[65];
ALGEBRAIC[85] =  CONSTANTS[132]*STATES[35]+ ALGEBRAIC[63]*STATES[41];
RATES[40] = ALGEBRAIC[85] - ALGEBRAIC[74];
ALGEBRAIC[37] =  3.00000*ALGEBRAIC[5];
ALGEBRAIC[84] =  (ALGEBRAIC[37]+ALGEBRAIC[62]+CONSTANTS[133])*STATES[36];
ALGEBRAIC[71] =  2.00000*ALGEBRAIC[16];
ALGEBRAIC[93] =  ALGEBRAIC[71]*STATES[37]+ CONSTANTS[138]*STATES[41]+ ALGEBRAIC[26]*STATES[35];
RATES[36] = ALGEBRAIC[93] - ALGEBRAIC[84];
ALGEBRAIC[38] = ( 3.00000*CONSTANTS[70]*ALGEBRAIC[5])/CONSTANTS[69];
ALGEBRAIC[86] =  (ALGEBRAIC[38]+CONSTANTS[138]+ALGEBRAIC[63])*STATES[41];
ALGEBRAIC[72] = ( 2.00000*CONSTANTS[65]*ALGEBRAIC[16])/CONSTANTS[66];
ALGEBRAIC[95] =  ALGEBRAIC[72]*STATES[42]+ CONSTANTS[133]*STATES[36]+ ALGEBRAIC[27]*STATES[40];
RATES[41] = ALGEBRAIC[95] - ALGEBRAIC[86];
ALGEBRAIC[45] =  2.00000*ALGEBRAIC[5];
ALGEBRAIC[94] =  (ALGEBRAIC[45]+ALGEBRAIC[71]+CONSTANTS[134])*STATES[37];
ALGEBRAIC[81] =  3.00000*ALGEBRAIC[16];
ALGEBRAIC[101] =  ALGEBRAIC[81]*STATES[38]+ CONSTANTS[139]*STATES[42]+ ALGEBRAIC[37]*STATES[36];
RATES[37] = ALGEBRAIC[101] - ALGEBRAIC[94];
ALGEBRAIC[46] = ( 2.00000*CONSTANTS[71]*ALGEBRAIC[5])/CONSTANTS[70];
ALGEBRAIC[96] =  (ALGEBRAIC[46]+CONSTANTS[139]+ALGEBRAIC[72])*STATES[42];
ALGEBRAIC[82] = ( 3.00000*CONSTANTS[66]*ALGEBRAIC[16])/CONSTANTS[67];
ALGEBRAIC[104] =  ALGEBRAIC[82]*STATES[43]+ CONSTANTS[134]*STATES[37]+ ALGEBRAIC[38]*STATES[41];
RATES[42] = ALGEBRAIC[104] - ALGEBRAIC[96];
ALGEBRAIC[53] = ALGEBRAIC[5];
ALGEBRAIC[102] =  (ALGEBRAIC[53]+ALGEBRAIC[81]+CONSTANTS[135])*STATES[38];
ALGEBRAIC[91] =  4.00000*ALGEBRAIC[16];
ALGEBRAIC[111] =  ALGEBRAIC[91]*STATES[39]+ CONSTANTS[140]*STATES[43]+ ALGEBRAIC[45]*STATES[37];
RATES[38] = ALGEBRAIC[111] - ALGEBRAIC[102];
ALGEBRAIC[103] =  (ALGEBRAIC[91]+CONSTANTS[136])*STATES[39];
ALGEBRAIC[112] =  ALGEBRAIC[53]*STATES[38]+ CONSTANTS[141]*STATES[44];
RATES[39] = ALGEBRAIC[112] - ALGEBRAIC[103];
ALGEBRAIC[54] = ( CONSTANTS[72]*ALGEBRAIC[5])/CONSTANTS[71];
ALGEBRAIC[105] =  (ALGEBRAIC[54]+CONSTANTS[140]+ALGEBRAIC[82])*STATES[43];
ALGEBRAIC[92] = ( 4.00000*CONSTANTS[67]*ALGEBRAIC[16])/CONSTANTS[68];
ALGEBRAIC[113] =  ALGEBRAIC[92]*STATES[44]+ CONSTANTS[135]*STATES[38]+ ALGEBRAIC[46]*STATES[42];
RATES[43] = ALGEBRAIC[113] - ALGEBRAIC[105];
ALGEBRAIC[106] =  (CONSTANTS[141]+ALGEBRAIC[92])*STATES[44];
ALGEBRAIC[114] =  CONSTANTS[136]*STATES[39]+ ALGEBRAIC[54]*STATES[43];
RATES[44] = ALGEBRAIC[114] - ALGEBRAIC[106];
ALGEBRAIC[4] =  CONSTANTS[49]*exp( CONSTANTS[50]*STATES[6]);
ALGEBRAIC[44] =  4.00000*ALGEBRAIC[4];
ALGEBRAIC[36] =  CONSTANTS[55]*exp( CONSTANTS[56]*STATES[6]);
ALGEBRAIC[138] = ALGEBRAIC[36];
ALGEBRAIC[152] =  (ALGEBRAIC[44]+ALGEBRAIC[138])*STATES[25];
ALGEBRAIC[15] =  CONSTANTS[51]*exp( - CONSTANTS[52]*STATES[6]);
ALGEBRAIC[118] = ALGEBRAIC[15];
ALGEBRAIC[25] =  CONSTANTS[53]*exp( - CONSTANTS[54]*STATES[6]);
ALGEBRAIC[145] = ALGEBRAIC[25];
ALGEBRAIC[160] =  ALGEBRAIC[118]*STATES[26]+ ALGEBRAIC[145]*STATES[30];
RATES[25] = ALGEBRAIC[160] - ALGEBRAIC[152];
ALGEBRAIC[80] =  4.00000*CONSTANTS[61]*ALGEBRAIC[4];
ALGEBRAIC[153] =  (ALGEBRAIC[145]+ALGEBRAIC[80])*STATES[30];
ALGEBRAIC[134] = ALGEBRAIC[15]/CONSTANTS[57];
ALGEBRAIC[162] =  ALGEBRAIC[138]*STATES[25]+ ALGEBRAIC[134]*STATES[31];
RATES[30] = ALGEBRAIC[162] - ALGEBRAIC[153];
ALGEBRAIC[52] =  3.00000*ALGEBRAIC[4];
ALGEBRAIC[143] =  CONSTANTS[57]*ALGEBRAIC[36];
ALGEBRAIC[161] =  (ALGEBRAIC[52]+ALGEBRAIC[118]+ALGEBRAIC[143])*STATES[26];
ALGEBRAIC[122] =  2.00000*ALGEBRAIC[15];
ALGEBRAIC[151] = ALGEBRAIC[25]/CONSTANTS[61];
ALGEBRAIC[169] =  ALGEBRAIC[122]*STATES[27]+ ALGEBRAIC[151]*STATES[31]+ ALGEBRAIC[44]*STATES[25];
RATES[26] = ALGEBRAIC[169] - ALGEBRAIC[161];
ALGEBRAIC[90] = ( 3.00000*CONSTANTS[62]*ALGEBRAIC[4])/CONSTANTS[61];
ALGEBRAIC[163] =  (ALGEBRAIC[90]+ALGEBRAIC[151]+ALGEBRAIC[134])*STATES[31];
ALGEBRAIC[139] = ( 2.00000*CONSTANTS[57]*ALGEBRAIC[15])/CONSTANTS[58];
ALGEBRAIC[171] =  ALGEBRAIC[139]*STATES[32]+ ALGEBRAIC[143]*STATES[26]+ ALGEBRAIC[80]*STATES[30];
RATES[31] = ALGEBRAIC[171] - ALGEBRAIC[163];
ALGEBRAIC[7] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 114007./CONSTANTS[145]+224.114/CONSTANTS[2]+( 0.286374*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[107] =  4.00000*ALGEBRAIC[7];
ALGEBRAIC[18] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 272470./CONSTANTS[145]+708.146/CONSTANTS[2]+( - 2.28528*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[173] = ALGEBRAIC[18];
RATES[50] =  - (CONSTANTS[165]+ALGEBRAIC[107])*STATES[50]+ ALGEBRAIC[173]*STATES[51]+ CONSTANTS[164]*STATES[55];
ALGEBRAIC[61] =  2.00000*ALGEBRAIC[4];
ALGEBRAIC[149] =  CONSTANTS[58]*ALGEBRAIC[36];
ALGEBRAIC[170] =  (ALGEBRAIC[61]+ALGEBRAIC[122]+ALGEBRAIC[149])*STATES[27];
ALGEBRAIC[126] =  3.00000*ALGEBRAIC[15];
ALGEBRAIC[159] = ALGEBRAIC[25]/CONSTANTS[62];
ALGEBRAIC[178] =  ALGEBRAIC[126]*STATES[28]+ ALGEBRAIC[159]*STATES[32]+ ALGEBRAIC[52]*STATES[26];
RATES[27] = ALGEBRAIC[178] - ALGEBRAIC[170];
ALGEBRAIC[100] = ( 2.00000*CONSTANTS[63]*ALGEBRAIC[4])/CONSTANTS[62];
ALGEBRAIC[172] =  (ALGEBRAIC[100]+ALGEBRAIC[159]+ALGEBRAIC[139])*STATES[32];
ALGEBRAIC[144] = ( 3.00000*CONSTANTS[58]*ALGEBRAIC[15])/CONSTANTS[59];
ALGEBRAIC[181] =  ALGEBRAIC[144]*STATES[33]+ ALGEBRAIC[149]*STATES[27]+ ALGEBRAIC[90]*STATES[31];
RATES[32] = ALGEBRAIC[181] - ALGEBRAIC[172];
ALGEBRAIC[115] =  3.00000*ALGEBRAIC[7];
ALGEBRAIC[184] =  2.00000*ALGEBRAIC[18];
RATES[51] = ( ALGEBRAIC[107]*STATES[50] -  (ALGEBRAIC[173]+ALGEBRAIC[115]+CONSTANTS[167])*STATES[51])+ ALGEBRAIC[184]*STATES[52]+ CONSTANTS[166]*STATES[56];
ALGEBRAIC[2] =  4.00000*1.20000*0.416000*exp( 0.0120000*(STATES[6] - 35.0000));
ALGEBRAIC[51] =  4.00000*ALGEBRAIC[2];
ALGEBRAIC[43] =  0.600000*0.0923300*STATES[3];
ALGEBRAIC[166] = ALGEBRAIC[43];
ALGEBRAIC[177] =  (ALGEBRAIC[51]+ALGEBRAIC[166])*STATES[13];
ALGEBRAIC[13] =  4.00000*0.450000*0.0490000*exp( - 0.0650000*(STATES[6] - 22.0000));
ALGEBRAIC[121] = ALGEBRAIC[13];
ALGEBRAIC[189] =  ALGEBRAIC[121]*STATES[14]+ CONSTANTS[131]*STATES[18];
RATES[13] = ALGEBRAIC[189] - ALGEBRAIC[177];
ALGEBRAIC[24] =  CONSTANTS[48]*ALGEBRAIC[2];
ALGEBRAIC[89] =  4.00000*ALGEBRAIC[24];
ALGEBRAIC[176] =  (ALGEBRAIC[89]+CONSTANTS[131])*STATES[18];
ALGEBRAIC[35] = ALGEBRAIC[13]/CONSTANTS[47];
ALGEBRAIC[137] = ALGEBRAIC[35];
ALGEBRAIC[187] =  ALGEBRAIC[137]*STATES[19]+ ALGEBRAIC[166]*STATES[13];
RATES[18] = ALGEBRAIC[187] - ALGEBRAIC[176];
ALGEBRAIC[70] = ALGEBRAIC[4];
ALGEBRAIC[157] =  CONSTANTS[59]*ALGEBRAIC[36];
ALGEBRAIC[179] =  (ALGEBRAIC[70]+ALGEBRAIC[126]+ALGEBRAIC[157])*STATES[28];
ALGEBRAIC[130] =  4.00000*ALGEBRAIC[15];
ALGEBRAIC[167] = ALGEBRAIC[25]/CONSTANTS[63];
ALGEBRAIC[191] =  ALGEBRAIC[130]*STATES[29]+ ALGEBRAIC[167]*STATES[33]+ ALGEBRAIC[61]*STATES[27];
RATES[28] = ALGEBRAIC[191] - ALGEBRAIC[179];
ALGEBRAIC[158] =  CONSTANTS[60]*ALGEBRAIC[36];
ALGEBRAIC[180] =  (ALGEBRAIC[130]+ALGEBRAIC[158])*STATES[29];
ALGEBRAIC[168] = ALGEBRAIC[25]/CONSTANTS[64];
ALGEBRAIC[192] =  ALGEBRAIC[70]*STATES[28]+ ALGEBRAIC[168]*STATES[34];
RATES[29] = ALGEBRAIC[192] - ALGEBRAIC[180];
ALGEBRAIC[110] = ( CONSTANTS[64]*ALGEBRAIC[4])/CONSTANTS[63];
ALGEBRAIC[182] =  (ALGEBRAIC[110]+ALGEBRAIC[167]+ALGEBRAIC[144])*STATES[33];
ALGEBRAIC[150] = ( 4.00000*CONSTANTS[59]*ALGEBRAIC[15])/CONSTANTS[60];
ALGEBRAIC[193] =  ALGEBRAIC[150]*STATES[34]+ ALGEBRAIC[157]*STATES[28]+ ALGEBRAIC[100]*STATES[32];
RATES[33] = ALGEBRAIC[193] - ALGEBRAIC[182];
ALGEBRAIC[183] =  (ALGEBRAIC[168]+ALGEBRAIC[150])*STATES[34];
ALGEBRAIC[194] =  ALGEBRAIC[158]*STATES[29]+ ALGEBRAIC[110]*STATES[33];
RATES[34] = ALGEBRAIC[194] - ALGEBRAIC[183];
ALGEBRAIC[119] =  2.00000*ALGEBRAIC[7];
ALGEBRAIC[195] =  3.00000*ALGEBRAIC[18];
RATES[52] = ( ALGEBRAIC[115]*STATES[51] -  (ALGEBRAIC[184]+ALGEBRAIC[119]+CONSTANTS[169])*STATES[52])+ ALGEBRAIC[195]*STATES[53]+ CONSTANTS[168]*STATES[57];
ALGEBRAIC[60] =  3.00000*ALGEBRAIC[2];
ALGEBRAIC[175] =  CONSTANTS[48]*ALGEBRAIC[166];
ALGEBRAIC[190] =  (ALGEBRAIC[121]+ALGEBRAIC[60]+ALGEBRAIC[175])*STATES[14];
ALGEBRAIC[125] =  2.00000*ALGEBRAIC[13];
ALGEBRAIC[200] =  ALGEBRAIC[51]*STATES[13]+ ALGEBRAIC[125]*STATES[15]+ CONSTANTS[147]*STATES[19];
RATES[14] = ALGEBRAIC[200] - ALGEBRAIC[190];
ALGEBRAIC[99] =  3.00000*ALGEBRAIC[24];
ALGEBRAIC[188] =  (ALGEBRAIC[137]+ALGEBRAIC[99]+CONSTANTS[147])*STATES[19];
ALGEBRAIC[142] =  2.00000*ALGEBRAIC[35];
ALGEBRAIC[198] =  ALGEBRAIC[89]*STATES[18]+ ALGEBRAIC[142]*STATES[20]+ ALGEBRAIC[175]*STATES[14];
RATES[19] = ALGEBRAIC[198] - ALGEBRAIC[188];
ALGEBRAIC[123] = ALGEBRAIC[7];
ALGEBRAIC[202] =  4.00000*ALGEBRAIC[18];
RATES[53] = ( ALGEBRAIC[119]*STATES[52] -  (ALGEBRAIC[195]+ALGEBRAIC[123]+CONSTANTS[171])*STATES[53])+ ALGEBRAIC[202]*STATES[54]+ CONSTANTS[170]*STATES[58];
ALGEBRAIC[69] =  2.00000*ALGEBRAIC[2];
ALGEBRAIC[186] =  CONSTANTS[48]*ALGEBRAIC[175];
ALGEBRAIC[201] =  (ALGEBRAIC[125]+ALGEBRAIC[69]+ALGEBRAIC[186])*STATES[15];
ALGEBRAIC[129] =  3.00000*ALGEBRAIC[13];
ALGEBRAIC[207] =  ALGEBRAIC[60]*STATES[14]+ ALGEBRAIC[129]*STATES[16]+ CONSTANTS[150]*STATES[20];
RATES[15] = ALGEBRAIC[207] - ALGEBRAIC[201];
ALGEBRAIC[109] =  2.00000*ALGEBRAIC[24];
ALGEBRAIC[199] =  (ALGEBRAIC[142]+ALGEBRAIC[109]+CONSTANTS[150])*STATES[20];
ALGEBRAIC[148] =  3.00000*ALGEBRAIC[35];
ALGEBRAIC[205] =  ALGEBRAIC[99]*STATES[19]+ ALGEBRAIC[148]*STATES[21]+ ALGEBRAIC[186]*STATES[15];
RATES[20] = ALGEBRAIC[205] - ALGEBRAIC[199];
ALGEBRAIC[79] = ALGEBRAIC[2];
ALGEBRAIC[197] =  CONSTANTS[48]*ALGEBRAIC[186];
ALGEBRAIC[208] =  (ALGEBRAIC[129]+ALGEBRAIC[79]+ALGEBRAIC[197])*STATES[16];
ALGEBRAIC[133] =  4.00000*ALGEBRAIC[13];
ALGEBRAIC[213] =  ALGEBRAIC[69]*STATES[15]+ ALGEBRAIC[133]*STATES[17]+ CONSTANTS[156]*STATES[21];
RATES[16] = ALGEBRAIC[213] - ALGEBRAIC[208];
ALGEBRAIC[117] = ALGEBRAIC[24];
ALGEBRAIC[206] =  (ALGEBRAIC[148]+ALGEBRAIC[117]+CONSTANTS[156])*STATES[21];
ALGEBRAIC[156] =  4.00000*ALGEBRAIC[35];
ALGEBRAIC[211] =  ALGEBRAIC[109]*STATES[20]+ ALGEBRAIC[156]*STATES[22]+ ALGEBRAIC[197]*STATES[16];
RATES[21] = ALGEBRAIC[211] - ALGEBRAIC[206];
ALGEBRAIC[28] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 196337./CONSTANTS[145]+529.952/CONSTANTS[2]+( 2.78085*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[127] = ALGEBRAIC[28];
ALGEBRAIC[87] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 147814./CONSTANTS[145]+338.915/CONSTANTS[2]+( 2.13600*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[164] = ALGEBRAIC[87];
ALGEBRAIC[39] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 133690./CONSTANTS[145]+229.205/CONSTANTS[2]+( - 1.55804*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[209] = ALGEBRAIC[39];
ALGEBRAIC[97] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 121322./CONSTANTS[145]+193.265/CONSTANTS[2]+( - 1.74290*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[216] = ALGEBRAIC[97];
RATES[54] = ( ALGEBRAIC[123]*STATES[53] -  (ALGEBRAIC[202]+ALGEBRAIC[127]+ALGEBRAIC[164]+CONSTANTS[173])*STATES[54])+ ALGEBRAIC[209]*STATES[45]+ ALGEBRAIC[216]*STATES[46]+ CONSTANTS[172]*STATES[59];
RATES[46] = ( ALGEBRAIC[164]*STATES[54]+ CONSTANTS[158]*STATES[45]) -  (ALGEBRAIC[216]+CONSTANTS[159])*STATES[46];
ALGEBRAIC[131] =  ALGEBRAIC[107]*CONSTANTS[162];
ALGEBRAIC[215] = ALGEBRAIC[173]/CONSTANTS[162];
RATES[55] = ( CONSTANTS[165]*STATES[50] -  (CONSTANTS[164]+ALGEBRAIC[131])*STATES[55])+ ALGEBRAIC[215]*STATES[56];
ALGEBRAIC[204] =  CONSTANTS[48]*ALGEBRAIC[197];
ALGEBRAIC[214] =  (ALGEBRAIC[133]+CONSTANTS[45]+ALGEBRAIC[204])*STATES[17];
ALGEBRAIC[219] =  ALGEBRAIC[79]*STATES[16]+ CONSTANTS[46]*STATES[23]+ CONSTANTS[161]*STATES[22];
RATES[17] = ALGEBRAIC[219] - ALGEBRAIC[214];
ALGEBRAIC[212] =  (ALGEBRAIC[156]+CONSTANTS[161])*STATES[22];
ALGEBRAIC[218] =  ALGEBRAIC[117]*STATES[21]+ ALGEBRAIC[204]*STATES[17];
RATES[22] = ALGEBRAIC[218] - ALGEBRAIC[212];
ALGEBRAIC[135] =  ALGEBRAIC[115]*CONSTANTS[162];
ALGEBRAIC[220] = ALGEBRAIC[184]/CONSTANTS[162];
RATES[56] = (( CONSTANTS[167]*STATES[51]+ ALGEBRAIC[131]*STATES[55]) -  (ALGEBRAIC[215]+CONSTANTS[166]+ALGEBRAIC[135])*STATES[56])+ ALGEBRAIC[220]*STATES[57];
ALGEBRAIC[140] =  ALGEBRAIC[119]*CONSTANTS[162];
ALGEBRAIC[222] = ALGEBRAIC[195]/CONSTANTS[162];
RATES[57] = (( CONSTANTS[169]*STATES[52]+ ALGEBRAIC[135]*STATES[56]) -  (ALGEBRAIC[140]+CONSTANTS[168]+ALGEBRAIC[220])*STATES[57])+ ALGEBRAIC[222]*STATES[58];
ALGEBRAIC[146] =  ALGEBRAIC[123]*CONSTANTS[162];
ALGEBRAIC[224] = ALGEBRAIC[202]/CONSTANTS[162];
RATES[58] = (( CONSTANTS[171]*STATES[53]+ ALGEBRAIC[140]*STATES[57]) -  (ALGEBRAIC[222]+CONSTANTS[170]+ALGEBRAIC[146])*STATES[58])+ ALGEBRAIC[224]*STATES[59];
ALGEBRAIC[64] =  CONSTANTS[149]*CONSTANTS[125]*exp(116431./CONSTANTS[145]+- 578.317/CONSTANTS[2]+( 0.764126*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[154] = ALGEBRAIC[64];
ALGEBRAIC[75] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 55700.7/CONSTANTS[145]+- 130.639/CONSTANTS[2]+( - 3.64982*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[226] = ALGEBRAIC[75];
RATES[59] = (( CONSTANTS[173]*STATES[54]+ ALGEBRAIC[146]*STATES[58]) -  (ALGEBRAIC[224]+CONSTANTS[172]+ALGEBRAIC[154])*STATES[59])+ ALGEBRAIC[226]*STATES[60];
ALGEBRAIC[55] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 97657.8/CONSTANTS[145]+1.51000/CONSTANTS[2]+( 0.0684862*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[228] = ALGEBRAIC[55];
ALGEBRAIC[47] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 62123.1/CONSTANTS[145]+39.2950/CONSTANTS[2]+( 0.288816*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[230] = ALGEBRAIC[47];
RATES[45] = ( ALGEBRAIC[127]*STATES[54] -  (ALGEBRAIC[209]+CONSTANTS[158]+ALGEBRAIC[230])*STATES[45])+ CONSTANTS[159]*STATES[46]+ ALGEBRAIC[228]*STATES[60];
RATES[60] = ( ALGEBRAIC[230]*STATES[45]+ ALGEBRAIC[154]*STATES[59]) -  (ALGEBRAIC[226]+ALGEBRAIC[228])*STATES[60];
ALGEBRAIC[116] =  CONSTANTS[115]*log(CONSTANTS[9]/STATES[0]);
ALGEBRAIC[128] =  CONSTANTS[82]*(STATES[45]+STATES[46])*(STATES[6] - ALGEBRAIC[116]);
ALGEBRAIC[225] =  CONSTANTS[83]*(STATES[6] - ALGEBRAIC[116]);
ALGEBRAIC[231] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[237] =  exp( CONSTANTS[90]*ALGEBRAIC[231])*pow(STATES[0], 3.00000)*CONSTANTS[10];
ALGEBRAIC[238] =  exp( (CONSTANTS[90] - 1.00000)*ALGEBRAIC[231])*pow(CONSTANTS[9], 3.00000)*STATES[2];
ALGEBRAIC[239] = 1.00000+ CONSTANTS[89]*exp( (CONSTANTS[90] - 1.00000)*ALGEBRAIC[231]);
ALGEBRAIC[240] = ( CONSTANTS[86]*(ALGEBRAIC[237] - ALGEBRAIC[238]))/( CONSTANTS[160]*ALGEBRAIC[239]*CONSTANTS[163]);
ALGEBRAIC[232] = 1.00000+ 0.124500*exp( - 0.100000*ALGEBRAIC[231]);
ALGEBRAIC[233] =  0.0365000*CONSTANTS[151]*exp( - 1.33000*ALGEBRAIC[231]);
ALGEBRAIC[234] = 1.00000/(ALGEBRAIC[232]+ALGEBRAIC[233]);
ALGEBRAIC[235] = 1.00000+pow(CONSTANTS[92]/STATES[0], 1.50000);
ALGEBRAIC[236] = ( CONSTANTS[91]*ALGEBRAIC[234]*CONSTANTS[155])/ALGEBRAIC[235];
ALGEBRAIC[147] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[155] =  1000.00*CONSTANTS[0]*ALGEBRAIC[147];
ALGEBRAIC[174] = exp(ALGEBRAIC[147]) - 1.00000;
ALGEBRAIC[196] =  STATES[0]*exp(ALGEBRAIC[147]) - CONSTANTS[9];
ALGEBRAIC[203] = ( 0.0200000*CONSTANTS[121]*STATES[39]*ALGEBRAIC[155]*ALGEBRAIC[196])/ALGEBRAIC[174];
RATES[0] =  - (ALGEBRAIC[128]+ALGEBRAIC[225]+ 3.00000*(ALGEBRAIC[240]+ALGEBRAIC[236])+ALGEBRAIC[203])*CONSTANTS[143];
ALGEBRAIC[124] =  0.500000*CONSTANTS[115]*log(CONSTANTS[10]/STATES[2]);
ALGEBRAIC[227] =  CONSTANTS[96]*(STATES[6] - ALGEBRAIC[124]);
ALGEBRAIC[229] = ( CONSTANTS[94]*STATES[2])/(CONSTANTS[95]+STATES[2]);
ALGEBRAIC[58] = (STATES[3] - STATES[2])/CONSTANTS[24];
ALGEBRAIC[88] =  CONSTANTS[25]*ALGEBRAIC[67]+ CONSTANTS[26]*ALGEBRAIC[78];
ALGEBRAIC[108] = 1.00000/(1.00000+( CONSTANTS[31]*CONSTANTS[34])/pow(STATES[2]+CONSTANTS[34], 2.00000)+( CONSTANTS[33]*CONSTANTS[36])/pow(STATES[2]+CONSTANTS[36], 2.00000));
RATES[2] =  ALGEBRAIC[108]*(((ALGEBRAIC[58] - ALGEBRAIC[42]) - ALGEBRAIC[88]) -  ((ALGEBRAIC[227] -  2.00000*ALGEBRAIC[240])+ALGEBRAIC[229])*0.500000*CONSTANTS[143]);
ALGEBRAIC[241] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[242] =  1000.00*CONSTANTS[0]*ALGEBRAIC[241];
ALGEBRAIC[243] =  0.00100000*exp( 2.00000*ALGEBRAIC[241]) -  CONSTANTS[10]*0.341000;
ALGEBRAIC[244] = exp( 2.00000*ALGEBRAIC[241]) - 1.00000;
ALGEBRAIC[245] = ( CONSTANTS[144]*4.00000*ALGEBRAIC[242]*ALGEBRAIC[243])/ALGEBRAIC[244];
ALGEBRAIC[246] =  ALGEBRAIC[245]*STATES[24]*STATES[23];
ALGEBRAIC[98] = 1.00000/(1.00000+( CONSTANTS[31]*CONSTANTS[34])/pow(STATES[3]+CONSTANTS[34], 2.00000)+( CONSTANTS[33]*CONSTANTS[36])/pow(STATES[3]+CONSTANTS[36], 2.00000));
RATES[3] =  ALGEBRAIC[98]*((( ALGEBRAIC[50]*CONSTANTS[5])/CONSTANTS[7] - ( ALGEBRAIC[58]*CONSTANTS[4])/CONSTANTS[7]) -  ALGEBRAIC[246]*CONSTANTS[148]);
ALGEBRAIC[120] =  CONSTANTS[115]*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[132] =  CONSTANTS[79]*CONSTANTS[122]*STATES[47]*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[136] =  CONSTANTS[80]*(STATES[48]+STATES[49])*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[221] = 1.00000/(0.940000+exp( (1.26000/CONSTANTS[115])*(STATES[6] - ALGEBRAIC[120])));
ALGEBRAIC[223] =  CONSTANTS[81]* pow((CONSTANTS[8]/1.00000), 1.0 / 2)*ALGEBRAIC[221]*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[247] = (ALGEBRAIC[245]>=0.00000 ? 0.00000 : ALGEBRAIC[245]);
ALGEBRAIC[248] = CONSTANTS[124]/(1.00000+ALGEBRAIC[247]/CONSTANTS[97]);
ALGEBRAIC[249] =  STATES[1]*exp(ALGEBRAIC[241]) - CONSTANTS[8];
ALGEBRAIC[250] = exp(ALGEBRAIC[241]) - 1.00000;
ALGEBRAIC[251] = ( ALGEBRAIC[248]*STATES[23]*STATES[24]*ALGEBRAIC[242]*ALGEBRAIC[249])/ALGEBRAIC[250];
ALGEBRAIC[0] =  floor(VOI/CONSTANTS[11])*CONSTANTS[11];
ALGEBRAIC[11] = (VOI - ALGEBRAIC[0]>=CONSTANTS[14]&&VOI - ALGEBRAIC[0]<=CONSTANTS[14]+CONSTANTS[12] ? CONSTANTS[13] : 0.00000);
ALGEBRAIC[141] =  CONSTANTS[120]*STATES[29]*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[165] =  STATES[1]*exp(ALGEBRAIC[147]) - CONSTANTS[8];
ALGEBRAIC[185] = ( CONSTANTS[121]*STATES[39]*ALGEBRAIC[155]*ALGEBRAIC[165])/ALGEBRAIC[174];
RATES[1] =  - (((ALGEBRAIC[132]+ALGEBRAIC[136]+ALGEBRAIC[223]+ALGEBRAIC[251]+ALGEBRAIC[11]) -  2.00000*ALGEBRAIC[236])+ALGEBRAIC[141]+ALGEBRAIC[185])*CONSTANTS[143];
ALGEBRAIC[210] = ALGEBRAIC[185]+ALGEBRAIC[203];
ALGEBRAIC[217] = ALGEBRAIC[141]+ALGEBRAIC[210];
ALGEBRAIC[252] = ALGEBRAIC[128]+ALGEBRAIC[246]+ALGEBRAIC[251]+ALGEBRAIC[132]+ALGEBRAIC[136]+ALGEBRAIC[223]+ALGEBRAIC[240]+ALGEBRAIC[236]+ALGEBRAIC[217]+ALGEBRAIC[229]+ALGEBRAIC[227]+ALGEBRAIC[225]+ALGEBRAIC[11];
RATES[6] = - ALGEBRAIC[252];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] =  - CONSTANTS[37]*pow( STATES[3]*1000.00, CONSTANTS[43])*STATES[11]+ CONSTANTS[38]*STATES[7];
ALGEBRAIC[12] =  CONSTANTS[39]*pow( STATES[3]*1000.00, CONSTANTS[44])*STATES[7] -  CONSTANTS[40]*STATES[8];
ALGEBRAIC[3] = (CONSTANTS[117]/(1.00000+exp((STATES[6]+28.5000)/7.80000))+1.00000) - CONSTANTS[117];
ALGEBRAIC[14] = 1.00000/(0.00653000/(0.500000+exp(STATES[6]/- 7.10000))+ 0.00512000*exp(STATES[6]/- 39.8000));
ALGEBRAIC[8] =  CONSTANTS[100]*CONSTANTS[101]*exp( CONSTANTS[102]*STATES[6]);
ALGEBRAIC[19] =  CONSTANTS[100]*CONSTANTS[103]*exp( CONSTANTS[104]*STATES[6]);
ALGEBRAIC[23] =  CONSTANTS[41]*STATES[7] -  CONSTANTS[42]*STATES[12];
ALGEBRAIC[32] =  0.216256*exp( - 1.88912e-05*STATES[6]);
ALGEBRAIC[21] =  0.00700807*exp( - 0.149998*STATES[6]);
ALGEBRAIC[10] =  0.00767254*exp( 0.0866295*STATES[6]);
ALGEBRAIC[31] =  0.00379738*exp( - 0.0142567*STATES[6]);
ALGEBRAIC[34] = - (ALGEBRAIC[1]+ALGEBRAIC[12]+ALGEBRAIC[23]);
ALGEBRAIC[30] =  ALGEBRAIC[8]*STATES[61]+ CONSTANTS[128]*STATES[63];
ALGEBRAIC[41] =  (ALGEBRAIC[19]+CONSTANTS[127])*STATES[62];
ALGEBRAIC[9] =  CONSTANTS[100]*CONSTANTS[105]*exp( CONSTANTS[106]*STATES[6]);
ALGEBRAIC[40] =  CONSTANTS[100]*CONSTANTS[111]*exp( CONSTANTS[112]*STATES[6]);
ALGEBRAIC[49] =  ALGEBRAIC[9]*STATES[63]+ ALGEBRAIC[40]*STATES[64];
ALGEBRAIC[20] =  CONSTANTS[100]*CONSTANTS[107]*exp( CONSTANTS[108]*STATES[6]);
ALGEBRAIC[29] =  CONSTANTS[100]*CONSTANTS[109]*exp( CONSTANTS[110]*STATES[6]);
ALGEBRAIC[57] =  (ALGEBRAIC[20]+ALGEBRAIC[29])*STATES[47];
ALGEBRAIC[22] = pow(STATES[2]/CONSTANTS[15], CONSTANTS[16]);
ALGEBRAIC[33] = pow(STATES[5]/CONSTANTS[17], CONSTANTS[18]);
ALGEBRAIC[42] = ( CONSTANTS[19]*( CONSTANTS[20]*ALGEBRAIC[22] -  CONSTANTS[21]*ALGEBRAIC[33]))/(1.00000+ALGEBRAIC[22]+ALGEBRAIC[33]);
ALGEBRAIC[59] = (STATES[5] - STATES[4])/CONSTANTS[23];
ALGEBRAIC[50] =  CONSTANTS[22]*(STATES[7]+STATES[8])*(STATES[4] - STATES[3]);
ALGEBRAIC[68] = 1.00000/(1.00000+( CONSTANTS[32]*CONSTANTS[35])/pow(STATES[4]+CONSTANTS[35], 2.00000));
ALGEBRAIC[67] =  CONSTANTS[29]*STATES[2]*(1.00000 - STATES[9]) -  CONSTANTS[30]*STATES[9];
ALGEBRAIC[48] =  CONSTANTS[100]*CONSTANTS[113]*exp( CONSTANTS[114]*STATES[6]);
ALGEBRAIC[56] = ( ALGEBRAIC[20]*ALGEBRAIC[40]*ALGEBRAIC[48])/( ALGEBRAIC[9]*ALGEBRAIC[29]);
ALGEBRAIC[65] =  CONSTANTS[127]*STATES[62]+ ALGEBRAIC[20]*STATES[47]+ ALGEBRAIC[56]*STATES[64];
ALGEBRAIC[76] =  (ALGEBRAIC[48]+ALGEBRAIC[9]+CONSTANTS[128])*STATES[63];
ALGEBRAIC[66] =  ALGEBRAIC[48]*STATES[63]+ ALGEBRAIC[29]*STATES[47];
ALGEBRAIC[77] =  (ALGEBRAIC[56]+ALGEBRAIC[40])*STATES[64];
ALGEBRAIC[78] =  CONSTANTS[27]*STATES[2]*(1.00000 - STATES[10]) -  CONSTANTS[28]*STATES[10];
ALGEBRAIC[5] =  CONSTANTS[73]*exp( CONSTANTS[74]*STATES[6]);
ALGEBRAIC[26] =  4.00000*ALGEBRAIC[5];
ALGEBRAIC[73] =  (ALGEBRAIC[26]+CONSTANTS[132])*STATES[35];
ALGEBRAIC[16] =  CONSTANTS[75]*exp( - CONSTANTS[76]*STATES[6]);
ALGEBRAIC[62] = ALGEBRAIC[16];
ALGEBRAIC[83] =  ALGEBRAIC[62]*STATES[36]+ CONSTANTS[137]*STATES[40];
ALGEBRAIC[27] =  4.00000*CONSTANTS[69]*ALGEBRAIC[5];
ALGEBRAIC[74] =  (CONSTANTS[137]+ALGEBRAIC[27])*STATES[40];
ALGEBRAIC[63] = ALGEBRAIC[16]/CONSTANTS[65];
ALGEBRAIC[85] =  CONSTANTS[132]*STATES[35]+ ALGEBRAIC[63]*STATES[41];
ALGEBRAIC[37] =  3.00000*ALGEBRAIC[5];
ALGEBRAIC[84] =  (ALGEBRAIC[37]+ALGEBRAIC[62]+CONSTANTS[133])*STATES[36];
ALGEBRAIC[71] =  2.00000*ALGEBRAIC[16];
ALGEBRAIC[93] =  ALGEBRAIC[71]*STATES[37]+ CONSTANTS[138]*STATES[41]+ ALGEBRAIC[26]*STATES[35];
ALGEBRAIC[38] = ( 3.00000*CONSTANTS[70]*ALGEBRAIC[5])/CONSTANTS[69];
ALGEBRAIC[86] =  (ALGEBRAIC[38]+CONSTANTS[138]+ALGEBRAIC[63])*STATES[41];
ALGEBRAIC[72] = ( 2.00000*CONSTANTS[65]*ALGEBRAIC[16])/CONSTANTS[66];
ALGEBRAIC[95] =  ALGEBRAIC[72]*STATES[42]+ CONSTANTS[133]*STATES[36]+ ALGEBRAIC[27]*STATES[40];
ALGEBRAIC[45] =  2.00000*ALGEBRAIC[5];
ALGEBRAIC[94] =  (ALGEBRAIC[45]+ALGEBRAIC[71]+CONSTANTS[134])*STATES[37];
ALGEBRAIC[81] =  3.00000*ALGEBRAIC[16];
ALGEBRAIC[101] =  ALGEBRAIC[81]*STATES[38]+ CONSTANTS[139]*STATES[42]+ ALGEBRAIC[37]*STATES[36];
ALGEBRAIC[46] = ( 2.00000*CONSTANTS[71]*ALGEBRAIC[5])/CONSTANTS[70];
ALGEBRAIC[96] =  (ALGEBRAIC[46]+CONSTANTS[139]+ALGEBRAIC[72])*STATES[42];
ALGEBRAIC[82] = ( 3.00000*CONSTANTS[66]*ALGEBRAIC[16])/CONSTANTS[67];
ALGEBRAIC[104] =  ALGEBRAIC[82]*STATES[43]+ CONSTANTS[134]*STATES[37]+ ALGEBRAIC[38]*STATES[41];
ALGEBRAIC[53] = ALGEBRAIC[5];
ALGEBRAIC[102] =  (ALGEBRAIC[53]+ALGEBRAIC[81]+CONSTANTS[135])*STATES[38];
ALGEBRAIC[91] =  4.00000*ALGEBRAIC[16];
ALGEBRAIC[111] =  ALGEBRAIC[91]*STATES[39]+ CONSTANTS[140]*STATES[43]+ ALGEBRAIC[45]*STATES[37];
ALGEBRAIC[103] =  (ALGEBRAIC[91]+CONSTANTS[136])*STATES[39];
ALGEBRAIC[112] =  ALGEBRAIC[53]*STATES[38]+ CONSTANTS[141]*STATES[44];
ALGEBRAIC[54] = ( CONSTANTS[72]*ALGEBRAIC[5])/CONSTANTS[71];
ALGEBRAIC[105] =  (ALGEBRAIC[54]+CONSTANTS[140]+ALGEBRAIC[82])*STATES[43];
ALGEBRAIC[92] = ( 4.00000*CONSTANTS[67]*ALGEBRAIC[16])/CONSTANTS[68];
ALGEBRAIC[113] =  ALGEBRAIC[92]*STATES[44]+ CONSTANTS[135]*STATES[38]+ ALGEBRAIC[46]*STATES[42];
ALGEBRAIC[106] =  (CONSTANTS[141]+ALGEBRAIC[92])*STATES[44];
ALGEBRAIC[114] =  CONSTANTS[136]*STATES[39]+ ALGEBRAIC[54]*STATES[43];
ALGEBRAIC[4] =  CONSTANTS[49]*exp( CONSTANTS[50]*STATES[6]);
ALGEBRAIC[44] =  4.00000*ALGEBRAIC[4];
ALGEBRAIC[36] =  CONSTANTS[55]*exp( CONSTANTS[56]*STATES[6]);
ALGEBRAIC[138] = ALGEBRAIC[36];
ALGEBRAIC[152] =  (ALGEBRAIC[44]+ALGEBRAIC[138])*STATES[25];
ALGEBRAIC[15] =  CONSTANTS[51]*exp( - CONSTANTS[52]*STATES[6]);
ALGEBRAIC[118] = ALGEBRAIC[15];
ALGEBRAIC[25] =  CONSTANTS[53]*exp( - CONSTANTS[54]*STATES[6]);
ALGEBRAIC[145] = ALGEBRAIC[25];
ALGEBRAIC[160] =  ALGEBRAIC[118]*STATES[26]+ ALGEBRAIC[145]*STATES[30];
ALGEBRAIC[80] =  4.00000*CONSTANTS[61]*ALGEBRAIC[4];
ALGEBRAIC[153] =  (ALGEBRAIC[145]+ALGEBRAIC[80])*STATES[30];
ALGEBRAIC[134] = ALGEBRAIC[15]/CONSTANTS[57];
ALGEBRAIC[162] =  ALGEBRAIC[138]*STATES[25]+ ALGEBRAIC[134]*STATES[31];
ALGEBRAIC[52] =  3.00000*ALGEBRAIC[4];
ALGEBRAIC[143] =  CONSTANTS[57]*ALGEBRAIC[36];
ALGEBRAIC[161] =  (ALGEBRAIC[52]+ALGEBRAIC[118]+ALGEBRAIC[143])*STATES[26];
ALGEBRAIC[122] =  2.00000*ALGEBRAIC[15];
ALGEBRAIC[151] = ALGEBRAIC[25]/CONSTANTS[61];
ALGEBRAIC[169] =  ALGEBRAIC[122]*STATES[27]+ ALGEBRAIC[151]*STATES[31]+ ALGEBRAIC[44]*STATES[25];
ALGEBRAIC[90] = ( 3.00000*CONSTANTS[62]*ALGEBRAIC[4])/CONSTANTS[61];
ALGEBRAIC[163] =  (ALGEBRAIC[90]+ALGEBRAIC[151]+ALGEBRAIC[134])*STATES[31];
ALGEBRAIC[139] = ( 2.00000*CONSTANTS[57]*ALGEBRAIC[15])/CONSTANTS[58];
ALGEBRAIC[171] =  ALGEBRAIC[139]*STATES[32]+ ALGEBRAIC[143]*STATES[26]+ ALGEBRAIC[80]*STATES[30];
ALGEBRAIC[7] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 114007./CONSTANTS[145]+224.114/CONSTANTS[2]+( 0.286374*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[107] =  4.00000*ALGEBRAIC[7];
ALGEBRAIC[18] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 272470./CONSTANTS[145]+708.146/CONSTANTS[2]+( - 2.28528*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[173] = ALGEBRAIC[18];
ALGEBRAIC[61] =  2.00000*ALGEBRAIC[4];
ALGEBRAIC[149] =  CONSTANTS[58]*ALGEBRAIC[36];
ALGEBRAIC[170] =  (ALGEBRAIC[61]+ALGEBRAIC[122]+ALGEBRAIC[149])*STATES[27];
ALGEBRAIC[126] =  3.00000*ALGEBRAIC[15];
ALGEBRAIC[159] = ALGEBRAIC[25]/CONSTANTS[62];
ALGEBRAIC[178] =  ALGEBRAIC[126]*STATES[28]+ ALGEBRAIC[159]*STATES[32]+ ALGEBRAIC[52]*STATES[26];
ALGEBRAIC[100] = ( 2.00000*CONSTANTS[63]*ALGEBRAIC[4])/CONSTANTS[62];
ALGEBRAIC[172] =  (ALGEBRAIC[100]+ALGEBRAIC[159]+ALGEBRAIC[139])*STATES[32];
ALGEBRAIC[144] = ( 3.00000*CONSTANTS[58]*ALGEBRAIC[15])/CONSTANTS[59];
ALGEBRAIC[181] =  ALGEBRAIC[144]*STATES[33]+ ALGEBRAIC[149]*STATES[27]+ ALGEBRAIC[90]*STATES[31];
ALGEBRAIC[115] =  3.00000*ALGEBRAIC[7];
ALGEBRAIC[184] =  2.00000*ALGEBRAIC[18];
ALGEBRAIC[2] =  4.00000*1.20000*0.416000*exp( 0.0120000*(STATES[6] - 35.0000));
ALGEBRAIC[51] =  4.00000*ALGEBRAIC[2];
ALGEBRAIC[43] =  0.600000*0.0923300*STATES[3];
ALGEBRAIC[166] = ALGEBRAIC[43];
ALGEBRAIC[177] =  (ALGEBRAIC[51]+ALGEBRAIC[166])*STATES[13];
ALGEBRAIC[13] =  4.00000*0.450000*0.0490000*exp( - 0.0650000*(STATES[6] - 22.0000));
ALGEBRAIC[121] = ALGEBRAIC[13];
ALGEBRAIC[189] =  ALGEBRAIC[121]*STATES[14]+ CONSTANTS[131]*STATES[18];
ALGEBRAIC[24] =  CONSTANTS[48]*ALGEBRAIC[2];
ALGEBRAIC[89] =  4.00000*ALGEBRAIC[24];
ALGEBRAIC[176] =  (ALGEBRAIC[89]+CONSTANTS[131])*STATES[18];
ALGEBRAIC[35] = ALGEBRAIC[13]/CONSTANTS[47];
ALGEBRAIC[137] = ALGEBRAIC[35];
ALGEBRAIC[187] =  ALGEBRAIC[137]*STATES[19]+ ALGEBRAIC[166]*STATES[13];
ALGEBRAIC[70] = ALGEBRAIC[4];
ALGEBRAIC[157] =  CONSTANTS[59]*ALGEBRAIC[36];
ALGEBRAIC[179] =  (ALGEBRAIC[70]+ALGEBRAIC[126]+ALGEBRAIC[157])*STATES[28];
ALGEBRAIC[130] =  4.00000*ALGEBRAIC[15];
ALGEBRAIC[167] = ALGEBRAIC[25]/CONSTANTS[63];
ALGEBRAIC[191] =  ALGEBRAIC[130]*STATES[29]+ ALGEBRAIC[167]*STATES[33]+ ALGEBRAIC[61]*STATES[27];
ALGEBRAIC[158] =  CONSTANTS[60]*ALGEBRAIC[36];
ALGEBRAIC[180] =  (ALGEBRAIC[130]+ALGEBRAIC[158])*STATES[29];
ALGEBRAIC[168] = ALGEBRAIC[25]/CONSTANTS[64];
ALGEBRAIC[192] =  ALGEBRAIC[70]*STATES[28]+ ALGEBRAIC[168]*STATES[34];
ALGEBRAIC[110] = ( CONSTANTS[64]*ALGEBRAIC[4])/CONSTANTS[63];
ALGEBRAIC[182] =  (ALGEBRAIC[110]+ALGEBRAIC[167]+ALGEBRAIC[144])*STATES[33];
ALGEBRAIC[150] = ( 4.00000*CONSTANTS[59]*ALGEBRAIC[15])/CONSTANTS[60];
ALGEBRAIC[193] =  ALGEBRAIC[150]*STATES[34]+ ALGEBRAIC[157]*STATES[28]+ ALGEBRAIC[100]*STATES[32];
ALGEBRAIC[183] =  (ALGEBRAIC[168]+ALGEBRAIC[150])*STATES[34];
ALGEBRAIC[194] =  ALGEBRAIC[158]*STATES[29]+ ALGEBRAIC[110]*STATES[33];
ALGEBRAIC[119] =  2.00000*ALGEBRAIC[7];
ALGEBRAIC[195] =  3.00000*ALGEBRAIC[18];
ALGEBRAIC[60] =  3.00000*ALGEBRAIC[2];
ALGEBRAIC[175] =  CONSTANTS[48]*ALGEBRAIC[166];
ALGEBRAIC[190] =  (ALGEBRAIC[121]+ALGEBRAIC[60]+ALGEBRAIC[175])*STATES[14];
ALGEBRAIC[125] =  2.00000*ALGEBRAIC[13];
ALGEBRAIC[200] =  ALGEBRAIC[51]*STATES[13]+ ALGEBRAIC[125]*STATES[15]+ CONSTANTS[147]*STATES[19];
ALGEBRAIC[99] =  3.00000*ALGEBRAIC[24];
ALGEBRAIC[188] =  (ALGEBRAIC[137]+ALGEBRAIC[99]+CONSTANTS[147])*STATES[19];
ALGEBRAIC[142] =  2.00000*ALGEBRAIC[35];
ALGEBRAIC[198] =  ALGEBRAIC[89]*STATES[18]+ ALGEBRAIC[142]*STATES[20]+ ALGEBRAIC[175]*STATES[14];
ALGEBRAIC[123] = ALGEBRAIC[7];
ALGEBRAIC[202] =  4.00000*ALGEBRAIC[18];
ALGEBRAIC[69] =  2.00000*ALGEBRAIC[2];
ALGEBRAIC[186] =  CONSTANTS[48]*ALGEBRAIC[175];
ALGEBRAIC[201] =  (ALGEBRAIC[125]+ALGEBRAIC[69]+ALGEBRAIC[186])*STATES[15];
ALGEBRAIC[129] =  3.00000*ALGEBRAIC[13];
ALGEBRAIC[207] =  ALGEBRAIC[60]*STATES[14]+ ALGEBRAIC[129]*STATES[16]+ CONSTANTS[150]*STATES[20];
ALGEBRAIC[109] =  2.00000*ALGEBRAIC[24];
ALGEBRAIC[199] =  (ALGEBRAIC[142]+ALGEBRAIC[109]+CONSTANTS[150])*STATES[20];
ALGEBRAIC[148] =  3.00000*ALGEBRAIC[35];
ALGEBRAIC[205] =  ALGEBRAIC[99]*STATES[19]+ ALGEBRAIC[148]*STATES[21]+ ALGEBRAIC[186]*STATES[15];
ALGEBRAIC[79] = ALGEBRAIC[2];
ALGEBRAIC[197] =  CONSTANTS[48]*ALGEBRAIC[186];
ALGEBRAIC[208] =  (ALGEBRAIC[129]+ALGEBRAIC[79]+ALGEBRAIC[197])*STATES[16];
ALGEBRAIC[133] =  4.00000*ALGEBRAIC[13];
ALGEBRAIC[213] =  ALGEBRAIC[69]*STATES[15]+ ALGEBRAIC[133]*STATES[17]+ CONSTANTS[156]*STATES[21];
ALGEBRAIC[117] = ALGEBRAIC[24];
ALGEBRAIC[206] =  (ALGEBRAIC[148]+ALGEBRAIC[117]+CONSTANTS[156])*STATES[21];
ALGEBRAIC[156] =  4.00000*ALGEBRAIC[35];
ALGEBRAIC[211] =  ALGEBRAIC[109]*STATES[20]+ ALGEBRAIC[156]*STATES[22]+ ALGEBRAIC[197]*STATES[16];
ALGEBRAIC[28] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 196337./CONSTANTS[145]+529.952/CONSTANTS[2]+( 2.78085*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[127] = ALGEBRAIC[28];
ALGEBRAIC[87] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 147814./CONSTANTS[145]+338.915/CONSTANTS[2]+( 2.13600*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[164] = ALGEBRAIC[87];
ALGEBRAIC[39] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 133690./CONSTANTS[145]+229.205/CONSTANTS[2]+( - 1.55804*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[209] = ALGEBRAIC[39];
ALGEBRAIC[97] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 121322./CONSTANTS[145]+193.265/CONSTANTS[2]+( - 1.74290*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[216] = ALGEBRAIC[97];
ALGEBRAIC[131] =  ALGEBRAIC[107]*CONSTANTS[162];
ALGEBRAIC[215] = ALGEBRAIC[173]/CONSTANTS[162];
ALGEBRAIC[204] =  CONSTANTS[48]*ALGEBRAIC[197];
ALGEBRAIC[214] =  (ALGEBRAIC[133]+CONSTANTS[45]+ALGEBRAIC[204])*STATES[17];
ALGEBRAIC[219] =  ALGEBRAIC[79]*STATES[16]+ CONSTANTS[46]*STATES[23]+ CONSTANTS[161]*STATES[22];
ALGEBRAIC[212] =  (ALGEBRAIC[156]+CONSTANTS[161])*STATES[22];
ALGEBRAIC[218] =  ALGEBRAIC[117]*STATES[21]+ ALGEBRAIC[204]*STATES[17];
ALGEBRAIC[135] =  ALGEBRAIC[115]*CONSTANTS[162];
ALGEBRAIC[220] = ALGEBRAIC[184]/CONSTANTS[162];
ALGEBRAIC[140] =  ALGEBRAIC[119]*CONSTANTS[162];
ALGEBRAIC[222] = ALGEBRAIC[195]/CONSTANTS[162];
ALGEBRAIC[146] =  ALGEBRAIC[123]*CONSTANTS[162];
ALGEBRAIC[224] = ALGEBRAIC[202]/CONSTANTS[162];
ALGEBRAIC[64] =  CONSTANTS[149]*CONSTANTS[125]*exp(116431./CONSTANTS[145]+- 578.317/CONSTANTS[2]+( 0.764126*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[154] = ALGEBRAIC[64];
ALGEBRAIC[75] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 55700.7/CONSTANTS[145]+- 130.639/CONSTANTS[2]+( - 3.64982*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[226] = ALGEBRAIC[75];
ALGEBRAIC[55] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 97657.8/CONSTANTS[145]+1.51000/CONSTANTS[2]+( 0.0684862*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[228] = ALGEBRAIC[55];
ALGEBRAIC[47] =  CONSTANTS[149]*CONSTANTS[125]*exp(- 62123.1/CONSTANTS[145]+39.2950/CONSTANTS[2]+( 0.288816*STATES[6])/CONSTANTS[126]);
ALGEBRAIC[230] = ALGEBRAIC[47];
ALGEBRAIC[116] =  CONSTANTS[115]*log(CONSTANTS[9]/STATES[0]);
ALGEBRAIC[128] =  CONSTANTS[82]*(STATES[45]+STATES[46])*(STATES[6] - ALGEBRAIC[116]);
ALGEBRAIC[225] =  CONSTANTS[83]*(STATES[6] - ALGEBRAIC[116]);
ALGEBRAIC[231] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[237] =  exp( CONSTANTS[90]*ALGEBRAIC[231])*pow(STATES[0], 3.00000)*CONSTANTS[10];
ALGEBRAIC[238] =  exp( (CONSTANTS[90] - 1.00000)*ALGEBRAIC[231])*pow(CONSTANTS[9], 3.00000)*STATES[2];
ALGEBRAIC[239] = 1.00000+ CONSTANTS[89]*exp( (CONSTANTS[90] - 1.00000)*ALGEBRAIC[231]);
ALGEBRAIC[240] = ( CONSTANTS[86]*(ALGEBRAIC[237] - ALGEBRAIC[238]))/( CONSTANTS[160]*ALGEBRAIC[239]*CONSTANTS[163]);
ALGEBRAIC[232] = 1.00000+ 0.124500*exp( - 0.100000*ALGEBRAIC[231]);
ALGEBRAIC[233] =  0.0365000*CONSTANTS[151]*exp( - 1.33000*ALGEBRAIC[231]);
ALGEBRAIC[234] = 1.00000/(ALGEBRAIC[232]+ALGEBRAIC[233]);
ALGEBRAIC[235] = 1.00000+pow(CONSTANTS[92]/STATES[0], 1.50000);
ALGEBRAIC[236] = ( CONSTANTS[91]*ALGEBRAIC[234]*CONSTANTS[155])/ALGEBRAIC[235];
ALGEBRAIC[147] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[155] =  1000.00*CONSTANTS[0]*ALGEBRAIC[147];
ALGEBRAIC[174] = exp(ALGEBRAIC[147]) - 1.00000;
ALGEBRAIC[196] =  STATES[0]*exp(ALGEBRAIC[147]) - CONSTANTS[9];
ALGEBRAIC[203] = ( 0.0200000*CONSTANTS[121]*STATES[39]*ALGEBRAIC[155]*ALGEBRAIC[196])/ALGEBRAIC[174];
ALGEBRAIC[124] =  0.500000*CONSTANTS[115]*log(CONSTANTS[10]/STATES[2]);
ALGEBRAIC[227] =  CONSTANTS[96]*(STATES[6] - ALGEBRAIC[124]);
ALGEBRAIC[229] = ( CONSTANTS[94]*STATES[2])/(CONSTANTS[95]+STATES[2]);
ALGEBRAIC[58] = (STATES[3] - STATES[2])/CONSTANTS[24];
ALGEBRAIC[88] =  CONSTANTS[25]*ALGEBRAIC[67]+ CONSTANTS[26]*ALGEBRAIC[78];
ALGEBRAIC[108] = 1.00000/(1.00000+( CONSTANTS[31]*CONSTANTS[34])/pow(STATES[2]+CONSTANTS[34], 2.00000)+( CONSTANTS[33]*CONSTANTS[36])/pow(STATES[2]+CONSTANTS[36], 2.00000));
ALGEBRAIC[241] = STATES[6]/CONSTANTS[115];
ALGEBRAIC[242] =  1000.00*CONSTANTS[0]*ALGEBRAIC[241];
ALGEBRAIC[243] =  0.00100000*exp( 2.00000*ALGEBRAIC[241]) -  CONSTANTS[10]*0.341000;
ALGEBRAIC[244] = exp( 2.00000*ALGEBRAIC[241]) - 1.00000;
ALGEBRAIC[245] = ( CONSTANTS[144]*4.00000*ALGEBRAIC[242]*ALGEBRAIC[243])/ALGEBRAIC[244];
ALGEBRAIC[246] =  ALGEBRAIC[245]*STATES[24]*STATES[23];
ALGEBRAIC[98] = 1.00000/(1.00000+( CONSTANTS[31]*CONSTANTS[34])/pow(STATES[3]+CONSTANTS[34], 2.00000)+( CONSTANTS[33]*CONSTANTS[36])/pow(STATES[3]+CONSTANTS[36], 2.00000));
ALGEBRAIC[120] =  CONSTANTS[115]*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[132] =  CONSTANTS[79]*CONSTANTS[122]*STATES[47]*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[136] =  CONSTANTS[80]*(STATES[48]+STATES[49])*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[221] = 1.00000/(0.940000+exp( (1.26000/CONSTANTS[115])*(STATES[6] - ALGEBRAIC[120])));
ALGEBRAIC[223] =  CONSTANTS[81]* pow((CONSTANTS[8]/1.00000), 1.0 / 2)*ALGEBRAIC[221]*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[247] = (ALGEBRAIC[245]>=0.00000 ? 0.00000 : ALGEBRAIC[245]);
ALGEBRAIC[248] = CONSTANTS[124]/(1.00000+ALGEBRAIC[247]/CONSTANTS[97]);
ALGEBRAIC[249] =  STATES[1]*exp(ALGEBRAIC[241]) - CONSTANTS[8];
ALGEBRAIC[250] = exp(ALGEBRAIC[241]) - 1.00000;
ALGEBRAIC[251] = ( ALGEBRAIC[248]*STATES[23]*STATES[24]*ALGEBRAIC[242]*ALGEBRAIC[249])/ALGEBRAIC[250];
ALGEBRAIC[0] =  floor(VOI/CONSTANTS[11])*CONSTANTS[11];
ALGEBRAIC[11] = (VOI - ALGEBRAIC[0]>=CONSTANTS[14]&&VOI - ALGEBRAIC[0]<=CONSTANTS[14]+CONSTANTS[12] ? CONSTANTS[13] : 0.00000);
ALGEBRAIC[141] =  CONSTANTS[120]*STATES[29]*(STATES[6] - ALGEBRAIC[120]);
ALGEBRAIC[165] =  STATES[1]*exp(ALGEBRAIC[147]) - CONSTANTS[8];
ALGEBRAIC[185] = ( CONSTANTS[121]*STATES[39]*ALGEBRAIC[155]*ALGEBRAIC[165])/ALGEBRAIC[174];
ALGEBRAIC[210] = ALGEBRAIC[185]+ALGEBRAIC[203];
ALGEBRAIC[217] = ALGEBRAIC[141]+ALGEBRAIC[210];
ALGEBRAIC[252] = ALGEBRAIC[128]+ALGEBRAIC[246]+ALGEBRAIC[251]+ALGEBRAIC[132]+ALGEBRAIC[136]+ALGEBRAIC[223]+ALGEBRAIC[240]+ALGEBRAIC[236]+ALGEBRAIC[217]+ALGEBRAIC[229]+ALGEBRAIC[227]+ALGEBRAIC[225]+ALGEBRAIC[11];
ALGEBRAIC[6] = STATES[1]+ 0.0183300*STATES[0];
ALGEBRAIC[17] =  CONSTANTS[115]*log(CONSTANTS[142]/ALGEBRAIC[6]);
}