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 66 entries in the algebraic variable array.
   There are a total of 35 entries in each of the rate and state variable arrays.
   There are a total of 76 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (millisecond).
 * STATES[0] is V in component membrane (millivolt).
 * CONSTANTS[0] is Cm in component membrane (microF_per_cm2).
 * CONSTANTS[1] is Vmyo in component membrane (microlitre).
 * CONSTANTS[2] is VJSR in component membrane (microlitre).
 * CONSTANTS[3] is VNSR in component membrane (microlitre).
 * CONSTANTS[4] is Vss in component membrane (microlitre).
 * CONSTANTS[5] is Acap in component membrane (cm2).
 * CONSTANTS[6] is Ko in component membrane (micromolar).
 * CONSTANTS[7] is Nao in component membrane (micromolar).
 * CONSTANTS[8] is Cao in component membrane (micromolar).
 * CONSTANTS[9] is Clo in component membrane (micromolar).
 * CONSTANTS[10] is Cli in component membrane (micromolar).
 * CONSTANTS[11] is R in component membrane (joule_per_mole_kelvin).
 * CONSTANTS[12] is T in component membrane (kelvin).
 * CONSTANTS[13] is F in component membrane (coulomb_per_millimole).
 * ALGEBRAIC[0] is i_stim in component membrane (picoA_per_picoF).
 * CONSTANTS[14] is stim_start in component membrane (millisecond).
 * CONSTANTS[15] is stim_end in component membrane (millisecond).
 * CONSTANTS[16] is stim_period in component membrane (millisecond).
 * CONSTANTS[17] is stim_duration in component membrane (millisecond).
 * CONSTANTS[18] is stim_amplitude in component membrane (picoA_per_picoF).
 * ALGEBRAIC[43] is i_CaL in component L_type_calcium_current (picoA_per_picoF).
 * ALGEBRAIC[45] is i_CaT in component T_type_calcium_current (picoA_per_picoF).
 * ALGEBRAIC[47] is i_pCa in component calcium_pump_current (picoA_per_picoF).
 * ALGEBRAIC[49] is i_NaCa in component sodium_calcium_exchange_current (picoA_per_picoF).
 * ALGEBRAIC[51] is i_Cab in component calcium_background_current (picoA_per_picoF).
 * ALGEBRAIC[53] is i_Na in component fast_sodium_current (picoA_per_picoF).
 * ALGEBRAIC[54] is i_Nab in component sodium_background_current (picoA_per_picoF).
 * ALGEBRAIC[63] is i_NaK in component sodium_potassium_pump_current (picoA_per_picoF).
 * ALGEBRAIC[56] is i_Kto_f in component fast_transient_outward_potassium_current (picoA_per_picoF).
 * ALGEBRAIC[59] is i_K1 in component time_independent_potassium_current (picoA_per_picoF).
 * ALGEBRAIC[60] is i_Ks in component slow_delayed_rectifier_potassium_current (picoA_per_picoF).
 * ALGEBRAIC[57] is i_Kur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (picoA_per_picoF).
 * ALGEBRAIC[58] is i_Kss in component non_inactivating_steady_state_potassium_current (picoA_per_picoF).
 * ALGEBRAIC[65] is i_ClCa in component calcium_activated_chloride_current (picoA_per_picoF).
 * ALGEBRAIC[61] is i_Kr in component rapid_delayed_rectifier_potassium_current (picoA_per_picoF).
 * STATES[1] is Cai in component calcium_concentration (micromolar).
 * STATES[2] is Cass in component calcium_concentration (micromolar).
 * STATES[3] is CaJSR in component calcium_concentration (micromolar).
 * STATES[4] is CaNSR in component calcium_concentration (micromolar).
 * ALGEBRAIC[15] is Bi in component calcium_concentration (dimensionless).
 * ALGEBRAIC[30] is Bss in component calcium_concentration (dimensionless).
 * ALGEBRAIC[35] is BJSR in component calcium_concentration (dimensionless).
 * CONSTANTS[19] is CMDN_tot in component calcium_concentration (micromolar).
 * CONSTANTS[20] is CSQN_tot in component calcium_concentration (micromolar).
 * CONSTANTS[21] is Km_CMDN in component calcium_concentration (micromolar).
 * CONSTANTS[22] is Km_CSQN in component calcium_concentration (micromolar).
 * ALGEBRAIC[44] is J_leak in component calcium_fluxes (micromolar_per_millisecond).
 * ALGEBRAIC[39] is J_rel in component calcium_fluxes (micromolar_per_millisecond).
 * ALGEBRAIC[46] is J_up in component calcium_fluxes (micromolar_per_millisecond).
 * ALGEBRAIC[41] is J_tr in component calcium_fluxes (micromolar_per_millisecond).
 * ALGEBRAIC[48] is J_trpn in component calcium_fluxes (micromolar_per_millisecond).
 * ALGEBRAIC[42] is J_xfer in component calcium_fluxes (micromolar_per_millisecond).
 * CONSTANTS[23] is k_plus_htrpn in component calcium_fluxes (per_micromolar_millisecond).
 * CONSTANTS[24] is k_minus_htrpn in component calcium_fluxes (per_millisecond).
 * CONSTANTS[25] is k_plus_ltrpn in component calcium_fluxes (per_micromolar_millisecond).
 * CONSTANTS[26] is k_minus_ltrpn in component calcium_fluxes (per_millisecond).
 * STATES[5] is P_RyR in component calcium_fluxes (dimensionless).
 * CONSTANTS[27] is v1 in component calcium_fluxes (per_millisecond).
 * CONSTANTS[28] is tau_tr in component calcium_fluxes (millisecond).
 * CONSTANTS[29] is v2 in component calcium_fluxes (per_millisecond).
 * CONSTANTS[30] is tau_xfer in component calcium_fluxes (millisecond).
 * CONSTANTS[31] is v3 in component calcium_fluxes (micromolar_per_millisecond).
 * CONSTANTS[32] is Km_up in component calcium_fluxes (micromolar).
 * CONSTANTS[33] is LTRPN_tot in component calcium_buffering (micromolar).
 * CONSTANTS[34] is HTRPN_tot in component calcium_buffering (micromolar).
 * STATES[6] is LTRPN_Ca in component calcium_buffering (micromolar).
 * STATES[7] is HTRPN_Ca in component calcium_buffering (micromolar).
 * CONSTANTS[35] is i_CaL_max in component L_type_calcium_current (picoA_per_picoF).
 * STATES[8] is P_O1 in component ryanodine_receptors (dimensionless).
 * STATES[9] is P_O2 in component ryanodine_receptors (dimensionless).
 * ALGEBRAIC[1] is P_C1 in component ryanodine_receptors (dimensionless).
 * STATES[10] is P_C2 in component ryanodine_receptors (dimensionless).
 * CONSTANTS[36] is k_plus_a in component ryanodine_receptors (micromolar4_per_millisecond).
 * CONSTANTS[37] is k_minus_a in component ryanodine_receptors (per_millisecond).
 * CONSTANTS[38] is k_plus_b in component ryanodine_receptors (micromolar3_per_millisecond).
 * CONSTANTS[39] is k_minus_b in component ryanodine_receptors (per_millisecond).
 * CONSTANTS[40] is k_plus_c in component ryanodine_receptors (per_millisecond).
 * CONSTANTS[41] is k_minus_c in component ryanodine_receptors (per_millisecond).
 * CONSTANTS[42] is m in component ryanodine_receptors (dimensionless).
 * CONSTANTS[43] is n in component ryanodine_receptors (dimensionless).
 * CONSTANTS[44] is E_CaL in component L_type_calcium_current (millivolt).
 * CONSTANTS[45] is g_CaL in component L_type_calcium_current (milliS_per_microF).
 * STATES[11] is O in component L_type_calcium_current (dimensionless).
 * ALGEBRAIC[2] is C1 in component L_type_calcium_current (dimensionless).
 * STATES[12] is C2 in component L_type_calcium_current (dimensionless).
 * STATES[13] is C3 in component L_type_calcium_current (dimensionless).
 * STATES[14] is C4 in component L_type_calcium_current (dimensionless).
 * STATES[15] is I1 in component L_type_calcium_current (dimensionless).
 * STATES[16] is I2 in component L_type_calcium_current (dimensionless).
 * STATES[17] is I3 in component L_type_calcium_current (dimensionless).
 * ALGEBRAIC[16] is alpha in component L_type_calcium_current (per_millisecond).
 * ALGEBRAIC[31] is beta in component L_type_calcium_current (per_millisecond).
 * ALGEBRAIC[36] is gamma in component L_type_calcium_current (per_millisecond).
 * ALGEBRAIC[40] is Kpcf in component L_type_calcium_current (per_millisecond).
 * CONSTANTS[46] is Kpcb in component L_type_calcium_current (per_millisecond).
 * CONSTANTS[47] is Kpc_max in component L_type_calcium_current (per_millisecond).
 * CONSTANTS[48] is Kpc_half in component L_type_calcium_current (micromolar).
 * CONSTANTS[49] is E_CaT in component T_type_calcium_current (millivolt).
 * CONSTANTS[50] is g_CaT in component T_type_calcium_current (milliS_per_microF).
 * STATES[18] is b in component T_type_calcium_current_b_gate (dimensionless).
 * STATES[19] is g in component T_type_calcium_current_g_gate (dimensionless).
 * ALGEBRAIC[3] is b_infinity in component T_type_calcium_current_b_gate (dimensionless).
 * ALGEBRAIC[17] is tau_b in component T_type_calcium_current_b_gate (millisecond).
 * ALGEBRAIC[4] is g_infinity in component T_type_calcium_current_g_gate (dimensionless).
 * ALGEBRAIC[18] is tau_g in component T_type_calcium_current_g_gate (millisecond).
 * CONSTANTS[51] is i_pCa_max in component calcium_pump_current (picoA_per_picoF).
 * CONSTANTS[52] is Km_pCa in component calcium_pump_current (micromolar).
 * CONSTANTS[53] is k_NaCa in component sodium_calcium_exchange_current (picoA_per_picoF).
 * CONSTANTS[54] is K_mNa in component sodium_calcium_exchange_current (micromolar).
 * CONSTANTS[55] is K_mCa in component sodium_calcium_exchange_current (micromolar).
 * CONSTANTS[56] is k_sat in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[57] is eta in component sodium_calcium_exchange_current (dimensionless).
 * STATES[20] is Nai in component sodium_concentration (micromolar).
 * CONSTANTS[58] is g_Cab in component calcium_background_current (milliS_per_microF).
 * ALGEBRAIC[50] is E_CaN in component calcium_background_current (millivolt).
 * ALGEBRAIC[52] is E_Na in component fast_sodium_current (millivolt).
 * CONSTANTS[59] is g_Na in component fast_sodium_current (milliS_per_microF).
 * STATES[21] is m in component fast_sodium_current_m_gate (dimensionless).
 * STATES[22] is h in component fast_sodium_current_h_gate (dimensionless).
 * STATES[23] is j in component fast_sodium_current_j_gate (dimensionless).
 * ALGEBRAIC[5] is alpha_m in component fast_sodium_current_m_gate (per_millisecond).
 * ALGEBRAIC[19] is beta_m in component fast_sodium_current_m_gate (per_millisecond).
 * ALGEBRAIC[6] is alpha_h in component fast_sodium_current_h_gate (per_millisecond).
 * ALGEBRAIC[20] is beta_h in component fast_sodium_current_h_gate (per_millisecond).
 * ALGEBRAIC[7] is alpha_j in component fast_sodium_current_j_gate (per_millisecond).
 * ALGEBRAIC[21] is beta_j in component fast_sodium_current_j_gate (per_millisecond).
 * CONSTANTS[60] is g_Nab in component sodium_background_current (milliS_per_microF).
 * STATES[24] is Ki in component potassium_concentration (micromolar).
 * ALGEBRAIC[55] is E_K in component fast_transient_outward_potassium_current (millivolt).
 * CONSTANTS[61] is g_Kto_f in component fast_transient_outward_potassium_current (milliS_per_microF).
 * STATES[25] is ato_f in component fast_transient_outward_potassium_current (dimensionless).
 * STATES[26] is ito_f in component fast_transient_outward_potassium_current (dimensionless).
 * ALGEBRAIC[37] is ito_f_infinity in component fast_transient_outward_potassium_current (dimensionless).
 * ALGEBRAIC[8] is tau_ito_f in component fast_transient_outward_potassium_current (millisecond).
 * ALGEBRAIC[9] is alpha_a in component fast_transient_outward_potassium_current (per_millisecond).
 * ALGEBRAIC[23] is beta_a in component fast_transient_outward_potassium_current (per_millisecond).
 * ALGEBRAIC[22] is alpha_i in component fast_transient_outward_potassium_current (per_millisecond).
 * ALGEBRAIC[32] is beta_i in component fast_transient_outward_potassium_current (per_millisecond).
 * ALGEBRAIC[10] is ass in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless).
 * CONSTANTS[62] is g_Kur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (milliS_per_microF).
 * STATES[27] is aur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless).
 * STATES[28] is iur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless).
 * ALGEBRAIC[24] is tau_aur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (millisecond).
 * ALGEBRAIC[25] is tau_iur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (millisecond).
 * ALGEBRAIC[11] is iss in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless).
 * CONSTANTS[63] is g_Kss in component non_inactivating_steady_state_potassium_current (milliS_per_microF).
 * STATES[29] is aKss in component non_inactivating_steady_state_potassium_current (dimensionless).
 * ALGEBRAIC[26] is tau_Kss in component non_inactivating_steady_state_potassium_current (millisecond).
 * CONSTANTS[64] is g_K1 in component time_independent_potassium_current (milliS_per_microF).
 * CONSTANTS[65] is g_Ks in component slow_delayed_rectifier_potassium_current (milliS_per_microF).
 * STATES[30] is nKs in component slow_delayed_rectifier_potassium_current (dimensionless).
 * ALGEBRAIC[12] is alpha_n in component slow_delayed_rectifier_potassium_current (per_millisecond).
 * ALGEBRAIC[27] is beta_n in component slow_delayed_rectifier_potassium_current (per_millisecond).
 * CONSTANTS[66] is g_Kr in component rapid_delayed_rectifier_potassium_current (milliS_per_microF).
 * STATES[31] is O_K in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * STATES[32] is C_K1 in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * STATES[33] is C_K2 in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * ALGEBRAIC[13] is C_K0 in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * STATES[34] is I_K in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * ALGEBRAIC[28] is alpha_a0 in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * ALGEBRAIC[33] is beta_a0 in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * CONSTANTS[67] is kb in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * CONSTANTS[68] is kf in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * ALGEBRAIC[14] is alpha_a1 in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * ALGEBRAIC[29] is beta_a1 in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * ALGEBRAIC[34] is alpha_i in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * ALGEBRAIC[38] is beta_i in component rapid_delayed_rectifier_potassium_current (per_millisecond).
 * CONSTANTS[69] is i_NaK_max in component sodium_potassium_pump_current (picoA_per_picoF).
 * CONSTANTS[70] is Km_Nai in component sodium_potassium_pump_current (micromolar).
 * CONSTANTS[71] is Km_Ko in component sodium_potassium_pump_current (micromolar).
 * ALGEBRAIC[62] is f_NaK in component sodium_potassium_pump_current (dimensionless).
 * CONSTANTS[75] is sigma in component sodium_potassium_pump_current (dimensionless).
 * CONSTANTS[72] is P_ClCa in component calcium_activated_chloride_current (cm_per_second).
 * ALGEBRAIC[64] is f_ClCa in component calcium_activated_chloride_current (dimensionless).
 * CONSTANTS[73] is Km_Cl in component calcium_activated_chloride_current (micromolar).
 * CONSTANTS[74] is n_ClCa in component calcium_activated_chloride_current (dimensionless).
 * RATES[0] is d/dt V in component membrane (millivolt).
 * RATES[1] is d/dt Cai in component calcium_concentration (micromolar).
 * RATES[2] is d/dt Cass in component calcium_concentration (micromolar).
 * RATES[3] is d/dt CaJSR in component calcium_concentration (micromolar).
 * RATES[4] is d/dt CaNSR in component calcium_concentration (micromolar).
 * RATES[5] is d/dt P_RyR in component calcium_fluxes (dimensionless).
 * RATES[6] is d/dt LTRPN_Ca in component calcium_buffering (micromolar).
 * RATES[7] is d/dt HTRPN_Ca in component calcium_buffering (micromolar).
 * RATES[8] is d/dt P_O1 in component ryanodine_receptors (dimensionless).
 * RATES[9] is d/dt P_O2 in component ryanodine_receptors (dimensionless).
 * RATES[10] is d/dt P_C2 in component ryanodine_receptors (dimensionless).
 * RATES[11] is d/dt O in component L_type_calcium_current (dimensionless).
 * RATES[12] is d/dt C2 in component L_type_calcium_current (dimensionless).
 * RATES[13] is d/dt C3 in component L_type_calcium_current (dimensionless).
 * RATES[14] is d/dt C4 in component L_type_calcium_current (dimensionless).
 * RATES[15] is d/dt I1 in component L_type_calcium_current (dimensionless).
 * RATES[16] is d/dt I2 in component L_type_calcium_current (dimensionless).
 * RATES[17] is d/dt I3 in component L_type_calcium_current (dimensionless).
 * RATES[18] is d/dt b in component T_type_calcium_current_b_gate (dimensionless).
 * RATES[19] is d/dt g in component T_type_calcium_current_g_gate (dimensionless).
 * RATES[20] is d/dt Nai in component sodium_concentration (micromolar).
 * RATES[21] is d/dt m in component fast_sodium_current_m_gate (dimensionless).
 * RATES[22] is d/dt h in component fast_sodium_current_h_gate (dimensionless).
 * RATES[23] is d/dt j in component fast_sodium_current_j_gate (dimensionless).
 * RATES[24] is d/dt Ki in component potassium_concentration (micromolar).
 * RATES[25] is d/dt ato_f in component fast_transient_outward_potassium_current (dimensionless).
 * RATES[26] is d/dt ito_f in component fast_transient_outward_potassium_current (dimensionless).
 * RATES[27] is d/dt aur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless).
 * RATES[28] is d/dt iur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless).
 * RATES[29] is d/dt aKss in component non_inactivating_steady_state_potassium_current (dimensionless).
 * RATES[30] is d/dt nKs in component slow_delayed_rectifier_potassium_current (dimensionless).
 * RATES[33] is d/dt C_K2 in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * RATES[32] is d/dt C_K1 in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * RATES[31] is d/dt O_K in component rapid_delayed_rectifier_potassium_current (dimensionless).
 * RATES[34] is d/dt I_K in component rapid_delayed_rectifier_potassium_current (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -80.6475;
CONSTANTS[0] = 1;
CONSTANTS[1] = 2.2826e-6;
CONSTANTS[2] = 0.12e-8;
CONSTANTS[3] = 2.098e-7;
CONSTANTS[4] = 3.0734e-8;
CONSTANTS[5] = 1.5410e-5;
CONSTANTS[6] = 5400;
CONSTANTS[7] = 140000;
CONSTANTS[8] = 1000;
CONSTANTS[9] = 132000;
CONSTANTS[10] = 30000;
CONSTANTS[11] = 8.314;
CONSTANTS[12] = 298;
CONSTANTS[13] = 96.5;
CONSTANTS[14] = 20;
CONSTANTS[15] = 500000;
CONSTANTS[16] = 2000.0;
CONSTANTS[17] = 0.5;
CONSTANTS[18] = -80;
STATES[1] = 0.2049;
STATES[2] = 0.1867;
STATES[3] = 705.5417;
STATES[4] = 703.3063;
CONSTANTS[19] = 25;
CONSTANTS[20] = 15000;
CONSTANTS[21] = 0.238;
CONSTANTS[22] = 800;
CONSTANTS[23] = 0.00237;
CONSTANTS[24] = 3.2e-5;
CONSTANTS[25] = 0.0327;
CONSTANTS[26] = 0.0196;
STATES[5] = 2.8603e-14;
CONSTANTS[27] = 0.45;
CONSTANTS[28] = 20;
CONSTANTS[29] = 2.088e-5;
CONSTANTS[30] = 8;
CONSTANTS[31] = 0.09;
CONSTANTS[32] = 0.5;
CONSTANTS[33] = 35;
CONSTANTS[34] = 70;
STATES[6] = 8.9220;
STATES[7] = 66.0407;
CONSTANTS[35] = 7;
STATES[8] = 1.1e-3;
STATES[9] = 2.9799e-8;
STATES[10] = 0.1003;
CONSTANTS[36] = 0.00608;
CONSTANTS[37] = 0.07125;
CONSTANTS[38] = 0.00405;
CONSTANTS[39] = 0.965;
CONSTANTS[40] = 0.009;
CONSTANTS[41] = 0.0008;
CONSTANTS[42] = 3;
CONSTANTS[43] = 4;
CONSTANTS[44] = 63;
CONSTANTS[45] = 0.19019;
STATES[11] = 4.4776e-12;
STATES[12] = 1.6996e-4;
STATES[13] = 1.1754e-8;
STATES[14] = 6.1239e-10;
STATES[15] = 1.5450e-9;
STATES[16] = 6.4226e-8;
STATES[17] = 6.5201e-7;
CONSTANTS[46] = 0.0005;
CONSTANTS[47] = 0.23324;
CONSTANTS[48] = 20;
CONSTANTS[49] = 50;
CONSTANTS[50] = 0.055;
STATES[18] = 0.0047;
STATES[19] = 0.9020;
CONSTANTS[51] = 0.2;
CONSTANTS[52] = 0.5;
CONSTANTS[53] = 907.68;
CONSTANTS[54] = 87500;
CONSTANTS[55] = 1380;
CONSTANTS[56] = 0.1;
CONSTANTS[57] = 0.35;
STATES[20] = 2.1747e4;
CONSTANTS[58] = 0.00025;
CONSTANTS[59] = 10;
STATES[21] = 0.0032;
STATES[22] = 0.9606;
STATES[23] = 0.9744;
CONSTANTS[60] = 0.0026;
STATES[24] = 1.3645e5;
CONSTANTS[61] = 0.1017;
STATES[25] = 0.0032;
STATES[26] = 0.7530;
CONSTANTS[62] = 0.0048;
STATES[27] = 5.2471e-4;
STATES[28] = 0.9866;
CONSTANTS[63] = 0.015;
STATES[29] = 0.8500;
CONSTANTS[64] = 0.235;
CONSTANTS[65] = 0.046;
STATES[30] = 0.0026;
CONSTANTS[66] = 01.17;
STATES[31] = 2.1e-3;
STATES[32] = 1.1e-3;
STATES[33] = 8.3641e-4;
STATES[34] = 4.3522e-4;
CONSTANTS[67] = 0.036778;
CONSTANTS[68] = 0.023761;
CONSTANTS[69] = 0.88;
CONSTANTS[70] = 21000;
CONSTANTS[71] = 1500;
CONSTANTS[72] = 2.74e-7;
CONSTANTS[73] = 4;
CONSTANTS[74] = 3;
CONSTANTS[75] =  (1.00000/7.00000)*(exp(CONSTANTS[7]/67300.0) - 1.00000);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[6] =  CONSTANTS[25]*STATES[1]*(CONSTANTS[33] - STATES[6]) -  CONSTANTS[26]*STATES[6];
RATES[7] =  CONSTANTS[23]*STATES[1]*(CONSTANTS[34] - STATES[7]) -  CONSTANTS[24]*STATES[7];
RATES[9] =  CONSTANTS[38]*pow(STATES[2], CONSTANTS[42])*STATES[8] -  CONSTANTS[39]*STATES[9];
RATES[10] =  CONSTANTS[40]*STATES[8] -  CONSTANTS[41]*STATES[10];
ALGEBRAIC[1] = 1.00000 - (STATES[10]+STATES[8]+STATES[9]);
RATES[8] = ( CONSTANTS[36]*pow(STATES[2], CONSTANTS[43])*ALGEBRAIC[1]+ CONSTANTS[39]*STATES[9]+ CONSTANTS[41]*STATES[10]) - ( CONSTANTS[37]*STATES[8]+ CONSTANTS[38]*pow(STATES[2], CONSTANTS[42])*STATES[8]+ CONSTANTS[40]*STATES[8]);
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (STATES[0]+48.0000)/6.10000));
ALGEBRAIC[17] = 0.100000+5.40000/(1.00000+exp((STATES[0]+100.000)/6.60000));
RATES[18] = (ALGEBRAIC[3] - STATES[18])/ALGEBRAIC[17];
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+66.0000)/6.60000));
ALGEBRAIC[18] = 8.00000+32.0000/(1.00000+exp((STATES[0]+65.0000)/5.00000));
RATES[19] = (ALGEBRAIC[4] - STATES[19])/ALGEBRAIC[18];
ALGEBRAIC[5] = ( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( -0.100000*(STATES[0]+47.1300)));
ALGEBRAIC[19] =  0.0800000*exp(- STATES[0]/11.0000);
RATES[21] = ALGEBRAIC[5] -  (ALGEBRAIC[5]+ALGEBRAIC[19])*STATES[21];
ALGEBRAIC[6] = (STATES[0]<- 40.0000 ?  0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[20] = (STATES[0]<- 40.0000 ?  3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/- 11.1000))));
RATES[22] = ALGEBRAIC[6] -  (ALGEBRAIC[6]+ALGEBRAIC[20])*STATES[22];
ALGEBRAIC[7] = (STATES[0]<- 40.0000 ? ( - ( 127140.*exp( 0.244400*STATES[0])+ 3.47400e-05*exp( - 0.0439100*STATES[0]))*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[21] = (STATES[0]<- 40.0000 ? ( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( - 2.53500e-07*STATES[0]))/(1.00000+exp( - 0.100000*(STATES[0]+32.0000))));
RATES[23] = ALGEBRAIC[7] -  (ALGEBRAIC[7]+ALGEBRAIC[21])*STATES[23];
ALGEBRAIC[9] =  0.180640*exp( 0.0357700*(STATES[0]+30.0000));
ALGEBRAIC[23] =  0.395600*exp( - 0.0623700*(STATES[0]+30.0000));
RATES[25] =  ALGEBRAIC[9]*(1.00000 - STATES[25]) -  ALGEBRAIC[23]*STATES[25];
ALGEBRAIC[10] = 1.00000/(1.00000+exp(- (STATES[0]+22.5000)/7.70000));
ALGEBRAIC[24] =  0.493000*exp( - 0.0629000*STATES[0])+2.05800;
RATES[27] = (ALGEBRAIC[10] - STATES[27])/ALGEBRAIC[24];
ALGEBRAIC[25] = 1200.00 - 170.000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[11] = 1.00000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
RATES[28] = (ALGEBRAIC[11] - STATES[28])/ALGEBRAIC[25];
ALGEBRAIC[26] =  39.3000*exp( - 0.0862000*STATES[0])+13.1700;
RATES[29] = (ALGEBRAIC[10] - STATES[29])/ALGEBRAIC[26];
ALGEBRAIC[12] = ( 4.81333e-06*(STATES[0]+26.5000))/(1.00000 - exp( - 0.128000*(STATES[0]+26.5000)));
ALGEBRAIC[27] =  9.53333e-05*exp( - 0.0380000*(STATES[0]+26.5000));
RATES[30] =  ALGEBRAIC[12]*(1.00000 - STATES[30]) -  ALGEBRAIC[27]*STATES[30];
ALGEBRAIC[14] =  0.0137330*exp( 0.0381980*STATES[0]);
ALGEBRAIC[29] =  6.89000e-05*exp( - 0.0417800*STATES[0]);
RATES[33] = ( CONSTANTS[68]*STATES[32]+ ALGEBRAIC[29]*STATES[31]) - ( CONSTANTS[67]*STATES[33]+ ALGEBRAIC[14]*STATES[33]);
ALGEBRAIC[2] = 1.00000 - (STATES[11]+STATES[12]+STATES[12]+STATES[13]+STATES[14]+STATES[15]+STATES[16]+STATES[17]);
ALGEBRAIC[16] = ( 0.400000*exp((STATES[0]+12.0000)/10.0000)*((1.00000+ 0.700000*exp(- pow(STATES[0]+40.0000, 2.00000)/10.0000)) -  0.750000*exp(- pow(STATES[0]+20.0000, 2.00000)/400.000)))/(1.00000+ 0.120000*exp((STATES[0]+12.0000)/10.0000));
ALGEBRAIC[31] =  0.0500000*exp(- (STATES[0]+12.0000)/13.0000);
RATES[12] = ( 4.00000*ALGEBRAIC[16]*ALGEBRAIC[2]+ 2.00000*ALGEBRAIC[31]*STATES[13]) - ( ALGEBRAIC[31]*STATES[12]+ 3.00000*ALGEBRAIC[16]*STATES[12]);
RATES[13] = ( 3.00000*ALGEBRAIC[16]*STATES[12]+ 3.00000*ALGEBRAIC[31]*STATES[14]) - ( 2.00000*ALGEBRAIC[31]*STATES[13]+ 2.00000*ALGEBRAIC[16]*STATES[13]);
ALGEBRAIC[13] = 1.00000 - (STATES[32]+STATES[33]+STATES[31]+STATES[34]);
ALGEBRAIC[28] =  0.0223480*exp( 0.0117600*STATES[0]);
ALGEBRAIC[33] =  0.0470020*exp( - 0.0631000*STATES[0]);
RATES[32] = ( ALGEBRAIC[28]*ALGEBRAIC[13]+ CONSTANTS[67]*STATES[33]) - ( ALGEBRAIC[33]*STATES[32]+ CONSTANTS[68]*STATES[32]);
ALGEBRAIC[22] = ( 0.000152000*exp(- (STATES[0] - 3.81000)/15.7500))/( 0.00670830*exp(- (STATES[0]+132.050)/15.7500)+1.00000);
ALGEBRAIC[32] = ( 0.000950000*exp((STATES[0]+132.050)/15.7500))/( 0.0513350*exp((STATES[0]+132.050)/15.7500)+1.00000);
ALGEBRAIC[37] = ALGEBRAIC[22]/(ALGEBRAIC[22]+ALGEBRAIC[32]);
ALGEBRAIC[8] = 1.00000/(( 0.000152000*exp(- (STATES[0]+13.5000)/7.00000))/( 0.0670830*exp(- (STATES[0]+33.5000)/7.00000)+1.00000)+( 0.000950000*exp((STATES[0]+33.5000)/7.00000))/( 0.0513350*exp((STATES[0]+33.5000)/7.00000)+1.00000));
RATES[26] = (ALGEBRAIC[37] - STATES[26])/ALGEBRAIC[8];
ALGEBRAIC[34] =  0.0908210*exp( 0.0233910*(STATES[0]+5.00000));
ALGEBRAIC[38] =  0.00649700*exp( - 0.0326800*(STATES[0]+5.00000));
RATES[31] = ( ALGEBRAIC[14]*STATES[33]+ ALGEBRAIC[38]*STATES[34]) - ( ALGEBRAIC[29]*STATES[31]+ ALGEBRAIC[34]*STATES[31]);
RATES[34] =  ALGEBRAIC[34]*STATES[31] -  ALGEBRAIC[38]*STATES[34];
ALGEBRAIC[36] = ( CONSTANTS[47]*STATES[2])/(CONSTANTS[48]+STATES[2]);
ALGEBRAIC[40] =  13.0000*(1.00000 - exp(- pow(STATES[0]+14.5000, 2.00000)/100.000));
RATES[11] = ( ALGEBRAIC[16]*STATES[14]+ CONSTANTS[46]*STATES[15]+ 0.00100000*( ALGEBRAIC[16]*STATES[16] -  ALGEBRAIC[40]*STATES[11])) - ( 4.00000*ALGEBRAIC[31]*STATES[11]+ ALGEBRAIC[36]*STATES[11]);
RATES[14] = ( 2.00000*ALGEBRAIC[16]*STATES[13]+ 4.00000*ALGEBRAIC[31]*STATES[11]+ 0.0100000*( 4.00000*CONSTANTS[46]*ALGEBRAIC[31]*STATES[15] -  ALGEBRAIC[16]*ALGEBRAIC[36]*STATES[14])+ 0.00200000*( 4.00000*ALGEBRAIC[31]*STATES[16] -  ALGEBRAIC[40]*STATES[14])+ 4.00000*ALGEBRAIC[31]*CONSTANTS[46]*STATES[17]) - ( 3.00000*ALGEBRAIC[31]*STATES[14]+ ALGEBRAIC[16]*STATES[14]+ 1.00000*ALGEBRAIC[36]*ALGEBRAIC[40]*STATES[14]);
RATES[15] = ( ALGEBRAIC[36]*STATES[11]+ 0.00100000*( ALGEBRAIC[16]*STATES[17] -  ALGEBRAIC[40]*STATES[15])+ 0.0100000*( ALGEBRAIC[16]*ALGEBRAIC[36]*STATES[14] -  4.00000*ALGEBRAIC[31]*ALGEBRAIC[40]*STATES[15])) -  CONSTANTS[46]*STATES[15];
RATES[16] = ( 0.00100000*( ALGEBRAIC[40]*STATES[11] -  ALGEBRAIC[16]*STATES[16])+ CONSTANTS[46]*STATES[17]+ 0.00200000*( ALGEBRAIC[40]*STATES[14] -  4.00000*ALGEBRAIC[31]*STATES[16])) -  ALGEBRAIC[36]*STATES[16];
RATES[17] = ( 0.00100000*( ALGEBRAIC[40]*STATES[15] -  ALGEBRAIC[16]*STATES[17])+ ALGEBRAIC[36]*STATES[16]+ 1.00000*ALGEBRAIC[36]*ALGEBRAIC[40]*STATES[14]) - ( 4.00000*ALGEBRAIC[31]*CONSTANTS[46]*STATES[17]+ CONSTANTS[46]*STATES[17]);
ALGEBRAIC[35] = pow(1.00000+( CONSTANTS[20]*CONSTANTS[22])/pow(CONSTANTS[22]+STATES[3], 2.00000), - 1.00000);
ALGEBRAIC[39] =  CONSTANTS[27]*(STATES[8]+STATES[9])*(STATES[3] - STATES[2])*STATES[5];
ALGEBRAIC[42] = (STATES[2] - STATES[1])/CONSTANTS[30];
RATES[3] =  ALGEBRAIC[35]*(ALGEBRAIC[42] - ALGEBRAIC[39]);
ALGEBRAIC[43] =  CONSTANTS[45]*STATES[11]*(STATES[0] - CONSTANTS[44]);
RATES[5] =  - 0.0400000*STATES[5] -  (( 0.100000*ALGEBRAIC[43])/CONSTANTS[35])*exp(- pow(STATES[0] - 5.00000, 2.00000)/648.000);
ALGEBRAIC[44] =  CONSTANTS[29]*(STATES[4] - STATES[1]);
ALGEBRAIC[46] = ( CONSTANTS[31]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[32], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[41] = (STATES[4] - STATES[3])/CONSTANTS[28];
RATES[4] = ( (ALGEBRAIC[46] - ALGEBRAIC[44])*CONSTANTS[1])/CONSTANTS[3] - ( ALGEBRAIC[41]*CONSTANTS[2])/CONSTANTS[3];
ALGEBRAIC[15] = pow(1.00000+( CONSTANTS[19]*CONSTANTS[21])/pow(CONSTANTS[21]+STATES[1], 2.00000), - 1.00000);
ALGEBRAIC[48] = ( CONSTANTS[23]*STATES[1]*(CONSTANTS[34] - STATES[7])+ CONSTANTS[25]*STATES[1]*(CONSTANTS[33] - STATES[6])) - ( CONSTANTS[24]*STATES[7]+ CONSTANTS[26]*STATES[6]);
RATES[1] =  ALGEBRAIC[15]*((ALGEBRAIC[44]+ALGEBRAIC[42]) - (ALGEBRAIC[46]+ALGEBRAIC[48]));
ALGEBRAIC[45] =  CONSTANTS[50]*STATES[18]*STATES[19]*(STATES[0] - CONSTANTS[49]);
ALGEBRAIC[47] = ( CONSTANTS[51]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[52], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[49] =  (( (( (( CONSTANTS[53]*1.00000)/(pow(CONSTANTS[54], 3.00000)+pow(CONSTANTS[7], 3.00000)))*1.00000)/(CONSTANTS[55]+CONSTANTS[8]))*1.00000)/(1.00000+ CONSTANTS[56]*exp(( (CONSTANTS[57] - 1.00000)*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))))*( exp(( CONSTANTS[57]*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))*pow(STATES[20], 3.00000)*CONSTANTS[8] -  exp(( (CONSTANTS[57] - 1.00000)*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))*pow(CONSTANTS[7], 3.00000)*STATES[1]);
ALGEBRAIC[50] =  (( CONSTANTS[11]*CONSTANTS[12])/( 2.00000*CONSTANTS[13]))*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[51] =  CONSTANTS[58]*(STATES[0] - ALGEBRAIC[50]);
ALGEBRAIC[30] = pow(1.00000+( CONSTANTS[19]*CONSTANTS[21])/pow(CONSTANTS[21]+STATES[2], 2.00000), - 1.00000);
RATES[2] =  ALGEBRAIC[30]*(( ALGEBRAIC[39]*CONSTANTS[2])/CONSTANTS[4] - (( ALGEBRAIC[42]*CONSTANTS[1])/CONSTANTS[4]+ ((ALGEBRAIC[43]+ALGEBRAIC[51]+ALGEBRAIC[47]+ALGEBRAIC[45]) -  2.00000*ALGEBRAIC[49])*(( CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[4]*CONSTANTS[13]))));
ALGEBRAIC[52] =  (( CONSTANTS[11]*CONSTANTS[12])/CONSTANTS[13])*log(CONSTANTS[7]/STATES[20]);
ALGEBRAIC[53] =  CONSTANTS[59]*pow(STATES[21], 3.00000)*STATES[22]*STATES[23]*(STATES[0] - ALGEBRAIC[52]);
ALGEBRAIC[54] =  CONSTANTS[60]*(STATES[0] - ALGEBRAIC[52]);
ALGEBRAIC[62] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))+ 0.0365000*CONSTANTS[75]*exp(( - STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12])));
ALGEBRAIC[63] = ( (( CONSTANTS[69]*ALGEBRAIC[62]*1.00000)/(1.00000+pow(CONSTANTS[70]/STATES[20], 1.50000)))*CONSTANTS[6])/(CONSTANTS[6]+CONSTANTS[71]);
RATES[20] = ( - (ALGEBRAIC[53]+ALGEBRAIC[54]+ 3.00000*ALGEBRAIC[63]+ 3.00000*ALGEBRAIC[49])*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[13]);
ALGEBRAIC[55] =  (( CONSTANTS[11]*CONSTANTS[12])/CONSTANTS[13])*log(CONSTANTS[6]/STATES[24]);
ALGEBRAIC[56] =  CONSTANTS[61]*pow(STATES[25], 6.50000)*STATES[26]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[59] =  CONSTANTS[64]*(CONSTANTS[6]/(CONSTANTS[6]+210.000))*((STATES[0] - ALGEBRAIC[55])/(1.00000+exp( 0.0896000*(STATES[0] - ALGEBRAIC[55]))));
ALGEBRAIC[60] =  CONSTANTS[65]*pow(STATES[30], 2.00000)*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[57] =  CONSTANTS[62]*STATES[27]*STATES[28]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[58] =  CONSTANTS[63]*STATES[29]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[61] =  CONSTANTS[66]*STATES[31]*(STATES[0] - ALGEBRAIC[55]);
RATES[24] = ( - ((ALGEBRAIC[56]+ALGEBRAIC[59]+ALGEBRAIC[60]+ALGEBRAIC[58]+ALGEBRAIC[57]+ALGEBRAIC[61]) -  2.00000*ALGEBRAIC[63])*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[13]);
ALGEBRAIC[0] = (VOI>=CONSTANTS[14]&&VOI<=CONSTANTS[15]&&(VOI - CONSTANTS[14]) -  floor((VOI - CONSTANTS[14])/CONSTANTS[16])*CONSTANTS[16]<=CONSTANTS[17] ? CONSTANTS[18] : 0.00000);
ALGEBRAIC[64] = pow(STATES[2]/(CONSTANTS[73]+STATES[2]), CONSTANTS[74]);
ALGEBRAIC[65] =  1.00000*CONSTANTS[72]*ALGEBRAIC[64]*(( STATES[0]*pow(CONSTANTS[13], 2.00000))/( CONSTANTS[11]*CONSTANTS[12]))*(( CONSTANTS[9]*exp(( STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12])) - CONSTANTS[10])/(exp(( STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12])) - 1.00000));
RATES[0] = - (ALGEBRAIC[43]+ALGEBRAIC[45]+ALGEBRAIC[47]+ALGEBRAIC[49]+ALGEBRAIC[51]+ALGEBRAIC[53]+ALGEBRAIC[54]+ALGEBRAIC[63]+ALGEBRAIC[56]+ALGEBRAIC[59]+ALGEBRAIC[60]+ALGEBRAIC[57]+ALGEBRAIC[58]+ALGEBRAIC[61]+ALGEBRAIC[65]+ALGEBRAIC[0]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = 1.00000 - (STATES[10]+STATES[8]+STATES[9]);
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (STATES[0]+48.0000)/6.10000));
ALGEBRAIC[17] = 0.100000+5.40000/(1.00000+exp((STATES[0]+100.000)/6.60000));
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+66.0000)/6.60000));
ALGEBRAIC[18] = 8.00000+32.0000/(1.00000+exp((STATES[0]+65.0000)/5.00000));
ALGEBRAIC[5] = ( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( -0.100000*(STATES[0]+47.1300)));
ALGEBRAIC[19] =  0.0800000*exp(- STATES[0]/11.0000);
ALGEBRAIC[6] = (STATES[0]<- 40.0000 ?  0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[20] = (STATES[0]<- 40.0000 ?  3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/- 11.1000))));
ALGEBRAIC[7] = (STATES[0]<- 40.0000 ? ( - ( 127140.*exp( 0.244400*STATES[0])+ 3.47400e-05*exp( - 0.0439100*STATES[0]))*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[21] = (STATES[0]<- 40.0000 ? ( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( - 2.53500e-07*STATES[0]))/(1.00000+exp( - 0.100000*(STATES[0]+32.0000))));
ALGEBRAIC[9] =  0.180640*exp( 0.0357700*(STATES[0]+30.0000));
ALGEBRAIC[23] =  0.395600*exp( - 0.0623700*(STATES[0]+30.0000));
ALGEBRAIC[10] = 1.00000/(1.00000+exp(- (STATES[0]+22.5000)/7.70000));
ALGEBRAIC[24] =  0.493000*exp( - 0.0629000*STATES[0])+2.05800;
ALGEBRAIC[25] = 1200.00 - 170.000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[11] = 1.00000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[26] =  39.3000*exp( - 0.0862000*STATES[0])+13.1700;
ALGEBRAIC[12] = ( 4.81333e-06*(STATES[0]+26.5000))/(1.00000 - exp( - 0.128000*(STATES[0]+26.5000)));
ALGEBRAIC[27] =  9.53333e-05*exp( - 0.0380000*(STATES[0]+26.5000));
ALGEBRAIC[14] =  0.0137330*exp( 0.0381980*STATES[0]);
ALGEBRAIC[29] =  6.89000e-05*exp( - 0.0417800*STATES[0]);
ALGEBRAIC[2] = 1.00000 - (STATES[11]+STATES[12]+STATES[12]+STATES[13]+STATES[14]+STATES[15]+STATES[16]+STATES[17]);
ALGEBRAIC[16] = ( 0.400000*exp((STATES[0]+12.0000)/10.0000)*((1.00000+ 0.700000*exp(- pow(STATES[0]+40.0000, 2.00000)/10.0000)) -  0.750000*exp(- pow(STATES[0]+20.0000, 2.00000)/400.000)))/(1.00000+ 0.120000*exp((STATES[0]+12.0000)/10.0000));
ALGEBRAIC[31] =  0.0500000*exp(- (STATES[0]+12.0000)/13.0000);
ALGEBRAIC[13] = 1.00000 - (STATES[32]+STATES[33]+STATES[31]+STATES[34]);
ALGEBRAIC[28] =  0.0223480*exp( 0.0117600*STATES[0]);
ALGEBRAIC[33] =  0.0470020*exp( - 0.0631000*STATES[0]);
ALGEBRAIC[22] = ( 0.000152000*exp(- (STATES[0] - 3.81000)/15.7500))/( 0.00670830*exp(- (STATES[0]+132.050)/15.7500)+1.00000);
ALGEBRAIC[32] = ( 0.000950000*exp((STATES[0]+132.050)/15.7500))/( 0.0513350*exp((STATES[0]+132.050)/15.7500)+1.00000);
ALGEBRAIC[37] = ALGEBRAIC[22]/(ALGEBRAIC[22]+ALGEBRAIC[32]);
ALGEBRAIC[8] = 1.00000/(( 0.000152000*exp(- (STATES[0]+13.5000)/7.00000))/( 0.0670830*exp(- (STATES[0]+33.5000)/7.00000)+1.00000)+( 0.000950000*exp((STATES[0]+33.5000)/7.00000))/( 0.0513350*exp((STATES[0]+33.5000)/7.00000)+1.00000));
ALGEBRAIC[34] =  0.0908210*exp( 0.0233910*(STATES[0]+5.00000));
ALGEBRAIC[38] =  0.00649700*exp( - 0.0326800*(STATES[0]+5.00000));
ALGEBRAIC[36] = ( CONSTANTS[47]*STATES[2])/(CONSTANTS[48]+STATES[2]);
ALGEBRAIC[40] =  13.0000*(1.00000 - exp(- pow(STATES[0]+14.5000, 2.00000)/100.000));
ALGEBRAIC[35] = pow(1.00000+( CONSTANTS[20]*CONSTANTS[22])/pow(CONSTANTS[22]+STATES[3], 2.00000), - 1.00000);
ALGEBRAIC[39] =  CONSTANTS[27]*(STATES[8]+STATES[9])*(STATES[3] - STATES[2])*STATES[5];
ALGEBRAIC[42] = (STATES[2] - STATES[1])/CONSTANTS[30];
ALGEBRAIC[43] =  CONSTANTS[45]*STATES[11]*(STATES[0] - CONSTANTS[44]);
ALGEBRAIC[44] =  CONSTANTS[29]*(STATES[4] - STATES[1]);
ALGEBRAIC[46] = ( CONSTANTS[31]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[32], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[41] = (STATES[4] - STATES[3])/CONSTANTS[28];
ALGEBRAIC[15] = pow(1.00000+( CONSTANTS[19]*CONSTANTS[21])/pow(CONSTANTS[21]+STATES[1], 2.00000), - 1.00000);
ALGEBRAIC[48] = ( CONSTANTS[23]*STATES[1]*(CONSTANTS[34] - STATES[7])+ CONSTANTS[25]*STATES[1]*(CONSTANTS[33] - STATES[6])) - ( CONSTANTS[24]*STATES[7]+ CONSTANTS[26]*STATES[6]);
ALGEBRAIC[45] =  CONSTANTS[50]*STATES[18]*STATES[19]*(STATES[0] - CONSTANTS[49]);
ALGEBRAIC[47] = ( CONSTANTS[51]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[52], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[49] =  (( (( (( CONSTANTS[53]*1.00000)/(pow(CONSTANTS[54], 3.00000)+pow(CONSTANTS[7], 3.00000)))*1.00000)/(CONSTANTS[55]+CONSTANTS[8]))*1.00000)/(1.00000+ CONSTANTS[56]*exp(( (CONSTANTS[57] - 1.00000)*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))))*( exp(( CONSTANTS[57]*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))*pow(STATES[20], 3.00000)*CONSTANTS[8] -  exp(( (CONSTANTS[57] - 1.00000)*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))*pow(CONSTANTS[7], 3.00000)*STATES[1]);
ALGEBRAIC[50] =  (( CONSTANTS[11]*CONSTANTS[12])/( 2.00000*CONSTANTS[13]))*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[51] =  CONSTANTS[58]*(STATES[0] - ALGEBRAIC[50]);
ALGEBRAIC[30] = pow(1.00000+( CONSTANTS[19]*CONSTANTS[21])/pow(CONSTANTS[21]+STATES[2], 2.00000), - 1.00000);
ALGEBRAIC[52] =  (( CONSTANTS[11]*CONSTANTS[12])/CONSTANTS[13])*log(CONSTANTS[7]/STATES[20]);
ALGEBRAIC[53] =  CONSTANTS[59]*pow(STATES[21], 3.00000)*STATES[22]*STATES[23]*(STATES[0] - ALGEBRAIC[52]);
ALGEBRAIC[54] =  CONSTANTS[60]*(STATES[0] - ALGEBRAIC[52]);
ALGEBRAIC[62] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12]))+ 0.0365000*CONSTANTS[75]*exp(( - STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12])));
ALGEBRAIC[63] = ( (( CONSTANTS[69]*ALGEBRAIC[62]*1.00000)/(1.00000+pow(CONSTANTS[70]/STATES[20], 1.50000)))*CONSTANTS[6])/(CONSTANTS[6]+CONSTANTS[71]);
ALGEBRAIC[55] =  (( CONSTANTS[11]*CONSTANTS[12])/CONSTANTS[13])*log(CONSTANTS[6]/STATES[24]);
ALGEBRAIC[56] =  CONSTANTS[61]*pow(STATES[25], 6.50000)*STATES[26]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[59] =  CONSTANTS[64]*(CONSTANTS[6]/(CONSTANTS[6]+210.000))*((STATES[0] - ALGEBRAIC[55])/(1.00000+exp( 0.0896000*(STATES[0] - ALGEBRAIC[55]))));
ALGEBRAIC[60] =  CONSTANTS[65]*pow(STATES[30], 2.00000)*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[57] =  CONSTANTS[62]*STATES[27]*STATES[28]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[58] =  CONSTANTS[63]*STATES[29]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[61] =  CONSTANTS[66]*STATES[31]*(STATES[0] - ALGEBRAIC[55]);
ALGEBRAIC[0] = (VOI>=CONSTANTS[14]&&VOI<=CONSTANTS[15]&&(VOI - CONSTANTS[14]) -  floor((VOI - CONSTANTS[14])/CONSTANTS[16])*CONSTANTS[16]<=CONSTANTS[17] ? CONSTANTS[18] : 0.00000);
ALGEBRAIC[64] = pow(STATES[2]/(CONSTANTS[73]+STATES[2]), CONSTANTS[74]);
ALGEBRAIC[65] =  1.00000*CONSTANTS[72]*ALGEBRAIC[64]*(( STATES[0]*pow(CONSTANTS[13], 2.00000))/( CONSTANTS[11]*CONSTANTS[12]))*(( CONSTANTS[9]*exp(( STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12])) - CONSTANTS[10])/(exp(( STATES[0]*CONSTANTS[13])/( CONSTANTS[11]*CONSTANTS[12])) - 1.00000));
}