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 75 entries in the algebraic variable array.
   There are a total of 27 entries in each of the rate and state variable arrays.
   There are a total of 66 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (millisecond).
 * STATES[0] is V in component membrane (millivolt).
 * CONSTANTS[0] is R in component membrane (attojoule_per_millimole_kelvin).
 * CONSTANTS[1] is T in component membrane (kelvin).
 * CONSTANTS[2] is F in component membrane (femtocoulomb_per_millimole).
 * CONSTANTS[3] is Cm in component membrane (picoF).
 * ALGEBRAIC[39] is i_CaT in component T_type_calcium_channel_current (picoA).
 * ALGEBRAIC[34] is i_CaL in component L_type_calcium_channel_current (picoA).
 * ALGEBRAIC[40] is i_Kr in component rapidly_activating_delayed_rectifier_potassium_current (picoA).
 * ALGEBRAIC[42] is i_Ks in component slowly_activating_delayed_rectifier_potassium_current (picoA).
 * ALGEBRAIC[43] is i_to in component AP_sensitive_currents (picoA).
 * ALGEBRAIC[44] is i_sus in component AP_sensitive_currents (picoA).
 * ALGEBRAIC[47] is i_h in component hyperpolarisation_activated_current (picoA).
 * ALGEBRAIC[48] is i_st in component sustained_inward_current (picoA).
 * ALGEBRAIC[49] is i_b_Na in component sodium_dependent_background_current (picoA).
 * ALGEBRAIC[65] is i_NaCa in component sodium_calcium_exchange_current (picoA).
 * ALGEBRAIC[51] is i_NaK in component sodium_potassium_pump_current (picoA).
 * ALGEBRAIC[50] is i_K_ACh in component background_muscarinic_potassium_channel_current (picoA).
 * ALGEBRAIC[0] is E_Na in component reversal_potentials (millivolt).
 * ALGEBRAIC[29] is E_K in component reversal_potentials (millivolt).
 * STATES[1] is Nai in component intracellular_ion_concentrations (millimolar).
 * CONSTANTS[4] is Nao in component intracellular_ion_concentrations (millimolar).
 * STATES[2] is Ki in component intracellular_ion_concentrations (millimolar).
 * CONSTANTS[5] is Ko in component intracellular_ion_concentrations (millimolar).
 * CONSTANTS[6] is g_CaL in component L_type_calcium_channel_current (nanoS).
 * CONSTANTS[7] is E_CaL in component L_type_calcium_channel_current (millivolt).
 * STATES[3] is Ca_sub in component intracellular_ion_concentrations (millimolar).
 * STATES[4] is d in component L_type_calcium_channel_current_d_gate (dimensionless).
 * STATES[5] is f in component L_type_calcium_channel_current_f_gate (dimensionless).
 * STATES[6] is fCa in component L_type_calcium_channel_current_fCa_gate (dimensionless).
 * ALGEBRAIC[1] is d_infinity in component L_type_calcium_channel_current_d_gate (dimensionless).
 * ALGEBRAIC[35] is tau_d in component L_type_calcium_channel_current_d_gate (millisecond).
 * ALGEBRAIC[15] is alpha_d in component L_type_calcium_channel_current_d_gate (per_millisecond).
 * ALGEBRAIC[30] is beta_d in component L_type_calcium_channel_current_d_gate (per_millisecond).
 * ALGEBRAIC[2] is f_infinity in component L_type_calcium_channel_current_f_gate (dimensionless).
 * ALGEBRAIC[16] is tau_f in component L_type_calcium_channel_current_f_gate (millisecond).
 * CONSTANTS[62] is alpha_fCa in component L_type_calcium_channel_current_fCa_gate (per_millisecond).
 * CONSTANTS[8] is beta_fCa in component L_type_calcium_channel_current_fCa_gate (per_millimolar_millisecond).
 * ALGEBRAIC[3] is fCa_infinity in component L_type_calcium_channel_current_fCa_gate (dimensionless).
 * ALGEBRAIC[17] is tau_fCa in component L_type_calcium_channel_current_fCa_gate (millisecond).
 * CONSTANTS[9] is Km_fCa in component L_type_calcium_channel_current_fCa_gate (millimolar).
 * CONSTANTS[10] is g_CaT in component T_type_calcium_channel_current (nanoS).
 * CONSTANTS[11] is E_CaT in component T_type_calcium_channel_current (millivolt).
 * STATES[7] is d in component T_type_calcium_channel_current_d_gate (dimensionless).
 * STATES[8] is f in component T_type_calcium_channel_current_f_gate (dimensionless).
 * ALGEBRAIC[4] is d_infinity in component T_type_calcium_channel_current_d_gate (dimensionless).
 * ALGEBRAIC[18] is tau_d in component T_type_calcium_channel_current_d_gate (millisecond).
 * ALGEBRAIC[5] is f_infinity in component T_type_calcium_channel_current_f_gate (dimensionless).
 * ALGEBRAIC[19] is tau_f in component T_type_calcium_channel_current_f_gate (millisecond).
 * CONSTANTS[63] is g_Kr in component rapidly_activating_delayed_rectifier_potassium_current (nanoS).
 * STATES[9] is paS in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
 * STATES[10] is paF in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
 * STATES[11] is piy in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (dimensionless).
 * ALGEBRAIC[6] is pa_infinity in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
 * ALGEBRAIC[20] is tau_paS in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (millisecond).
 * ALGEBRAIC[21] is tau_paF in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (millisecond).
 * ALGEBRAIC[7] is pi_infinity in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (dimensionless).
 * ALGEBRAIC[22] is tau_pi in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (millisecond).
 * CONSTANTS[12] is g_Ks in component slowly_activating_delayed_rectifier_potassium_current (nanoS).
 * ALGEBRAIC[41] is E_Ks in component slowly_activating_delayed_rectifier_potassium_current (millivolt).
 * STATES[12] is n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (dimensionless).
 * ALGEBRAIC[31] is n_infinity in component slowly_activating_delayed_rectifier_potassium_current_n_gate (dimensionless).
 * ALGEBRAIC[36] is tau_n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (millisecond).
 * ALGEBRAIC[8] is alpha_n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (per_millisecond).
 * ALGEBRAIC[23] is beta_n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (per_millisecond).
 * CONSTANTS[13] is g_to in component AP_sensitive_currents (nanoS).
 * CONSTANTS[14] is g_sus in component AP_sensitive_currents (nanoS).
 * STATES[13] is q in component AP_sensitive_currents_q_gate (dimensionless).
 * STATES[14] is r in component AP_sensitive_currents_r_gate (dimensionless).
 * ALGEBRAIC[9] is q_infinity in component AP_sensitive_currents_q_gate (dimensionless).
 * ALGEBRAIC[24] is tau_q in component AP_sensitive_currents_q_gate (millisecond).
 * ALGEBRAIC[10] is r_infinity in component AP_sensitive_currents_r_gate (dimensionless).
 * ALGEBRAIC[25] is tau_r in component AP_sensitive_currents_r_gate (millisecond).
 * ALGEBRAIC[45] is i_h_Na in component hyperpolarisation_activated_current (picoA).
 * ALGEBRAIC[46] is i_h_K in component hyperpolarisation_activated_current (picoA).
 * CONSTANTS[15] is g_h_Na in component hyperpolarisation_activated_current (nanoS).
 * CONSTANTS[16] is g_h_K in component hyperpolarisation_activated_current (nanoS).
 * STATES[15] is y in component hyperpolarisation_activated_current_y_gate (dimensionless).
 * ALGEBRAIC[11] is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless).
 * ALGEBRAIC[26] is tau_y in component hyperpolarisation_activated_current_y_gate (millisecond).
 * CONSTANTS[17] is g_st in component sustained_inward_current (nanoS).
 * CONSTANTS[18] is E_st in component sustained_inward_current (millivolt).
 * STATES[16] is qa in component sustained_inward_current_qa_gate (dimensionless).
 * STATES[17] is qi in component sustained_inward_current_qi_gate (dimensionless).
 * ALGEBRAIC[12] is qa_infinity in component sustained_inward_current_qa_gate (dimensionless).
 * ALGEBRAIC[37] is tau_qa in component sustained_inward_current_qa_gate (millisecond).
 * ALGEBRAIC[27] is alpha_qa in component sustained_inward_current_qa_gate (per_millisecond).
 * ALGEBRAIC[32] is beta_qa in component sustained_inward_current_qa_gate (per_millisecond).
 * ALGEBRAIC[33] is qi_infinity in component sustained_inward_current_qi_gate (dimensionless).
 * ALGEBRAIC[38] is tau_qi in component sustained_inward_current_qi_gate (millisecond).
 * ALGEBRAIC[13] is alpha_qi in component sustained_inward_current_qi_gate (per_millisecond).
 * ALGEBRAIC[28] is beta_qi in component sustained_inward_current_qi_gate (per_millisecond).
 * CONSTANTS[19] is g_b_Na in component sodium_dependent_background_current (nanoS).
 * CONSTANTS[64] is g_K_ACh in component background_muscarinic_potassium_channel_current (picoA).
 * CONSTANTS[20] is Km_Kp in component sodium_potassium_pump_current (millimolar).
 * CONSTANTS[21] is Km_Nap in component sodium_potassium_pump_current (millimolar).
 * CONSTANTS[22] is i_NaK_max in component sodium_potassium_pump_current (picoA).
 * CONSTANTS[23] is kNaCa in component sodium_calcium_exchange_current (picoA).
 * ALGEBRAIC[62] is x1 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[58] is x2 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[63] is x3 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[64] is x4 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[53] is k41 in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[65] is k34 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[61] is k23 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[60] is k21 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[57] is k32 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[52] is k43 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[55] is k12 in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[56] is k14 in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[24] is Qci in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[25] is Qn in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[26] is Qco in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[27] is K3ni in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[28] is Kci in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[29] is K1ni in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[30] is K2ni in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[31] is Kcni in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[32] is K3no in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[33] is K1no in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[34] is K2no in component sodium_calcium_exchange_current (millimolar).
 * CONSTANTS[35] is Kco in component sodium_calcium_exchange_current (millimolar).
 * ALGEBRAIC[59] is do in component sodium_calcium_exchange_current (dimensionless).
 * ALGEBRAIC[54] is di in component sodium_calcium_exchange_current (dimensionless).
 * CONSTANTS[36] is Cao in component intracellular_ion_concentrations (millimolar).
 * ALGEBRAIC[66] is j_Ca_dif in component intracellular_calcium_dynamics (millimolar_per_millisecond).
 * ALGEBRAIC[67] is j_rel in component intracellular_calcium_dynamics (millimolar_per_millisecond).
 * ALGEBRAIC[68] is j_up in component intracellular_calcium_dynamics (millimolar_per_millisecond).
 * ALGEBRAIC[70] is j_tr in component intracellular_calcium_dynamics (millimolar_per_millisecond).
 * CONSTANTS[37] is tau_dif_Ca in component intracellular_calcium_dynamics (millisecond).
 * CONSTANTS[38] is tau_tr in component intracellular_calcium_dynamics (millisecond).
 * CONSTANTS[39] is K_rel in component intracellular_calcium_dynamics (millimolar).
 * CONSTANTS[40] is P_up in component intracellular_calcium_dynamics (millimolar_per_millisecond).
 * CONSTANTS[41] is P_rel in component intracellular_calcium_dynamics (per_millisecond).
 * CONSTANTS[42] is K_up in component intracellular_calcium_dynamics (millimolar).
 * STATES[18] is Ca_up in component intracellular_ion_concentrations (millimolar).
 * STATES[19] is Cai in component intracellular_ion_concentrations (millimolar).
 * STATES[20] is Ca_rel in component intracellular_ion_concentrations (millimolar).
 * CONSTANTS[43] is Mgi in component intracellular_ion_concentrations (millimolar).
 * CONSTANTS[44] is V_i in component intracellular_ion_concentrations (litre).
 * CONSTANTS[45] is V_rel in component intracellular_ion_concentrations (litre).
 * CONSTANTS[46] is V_up in component intracellular_ion_concentrations (litre).
 * CONSTANTS[47] is V_sub in component intracellular_ion_concentrations (litre).
 * CONSTANTS[48] is TMC_tot in component calcium_buffering (dimensionless).
 * CONSTANTS[49] is CM_tot in component calcium_buffering (dimensionless).
 * CONSTANTS[50] is TC_tot in component calcium_buffering (dimensionless).
 * CONSTANTS[51] is CQ_tot in component calcium_buffering (dimensionless).
 * ALGEBRAIC[72] is delta_fTMC in component calcium_buffering (millimolar_per_millisecond).
 * ALGEBRAIC[74] is delta_fCMi in component calcium_buffering (millimolar_per_millisecond).
 * ALGEBRAIC[69] is delta_fCMs in component calcium_buffering (millimolar_per_millisecond).
 * ALGEBRAIC[71] is delta_fTC in component calcium_buffering (millimolar_per_millisecond).
 * ALGEBRAIC[73] is delta_fCQ in component calcium_buffering (millimolar_per_millisecond).
 * ALGEBRAIC[14] is delta_fTMM in component calcium_buffering (millimolar_per_millisecond).
 * STATES[21] is fTMM in component calcium_buffering (millimolar).
 * STATES[22] is fCMi in component calcium_buffering (millimolar).
 * STATES[23] is fCMs in component calcium_buffering (millimolar).
 * STATES[24] is fTC in component calcium_buffering (millimolar).
 * STATES[25] is fTMC in component calcium_buffering (millimolar).
 * STATES[26] is fCQ in component calcium_buffering (millimolar).
 * CONSTANTS[52] is kf_TC in component calcium_buffering (per_millimolar_millisecond).
 * CONSTANTS[53] is kf_TMM in component calcium_buffering (per_millimolar_millisecond).
 * CONSTANTS[54] is kf_TMC in component calcium_buffering (per_millimolar_millisecond).
 * CONSTANTS[55] is kf_CM in component calcium_buffering (per_millimolar_millisecond).
 * CONSTANTS[56] is kf_CQ in component calcium_buffering (per_millimolar_millisecond).
 * CONSTANTS[57] is kb_TC in component calcium_buffering (per_millisecond).
 * CONSTANTS[58] is kb_TMC in component calcium_buffering (per_millisecond).
 * CONSTANTS[59] is kb_TMM in component calcium_buffering (per_millisecond).
 * CONSTANTS[60] is kb_CM in component calcium_buffering (per_millisecond).
 * CONSTANTS[61] is kb_CQ in component calcium_buffering (per_millisecond).
 * RATES[0] is d/dt V in component membrane (millivolt).
 * RATES[4] is d/dt d in component L_type_calcium_channel_current_d_gate (dimensionless).
 * RATES[5] is d/dt f in component L_type_calcium_channel_current_f_gate (dimensionless).
 * RATES[6] is d/dt fCa in component L_type_calcium_channel_current_fCa_gate (dimensionless).
 * RATES[7] is d/dt d in component T_type_calcium_channel_current_d_gate (dimensionless).
 * RATES[8] is d/dt f in component T_type_calcium_channel_current_f_gate (dimensionless).
 * RATES[9] is d/dt paS in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
 * RATES[10] is d/dt paF in component rapidly_activating_delayed_rectifier_potassium_current_pa_gate (dimensionless).
 * RATES[11] is d/dt piy in component rapidly_activating_delayed_rectifier_potassium_current_pi_gate (dimensionless).
 * RATES[12] is d/dt n in component slowly_activating_delayed_rectifier_potassium_current_n_gate (dimensionless).
 * RATES[13] is d/dt q in component AP_sensitive_currents_q_gate (dimensionless).
 * RATES[14] is d/dt r in component AP_sensitive_currents_r_gate (dimensionless).
 * RATES[15] is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless).
 * RATES[16] is d/dt qa in component sustained_inward_current_qa_gate (dimensionless).
 * RATES[17] is d/dt qi in component sustained_inward_current_qi_gate (dimensionless).
 * RATES[1] is d/dt Nai in component intracellular_ion_concentrations (millimolar).
 * RATES[2] is d/dt Ki in component intracellular_ion_concentrations (millimolar).
 * RATES[19] is d/dt Cai in component intracellular_ion_concentrations (millimolar).
 * RATES[3] is d/dt Ca_sub in component intracellular_ion_concentrations (millimolar).
 * RATES[18] is d/dt Ca_up in component intracellular_ion_concentrations (millimolar).
 * RATES[20] is d/dt Ca_rel in component intracellular_ion_concentrations (millimolar).
 * RATES[24] is d/dt fTC in component calcium_buffering (millimolar).
 * RATES[25] is d/dt fTMC in component calcium_buffering (millimolar).
 * RATES[21] is d/dt fTMM in component calcium_buffering (millimolar).
 * RATES[22] is d/dt fCMi in component calcium_buffering (millimolar).
 * RATES[23] is d/dt fCMs in component calcium_buffering (millimolar).
 * RATES[26] is d/dt fCQ in component calcium_buffering (millimolar).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -58.600291137693;
