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 81 entries in the algebraic variable array.
There are a total of 36 entries in each of the rate and state variable arrays.
There are a total of 77 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 R in component membrane (joule_per_mole_kelvin).
* CONSTANTS[10] is T in component membrane (kelvin).
* CONSTANTS[11] is F in component membrane (coulomb_per_millimole).
* ALGEBRAIC[58] is i_CaL in component L_type_calcium_current (picoA_per_picoF).
* ALGEBRAIC[59] is i_pCa in component calcium_pump_current (picoA_per_picoF).
* ALGEBRAIC[60] is i_NaCa in component sodium_calcium_exchange_current (picoA_per_picoF).
* ALGEBRAIC[64] is i_Cab in component calcium_background_current (picoA_per_picoF).
* ALGEBRAIC[67] is i_Na in component fast_sodium_current (picoA_per_picoF).
* ALGEBRAIC[68] is i_Nab in component sodium_background_current (picoA_per_picoF).
* ALGEBRAIC[78] is i_NaK in component sodium_potassium_pump_current (picoA_per_picoF).
* ALGEBRAIC[70] is i_Kto_f in component fast_transient_outward_K_I (picoA_per_picoF).
* ALGEBRAIC[71] is i_Kto_s in component slow_transient_outward_K_I (picoA_per_picoF).
* ALGEBRAIC[72] is i_K1 in component time_independent_K_I (picoA_per_picoF).
* ALGEBRAIC[73] is i_Ks in component slow_delayed_rectifier_K_I (picoA_per_picoF).
* ALGEBRAIC[74] is i_Kur in component ultra_rapidly_activating_delayed_rectifier_K_I (picoA_per_picoF).
* ALGEBRAIC[75] is i_Kss in component non_inactivating_steady_state_K_I (picoA_per_picoF).
* ALGEBRAIC[80] is i_ClCa in component calcium_activated_chloride_current (picoA_per_picoF).
* ALGEBRAIC[76] is i_Kr in component rapid_delayed_rectifier_K_I (picoA_per_picoF).
* CONSTANTS[12] is stim_offset in component membrane (millisecond).
* CONSTANTS[13] is stim_period in component membrane (millisecond).
* CONSTANTS[14] is stim_duration in component membrane (millisecond).
* CONSTANTS[15] is stim_amplitude in component membrane (picoA_per_picoF).
* ALGEBRAIC[12] is i_Stim in component membrane (picoA_per_picoF).
* ALGEBRAIC[0] is past in component membrane (millisecond).
* 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[26] is Bi in component calcium_concentration (dimensionless).
* ALGEBRAIC[31] is Bss in component calcium_concentration (dimensionless).
* ALGEBRAIC[35] is BJSR in component calcium_concentration (dimensionless).
* CONSTANTS[16] is Bmax in component calcium_concentration (micromolar).
* CONSTANTS[17] is CSQN_tot in component calcium_concentration (micromolar).
* CONSTANTS[18] is Kd in component calcium_concentration (micromolar).
* CONSTANTS[19] is Km_CSQN in component calcium_concentration (micromolar).
* ALGEBRAIC[44] is J_leak in component calcium_fluxes (micromolar_per_millisecond).
* ALGEBRAIC[38] is J_rel in component calcium_fluxes (micromolar_per_millisecond).
* ALGEBRAIC[52] is J_serca in component calcium_fluxes (micromolar_per_millisecond).
* ALGEBRAIC[40] is J_tr in component calcium_fluxes (micromolar_per_millisecond).
* ALGEBRAIC[42] is J_xfer in component calcium_fluxes (micromolar_per_millisecond).
* STATES[5] is P_RyR in component calcium_fluxes (dimensionless).
* CONSTANTS[20] is v1 in component calcium_fluxes (per_millisecond).
* CONSTANTS[21] is tau_tr in component calcium_fluxes (millisecond).
* CONSTANTS[22] is v2 in component calcium_fluxes (per_millisecond).
* CONSTANTS[23] is tau_xfer in component calcium_fluxes (millisecond).
* ALGEBRAIC[46] is CaMKb in component calcium_fluxes (dimensionless).
* STATES[6] is CaMKt in component calcium_fluxes (dimensionless).
* ALGEBRAIC[48] is CaMKa in component calcium_fluxes (dimensionless).
* CONSTANTS[24] is on_rate in component calcium_fluxes (per_millisecond).
* CONSTANTS[25] is off_rate in component calcium_fluxes (per_millisecond).
* ALGEBRAIC[50] is vmup in component calcium_fluxes (micromolar_per_millisecond).
* CONSTANTS[26] is Km_up in component calcium_fluxes (micromolar).
* CONSTANTS[27] is i_CaL_max in component L_type_calcium_current (picoA_per_picoF).
* STATES[7] is P_O1 in component ryanodine_receptors (dimensionless).
* STATES[8] is P_O2 in component ryanodine_receptors (dimensionless).
* CONSTANTS[28] is vmup_init in component calcium_fluxes (micromolar_per_millisecond).
* CONSTANTS[29] is P_ryr_const1 in component calcium_fluxes (per_millisecond).
* CONSTANTS[30] is P_ryr_const2 in component calcium_fluxes (per_millisecond).
* ALGEBRAIC[1] is P_C1 in component ryanodine_receptors (dimensionless).
* STATES[9] is P_C2 in component ryanodine_receptors (dimensionless).
* CONSTANTS[31] is k_plus_a in component ryanodine_receptors (micromolar4_per_millisecond).
* CONSTANTS[32] is k_minus_a in component ryanodine_receptors (per_millisecond).
* CONSTANTS[33] is k_plus_b in component ryanodine_receptors (micromolar3_per_millisecond).
* CONSTANTS[34] is k_minus_b in component ryanodine_receptors (per_millisecond).
* CONSTANTS[35] is k_plus_c in component ryanodine_receptors (per_millisecond).
* CONSTANTS[36] is k_minus_c in component ryanodine_receptors (per_millisecond).
* CONSTANTS[37] is m in component ryanodine_receptors (dimensionless).
* CONSTANTS[38] is n in component ryanodine_receptors (dimensionless).
* CONSTANTS[39] is P_CaL in component L_type_calcium_current (per_millisecond).
* STATES[10] is O in component L_type_calcium_current (dimensionless).
* ALGEBRAIC[36] is C in component L_type_calcium_current (dimensionless).
* STATES[11] is I in component L_type_calcium_current (dimensionless).
* STATES[12] is y_gate in component L_type_calcium_current (dimensionless).
* ALGEBRAIC[2] is y_gate_inf in component L_type_calcium_current (dimensionless).
* ALGEBRAIC[13] is y_gate_tau in component L_type_calcium_current (millisecond).
* ALGEBRAIC[14] is alpha_p in component L_type_calcium_current (per_millisecond).
* CONSTANTS[74] is alpha_m in component L_type_calcium_current (per_millisecond).
* ALGEBRAIC[27] is epsilon_p in component L_type_calcium_current (per_micromolar_millisecond).
* ALGEBRAIC[32] is epsilon_m in component L_type_calcium_current (per_millisecond).
* CONSTANTS[40] is V_L in component L_type_calcium_current (millivolt).
* CONSTANTS[41] is delta_V_L in component L_type_calcium_current (millivolt).
* CONSTANTS[42] is t_L in component L_type_calcium_current (millisecond).
* CONSTANTS[43] is phi_L in component L_type_calcium_current (dimensionless).
* CONSTANTS[44] is a in component L_type_calcium_current (dimensionless).
* CONSTANTS[45] is b in component L_type_calcium_current (dimensionless).
* CONSTANTS[46] is tau_L in component L_type_calcium_current (millisecond).
* CONSTANTS[47] is K_L in component L_type_calcium_current (micromolar).
* ALGEBRAIC[3] is expVL in component L_type_calcium_current (dimensionless).
* ALGEBRAIC[54] is FVRT in component L_type_calcium_current (dimensionless).
* ALGEBRAIC[56] is FVRT_Ca in component L_type_calcium_current (dimensionless).
* CONSTANTS[48] is const5 in component L_type_calcium_current (millivolt).
* CONSTANTS[49] is i_pCa_max in component calcium_pump_current (picoA_per_picoF).
* CONSTANTS[50] is Km_pCa in component calcium_pump_current (micromolar).
* ALGEBRAIC[61] is J_pCa in component calcium_pump_current (micromolar_per_millisecond).
* ALGEBRAIC[63] is J_ncx in component sodium_calcium_exchange_current (micromolar_per_millisecond).
* CONSTANTS[51] is k_NaCa in component sodium_calcium_exchange_current (picoA_per_picoF).
* CONSTANTS[52] is K_mNa in component sodium_calcium_exchange_current (micromolar).
* CONSTANTS[53] is K_mCa in component sodium_calcium_exchange_current (micromolar).
* CONSTANTS[54] is k_sat in component sodium_calcium_exchange_current (dimensionless).
* CONSTANTS[55] is eta in component sodium_calcium_exchange_current (dimensionless).
* STATES[13] is Nai in component sodium_concentration (micromolar).
* CONSTANTS[56] is g_Cab in component calcium_background_current (milliS_per_microF).
* ALGEBRAIC[62] is E_CaN in component calcium_background_current (millivolt).
* ALGEBRAIC[66] is J_Cab in component calcium_background_current (micromolar_per_millisecond).
* ALGEBRAIC[65] is E_Na in component fast_sodium_current (millivolt).
* CONSTANTS[57] is g_Na in component fast_sodium_current (milliS_per_microF).
* STATES[14] is O_Na in component fast_sodium_current (dimensionless).
* STATES[15] is C_Na1 in component fast_sodium_current (dimensionless).
* STATES[16] is C_Na2 in component fast_sodium_current (dimensionless).
* ALGEBRAIC[4] is C_Na3 in component fast_sodium_current (dimensionless).
* STATES[17] is I1_Na in component fast_sodium_current (dimensionless).
* STATES[18] is I2_Na in component fast_sodium_current (dimensionless).
* STATES[19] is IF_Na in component fast_sodium_current (dimensionless).
* STATES[20] is IC_Na2 in component fast_sodium_current (dimensionless).
* STATES[21] is IC_Na3 in component fast_sodium_current (dimensionless).
* ALGEBRAIC[15] is alpha_Na11 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[37] is beta_Na11 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[28] is alpha_Na12 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[39] is beta_Na12 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[33] is alpha_Na13 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[41] is beta_Na13 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[43] is alpha_Na3 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[45] is beta_Na3 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[47] is alpha_Na2 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[49] is beta_Na2 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[51] is alpha_Na4 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[53] is beta_Na4 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[55] is alpha_Na5 in component fast_sodium_current (per_millisecond).
* ALGEBRAIC[57] is beta_Na5 in component fast_sodium_current (per_millisecond).
* STATES[22] is Ki in component potassium_concentration (micromolar).
* CONSTANTS[58] is g_Nab in component sodium_background_current (milliS_per_microF).
* ALGEBRAIC[69] is E_K in component fast_transient_outward_K_I (millivolt).
* CONSTANTS[59] is g_Kto_f in component fast_transient_outward_K_I (milliS_per_microF).
* STATES[23] is ato_f in component fast_transient_outward_K_I (dimensionless).
* STATES[24] is ito_f in component fast_transient_outward_K_I (dimensionless).
* ALGEBRAIC[5] is alpha_a in component fast_transient_outward_K_I (per_millisecond).
* ALGEBRAIC[16] is beta_a in component fast_transient_outward_K_I (per_millisecond).
* ALGEBRAIC[6] is alpha_i in component fast_transient_outward_K_I (per_millisecond).
* ALGEBRAIC[17] is beta_i in component fast_transient_outward_K_I (per_millisecond).
* ALGEBRAIC[7] is ass in component slow_transient_outward_K_I (dimensionless).
* ALGEBRAIC[8] is iss in component slow_transient_outward_K_I (dimensionless).
* CONSTANTS[60] is g_Kto_s in component slow_transient_outward_K_I (milliS_per_microF).
* STATES[25] is ato_s in component slow_transient_outward_K_I (dimensionless).
* STATES[26] is ito_s in component slow_transient_outward_K_I (dimensionless).
* ALGEBRAIC[18] is tau_ta_s in component slow_transient_outward_K_I (millisecond).
* ALGEBRAIC[19] is tau_ti_s in component slow_transient_outward_K_I (millisecond).
* CONSTANTS[61] is g_K1 in component time_independent_K_I (milliS_per_microF).
* CONSTANTS[62] is g_Ks in component slow_delayed_rectifier_K_I (milliS_per_microF).
* STATES[27] is nKs in component slow_delayed_rectifier_K_I (dimensionless).
* ALGEBRAIC[9] is alpha_n in component slow_delayed_rectifier_K_I (per_millisecond).
* ALGEBRAIC[20] is beta_n in component slow_delayed_rectifier_K_I (per_millisecond).
* CONSTANTS[63] is g_Kur in component ultra_rapidly_activating_delayed_rectifier_K_I (milliS_per_microF).
* STATES[28] is aur in component ultra_rapidly_activating_delayed_rectifier_K_I (dimensionless).
* STATES[29] is iur in component ultra_rapidly_activating_delayed_rectifier_K_I (dimensionless).
* ALGEBRAIC[21] is tau_aur in component ultra_rapidly_activating_delayed_rectifier_K_I (millisecond).
* ALGEBRAIC[22] is tau_iur in component ultra_rapidly_activating_delayed_rectifier_K_I (millisecond).
* CONSTANTS[64] is g_Kss in component non_inactivating_steady_state_K_I (milliS_per_microF).
* STATES[30] is aKss in component non_inactivating_steady_state_K_I (dimensionless).
* STATES[31] is iKss in component non_inactivating_steady_state_K_I (dimensionless).
* ALGEBRAIC[23] is tau_Kss in component non_inactivating_steady_state_K_I (millisecond).
* CONSTANTS[65] is g_Kr in component rapid_delayed_rectifier_K_I (milliS_per_microF).
* STATES[32] is O_K in component rapid_delayed_rectifier_K_I (dimensionless).
* STATES[33] is C_K1 in component rapid_delayed_rectifier_K_I (dimensionless).
* STATES[34] is C_K2 in component rapid_delayed_rectifier_K_I (dimensionless).
* ALGEBRAIC[10] is C_K0 in component rapid_delayed_rectifier_K_I (dimensionless).
* STATES[35] is I_K in component rapid_delayed_rectifier_K_I (dimensionless).
* ALGEBRAIC[24] is alpha_a0 in component rapid_delayed_rectifier_K_I (per_millisecond).
* ALGEBRAIC[29] is beta_a0 in component rapid_delayed_rectifier_K_I (per_millisecond).
* CONSTANTS[66] is kb in component rapid_delayed_rectifier_K_I (per_millisecond).
* CONSTANTS[67] is kf in component rapid_delayed_rectifier_K_I (per_millisecond).
* ALGEBRAIC[11] is alpha_a1 in component rapid_delayed_rectifier_K_I (per_millisecond).
* ALGEBRAIC[25] is beta_a1 in component rapid_delayed_rectifier_K_I (per_millisecond).
* ALGEBRAIC[30] is alpha_i in component rapid_delayed_rectifier_K_I (per_millisecond).
* ALGEBRAIC[34] is beta_i in component rapid_delayed_rectifier_K_I (per_millisecond).
* CONSTANTS[68] is i_NaK_max in component sodium_potassium_pump_current (picoA_per_picoF).
* CONSTANTS[69] is Km_Nai in component sodium_potassium_pump_current (micromolar).
* CONSTANTS[70] is Km_Ko in component sodium_potassium_pump_current (micromolar).
* ALGEBRAIC[77] is f_NaK in component sodium_potassium_pump_current (dimensionless).
* CONSTANTS[75] is sigma in component sodium_potassium_pump_current (dimensionless).
* CONSTANTS[71] is g_ClCa in component calcium_activated_chloride_current (milliS_per_microF).
* ALGEBRAIC[79] is O_ClCa in component calcium_activated_chloride_current (dimensionless).
* CONSTANTS[72] is E_Cl in component calcium_activated_chloride_current (millivolt).
* CONSTANTS[73] is Km_Cl in component calcium_activated_chloride_current (micromolar).
* 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[6] is d/dt CaMKt in component calcium_fluxes (dimensionless).
* RATES[5] is d/dt P_RyR in component calcium_fluxes (dimensionless).
* RATES[7] is d/dt P_O1 in component ryanodine_receptors (dimensionless).
* RATES[8] is d/dt P_O2 in component ryanodine_receptors (dimensionless).
* RATES[9] is d/dt P_C2 in component ryanodine_receptors (dimensionless).
* RATES[10] is d/dt O in component L_type_calcium_current (dimensionless).
* RATES[11] is d/dt I in component L_type_calcium_current (dimensionless).
* RATES[12] is d/dt y_gate in component L_type_calcium_current (dimensionless).
* RATES[13] is d/dt Nai in component sodium_concentration (micromolar).
* RATES[16] is d/dt C_Na2 in component fast_sodium_current (dimensionless).
* RATES[15] is d/dt C_Na1 in component fast_sodium_current (dimensionless).
* RATES[14] is d/dt O_Na in component fast_sodium_current (dimensionless).
* RATES[19] is d/dt IF_Na in component fast_sodium_current (dimensionless).
* RATES[17] is d/dt I1_Na in component fast_sodium_current (dimensionless).
* RATES[18] is d/dt I2_Na in component fast_sodium_current (dimensionless).
* RATES[20] is d/dt IC_Na2 in component fast_sodium_current (dimensionless).
* RATES[21] is d/dt IC_Na3 in component fast_sodium_current (dimensionless).
* RATES[22] is d/dt Ki in component potassium_concentration (micromolar).
* RATES[23] is d/dt ato_f in component fast_transient_outward_K_I (dimensionless).
* RATES[24] is d/dt ito_f in component fast_transient_outward_K_I (dimensionless).
* RATES[25] is d/dt ato_s in component slow_transient_outward_K_I (dimensionless).
* RATES[26] is d/dt ito_s in component slow_transient_outward_K_I (dimensionless).
* RATES[27] is d/dt nKs in component slow_delayed_rectifier_K_I (dimensionless).
* RATES[28] is d/dt aur in component ultra_rapidly_activating_delayed_rectifier_K_I (dimensionless).
* RATES[29] is d/dt iur in component ultra_rapidly_activating_delayed_rectifier_K_I (dimensionless).
* RATES[30] is d/dt aKss in component non_inactivating_steady_state_K_I (dimensionless).
* RATES[31] is d/dt iKss in component non_inactivating_steady_state_K_I (dimensionless).
* RATES[34] is d/dt C_K2 in component rapid_delayed_rectifier_K_I (dimensionless).
* RATES[33] is d/dt C_K1 in component rapid_delayed_rectifier_K_I (dimensionless).
* RATES[32] is d/dt O_K in component rapid_delayed_rectifier_K_I (dimensionless).
* RATES[35] is d/dt I_K in component rapid_delayed_rectifier_K_I (dimensionless).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -85.64004;
CONSTANTS[0] = 1;
CONSTANTS[1] = 2.2e-5;
CONSTANTS[2] = 1.022e-7;
CONSTANTS[3] = 1.786e-6;
CONSTANTS[4] = 1.264e-9;
CONSTANTS[5] = 0.00013866;
CONSTANTS[6] = 5400;
CONSTANTS[7] = 134000;
CONSTANTS[8] = 1400;
CONSTANTS[9] = 8.314;
CONSTANTS[10] = 308;
CONSTANTS[11] = 96.5;
CONSTANTS[12] = 100;
CONSTANTS[13] = 333.3333333;
CONSTANTS[14] = 3;
CONSTANTS[15] = -15;
STATES[1] = 0.1040595;
STATES[2] = 0.1043777;
STATES[3] = 730.0589;
STATES[4] = 841.106;
CONSTANTS[16] = 109;
CONSTANTS[17] = 15000;
CONSTANTS[18] = 0.6;
CONSTANTS[19] = 800;
STATES[5] = 2.290355e-9;
CONSTANTS[20] = 4.5;
CONSTANTS[21] = 20;
CONSTANTS[22] = 1e-5;
CONSTANTS[23] = 8;
STATES[6] = 0.08989079;
CONSTANTS[24] = 0.05;
CONSTANTS[25] = 0.00068;
CONSTANTS[26] = 0.412;
CONSTANTS[27] = 4;
STATES[7] = 0.003825599;
STATES[8] = 1.835831e-8;
CONSTANTS[28] = 0.3158;
CONSTANTS[29] = -0.09;
CONSTANTS[30] = -0.225;
STATES[9] = 0.3797679;
CONSTANTS[31] = 0.006075;
CONSTANTS[32] = 0.07125;
CONSTANTS[33] = 0.00405;
CONSTANTS[34] = 0.965;
CONSTANTS[35] = 0.009;
CONSTANTS[36] = 0.0008;
CONSTANTS[37] = 3;
CONSTANTS[38] = 4;
CONSTANTS[39] = 19.1078;
STATES[10] = 4.373318e-6;
STATES[11] = 0.009171979;
STATES[12] = 0.8876797;
CONSTANTS[40] = -2;
CONSTANTS[41] = 7.0671;
CONSTANTS[42] = 1.1683;
CONSTANTS[43] = 1.6411;
CONSTANTS[44] = 0.07;
CONSTANTS[45] = 14;
CONSTANTS[46] = 972.9715;
CONSTANTS[47] = 0.0964;
CONSTANTS[48] = 6.6755;
CONSTANTS[49] = 0.9;
CONSTANTS[50] = 0.4;
CONSTANTS[51] = 772.8991;
CONSTANTS[52] = 86500;
CONSTANTS[53] = 1380;
CONSTANTS[54] = 0.1;
CONSTANTS[55] = 0.35;
STATES[13] = 16522.45;
CONSTANTS[56] = 0.00088;
CONSTANTS[57] = 13;
STATES[14] = 2.639399e-7;
STATES[15] = 0.0001581035;
STATES[16] = 0.01702105;
STATES[17] = 0.00001799179;
STATES[18] = 0.000005460299;
STATES[19] = 0.0000556206;
STATES[20] = 0.005985434;
STATES[21] = 0.2543133;
STATES[22] = 141474;
CONSTANTS[58] = 0.0026;
CONSTANTS[59] = 0.4067;
STATES[23] = 0.001937245;
STATES[24] = 0.9999985;
CONSTANTS[60] = 0;
STATES[25] = 0.02000568;
STATES[26] = 0.9308568;
CONSTANTS[61] = 0.2938;
CONSTANTS[62] = 0.00575;
STATES[27] = 0.002206261;
CONSTANTS[63] = 0.16;
STATES[28] = 0.02000568;
STATES[29] = 0.9822006;
CONSTANTS[64] = 0.05;
STATES[30] = 0.8883113;
STATES[31] = 1;
CONSTANTS[65] = 0.078;
STATES[32] = 0.0004858865;
STATES[33] = 0.0007799137;
STATES[34] = 0.0005301217;
STATES[35] = 0.00007519518;
CONSTANTS[66] = 0.036778;
CONSTANTS[67] = 0.023761;
CONSTANTS[68] = 1.66;
CONSTANTS[69] = 21000;
CONSTANTS[70] = 1500;
CONSTANTS[71] = 10;
CONSTANTS[72] = -40;
CONSTANTS[73] = 10;
CONSTANTS[74] = CONSTANTS[43]/CONSTANTS[42];
CONSTANTS[75] = (1.00000/7.00000)*(exp(CONSTANTS[7]/67300.0) - 1.00000);
CONSTANTS[76] = 0.00000;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[31] = CONSTANTS[76];
RATES[8] = CONSTANTS[33]*pow(STATES[2], CONSTANTS[37])*STATES[7] - CONSTANTS[34]*STATES[8];
RATES[9] = CONSTANTS[35]*STATES[7] - CONSTANTS[36]*STATES[9];
ALGEBRAIC[1] = 1.00000 - (STATES[9]+STATES[7]+STATES[8]);
RATES[7] = ( CONSTANTS[31]*pow(STATES[2], CONSTANTS[38])*ALGEBRAIC[1]+ CONSTANTS[34]*STATES[8]+ CONSTANTS[36]*STATES[9]) - ( CONSTANTS[32]*STATES[7]+ CONSTANTS[33]*pow(STATES[2], CONSTANTS[37])*STATES[7]+ CONSTANTS[35]*STATES[7]);
ALGEBRAIC[2] = 1.00000/(1.00000+exp((STATES[0]+16.6577)/CONSTANTS[48]))+0.100000/(1.00000+exp((- STATES[0]+40.0000)/6.00000));
ALGEBRAIC[13] = 20.0000+600.000/(1.00000+exp((STATES[0]+30.0000)/9.60000));
RATES[12] = (ALGEBRAIC[2] - STATES[12])/ALGEBRAIC[13];
ALGEBRAIC[5] = 0.180640*exp( 0.0357700*(STATES[0]+30.0000));
ALGEBRAIC[16] = 0.395600*exp( - 0.0623700*(STATES[0]+30.0000));
RATES[23] = ALGEBRAIC[5]*(1.00000 - STATES[23]) - ALGEBRAIC[16]*STATES[23];
ALGEBRAIC[6] = ( 0.000152000*exp(- (STATES[0]+13.5000)/7.00000))/( 0.00670830*exp(- (STATES[0]+33.5000)/7.00000)+1.00000);
ALGEBRAIC[17] = ( 0.000950000*exp((STATES[0]+33.5000)/7.00000))/( 0.0513350*exp((STATES[0]+33.5000)/7.00000)+1.00000);
RATES[24] = ALGEBRAIC[6]*(1.00000 - STATES[24]) - ALGEBRAIC[17]*STATES[24];
ALGEBRAIC[7] = 1.00000/(1.00000+exp(- (STATES[0]+22.5000)/7.70000));
ALGEBRAIC[18] = 0.493000*exp( - 0.0629000*STATES[0])+2.05800;
RATES[25] = (ALGEBRAIC[7] - STATES[25])/ALGEBRAIC[18];
ALGEBRAIC[8] = 1.00000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[19] = 270.000+1050.00/(1.00000+exp((STATES[0]+45.2000)/5.70000));
RATES[26] = (ALGEBRAIC[8] - STATES[26])/ALGEBRAIC[19];
ALGEBRAIC[9] = ( 4.81333e-06*(STATES[0]+26.5000))/(1.00000 - exp( - 0.128000*(STATES[0]+26.5000)));
ALGEBRAIC[20] = 9.53333e-05*exp( - 0.0380000*(STATES[0]+26.5000));
RATES[27] = ALGEBRAIC[9]*(1.00000 - STATES[27]) - ALGEBRAIC[20]*STATES[27];
ALGEBRAIC[21] = 0.493000*exp( - 0.0629000*STATES[0])+2.05800;
RATES[28] = (ALGEBRAIC[7] - STATES[28])/ALGEBRAIC[21];
ALGEBRAIC[22] = 1200.00 - 170.000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
RATES[29] = (ALGEBRAIC[8] - STATES[29])/ALGEBRAIC[22];
ALGEBRAIC[23] = 39.3000*exp( - 0.0862000*STATES[0])+13.1700;
RATES[30] = (ALGEBRAIC[7] - STATES[30])/ALGEBRAIC[23];
ALGEBRAIC[11] = 0.0137330*exp( 0.0381980*STATES[0]);
ALGEBRAIC[25] = 6.89000e-05*exp( - 0.0417800*STATES[0]);
RATES[34] = ( CONSTANTS[67]*STATES[33]+ ALGEBRAIC[25]*STATES[32]) - ( CONSTANTS[66]*STATES[34]+ ALGEBRAIC[11]*STATES[34]);
ALGEBRAIC[10] = 1.00000 - (STATES[33]+STATES[34]+STATES[32]+STATES[35]);
ALGEBRAIC[24] = 0.0223480*exp( 0.0117600*STATES[0]);
ALGEBRAIC[29] = 0.0470020*exp( - 0.0631000*STATES[0]);
RATES[33] = ( ALGEBRAIC[24]*ALGEBRAIC[10]+ CONSTANTS[66]*STATES[34]) - ( ALGEBRAIC[29]*STATES[33]+ CONSTANTS[67]*STATES[33]);
ALGEBRAIC[30] = 0.0908210*exp( 0.0233910*(STATES[0]+5.00000));
ALGEBRAIC[34] = 0.00649700*exp( - 0.0326800*(STATES[0]+5.00000));
RATES[32] = ( ALGEBRAIC[11]*STATES[34]+ ALGEBRAIC[34]*STATES[35]) - ( ALGEBRAIC[25]*STATES[32]+ ALGEBRAIC[30]*STATES[32]);
RATES[35] = ALGEBRAIC[30]*STATES[32] - ALGEBRAIC[34]*STATES[35];
ALGEBRAIC[36] = (1.00000 - STATES[10]) - STATES[11];
ALGEBRAIC[3] = exp((STATES[0] - CONSTANTS[40])/CONSTANTS[41]);
ALGEBRAIC[14] = ALGEBRAIC[3]/( CONSTANTS[42]*(ALGEBRAIC[3]+1.00000));
RATES[10] = ALGEBRAIC[14]*ALGEBRAIC[36] - CONSTANTS[74]*STATES[10];
ALGEBRAIC[27] = (ALGEBRAIC[3]+CONSTANTS[44])/( CONSTANTS[46]*CONSTANTS[47]*(ALGEBRAIC[3]+1.00000));
ALGEBRAIC[32] = ( CONSTANTS[45]*(ALGEBRAIC[3]+CONSTANTS[44]))/( CONSTANTS[46]*( CONSTANTS[45]*ALGEBRAIC[3]+CONSTANTS[44]));
RATES[11] = ALGEBRAIC[27]*STATES[2]*ALGEBRAIC[36] - ALGEBRAIC[32]*STATES[11];
ALGEBRAIC[35] = pow(1.00000+( CONSTANTS[17]*CONSTANTS[19])/pow(CONSTANTS[19]+STATES[3], 2.00000), - 1.00000);
ALGEBRAIC[38] = CONSTANTS[20]*(STATES[7]+STATES[8])*(STATES[3] - STATES[2])*STATES[5];
ALGEBRAIC[40] = (STATES[4] - STATES[3])/CONSTANTS[21];
RATES[3] = ALGEBRAIC[35]*(ALGEBRAIC[40] - ALGEBRAIC[38]);
ALGEBRAIC[4] = 1.00000 - (STATES[14]+STATES[15]+STATES[16]+STATES[19]+STATES[17]+STATES[18]+STATES[20]+STATES[21]);
ALGEBRAIC[15] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/17.0000)+ 0.200000*exp(- (STATES[0]+2.50000)/150.000));
ALGEBRAIC[37] = 0.191700*exp(- (STATES[0]+2.50000)/20.3000);
ALGEBRAIC[28] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/15.0000)+ 0.230000*exp(- (STATES[0]+2.50000)/150.000));
ALGEBRAIC[39] = 0.200000*exp(- (STATES[0] - 2.50000)/20.3000);
ALGEBRAIC[43] = 7.00000e-07*exp(- (STATES[0]+7.00000)/7.70000);
ALGEBRAIC[45] = 0.00840000+ 2.00000e-05*(STATES[0]+7.00000);
RATES[16] = ( ALGEBRAIC[15]*ALGEBRAIC[4]+ ALGEBRAIC[39]*STATES[15]+ ALGEBRAIC[43]*STATES[20]) - ( ALGEBRAIC[37]*STATES[16]+ ALGEBRAIC[28]*STATES[16]+ ALGEBRAIC[45]*STATES[16]);
ALGEBRAIC[33] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/12.0000)+ 0.250000*exp(- (STATES[0]+2.50000)/150.000));
ALGEBRAIC[41] = 0.220000*exp(- (STATES[0] - 7.50000)/20.3000);
RATES[15] = ( ALGEBRAIC[28]*STATES[16]+ ALGEBRAIC[41]*STATES[14]+ ALGEBRAIC[43]*STATES[19]) - ( ALGEBRAIC[39]*STATES[15]+ ALGEBRAIC[33]*STATES[15]+ ALGEBRAIC[45]*STATES[15]);
RATES[20] = ( ALGEBRAIC[15]*STATES[21]+ ALGEBRAIC[39]*STATES[19]+ ALGEBRAIC[45]*STATES[20]) - ( ALGEBRAIC[37]*STATES[20]+ ALGEBRAIC[28]*STATES[20]+ ALGEBRAIC[43]*STATES[20]);
RATES[21] = ( ALGEBRAIC[37]*STATES[20]+ ALGEBRAIC[45]*ALGEBRAIC[4]) - ( ALGEBRAIC[15]*STATES[21]+ ALGEBRAIC[43]*STATES[21]);
ALGEBRAIC[46] = ( 0.0500000*(1.00000 - STATES[6])*1.00000)/(1.00000+0.700000/STATES[2]);
RATES[6] = CONSTANTS[24]*ALGEBRAIC[46]*(ALGEBRAIC[46]+STATES[6]) - CONSTANTS[25]*STATES[6];
ALGEBRAIC[47] = 1.00000/( 0.188495*exp(- (STATES[0]+7.00000)/16.6000)+0.393956);
ALGEBRAIC[49] = ( ALGEBRAIC[33]*ALGEBRAIC[47]*ALGEBRAIC[43])/( ALGEBRAIC[41]*ALGEBRAIC[45]);
RATES[14] = ( ALGEBRAIC[33]*STATES[15]+ ALGEBRAIC[49]*STATES[19]) - ( ALGEBRAIC[41]*STATES[14]+ ALGEBRAIC[47]*STATES[14]);
ALGEBRAIC[44] = CONSTANTS[22]*(STATES[4] - STATES[1]);
ALGEBRAIC[48] = ALGEBRAIC[46]+STATES[6];
ALGEBRAIC[50] = (( 3.15120*pow(ALGEBRAIC[48], 2.20620))/(pow(0.158800, 2.20620)+pow(ALGEBRAIC[48], 2.20620))+1.00000)*CONSTANTS[28];
ALGEBRAIC[52] = ( ALGEBRAIC[50]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[26], 2.00000)+pow(STATES[1], 2.00000));
RATES[4] = ( (ALGEBRAIC[52] - ALGEBRAIC[44])*CONSTANTS[1])/CONSTANTS[3] - ( ALGEBRAIC[40]*CONSTANTS[2])/CONSTANTS[3];
ALGEBRAIC[51] = ALGEBRAIC[47]/1000.00;
ALGEBRAIC[53] = ALGEBRAIC[43];
RATES[19] = ( ALGEBRAIC[47]*STATES[14]+ ALGEBRAIC[45]*STATES[15]+ ALGEBRAIC[53]*STATES[17]+ ALGEBRAIC[28]*STATES[20]) - ( ALGEBRAIC[49]*STATES[19]+ ALGEBRAIC[43]*STATES[19]+ ALGEBRAIC[51]*STATES[19]+ ALGEBRAIC[39]*STATES[19]);
ALGEBRAIC[55] = ALGEBRAIC[47]/95000.0;
ALGEBRAIC[57] = ALGEBRAIC[43]/50.0000;
RATES[17] = ( ALGEBRAIC[51]*STATES[19]+ ALGEBRAIC[57]*STATES[18]) - ( ALGEBRAIC[53]*STATES[17]+ ALGEBRAIC[55]*STATES[17]);
RATES[18] = ALGEBRAIC[55]*STATES[17] - ALGEBRAIC[57]*STATES[18];
ALGEBRAIC[54] = ( CONSTANTS[11]*STATES[0])/( CONSTANTS[9]*CONSTANTS[10]);
ALGEBRAIC[56] = 2.00000*ALGEBRAIC[54];
ALGEBRAIC[58] = (fabs(ALGEBRAIC[56])>1.00000e-05 ? (( (( - CONSTANTS[39]*2.00000*CONSTANTS[4]*CONSTANTS[11])/( CONSTANTS[5]*CONSTANTS[0]))*STATES[10]*STATES[12]*ALGEBRAIC[56])/(1.00000 - exp(- ALGEBRAIC[56])))*( CONSTANTS[8]*exp(- ALGEBRAIC[56]) - STATES[2]) : (( (( - CONSTANTS[39]*2.00000*CONSTANTS[4]*CONSTANTS[11])/( CONSTANTS[5]*CONSTANTS[0]))*STATES[10]*STATES[12]*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*( CONSTANTS[8]*exp(- 1.00000e-05) - STATES[2]));
ALGEBRAIC[31] = pow(1.00000+( CONSTANTS[16]*CONSTANTS[18])/pow(CONSTANTS[18]+STATES[2], 2.00000), - 1.00000);
ALGEBRAIC[42] = (STATES[2] - STATES[1])/CONSTANTS[23];
RATES[2] = ALGEBRAIC[31]*(( ALGEBRAIC[38]*CONSTANTS[2])/CONSTANTS[4] - (( ALGEBRAIC[42]*CONSTANTS[1])/CONSTANTS[4]+( ALGEBRAIC[58]*CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[4]*CONSTANTS[11])));
RATES[5] = CONSTANTS[29]*STATES[5]+ (( CONSTANTS[30]*ALGEBRAIC[58])/CONSTANTS[27])*exp(- pow(STATES[0] - 5.00000, 2.00000)/648.000);
ALGEBRAIC[59] = ( CONSTANTS[49]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[50], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[60] = (( (( (( CONSTANTS[51]*1.00000)/(pow(CONSTANTS[52], 3.00000)+pow(CONSTANTS[7], 3.00000)))*1.00000)/(CONSTANTS[53]+CONSTANTS[8]))*1.00000)/(1.00000+ CONSTANTS[54]*exp(( (CONSTANTS[55] - 1.00000)*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))))*( exp(( CONSTANTS[55]*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))*pow(STATES[13], 3.00000)*CONSTANTS[8] - exp(( (CONSTANTS[55] - 1.00000)*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))*pow(CONSTANTS[7], 3.00000)*STATES[1]);
ALGEBRAIC[62] = (( CONSTANTS[9]*CONSTANTS[10])/( 2.00000*CONSTANTS[11]))*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[64] = CONSTANTS[56]*(STATES[0] - ALGEBRAIC[62]);
ALGEBRAIC[26] = pow(1.00000+( CONSTANTS[16]*CONSTANTS[18])/pow(CONSTANTS[18]+STATES[1], 2.00000), - 1.00000);
RATES[1] = ALGEBRAIC[26]*(((ALGEBRAIC[44]+ALGEBRAIC[42]) - ( ((ALGEBRAIC[64]+ALGEBRAIC[59]) - 2.00000*ALGEBRAIC[60])*CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[1]*CONSTANTS[11])) - ALGEBRAIC[52]);
ALGEBRAIC[65] = (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(( 0.900000*CONSTANTS[7]+ 0.100000*CONSTANTS[6])/( 0.900000*STATES[13]+ 0.100000*STATES[22]));
ALGEBRAIC[67] = CONSTANTS[57]*STATES[14]*(STATES[0] - ALGEBRAIC[65]);
ALGEBRAIC[68] = CONSTANTS[58]*(STATES[0] - ALGEBRAIC[65]);
ALGEBRAIC[77] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))+ 0.0365000*CONSTANTS[75]*exp(( - STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10])));
ALGEBRAIC[78] = ( (( CONSTANTS[68]*ALGEBRAIC[77]*1.00000)/(1.00000+pow(CONSTANTS[69]/STATES[13], 1.50000)))*CONSTANTS[6])/(CONSTANTS[6]+CONSTANTS[70]);
RATES[13] = ( - (ALGEBRAIC[67]+ALGEBRAIC[68]+ 3.00000*ALGEBRAIC[78]+ 3.00000*ALGEBRAIC[60])*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[11]);
ALGEBRAIC[69] = (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(CONSTANTS[6]/STATES[22]);
ALGEBRAIC[70] = CONSTANTS[59]*pow(STATES[23], 3.00000)*STATES[24]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[71] = CONSTANTS[60]*STATES[25]*STATES[26]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[72] = ( (( CONSTANTS[61]*CONSTANTS[6])/(CONSTANTS[6]+210.000))*(STATES[0] - ALGEBRAIC[69]))/(1.00000+exp( 0.0896000*(STATES[0] - ALGEBRAIC[69])));
ALGEBRAIC[73] = CONSTANTS[62]*pow(STATES[27], 2.00000)*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[74] = CONSTANTS[63]*STATES[28]*STATES[29]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[75] = CONSTANTS[64]*STATES[30]*STATES[31]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[76] = CONSTANTS[65]*STATES[32]*(STATES[0] - (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(( 0.980000*CONSTANTS[6]+ 0.0200000*CONSTANTS[7])/( 0.980000*STATES[22]+ 0.0200000*STATES[13])));
ALGEBRAIC[0] = floor(VOI/CONSTANTS[13])*CONSTANTS[13];
ALGEBRAIC[12] = (VOI - ALGEBRAIC[0]>=CONSTANTS[12]&&VOI - ALGEBRAIC[0]<=CONSTANTS[12]+CONSTANTS[14] ? CONSTANTS[15] : 0.00000);
RATES[22] = ( - ((ALGEBRAIC[12]+ALGEBRAIC[70]+ALGEBRAIC[71]+ALGEBRAIC[72]+ALGEBRAIC[73]+ALGEBRAIC[75]+ALGEBRAIC[74]+ALGEBRAIC[76]) - 2.00000*ALGEBRAIC[78])*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[11]);
ALGEBRAIC[79] = 0.200000/(1.00000+exp(- (STATES[0] - 46.7000)/7.80000));
ALGEBRAIC[80] = (( CONSTANTS[71]*ALGEBRAIC[79]*STATES[1])/(STATES[1]+CONSTANTS[73]))*(STATES[0] - CONSTANTS[72]);
RATES[0] = - (ALGEBRAIC[58]+ALGEBRAIC[59]+ALGEBRAIC[60]+ALGEBRAIC[64]+ALGEBRAIC[67]+ALGEBRAIC[68]+ALGEBRAIC[78]+ALGEBRAIC[70]+ALGEBRAIC[71]+ALGEBRAIC[72]+ALGEBRAIC[73]+ALGEBRAIC[74]+ALGEBRAIC[75]+ALGEBRAIC[76]+ALGEBRAIC[80]+ALGEBRAIC[12]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] = 1.00000 - (STATES[9]+STATES[7]+STATES[8]);
ALGEBRAIC[2] = 1.00000/(1.00000+exp((STATES[0]+16.6577)/CONSTANTS[48]))+0.100000/(1.00000+exp((- STATES[0]+40.0000)/6.00000));
ALGEBRAIC[13] = 20.0000+600.000/(1.00000+exp((STATES[0]+30.0000)/9.60000));
ALGEBRAIC[5] = 0.180640*exp( 0.0357700*(STATES[0]+30.0000));
ALGEBRAIC[16] = 0.395600*exp( - 0.0623700*(STATES[0]+30.0000));
ALGEBRAIC[6] = ( 0.000152000*exp(- (STATES[0]+13.5000)/7.00000))/( 0.00670830*exp(- (STATES[0]+33.5000)/7.00000)+1.00000);
ALGEBRAIC[17] = ( 0.000950000*exp((STATES[0]+33.5000)/7.00000))/( 0.0513350*exp((STATES[0]+33.5000)/7.00000)+1.00000);
ALGEBRAIC[7] = 1.00000/(1.00000+exp(- (STATES[0]+22.5000)/7.70000));
ALGEBRAIC[18] = 0.493000*exp( - 0.0629000*STATES[0])+2.05800;
ALGEBRAIC[8] = 1.00000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[19] = 270.000+1050.00/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[9] = ( 4.81333e-06*(STATES[0]+26.5000))/(1.00000 - exp( - 0.128000*(STATES[0]+26.5000)));
ALGEBRAIC[20] = 9.53333e-05*exp( - 0.0380000*(STATES[0]+26.5000));
ALGEBRAIC[21] = 0.493000*exp( - 0.0629000*STATES[0])+2.05800;
ALGEBRAIC[22] = 1200.00 - 170.000/(1.00000+exp((STATES[0]+45.2000)/5.70000));
ALGEBRAIC[23] = 39.3000*exp( - 0.0862000*STATES[0])+13.1700;
ALGEBRAIC[11] = 0.0137330*exp( 0.0381980*STATES[0]);
ALGEBRAIC[25] = 6.89000e-05*exp( - 0.0417800*STATES[0]);
ALGEBRAIC[10] = 1.00000 - (STATES[33]+STATES[34]+STATES[32]+STATES[35]);
ALGEBRAIC[24] = 0.0223480*exp( 0.0117600*STATES[0]);
ALGEBRAIC[29] = 0.0470020*exp( - 0.0631000*STATES[0]);
ALGEBRAIC[30] = 0.0908210*exp( 0.0233910*(STATES[0]+5.00000));
ALGEBRAIC[34] = 0.00649700*exp( - 0.0326800*(STATES[0]+5.00000));
ALGEBRAIC[36] = (1.00000 - STATES[10]) - STATES[11];
ALGEBRAIC[3] = exp((STATES[0] - CONSTANTS[40])/CONSTANTS[41]);
ALGEBRAIC[14] = ALGEBRAIC[3]/( CONSTANTS[42]*(ALGEBRAIC[3]+1.00000));
ALGEBRAIC[27] = (ALGEBRAIC[3]+CONSTANTS[44])/( CONSTANTS[46]*CONSTANTS[47]*(ALGEBRAIC[3]+1.00000));
ALGEBRAIC[32] = ( CONSTANTS[45]*(ALGEBRAIC[3]+CONSTANTS[44]))/( CONSTANTS[46]*( CONSTANTS[45]*ALGEBRAIC[3]+CONSTANTS[44]));
ALGEBRAIC[35] = pow(1.00000+( CONSTANTS[17]*CONSTANTS[19])/pow(CONSTANTS[19]+STATES[3], 2.00000), - 1.00000);
ALGEBRAIC[38] = CONSTANTS[20]*(STATES[7]+STATES[8])*(STATES[3] - STATES[2])*STATES[5];
ALGEBRAIC[40] = (STATES[4] - STATES[3])/CONSTANTS[21];
ALGEBRAIC[4] = 1.00000 - (STATES[14]+STATES[15]+STATES[16]+STATES[19]+STATES[17]+STATES[18]+STATES[20]+STATES[21]);
ALGEBRAIC[15] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/17.0000)+ 0.200000*exp(- (STATES[0]+2.50000)/150.000));
ALGEBRAIC[37] = 0.191700*exp(- (STATES[0]+2.50000)/20.3000);
ALGEBRAIC[28] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/15.0000)+ 0.230000*exp(- (STATES[0]+2.50000)/150.000));
ALGEBRAIC[39] = 0.200000*exp(- (STATES[0] - 2.50000)/20.3000);
ALGEBRAIC[43] = 7.00000e-07*exp(- (STATES[0]+7.00000)/7.70000);
ALGEBRAIC[45] = 0.00840000+ 2.00000e-05*(STATES[0]+7.00000);
ALGEBRAIC[33] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/12.0000)+ 0.250000*exp(- (STATES[0]+2.50000)/150.000));
ALGEBRAIC[41] = 0.220000*exp(- (STATES[0] - 7.50000)/20.3000);
ALGEBRAIC[46] = ( 0.0500000*(1.00000 - STATES[6])*1.00000)/(1.00000+0.700000/STATES[2]);
ALGEBRAIC[47] = 1.00000/( 0.188495*exp(- (STATES[0]+7.00000)/16.6000)+0.393956);
ALGEBRAIC[49] = ( ALGEBRAIC[33]*ALGEBRAIC[47]*ALGEBRAIC[43])/( ALGEBRAIC[41]*ALGEBRAIC[45]);
ALGEBRAIC[44] = CONSTANTS[22]*(STATES[4] - STATES[1]);
ALGEBRAIC[48] = ALGEBRAIC[46]+STATES[6];
ALGEBRAIC[50] = (( 3.15120*pow(ALGEBRAIC[48], 2.20620))/(pow(0.158800, 2.20620)+pow(ALGEBRAIC[48], 2.20620))+1.00000)*CONSTANTS[28];
ALGEBRAIC[52] = ( ALGEBRAIC[50]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[26], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[51] = ALGEBRAIC[47]/1000.00;
ALGEBRAIC[53] = ALGEBRAIC[43];
ALGEBRAIC[55] = ALGEBRAIC[47]/95000.0;
ALGEBRAIC[57] = ALGEBRAIC[43]/50.0000;
ALGEBRAIC[54] = ( CONSTANTS[11]*STATES[0])/( CONSTANTS[9]*CONSTANTS[10]);
ALGEBRAIC[56] = 2.00000*ALGEBRAIC[54];
ALGEBRAIC[58] = (fabs(ALGEBRAIC[56])>1.00000e-05 ? (( (( - CONSTANTS[39]*2.00000*CONSTANTS[4]*CONSTANTS[11])/( CONSTANTS[5]*CONSTANTS[0]))*STATES[10]*STATES[12]*ALGEBRAIC[56])/(1.00000 - exp(- ALGEBRAIC[56])))*( CONSTANTS[8]*exp(- ALGEBRAIC[56]) - STATES[2]) : (( (( - CONSTANTS[39]*2.00000*CONSTANTS[4]*CONSTANTS[11])/( CONSTANTS[5]*CONSTANTS[0]))*STATES[10]*STATES[12]*1.00000e-05)/(1.00000 - exp(- 1.00000e-05)))*( CONSTANTS[8]*exp(- 1.00000e-05) - STATES[2]));
ALGEBRAIC[31] = pow(1.00000+( CONSTANTS[16]*CONSTANTS[18])/pow(CONSTANTS[18]+STATES[2], 2.00000), - 1.00000);
ALGEBRAIC[42] = (STATES[2] - STATES[1])/CONSTANTS[23];
ALGEBRAIC[59] = ( CONSTANTS[49]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[50], 2.00000)+pow(STATES[1], 2.00000));
ALGEBRAIC[60] = (( (( (( CONSTANTS[51]*1.00000)/(pow(CONSTANTS[52], 3.00000)+pow(CONSTANTS[7], 3.00000)))*1.00000)/(CONSTANTS[53]+CONSTANTS[8]))*1.00000)/(1.00000+ CONSTANTS[54]*exp(( (CONSTANTS[55] - 1.00000)*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))))*( exp(( CONSTANTS[55]*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))*pow(STATES[13], 3.00000)*CONSTANTS[8] - exp(( (CONSTANTS[55] - 1.00000)*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))*pow(CONSTANTS[7], 3.00000)*STATES[1]);
ALGEBRAIC[62] = (( CONSTANTS[9]*CONSTANTS[10])/( 2.00000*CONSTANTS[11]))*log(CONSTANTS[8]/STATES[1]);
ALGEBRAIC[64] = CONSTANTS[56]*(STATES[0] - ALGEBRAIC[62]);
ALGEBRAIC[26] = pow(1.00000+( CONSTANTS[16]*CONSTANTS[18])/pow(CONSTANTS[18]+STATES[1], 2.00000), - 1.00000);
ALGEBRAIC[65] = (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(( 0.900000*CONSTANTS[7]+ 0.100000*CONSTANTS[6])/( 0.900000*STATES[13]+ 0.100000*STATES[22]));
ALGEBRAIC[67] = CONSTANTS[57]*STATES[14]*(STATES[0] - ALGEBRAIC[65]);
ALGEBRAIC[68] = CONSTANTS[58]*(STATES[0] - ALGEBRAIC[65]);
ALGEBRAIC[77] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))+ 0.0365000*CONSTANTS[75]*exp(( - STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10])));
ALGEBRAIC[78] = ( (( CONSTANTS[68]*ALGEBRAIC[77]*1.00000)/(1.00000+pow(CONSTANTS[69]/STATES[13], 1.50000)))*CONSTANTS[6])/(CONSTANTS[6]+CONSTANTS[70]);
ALGEBRAIC[69] = (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(CONSTANTS[6]/STATES[22]);
ALGEBRAIC[70] = CONSTANTS[59]*pow(STATES[23], 3.00000)*STATES[24]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[71] = CONSTANTS[60]*STATES[25]*STATES[26]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[72] = ( (( CONSTANTS[61]*CONSTANTS[6])/(CONSTANTS[6]+210.000))*(STATES[0] - ALGEBRAIC[69]))/(1.00000+exp( 0.0896000*(STATES[0] - ALGEBRAIC[69])));
ALGEBRAIC[73] = CONSTANTS[62]*pow(STATES[27], 2.00000)*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[74] = CONSTANTS[63]*STATES[28]*STATES[29]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[75] = CONSTANTS[64]*STATES[30]*STATES[31]*(STATES[0] - ALGEBRAIC[69]);
ALGEBRAIC[76] = CONSTANTS[65]*STATES[32]*(STATES[0] - (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(( 0.980000*CONSTANTS[6]+ 0.0200000*CONSTANTS[7])/( 0.980000*STATES[22]+ 0.0200000*STATES[13])));
ALGEBRAIC[0] = floor(VOI/CONSTANTS[13])*CONSTANTS[13];
ALGEBRAIC[12] = (VOI - ALGEBRAIC[0]>=CONSTANTS[12]&&VOI - ALGEBRAIC[0]<=CONSTANTS[12]+CONSTANTS[14] ? CONSTANTS[15] : 0.00000);
ALGEBRAIC[79] = 0.200000/(1.00000+exp(- (STATES[0] - 46.7000)/7.80000));
ALGEBRAIC[80] = (( CONSTANTS[71]*ALGEBRAIC[79]*STATES[1])/(STATES[1]+CONSTANTS[73]))*(STATES[0] - CONSTANTS[72]);
ALGEBRAIC[61] = ( - ALGEBRAIC[59]*CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[1]*CONSTANTS[11]);
ALGEBRAIC[63] = ( ALGEBRAIC[60]*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[11]);
ALGEBRAIC[66] = ( - ALGEBRAIC[64]*CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[1]*CONSTANTS[11]);
}