CONSTANTS[0] = 8314400000000000;
CONSTANTS[1] = 310.15;
CONSTANTS[2] = 96485000000000000;
CONSTANTS[3] = 32;
STATES[1] = 9.438646305915;
CONSTANTS[4] = 140;
STATES[2] = 139.984146485614;
CONSTANTS[5] = 5.4;
CONSTANTS[6] = 0.58;
CONSTANTS[7] = 45;
STATES[3] = 0.00019074741;
STATES[4] = 0.000602055134;
STATES[5] = 0.626999773853;
STATES[6] = 0.589580408056;
CONSTANTS[8] = 60;
CONSTANTS[9] = 0.00035;
CONSTANTS[10] = 0.458;
CONSTANTS[11] = 45;
STATES[7] = 0.004571884917;
STATES[8] = 0.249637570396;
STATES[9] = 0.629323128348;
STATES[10] = 0.3493633709533;
STATES[11] = 0.852396631172;
CONSTANTS[12] = 0.0259;
STATES[12] = 0.054409723782;
CONSTANTS[13] = 0.18;
CONSTANTS[14] = 0.02;
STATES[13] = 0.531446952485;
STATES[14] = 0.005550489445;
CONSTANTS[15] = 0.1437375;
CONSTANTS[16] = 0.2312625;
STATES[15] = 0.067156687129;
CONSTANTS[17] = 0.015;
CONSTANTS[18] = 37.4;
STATES[16] = 0.426018100136;
STATES[17] = 0.333330378068;
CONSTANTS[19] = 0.0054;
CONSTANTS[20] = 1.4;
CONSTANTS[21] = 14;
CONSTANTS[22] = 3.6;
CONSTANTS[23] = 125;
CONSTANTS[24] = 0.1369;
CONSTANTS[25] = 0.4315;
CONSTANTS[26] = 0;
CONSTANTS[27] = 26.44;
CONSTANTS[28] = 0.0207;
CONSTANTS[29] = 395.3;
CONSTANTS[30] = 2.289;
CONSTANTS[31] = 26.44;
CONSTANTS[32] = 4.663;
CONSTANTS[33] = 1628;
CONSTANTS[34] = 561.4;
CONSTANTS[35] = 3.663;
CONSTANTS[36] = 2;
CONSTANTS[37] = 0.04;
CONSTANTS[38] = 60;
CONSTANTS[39] = 0.0012;
CONSTANTS[40] = 0.005;
CONSTANTS[41] = 0.5;
CONSTANTS[42] = 0.0006;
STATES[18] = 1.462338380106;
STATES[19] = 0.000312494921;
STATES[20] = 0.296742023718;
CONSTANTS[43] = 2.5;
CONSTANTS[44] = 0.0000000000015835;
CONSTANTS[45] = 0.0000000000000042223;
CONSTANTS[46] = 0.000000000000040816;
CONSTANTS[47] = 0.000000000000035098;
CONSTANTS[48] = 0.062;
CONSTANTS[49] = 0.045;
CONSTANTS[50] = 0.031;
CONSTANTS[51] = 10;
STATES[21] = 0.350600895635;
STATES[22] = 0.116947220413;
STATES[23] = 0.074631965653;
STATES[24] = 0.059206293446;
STATES[25] = 0.602955114871;
STATES[26] = 0.260317260703;
CONSTANTS[52] = 88.8;
CONSTANTS[53] = 2.277;
CONSTANTS[54] = 227.7;
CONSTANTS[55] = 227.7;
CONSTANTS[56] = 0.534;
CONSTANTS[57] = 0.446;
CONSTANTS[58] = 0.00751;
CONSTANTS[59] = 0.751;
CONSTANTS[60] = 0.542;
CONSTANTS[61] = 0.445;
CONSTANTS[62] =  CONSTANTS[9]*CONSTANTS[8];
CONSTANTS[63] =  0.0250000*pow(CONSTANTS[5]/1.00000, 0.590000);
CONSTANTS[64] =  0.00110000*pow(CONSTANTS[5]/1.00000, 0.410000);
CONSTANTS[65] = CONSTANTS[4]/(CONSTANTS[32]+CONSTANTS[4]);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[14] =  CONSTANTS[53]*CONSTANTS[43]*(1.00000 - (STATES[25]+STATES[21])) -  CONSTANTS[59]*STATES[21];
RATES[21] = ALGEBRAIC[14];
ALGEBRAIC[2] = 1.00000/(1.00000+exp((STATES[0]+30.0000)/5.00000));
ALGEBRAIC[16] = 44.3000+ 257.100*exp(- pow((STATES[0]+32.5000)/13.9000, 2.00000));
RATES[5] = (ALGEBRAIC[2] - STATES[5])/ALGEBRAIC[16];
ALGEBRAIC[3] = CONSTANTS[9]/(CONSTANTS[9]+STATES[3]);
ALGEBRAIC[17] = ALGEBRAIC[3]/CONSTANTS[62];
RATES[6] = (ALGEBRAIC[3] - STATES[6])/ALGEBRAIC[17];
ALGEBRAIC[4] = 1.00000/(1.00000+exp(- (STATES[0]+26.3000)/6.00000));
ALGEBRAIC[18] = 1.00000/( 1.06800*exp((STATES[0]+26.3000)/30.0000)+ 1.06800*exp(- (STATES[0]+26.3000)/30.0000));
RATES[7] = (ALGEBRAIC[4] - STATES[7])/ALGEBRAIC[18];
ALGEBRAIC[5] = 1.00000/(1.00000+exp((STATES[0]+61.7000)/5.60000));
ALGEBRAIC[19] = 1.00000/( 0.0153000*exp(- (STATES[0]+61.7000)/83.3000)+ 0.0150000*exp((STATES[0]+61.7000)/15.3800));
RATES[8] = (ALGEBRAIC[5] - STATES[8])/ALGEBRAIC[19];
ALGEBRAIC[6] = 1.00000/(1.00000+exp(- (STATES[0]+23.2000)/10.6000));
ALGEBRAIC[20] = 0.846554/( 0.00420000*exp(STATES[0]/17.0000)+ 0.000150000*exp(- STATES[0]/21.6000));
RATES[9] = (ALGEBRAIC[6] - STATES[9])/ALGEBRAIC[20];
ALGEBRAIC[21] = 0.846554/( 0.0372000*exp(STATES[0]/15.9000)+ 0.000960000*exp(- STATES[0]/22.5000));
RATES[10] = (ALGEBRAIC[6] - STATES[10])/ALGEBRAIC[21];
ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[0]+28.6000)/17.1000));
ALGEBRAIC[22] = 1.00000/( 0.100000*exp(- STATES[0]/54.6450)+ 0.656000*exp(STATES[0]/106.157));
RATES[11] = (ALGEBRAIC[7] - STATES[11])/ALGEBRAIC[22];
ALGEBRAIC[9] = 1.00000/(1.00000+exp((STATES[0]+49.0000)/13.0000));
ALGEBRAIC[24] =  0.600000*(65.1700/( 0.570000*exp( - 0.0800000*(STATES[0]+44.0000))+ 0.0650000*exp( 0.100000*(STATES[0]+45.9300)))+10.1000);
RATES[13] = (ALGEBRAIC[9] - STATES[13])/ALGEBRAIC[24];
ALGEBRAIC[10] = 1.00000/(1.00000+exp(- (STATES[0] - 19.3000)/15.0000));
ALGEBRAIC[25] =  0.660000*1.40000*(15.5900/( 1.03700*exp( 0.0900000*(STATES[0]+30.6100))+ 0.369000*exp( - 0.120000*(STATES[0]+23.8400)))+2.98000);
RATES[14] = (ALGEBRAIC[10] - STATES[14])/ALGEBRAIC[25];
ALGEBRAIC[11] = 1.00000/(1.00000+exp((STATES[0]+64.0000)/13.5000));
ALGEBRAIC[26] = 0.716653/(exp(- (STATES[0]+386.900)/45.3020)+exp((STATES[0] - 73.0800)/19.2310));
RATES[15] = (ALGEBRAIC[11] - STATES[15])/ALGEBRAIC[26];
ALGEBRAIC[1] = 1.00000/(1.00000+exp(- (STATES[0]+14.1000)/6.00000));
ALGEBRAIC[15] = ( - 0.0283900*(STATES[0]+35.0000))/(exp(- (STATES[0]+35.0000)/2.50000) - 1.00000) - ( 0.0849000*STATES[0])/(exp(- STATES[0]/4.80800) - 1.00000);
ALGEBRAIC[30] = ( 0.0114300*(STATES[0] - 5.00000))/(exp((STATES[0] - 5.00000)/2.50000) - 1.00000);
ALGEBRAIC[35] = 1.00000/(ALGEBRAIC[15]+ALGEBRAIC[30]);
RATES[4] = (ALGEBRAIC[1] - STATES[4])/ALGEBRAIC[35];
ALGEBRAIC[8] = 0.0140000/(1.00000+exp(- (STATES[0] - 40.0000)/9.00000));
ALGEBRAIC[23] =  0.00100000*exp(- STATES[0]/45.0000);
ALGEBRAIC[31] = ALGEBRAIC[8]/(ALGEBRAIC[8]+ALGEBRAIC[23]);
ALGEBRAIC[36] = 1.00000/(ALGEBRAIC[8]+ALGEBRAIC[23]);
RATES[12] = (ALGEBRAIC[31] - STATES[12])/ALGEBRAIC[36];
ALGEBRAIC[12] = 1.00000/(1.00000+exp(- (STATES[0]+57.0000)/5.00000));
ALGEBRAIC[27] = 1.00000/( 0.150000*exp(- STATES[0]/11.0000)+ 0.200000*exp(- STATES[0]/700.000));
ALGEBRAIC[32] = 1.00000/( 16.0000*exp(STATES[0]/8.00000)+ 15.0000*exp(STATES[0]/50.0000));
ALGEBRAIC[37] = 1.00000/(ALGEBRAIC[27]+ALGEBRAIC[32]);
RATES[16] = (ALGEBRAIC[12] - STATES[16])/ALGEBRAIC[37];
ALGEBRAIC[13] = 1.00000/( 3100.00*exp(STATES[0]/13.0000)+ 700.000*exp(STATES[0]/70.0000));
ALGEBRAIC[28] = 1.00000/( 95.0000*exp(- STATES[0]/10.0000)+ 50.0000*exp(- STATES[0]/700.000))+0.000229000/(1.00000+exp(- STATES[0]/5.00000));
ALGEBRAIC[33] = ALGEBRAIC[13]/(ALGEBRAIC[13]+ALGEBRAIC[28]);
ALGEBRAIC[38] = 6.65000/(ALGEBRAIC[13]+ALGEBRAIC[28]);
RATES[17] = (ALGEBRAIC[33] - STATES[17])/ALGEBRAIC[38];
ALGEBRAIC[29] =  (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[5]/STATES[2]);
ALGEBRAIC[40] =  CONSTANTS[63]*(STATES[0] - ALGEBRAIC[29])*( 0.600000*STATES[10]+ 0.400000*STATES[9])*STATES[11];
ALGEBRAIC[41] =  (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log((CONSTANTS[5]+ 0.120000*CONSTANTS[4])/(STATES[2]+ 0.120000*STATES[1]));
ALGEBRAIC[42] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[41])*pow(STATES[12], 2.00000);
ALGEBRAIC[43] =  CONSTANTS[13]*(STATES[0] - ALGEBRAIC[29])*STATES[13]*STATES[14];
ALGEBRAIC[44] =  CONSTANTS[14]*(STATES[0] - ALGEBRAIC[29])*STATES[14];
ALGEBRAIC[0] =  (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[4]/STATES[1]);
ALGEBRAIC[51] =  CONSTANTS[22]*pow(1.00000+pow(CONSTANTS[20]/CONSTANTS[5], 1.20000), - 1.00000)*pow(1.00000+pow(CONSTANTS[21]/STATES[1], 1.30000), - 1.00000)*pow(1.00000+exp(- ((STATES[0] - ALGEBRAIC[0])+120.000)/30.0000), - 1.00000);
ALGEBRAIC[50] = ( CONSTANTS[64]*(STATES[2] -  CONSTANTS[5]*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))/1.00000;
ALGEBRAIC[46] =  CONSTANTS[16]*(STATES[0] - ALGEBRAIC[29])*pow(STATES[15], 2.00000);
RATES[2] = ( - (ALGEBRAIC[40]+ALGEBRAIC[42]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[46]+ALGEBRAIC[50]+ - 2.00000*ALGEBRAIC[51])*CONSTANTS[3])/( 1.00000*CONSTANTS[2]*(CONSTANTS[44]+CONSTANTS[47]));
ALGEBRAIC[39] =  CONSTANTS[10]*(STATES[0] - CONSTANTS[11])*STATES[7]*STATES[8];
ALGEBRAIC[34] =  CONSTANTS[6]*(STATES[0] - CONSTANTS[7])*STATES[4]*STATES[5]*STATES[6];
ALGEBRAIC[45] =  CONSTANTS[15]*(STATES[0] - ALGEBRAIC[0])*pow(STATES[15], 2.00000);
ALGEBRAIC[47] = ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[48] =  CONSTANTS[17]*(STATES[0] - CONSTANTS[18])*STATES[16]*STATES[17];
ALGEBRAIC[49] =  CONSTANTS[19]*(STATES[0] - ALGEBRAIC[0]);
ALGEBRAIC[53] = exp(( - CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1]));
ALGEBRAIC[59] = 1.00000+ (CONSTANTS[36]/CONSTANTS[35])*(1.00000+exp(( CONSTANTS[26]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])))+ (CONSTANTS[4]/CONSTANTS[33])*(1.00000+CONSTANTS[4]/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[32]);
ALGEBRAIC[61] = ( (( (CONSTANTS[4]/CONSTANTS[33])*CONSTANTS[4])/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[32])*exp(( - CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[59];
ALGEBRAIC[60] = ( (CONSTANTS[36]/CONSTANTS[35])*exp(( CONSTANTS[26]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[59];
ALGEBRAIC[57] = exp(( CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1]));
ALGEBRAIC[52] = STATES[1]/(CONSTANTS[27]+STATES[1]);
ALGEBRAIC[62] =  ALGEBRAIC[53]*CONSTANTS[65]*(ALGEBRAIC[61]+ALGEBRAIC[60])+ ALGEBRAIC[60]*ALGEBRAIC[57]*(ALGEBRAIC[52]+ALGEBRAIC[53]);
ALGEBRAIC[54] = 1.00000+ (STATES[3]/CONSTANTS[28])*(1.00000+exp(( - CONSTANTS[24]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+STATES[1]/CONSTANTS[31])+ (STATES[1]/CONSTANTS[29])*(1.00000+ (STATES[1]/CONSTANTS[30])*(1.00000+STATES[1]/CONSTANTS[27]));
ALGEBRAIC[55] = ( (STATES[3]/CONSTANTS[28])*exp(( - CONSTANTS[24]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[54];
ALGEBRAIC[56] = ( (( (STATES[1]/CONSTANTS[29])*STATES[1])/CONSTANTS[30])*(1.00000+STATES[1]/CONSTANTS[27])*exp(( CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[54];
ALGEBRAIC[58] =  ALGEBRAIC[57]*ALGEBRAIC[52]*(ALGEBRAIC[56]+ALGEBRAIC[55])+ ALGEBRAIC[53]*ALGEBRAIC[55]*(CONSTANTS[65]+ALGEBRAIC[57]);
ALGEBRAIC[63] =  ALGEBRAIC[56]*ALGEBRAIC[52]*(ALGEBRAIC[61]+ALGEBRAIC[60])+ ALGEBRAIC[55]*ALGEBRAIC[61]*(ALGEBRAIC[52]+ALGEBRAIC[53]);
ALGEBRAIC[64] =  ALGEBRAIC[61]*CONSTANTS[65]*(ALGEBRAIC[56]+ALGEBRAIC[55])+ ALGEBRAIC[56]*ALGEBRAIC[60]*(CONSTANTS[65]+ALGEBRAIC[57]);
ALGEBRAIC[65] = ( CONSTANTS[23]*( ALGEBRAIC[58]*ALGEBRAIC[60] -  ALGEBRAIC[62]*ALGEBRAIC[55]))/(ALGEBRAIC[62]+ALGEBRAIC[58]+ALGEBRAIC[63]+ALGEBRAIC[64]);
RATES[0] = - (ALGEBRAIC[34]+ALGEBRAIC[39]+ALGEBRAIC[40]+ALGEBRAIC[42]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[47]+ALGEBRAIC[48]+ALGEBRAIC[49]+ALGEBRAIC[50]+ALGEBRAIC[51]+ALGEBRAIC[65])/1.00000;
RATES[1] = ( - (ALGEBRAIC[45]+ALGEBRAIC[48]+ALGEBRAIC[49]+ 3.00000*ALGEBRAIC[51]+ 3.00000*ALGEBRAIC[65])*CONSTANTS[3])/( 1.00000*CONSTANTS[2]*(CONSTANTS[44]+CONSTANTS[47]));
ALGEBRAIC[66] = (STATES[3] - STATES[19])/CONSTANTS[37];
ALGEBRAIC[67] = ( CONSTANTS[41]*(STATES[20] - STATES[3]))/(1.00000+pow(CONSTANTS[39]/STATES[3], 2.00000));
ALGEBRAIC[69] =  CONSTANTS[55]*STATES[3]*(1.00000 - STATES[23]) -  CONSTANTS[60]*STATES[23];
RATES[3] = (( - ((ALGEBRAIC[34]+ALGEBRAIC[39]) -  2.00000*ALGEBRAIC[65])*CONSTANTS[3])/( 1.00000*2.00000*CONSTANTS[2])+ ALGEBRAIC[67]*CONSTANTS[45])/CONSTANTS[47] - (ALGEBRAIC[66]+ CONSTANTS[49]*ALGEBRAIC[69]);
RATES[23] = ALGEBRAIC[69];
ALGEBRAIC[68] = CONSTANTS[40]/(1.00000+CONSTANTS[42]/STATES[19]);
ALGEBRAIC[70] = (STATES[18] - STATES[20])/CONSTANTS[38];
RATES[18] = ALGEBRAIC[68] - ( ALGEBRAIC[70]*CONSTANTS[45])/CONSTANTS[46];
ALGEBRAIC[71] =  CONSTANTS[52]*STATES[19]*(1.00000 - STATES[24]) -  CONSTANTS[57]*STATES[24];
RATES[24] = ALGEBRAIC[71];
ALGEBRAIC[73] =  CONSTANTS[56]*STATES[20]*(1.00000 - STATES[26]) -  CONSTANTS[61]*STATES[26];
RATES[20] = ALGEBRAIC[70] - (ALGEBRAIC[67]+ CONSTANTS[51]*ALGEBRAIC[73]);
ALGEBRAIC[72] =  CONSTANTS[54]*STATES[19]*(1.00000 - (STATES[25]+STATES[21])) -  CONSTANTS[58]*STATES[25];
RATES[25] = ALGEBRAIC[72];
RATES[26] = ALGEBRAIC[73];
ALGEBRAIC[74] =  CONSTANTS[55]*STATES[19]*(1.00000 - STATES[22]) -  CONSTANTS[60]*STATES[22];
RATES[19] = ( ALGEBRAIC[66]*CONSTANTS[47] -  ALGEBRAIC[68]*CONSTANTS[46])/CONSTANTS[44] - ( CONSTANTS[49]*ALGEBRAIC[74]+ CONSTANTS[50]*ALGEBRAIC[71]+ CONSTANTS[48]*ALGEBRAIC[72]);
RATES[22] = ALGEBRAIC[74];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[14] =  CONSTANTS[53]*CONSTANTS[43]*(1.00000 - (STATES[25]+STATES[21])) -  CONSTANTS[59]*STATES[21];
ALGEBRAIC[2] = 1.00000/(1.00000+exp((STATES[0]+30.0000)/5.00000));
ALGEBRAIC[16] = 44.3000+ 257.100*exp(- pow((STATES[0]+32.5000)/13.9000, 2.00000));
ALGEBRAIC[3] = CONSTANTS[9]/(CONSTANTS[9]+STATES[3]);
ALGEBRAIC[17] = ALGEBRAIC[3]/CONSTANTS[62];
ALGEBRAIC[4] = 1.00000/(1.00000+exp(- (STATES[0]+26.3000)/6.00000));
ALGEBRAIC[18] = 1.00000/( 1.06800*exp((STATES[0]+26.3000)/30.0000)+ 1.06800*exp(- (STATES[0]+26.3000)/30.0000));
ALGEBRAIC[5] = 1.00000/(1.00000+exp((STATES[0]+61.7000)/5.60000));
ALGEBRAIC[19] = 1.00000/( 0.0153000*exp(- (STATES[0]+61.7000)/83.3000)+ 0.0150000*exp((STATES[0]+61.7000)/15.3800));
ALGEBRAIC[6] = 1.00000/(1.00000+exp(- (STATES[0]+23.2000)/10.6000));
ALGEBRAIC[20] = 0.846554/( 0.00420000*exp(STATES[0]/17.0000)+ 0.000150000*exp(- STATES[0]/21.6000));
ALGEBRAIC[21] = 0.846554/( 0.0372000*exp(STATES[0]/15.9000)+ 0.000960000*exp(- STATES[0]/22.5000));
ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[0]+28.6000)/17.1000));
ALGEBRAIC[22] = 1.00000/( 0.100000*exp(- STATES[0]/54.6450)+ 0.656000*exp(STATES[0]/106.157));
ALGEBRAIC[9] = 1.00000/(1.00000+exp((STATES[0]+49.0000)/13.0000));
ALGEBRAIC[24] =  0.600000*(65.1700/( 0.570000*exp( - 0.0800000*(STATES[0]+44.0000))+ 0.0650000*exp( 0.100000*(STATES[0]+45.9300)))+10.1000);
ALGEBRAIC[10] = 1.00000/(1.00000+exp(- (STATES[0] - 19.3000)/15.0000));
ALGEBRAIC[25] =  0.660000*1.40000*(15.5900/( 1.03700*exp( 0.0900000*(STATES[0]+30.6100))+ 0.369000*exp( - 0.120000*(STATES[0]+23.8400)))+2.98000);
ALGEBRAIC[11] = 1.00000/(1.00000+exp((STATES[0]+64.0000)/13.5000));
ALGEBRAIC[26] = 0.716653/(exp(- (STATES[0]+386.900)/45.3020)+exp((STATES[0] - 73.0800)/19.2310));
ALGEBRAIC[1] = 1.00000/(1.00000+exp(- (STATES[0]+14.1000)/6.00000));
ALGEBRAIC[15] = ( - 0.0283900*(STATES[0]+35.0000))/(exp(- (STATES[0]+35.0000)/2.50000) - 1.00000) - ( 0.0849000*STATES[0])/(exp(- STATES[0]/4.80800) - 1.00000);
ALGEBRAIC[30] = ( 0.0114300*(STATES[0] - 5.00000))/(exp((STATES[0] - 5.00000)/2.50000) - 1.00000);
ALGEBRAIC[35] = 1.00000/(ALGEBRAIC[15]+ALGEBRAIC[30]);
ALGEBRAIC[8] = 0.0140000/(1.00000+exp(- (STATES[0] - 40.0000)/9.00000));
ALGEBRAIC[23] =  0.00100000*exp(- STATES[0]/45.0000);
ALGEBRAIC[31] = ALGEBRAIC[8]/(ALGEBRAIC[8]+ALGEBRAIC[23]);
ALGEBRAIC[36] = 1.00000/(ALGEBRAIC[8]+ALGEBRAIC[23]);
ALGEBRAIC[12] = 1.00000/(1.00000+exp(- (STATES[0]+57.0000)/5.00000));
ALGEBRAIC[27] = 1.00000/( 0.150000*exp(- STATES[0]/11.0000)+ 0.200000*exp(- STATES[0]/700.000));
ALGEBRAIC[32] = 1.00000/( 16.0000*exp(STATES[0]/8.00000)+ 15.0000*exp(STATES[0]/50.0000));
ALGEBRAIC[37] = 1.00000/(ALGEBRAIC[27]+ALGEBRAIC[32]);
ALGEBRAIC[13] = 1.00000/( 3100.00*exp(STATES[0]/13.0000)+ 700.000*exp(STATES[0]/70.0000));
ALGEBRAIC[28] = 1.00000/( 95.0000*exp(- STATES[0]/10.0000)+ 50.0000*exp(- STATES[0]/700.000))+0.000229000/(1.00000+exp(- STATES[0]/5.00000));
ALGEBRAIC[33] = ALGEBRAIC[13]/(ALGEBRAIC[13]+ALGEBRAIC[28]);
ALGEBRAIC[38] = 6.65000/(ALGEBRAIC[13]+ALGEBRAIC[28]);
ALGEBRAIC[29] =  (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[5]/STATES[2]);
ALGEBRAIC[40] =  CONSTANTS[63]*(STATES[0] - ALGEBRAIC[29])*( 0.600000*STATES[10]+ 0.400000*STATES[9])*STATES[11];
ALGEBRAIC[41] =  (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log((CONSTANTS[5]+ 0.120000*CONSTANTS[4])/(STATES[2]+ 0.120000*STATES[1]));
ALGEBRAIC[42] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[41])*pow(STATES[12], 2.00000);
ALGEBRAIC[43] =  CONSTANTS[13]*(STATES[0] - ALGEBRAIC[29])*STATES[13]*STATES[14];
ALGEBRAIC[44] =  CONSTANTS[14]*(STATES[0] - ALGEBRAIC[29])*STATES[14];
ALGEBRAIC[0] =  (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(CONSTANTS[4]/STATES[1]);
ALGEBRAIC[51] =  CONSTANTS[22]*pow(1.00000+pow(CONSTANTS[20]/CONSTANTS[5], 1.20000), - 1.00000)*pow(1.00000+pow(CONSTANTS[21]/STATES[1], 1.30000), - 1.00000)*pow(1.00000+exp(- ((STATES[0] - ALGEBRAIC[0])+120.000)/30.0000), - 1.00000);
ALGEBRAIC[50] = ( CONSTANTS[64]*(STATES[2] -  CONSTANTS[5]*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))/1.00000;
ALGEBRAIC[46] =  CONSTANTS[16]*(STATES[0] - ALGEBRAIC[29])*pow(STATES[15], 2.00000);
ALGEBRAIC[39] =  CONSTANTS[10]*(STATES[0] - CONSTANTS[11])*STATES[7]*STATES[8];
ALGEBRAIC[34] =  CONSTANTS[6]*(STATES[0] - CONSTANTS[7])*STATES[4]*STATES[5]*STATES[6];
ALGEBRAIC[45] =  CONSTANTS[15]*(STATES[0] - ALGEBRAIC[0])*pow(STATES[15], 2.00000);
ALGEBRAIC[47] = ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[48] =  CONSTANTS[17]*(STATES[0] - CONSTANTS[18])*STATES[16]*STATES[17];
ALGEBRAIC[49] =  CONSTANTS[19]*(STATES[0] - ALGEBRAIC[0]);
ALGEBRAIC[53] = exp(( - CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1]));
ALGEBRAIC[59] = 1.00000+ (CONSTANTS[36]/CONSTANTS[35])*(1.00000+exp(( CONSTANTS[26]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])))+ (CONSTANTS[4]/CONSTANTS[33])*(1.00000+CONSTANTS[4]/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[32]);
ALGEBRAIC[61] = ( (( (CONSTANTS[4]/CONSTANTS[33])*CONSTANTS[4])/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[32])*exp(( - CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[59];
ALGEBRAIC[60] = ( (CONSTANTS[36]/CONSTANTS[35])*exp(( CONSTANTS[26]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[59];
ALGEBRAIC[57] = exp(( CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1]));
ALGEBRAIC[52] = STATES[1]/(CONSTANTS[27]+STATES[1]);
ALGEBRAIC[62] =  ALGEBRAIC[53]*CONSTANTS[65]*(ALGEBRAIC[61]+ALGEBRAIC[60])+ ALGEBRAIC[60]*ALGEBRAIC[57]*(ALGEBRAIC[52]+ALGEBRAIC[53]);
ALGEBRAIC[54] = 1.00000+ (STATES[3]/CONSTANTS[28])*(1.00000+exp(( - CONSTANTS[24]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+STATES[1]/CONSTANTS[31])+ (STATES[1]/CONSTANTS[29])*(1.00000+ (STATES[1]/CONSTANTS[30])*(1.00000+STATES[1]/CONSTANTS[27]));
ALGEBRAIC[55] = ( (STATES[3]/CONSTANTS[28])*exp(( - CONSTANTS[24]*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[54];
ALGEBRAIC[56] = ( (( (STATES[1]/CONSTANTS[29])*STATES[1])/CONSTANTS[30])*(1.00000+STATES[1]/CONSTANTS[27])*exp(( CONSTANTS[25]*STATES[0]*CONSTANTS[2])/( 2.00000*CONSTANTS[0]*CONSTANTS[1])))/ALGEBRAIC[54];
ALGEBRAIC[58] =  ALGEBRAIC[57]*ALGEBRAIC[52]*(ALGEBRAIC[56]+ALGEBRAIC[55])+ ALGEBRAIC[53]*ALGEBRAIC[55]*(CONSTANTS[65]+ALGEBRAIC[57]);
ALGEBRAIC[63] =  ALGEBRAIC[56]*ALGEBRAIC[52]*(ALGEBRAIC[61]+ALGEBRAIC[60])+ ALGEBRAIC[55]*ALGEBRAIC[61]*(ALGEBRAIC[52]+ALGEBRAIC[53]);
ALGEBRAIC[64] =  ALGEBRAIC[61]*CONSTANTS[65]*(ALGEBRAIC[56]+ALGEBRAIC[55])+ ALGEBRAIC[56]*ALGEBRAIC[60]*(CONSTANTS[65]+ALGEBRAIC[57]);
ALGEBRAIC[65] = ( CONSTANTS[23]*( ALGEBRAIC[58]*ALGEBRAIC[60] -  ALGEBRAIC[62]*ALGEBRAIC[55]))/(ALGEBRAIC[62]+ALGEBRAIC[58]+ALGEBRAIC[63]+ALGEBRAIC[64]);
ALGEBRAIC[66] = (STATES[3] - STATES[19])/CONSTANTS[37];
ALGEBRAIC[67] = ( CONSTANTS[41]*(STATES[20] - STATES[3]))/(1.00000+pow(CONSTANTS[39]/STATES[3], 2.00000));
ALGEBRAIC[69] =  CONSTANTS[55]*STATES[3]*(1.00000 - STATES[23]) -  CONSTANTS[60]*STATES[23];
ALGEBRAIC[68] = CONSTANTS[40]/(1.00000+CONSTANTS[42]/STATES[19]);
ALGEBRAIC[70] = (STATES[18] - STATES[20])/CONSTANTS[38];
ALGEBRAIC[71] =  CONSTANTS[52]*STATES[19]*(1.00000 - STATES[24]) -  CONSTANTS[57]*STATES[24];
ALGEBRAIC[73] =  CONSTANTS[56]*STATES[20]*(1.00000 - STATES[26]) -  CONSTANTS[61]*STATES[26];
ALGEBRAIC[72] =  CONSTANTS[54]*STATES[19]*(1.00000 - (STATES[25]+STATES[21])) -  CONSTANTS[58]*STATES[25];
ALGEBRAIC[74] =  CONSTANTS[55]*STATES[19]*(1.00000 - STATES[22]) -  CONSTANTS[60]*STATES[22];
}