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 76 entries in the algebraic variable array.
There are a total of 29 entries in each of the rate and state variable arrays.
There are a total of 86 entries in the constant variable array.
*/
/*
* VOI is time in component environment (millisecond).
* STATES[0] is Vm in component Vm (millivolt).
* CONSTANTS[0] is Cm in component model_parameters (picoF).
* ALGEBRAIC[31] is i_CaT in component i_CaT (picoA).
* ALGEBRAIC[14] is i_CaL in component i_CaL (picoA).
* ALGEBRAIC[49] is i_f in component i_f (picoA).
* ALGEBRAIC[50] is i_st in component i_st (picoA).
* ALGEBRAIC[37] is i_Kr in component i_Kr (picoA).
* ALGEBRAIC[42] is i_Ks in component i_Ks (picoA).
* ALGEBRAIC[44] is i_to in component AP_sensitive_currents (picoA).
* ALGEBRAIC[46] is i_sus in component AP_sensitive_currents (picoA).
* ALGEBRAIC[52] is i_NaK in component i_NaK (picoA).
* ALGEBRAIC[66] is i_NaCa in component i_NaCa (picoA).
* ALGEBRAIC[53] is i_b_Ca in component i_b_Ca (picoA).
* ALGEBRAIC[51] is i_b_Na in component i_b_Na (picoA).
* CONSTANTS[75] is E_Na in component electric_potentials (millivolt).
* CONSTANTS[76] is E_K in component electric_potentials (millivolt).
* CONSTANTS[77] is E_Ks in component electric_potentials (millivolt).
* CONSTANTS[1] is R in component model_parameters (joule_per_kilomole_kelvin).
* CONSTANTS[2] is T in component model_parameters (kelvin).
* CONSTANTS[3] is F in component model_parameters (coulomb_per_mole).
* CONSTANTS[4] is Nai in component model_parameters (millimolar).
* CONSTANTS[5] is Nao in component model_parameters (millimolar).
* CONSTANTS[6] is Ki in component model_parameters (millimolar).
* CONSTANTS[7] is Ko in component model_parameters (millimolar).
* CONSTANTS[8] is E_CaL in component i_CaL (millivolt).
* CONSTANTS[9] is g_CaL in component i_CaL (nanoS_per_picoF).
* STATES[1] is Ca_sub in component calcium_dynamics (millimolar).
* STATES[2] is dL in component i_CaL_dL_gate (dimensionless).
* STATES[3] is fL in component i_CaL_fL_gate (dimensionless).
* STATES[4] is fCa in component i_CaL_fCa_gate (dimensionless).
* ALGEBRAIC[0] is dL_infinity in component i_CaL_dL_gate (dimensionless).
* ALGEBRAIC[45] is tau_dL in component i_CaL_dL_gate (millisecond).
* ALGEBRAIC[32] is alpha_dL in component i_CaL_dL_gate (per_millisecond).
* ALGEBRAIC[43] is beta_dL in component i_CaL_dL_gate (per_millisecond).
* ALGEBRAIC[16] is adVm in component i_CaL_dL_gate (millivolt).
* ALGEBRAIC[38] is bdVm in component i_CaL_dL_gate (millivolt).
* ALGEBRAIC[1] is fL_infinity in component i_CaL_fL_gate (dimensionless).
* ALGEBRAIC[17] is tau_fL in component i_CaL_fL_gate (millisecond).
* CONSTANTS[10] is alpha_fCa in component i_CaL_fCa_gate (per_millisecond).
* ALGEBRAIC[2] is fCa_infinity in component i_CaL_fCa_gate (dimensionless).
* ALGEBRAIC[18] is tau_fCa in component i_CaL_fCa_gate (millisecond).
* CONSTANTS[11] is Km_fCa in component i_CaL_fCa_gate (millimolar).
* CONSTANTS[12] is g_CaT in component i_CaT (nanoS_per_picoF).
* CONSTANTS[13] is E_CaT in component i_CaT (millivolt).
* STATES[5] is dT in component i_CaT_dT_gate (dimensionless).
* STATES[6] is fT in component i_CaT_fT_gate (dimensionless).
* ALGEBRAIC[3] is dT_infinity in component i_CaT_dT_gate (dimensionless).
* ALGEBRAIC[19] is tau_dT in component i_CaT_dT_gate (millisecond).
* ALGEBRAIC[4] is fT_infinity in component i_CaT_fT_gate (dimensionless).
* ALGEBRAIC[20] is tau_fT in component i_CaT_fT_gate (millisecond).
* CONSTANTS[14] is g_Kr in component i_Kr (nanoS_per_picoF).
* STATES[7] is paS in component i_Kr_pa_gate (dimensionless).
* STATES[8] is paF in component i_Kr_pa_gate (dimensionless).
* STATES[9] is pi_ in component i_Kr_pi_gate (dimensionless).
* ALGEBRAIC[5] is pa_infinity in component i_Kr_pa_gate (dimensionless).
* ALGEBRAIC[21] is tau_paS in component i_Kr_pa_gate (millisecond).
* ALGEBRAIC[22] is tau_paF in component i_Kr_pa_gate (millisecond).
* ALGEBRAIC[6] is pi_infinity in component i_Kr_pi_gate (dimensionless).
* ALGEBRAIC[23] is tau_pi in component i_Kr_pi_gate (millisecond).
* CONSTANTS[15] is g_Ks in component i_Ks (nanoS_per_picoF).
* STATES[10] is n in component i_Ks_n_gate (dimensionless).
* ALGEBRAIC[33] is n_infinity in component i_Ks_n_gate (dimensionless).
* ALGEBRAIC[39] is tau_n in component i_Ks_n_gate (millisecond).
* ALGEBRAIC[7] is alpha_n in component i_Ks_n_gate (per_millisecond).
* ALGEBRAIC[24] is beta_n in component i_Ks_n_gate (per_millisecond).
* CONSTANTS[16] is g_to in component AP_sensitive_currents (nanoS_per_picoF).
* CONSTANTS[17] is g_sus in component AP_sensitive_currents (nanoS_per_picoF).
* STATES[11] is q in component AP_sensitive_currents_q_gate (dimensionless).
* STATES[12] is r in component AP_sensitive_currents_r_gate (dimensionless).
* ALGEBRAIC[8] is q_infinity in component AP_sensitive_currents_q_gate (dimensionless).
* ALGEBRAIC[25] is tau_q in component AP_sensitive_currents_q_gate (millisecond).
* ALGEBRAIC[9] is r_infinity in component AP_sensitive_currents_r_gate (dimensionless).
* ALGEBRAIC[26] is tau_r in component AP_sensitive_currents_r_gate (millisecond).
* ALGEBRAIC[47] is i_f_Na in component i_f (picoA).
* ALGEBRAIC[48] is i_f_K in component i_f (picoA).
* CONSTANTS[18] is g_if in component i_f (nanoS_per_picoF).
* STATES[13] is y in component i_f_y_gate (dimensionless).
* ALGEBRAIC[10] is y_infinity in component i_f_y_gate (dimensionless).
* ALGEBRAIC[27] is tau_y in component i_f_y_gate (millisecond).
* CONSTANTS[19] is VIf_half in component i_f_y_gate (millivolt).
* CONSTANTS[20] is g_st in component i_st (nanoS_per_picoF).
* CONSTANTS[21] is E_st in component i_st (millivolt).
* STATES[14] is qa in component i_st_qa_gate (dimensionless).
* STATES[15] is qi in component i_st_qi_gate (dimensionless).
* ALGEBRAIC[11] is qa_infinity in component i_st_qa_gate (dimensionless).
* ALGEBRAIC[40] is tau_qa in component i_st_qa_gate (millisecond).
* ALGEBRAIC[28] is alpha_qa in component i_st_qa_gate (per_millisecond).
* ALGEBRAIC[34] is beta_qa in component i_st_qa_gate (per_millisecond).
* ALGEBRAIC[35] is qi_infinity in component i_st_qi_gate (dimensionless).
* ALGEBRAIC[41] is tau_qi in component i_st_qi_gate (millisecond).
* ALGEBRAIC[12] is alpha_qi in component i_st_qi_gate (per_millisecond).
* ALGEBRAIC[29] is beta_qi in component i_st_qi_gate (per_millisecond).
* CONSTANTS[22] is g_b_Na in component i_b_Na (nanoS_per_picoF).
* CONSTANTS[23] is Km_Kp in component i_NaK (millimolar).
* CONSTANTS[24] is Km_Nap in component i_NaK (millimolar).
* CONSTANTS[25] is i_NaK_max in component i_NaK (picoA_per_picoF).
* CONSTANTS[26] is g_b_Ca in component i_b_Ca (nanoS_per_picoF).
* CONSTANTS[27] is kNaCa in component i_NaCa (picoA_per_picoF).
* ALGEBRAIC[63] is x1 in component i_NaCa (dimensionless).
* ALGEBRAIC[59] is x2 in component i_NaCa (dimensionless).
* ALGEBRAIC[64] is x3 in component i_NaCa (dimensionless).
* ALGEBRAIC[65] is x4 in component i_NaCa (dimensionless).
* ALGEBRAIC[54] is k41 in component i_NaCa (dimensionless).
* CONSTANTS[80] is k34 in component i_NaCa (dimensionless).
* ALGEBRAIC[62] is k23 in component i_NaCa (dimensionless).
* ALGEBRAIC[61] is k21 in component i_NaCa (dimensionless).
* ALGEBRAIC[58] is k32 in component i_NaCa (dimensionless).
* CONSTANTS[79] is k43 in component i_NaCa (dimensionless).
* ALGEBRAIC[56] is k12 in component i_NaCa (dimensionless).
* ALGEBRAIC[57] is k14 in component i_NaCa (dimensionless).
* CONSTANTS[28] is Qci in component i_NaCa (dimensionless).
* CONSTANTS[29] is Qn in component i_NaCa (dimensionless).
* CONSTANTS[30] is Qco in component i_NaCa (dimensionless).
* CONSTANTS[31] is K3ni in component i_NaCa (millimolar).
* CONSTANTS[32] is Kci in component i_NaCa (millimolar).
* CONSTANTS[33] is K1ni in component i_NaCa (millimolar).
* CONSTANTS[34] is K2ni in component i_NaCa (millimolar).
* CONSTANTS[35] is Kcni in component i_NaCa (millimolar).
* CONSTANTS[36] is K3no in component i_NaCa (millimolar).
* CONSTANTS[37] is K1no in component i_NaCa (millimolar).
* CONSTANTS[38] is K2no in component i_NaCa (millimolar).
* CONSTANTS[39] is Kco in component i_NaCa (millimolar).
* CONSTANTS[78] is RTOnF in component i_NaCa (millivolt).
* ALGEBRAIC[60] is do in component i_NaCa (dimensionless).
* ALGEBRAIC[55] is di in component i_NaCa (dimensionless).
* CONSTANTS[40] is Cao in component model_parameters (millimolar).
* ALGEBRAIC[67] is j_SRCarel in component j_SRCarel (millimolar_per_millisecond).
* STATES[16] is R in component j_SRCarel (dimensionless).
* STATES[17] is O in component j_SRCarel (dimensionless).
* STATES[18] is I in component j_SRCarel (dimensionless).
* STATES[19] is RI in component j_SRCarel (dimensionless).
* CONSTANTS[41] is ks in component j_SRCarel (per_millisecond).
* CONSTANTS[42] is MaxSR in component j_SRCarel (dimensionless).
* CONSTANTS[43] is MinSR in component j_SRCarel (dimensionless).
* CONSTANTS[44] is EC50_SR in component j_SRCarel (millimolar).
* CONSTANTS[45] is HSR in component j_SRCarel (dimensionless).
* ALGEBRAIC[30] is koSRCa in component j_SRCarel (per_millimolar2_millisecond).
* ALGEBRAIC[36] is kiSRCa in component j_SRCarel (per_millimolar_millisecond).
* CONSTANTS[46] is koCa in component j_SRCarel (per_millimolar2_millisecond).
* CONSTANTS[47] is kiCa in component j_SRCarel (per_millimolar_millisecond).
* ALGEBRAIC[13] is kCaSR in component j_SRCarel (dimensionless).
* CONSTANTS[48] is kim in component j_SRCarel (per_millisecond).
* CONSTANTS[49] is kom in component j_SRCarel (per_millisecond).
* STATES[20] is Ca_jsr in component calcium_dynamics (millimolar).
* ALGEBRAIC[68] is j_Ca_dif in component intracellular_calcium_fluxes (millimolar_per_millisecond).
* ALGEBRAIC[70] is j_up in component intracellular_calcium_fluxes (millimolar_per_millisecond).
* ALGEBRAIC[72] is j_tr in component intracellular_calcium_fluxes (millimolar_per_millisecond).
* CONSTANTS[50] is tau_dif_Ca in component intracellular_calcium_fluxes (millisecond).
* CONSTANTS[51] is tau_tr in component intracellular_calcium_fluxes (millisecond).
* CONSTANTS[52] is P_up in component intracellular_calcium_fluxes (millimolar_per_millisecond).
* CONSTANTS[53] is K_up in component intracellular_calcium_fluxes (millimolar).
* STATES[21] is Ca_nsr in component calcium_dynamics (millimolar).
* STATES[22] is Cai in component calcium_dynamics (millimolar).
* CONSTANTS[54] is TC_tot in component calcium_buffering (millimolar).
* CONSTANTS[55] is TMC_tot in component calcium_buffering (millimolar).
* CONSTANTS[56] is CM_tot in component calcium_buffering (millimolar).
* CONSTANTS[57] is CQ_tot in component calcium_buffering (millimolar).
* ALGEBRAIC[71] is delta_fTC in component calcium_buffering (per_millisecond).
* ALGEBRAIC[73] is delta_fTMC in component calcium_buffering (per_millisecond).
* ALGEBRAIC[69] is delta_fCMs in component calcium_buffering (per_millisecond).
* ALGEBRAIC[75] is delta_fCMi in component calcium_buffering (per_millisecond).
* ALGEBRAIC[74] is delta_fCQ in component calcium_buffering (per_millisecond).
* ALGEBRAIC[15] is delta_fTMM in component calcium_buffering (per_millisecond).
* STATES[23] is fTMM in component calcium_buffering (dimensionless).
* STATES[24] is fCMi in component calcium_buffering (dimensionless).
* STATES[25] is fCMs in component calcium_buffering (dimensionless).
* STATES[26] is fTC in component calcium_buffering (dimensionless).
* STATES[27] is fTMC in component calcium_buffering (dimensionless).
* STATES[28] is fCQ in component calcium_buffering (dimensionless).
* CONSTANTS[58] is kf_TC in component calcium_buffering (per_millimolar_millisecond).
* CONSTANTS[59] is kf_TMM in component calcium_buffering (per_millimolar_millisecond).
* CONSTANTS[60] is kf_TMC in component calcium_buffering (per_millimolar_millisecond).
* CONSTANTS[61] is kf_CM in component calcium_buffering (per_millimolar_millisecond).
* CONSTANTS[62] is kf_CQ in component calcium_buffering (per_millimolar_millisecond).
* CONSTANTS[63] is kb_TC in component calcium_buffering (per_millisecond).
* CONSTANTS[64] is kb_TMC in component calcium_buffering (per_millisecond).
* CONSTANTS[65] is kb_TMM in component calcium_buffering (per_millisecond).
* CONSTANTS[66] is kb_CM in component calcium_buffering (per_millisecond).
* CONSTANTS[67] is kb_CQ in component calcium_buffering (per_millisecond).
* CONSTANTS[68] is Mgi in component model_parameters (millimolar).
* CONSTANTS[84] is V_i in component model_parameters (picolitre).
* CONSTANTS[83] is V_jsr in component model_parameters (picolitre).
* CONSTANTS[85] is V_nsr in component model_parameters (picolitre).
* CONSTANTS[82] is V_sub in component model_parameters (picolitre).
* CONSTANTS[81] is V_cell in component model_parameters (picolitre).
* CONSTANTS[69] is V_jsr_part in component model_parameters (dimensionless).
* CONSTANTS[70] is V_i_part in component model_parameters (dimensionless).
* CONSTANTS[71] is V_nsr_part in component model_parameters (dimensionless).
* CONSTANTS[72] is R_cell in component model_parameters (micrometre).
* CONSTANTS[73] is L_cell in component model_parameters (micrometre).
* CONSTANTS[74] is L_sub in component model_parameters (micrometre).
* RATES[0] is d/dt Vm in component Vm (millivolt).
* RATES[2] is d/dt dL in component i_CaL_dL_gate (dimensionless).
* RATES[3] is d/dt fL in component i_CaL_fL_gate (dimensionless).
* RATES[4] is d/dt fCa in component i_CaL_fCa_gate (dimensionless).
* RATES[5] is d/dt dT in component i_CaT_dT_gate (dimensionless).
* RATES[6] is d/dt fT in component i_CaT_fT_gate (dimensionless).
* RATES[7] is d/dt paS in component i_Kr_pa_gate (dimensionless).
* RATES[8] is d/dt paF in component i_Kr_pa_gate (dimensionless).
* RATES[9] is d/dt pi_ in component i_Kr_pi_gate (dimensionless).
* RATES[10] is d/dt n in component i_Ks_n_gate (dimensionless).
* RATES[11] is d/dt q in component AP_sensitive_currents_q_gate (dimensionless).
* RATES[12] is d/dt r in component AP_sensitive_currents_r_gate (dimensionless).
* RATES[13] is d/dt y in component i_f_y_gate (dimensionless).
* RATES[14] is d/dt qa in component i_st_qa_gate (dimensionless).
* RATES[15] is d/dt qi in component i_st_qi_gate (dimensionless).
* RATES[16] is d/dt R in component j_SRCarel (dimensionless).
* RATES[17] is d/dt O in component j_SRCarel (dimensionless).
* RATES[18] is d/dt I in component j_SRCarel (dimensionless).
* RATES[19] is d/dt RI in component j_SRCarel (dimensionless).
* RATES[26] is d/dt fTC in component calcium_buffering (dimensionless).
* RATES[27] is d/dt fTMC in component calcium_buffering (dimensionless).
* RATES[23] is d/dt fTMM in component calcium_buffering (dimensionless).
* RATES[24] is d/dt fCMi in component calcium_buffering (dimensionless).
* RATES[25] is d/dt fCMs in component calcium_buffering (dimensionless).
* RATES[28] is d/dt fCQ in component calcium_buffering (dimensionless).
* RATES[22] is d/dt Cai in component calcium_dynamics (millimolar).
* RATES[1] is d/dt Ca_sub in component calcium_dynamics (millimolar).
* RATES[21] is d/dt Ca_nsr in component calcium_dynamics (millimolar).
* RATES[20] is d/dt Ca_jsr in component calcium_dynamics (millimolar).
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -65;
CONSTANTS[0] = 32;
CONSTANTS[1] = 8314.4;
CONSTANTS[2] = 310.15;
CONSTANTS[3] = 96485;
CONSTANTS[4] = 10;
CONSTANTS[5] = 140;
CONSTANTS[6] = 140;
CONSTANTS[7] = 5.4;
CONSTANTS[8] = 45;
CONSTANTS[9] = 0.464;
STATES[1] = 0.000223;
STATES[2] = 0;
STATES[3] = 1;
STATES[4] = 1;
CONSTANTS[10] = 0.021;
CONSTANTS[11] = 0.00035;
CONSTANTS[12] = 0.1832;
CONSTANTS[13] = 45;
STATES[5] = 0;
STATES[6] = 1;
CONSTANTS[14] = 0.08113973;
STATES[7] = 0;
STATES[8] = 0;
STATES[9] = 1;
CONSTANTS[15] = 0.0259;
STATES[10] = 0;
CONSTANTS[16] = 0.252;
CONSTANTS[17] = 0.02;
STATES[11] = 1;
STATES[12] = 0;
CONSTANTS[18] = 0.15;
STATES[13] = 1;
CONSTANTS[19] = -64;
CONSTANTS[20] = 0.003;
CONSTANTS[21] = 37.4;
STATES[14] = 0;
STATES[15] = 1;
CONSTANTS[22] = 0.00486;
CONSTANTS[23] = 1.4;
CONSTANTS[24] = 14;
CONSTANTS[25] = 2.88;
CONSTANTS[26] = 0.0006;
CONSTANTS[27] = 187.5;
CONSTANTS[28] = 0.1369;
CONSTANTS[29] = 0.4315;
CONSTANTS[30] = 0;
CONSTANTS[31] = 26.44;
CONSTANTS[32] = 0.0207;
CONSTANTS[33] = 395.3;
CONSTANTS[34] = 2.289;
CONSTANTS[35] = 26.44;
CONSTANTS[36] = 4.663;
CONSTANTS[37] = 1628;
CONSTANTS[38] = 561.4;
CONSTANTS[39] = 3.663;
CONSTANTS[40] = 2;
STATES[16] = 0.7499955;
STATES[17] = 3.4e-6;
STATES[18] = 1.1e-6;
STATES[19] = 0.25;
CONSTANTS[41] = 250000;
CONSTANTS[42] = 15;
CONSTANTS[43] = 1;
CONSTANTS[44] = 0.45;
CONSTANTS[45] = 2.5;
CONSTANTS[46] = 10;
CONSTANTS[47] = 0.5;
CONSTANTS[48] = 0.005;
CONSTANTS[49] = 0.06;
STATES[20] = 0.029;
CONSTANTS[50] = 0.04;
CONSTANTS[51] = 40;
CONSTANTS[52] = 0.012;
CONSTANTS[53] = 0.0006;
STATES[21] = 1.35;
STATES[22] = 0.0001;
CONSTANTS[54] = 0.031;
CONSTANTS[55] = 0.062;
CONSTANTS[56] = 0.045;
CONSTANTS[57] = 10;
STATES[23] = 0.69;
STATES[24] = 0.042;
STATES[25] = 0.089;
STATES[26] = 0.02;
STATES[27] = 0.22;
STATES[28] = 0.032;
CONSTANTS[58] = 88.8;
CONSTANTS[59] = 2.277;
CONSTANTS[60] = 227.7;
CONSTANTS[61] = 227.7;
CONSTANTS[62] = 0.534;
CONSTANTS[63] = 0.446;
CONSTANTS[64] = 0.00751;
CONSTANTS[65] = 0.751;
CONSTANTS[66] = 0.542;
CONSTANTS[67] = 0.445;
CONSTANTS[68] = 2.5;
CONSTANTS[69] = 0.0012;
CONSTANTS[70] = 0.46;
CONSTANTS[71] = 0.0116;
CONSTANTS[72] = 4;
CONSTANTS[73] = 70;
CONSTANTS[74] = 0.02;
CONSTANTS[75] = (( CONSTANTS[1]*CONSTANTS[2])/CONSTANTS[3])*log(CONSTANTS[5]/CONSTANTS[4]);
CONSTANTS[76] = (( CONSTANTS[1]*CONSTANTS[2])/CONSTANTS[3])*log(CONSTANTS[7]/CONSTANTS[6]);
CONSTANTS[77] = (( CONSTANTS[1]*CONSTANTS[2])/CONSTANTS[3])*log((CONSTANTS[7]+ 0.120000*CONSTANTS[5])/(CONSTANTS[6]+ 0.120000*CONSTANTS[4]));
CONSTANTS[78] = ( CONSTANTS[1]*CONSTANTS[2])/CONSTANTS[3];
CONSTANTS[79] = CONSTANTS[4]/(CONSTANTS[31]+CONSTANTS[4]);
CONSTANTS[80] = CONSTANTS[5]/(CONSTANTS[36]+CONSTANTS[5]);
CONSTANTS[81] = 0.00100000* 3.14159265358979*pow(CONSTANTS[72], 2.00000)*CONSTANTS[73];
CONSTANTS[82] = 0.00100000*2.00000* 3.14159265358979*CONSTANTS[74]*(CONSTANTS[72] - CONSTANTS[74]/2.00000)*CONSTANTS[73];
CONSTANTS[83] = CONSTANTS[69]*CONSTANTS[81];
CONSTANTS[84] = CONSTANTS[70]*CONSTANTS[81] - CONSTANTS[82];
CONSTANTS[85] = CONSTANTS[71]*CONSTANTS[81];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[15] = CONSTANTS[59]*CONSTANTS[68]*(1.00000 - (STATES[27]+STATES[23])) - CONSTANTS[65]*STATES[23];
RATES[23] = ALGEBRAIC[15];
ALGEBRAIC[1] = 1.00000/(1.00000+exp((STATES[0]+35.0000)/7.30000));
ALGEBRAIC[17] = 44.3000+ 257.100*exp(- pow((STATES[0]+32.5000)/13.9000, 2.00000));
RATES[3] = (ALGEBRAIC[1] - STATES[3])/ALGEBRAIC[17];
ALGEBRAIC[2] = CONSTANTS[11]/(CONSTANTS[11]+STATES[1]);
ALGEBRAIC[18] = ALGEBRAIC[2]/CONSTANTS[10];
RATES[4] = (ALGEBRAIC[2] - STATES[4])/ALGEBRAIC[18];
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (STATES[0]+26.3000)/6.00000));
ALGEBRAIC[19] = 1.00000/( 1.06800*exp((STATES[0]+26.3000)/30.0000)+ 1.06800*exp(- (STATES[0]+26.3000)/30.0000));
RATES[5] = (ALGEBRAIC[3] - STATES[5])/ALGEBRAIC[19];
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+61.7000)/5.60000));
ALGEBRAIC[20] = 1.00000/( 0.0153000*exp(- (STATES[0]+61.7000)/83.3000)+ 0.0150000*exp((STATES[0]+61.7000)/15.3800));
RATES[6] = (ALGEBRAIC[4] - STATES[6])/ALGEBRAIC[20];
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0]+23.2000)/10.6000));
ALGEBRAIC[21] = 0.846554/( 0.00420000*exp(STATES[0]/17.0000)+ 0.000150000*exp(- STATES[0]/21.6000));
RATES[7] = (ALGEBRAIC[5] - STATES[7])/ALGEBRAIC[21];
ALGEBRAIC[22] = 0.846554/( 0.0372000*exp(STATES[0]/15.9000)+ 0.000960000*exp(- STATES[0]/22.5000));
RATES[8] = (ALGEBRAIC[5] - STATES[8])/ALGEBRAIC[22];
ALGEBRAIC[6] = 1.00000/(1.00000+exp((STATES[0]+28.6000)/17.1000));
ALGEBRAIC[23] = 1.00000/( 0.100000*exp(- STATES[0]/54.6450)+ 0.656000*exp(STATES[0]/106.157));
RATES[9] = (ALGEBRAIC[6] - STATES[9])/ALGEBRAIC[23];
ALGEBRAIC[8] = 1.00000/(1.00000+exp((STATES[0]+49.0000)/13.0000));
ALGEBRAIC[25] = 6.06000+39.1020/( 0.570000*exp( - 0.0800000*(STATES[0]+44.0000))+ 0.0650000*exp( 0.100000*(STATES[0]+45.9300)));
RATES[11] = (ALGEBRAIC[8] - STATES[11])/ALGEBRAIC[25];
ALGEBRAIC[9] = 1.00000/(1.00000+exp(- (STATES[0] - 19.3000)/15.0000));
ALGEBRAIC[26] = 2.75352+14.4052/( 1.03700*exp( 0.0900000*(STATES[0]+30.6100))+ 0.369000*exp( - 0.120000*(STATES[0]+23.8400)));
RATES[12] = (ALGEBRAIC[9] - STATES[12])/ALGEBRAIC[26];
ALGEBRAIC[10] = 1.00000/(1.00000+exp((STATES[0] - CONSTANTS[19])/13.5000));
ALGEBRAIC[27] = 0.716653/(exp(- (STATES[0]+386.900)/45.3020)+exp((STATES[0] - 73.0800)/19.2310));
RATES[13] = (ALGEBRAIC[10] - STATES[13])/ALGEBRAIC[27];
ALGEBRAIC[13] = CONSTANTS[42] - (CONSTANTS[42] - CONSTANTS[43])/(1.00000+pow(CONSTANTS[44]/STATES[20], CONSTANTS[45]));
ALGEBRAIC[30] = CONSTANTS[46]/ALGEBRAIC[13];
ALGEBRAIC[36] = CONSTANTS[47]*ALGEBRAIC[13];
RATES[16] = ( CONSTANTS[48]*STATES[19] - ALGEBRAIC[36]*STATES[1]*STATES[16]) - ( ALGEBRAIC[30]*pow(STATES[1], 2.00000)*STATES[16] - CONSTANTS[49]*STATES[17]);
RATES[17] = ( ALGEBRAIC[30]*pow(STATES[1], 2.00000)*STATES[16] - CONSTANTS[49]*STATES[17]) - ( ALGEBRAIC[36]*STATES[1]*STATES[17] - CONSTANTS[48]*STATES[18]);
RATES[18] = ( ALGEBRAIC[36]*STATES[1]*STATES[17] - CONSTANTS[48]*STATES[18]) - ( CONSTANTS[49]*STATES[18] - ALGEBRAIC[30]*pow(STATES[1], 2.00000)*STATES[19]);
RATES[19] = ( CONSTANTS[49]*STATES[18] - ALGEBRAIC[30]*pow(STATES[1], 2.00000)*STATES[19]) - ( CONSTANTS[48]*STATES[19] - ALGEBRAIC[36]*STATES[1]*STATES[16]);
ALGEBRAIC[7] = 0.0140000/(1.00000+exp(- (STATES[0] - 40.0000)/9.00000));
ALGEBRAIC[24] = 0.00100000*exp(- STATES[0]/45.0000);
ALGEBRAIC[33] = ALGEBRAIC[7]/(ALGEBRAIC[7]+ALGEBRAIC[24]);
ALGEBRAIC[39] = 1.00000/(ALGEBRAIC[7]+ALGEBRAIC[24]);
RATES[10] = (ALGEBRAIC[33] - STATES[10])/ALGEBRAIC[39];
ALGEBRAIC[11] = 1.00000/(1.00000+exp(- (STATES[0]+57.0000)/5.00000));
ALGEBRAIC[28] = 1.00000/( 0.150000*exp(- STATES[0]/11.0000)+ 0.200000*exp(- STATES[0]/700.000));
ALGEBRAIC[34] = 1.00000/( 16.0000*exp(STATES[0]/8.00000)+ 15.0000*exp(STATES[0]/50.0000));
ALGEBRAIC[40] = 1.00000/(ALGEBRAIC[28]+ALGEBRAIC[34]);
RATES[14] = (ALGEBRAIC[11] - STATES[14])/ALGEBRAIC[40];
ALGEBRAIC[12] = 1.00000/( 3100.00*exp(STATES[0]/13.0000)+ 700.000*exp(STATES[0]/70.0000));
ALGEBRAIC[29] = 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[35] = ALGEBRAIC[12]/(ALGEBRAIC[12]+ALGEBRAIC[29]);
ALGEBRAIC[41] = 6.65000/(ALGEBRAIC[12]+ALGEBRAIC[29]);
RATES[15] = (ALGEBRAIC[35] - STATES[15])/ALGEBRAIC[41];
ALGEBRAIC[0] = 1.00000/(1.00000+exp(- (STATES[0]+13.5000)/6.00000));
ALGEBRAIC[16] = (STATES[0]==- 35.0000 ? - 35.0000 : STATES[0]==0.00000 ? 1.00000e-05 : STATES[0]);
ALGEBRAIC[32] = ( - 0.0283900*(ALGEBRAIC[16]+35.0000))/(exp(- (ALGEBRAIC[16]+35.0000)/2.50000) - 1.00000) - ( 0.0849000*ALGEBRAIC[16])/(exp(- ALGEBRAIC[16]/4.80000) - 1.00000);
ALGEBRAIC[38] = (STATES[0]==5.00000 ? 5.00001 : STATES[0]);
ALGEBRAIC[43] = ( 0.0114300*(ALGEBRAIC[38] - 5.00000))/(exp((ALGEBRAIC[38] - 5.00000)/2.50000) - 1.00000);
ALGEBRAIC[45] = 1.00000/(ALGEBRAIC[32]+ALGEBRAIC[43]);
RATES[2] = (ALGEBRAIC[0] - STATES[2])/ALGEBRAIC[45];
ALGEBRAIC[31] = CONSTANTS[0]*CONSTANTS[12]*(STATES[0] - CONSTANTS[13])*STATES[5]*STATES[6];
ALGEBRAIC[14] = CONSTANTS[0]*CONSTANTS[9]*(STATES[0] - CONSTANTS[8])*STATES[2]*STATES[3]*STATES[4];
ALGEBRAIC[47] = CONSTANTS[0]*0.383300*CONSTANTS[18]*(STATES[0] - CONSTANTS[75])*pow(STATES[13], 2.00000);
ALGEBRAIC[48] = CONSTANTS[0]*0.616700*CONSTANTS[18]*(STATES[0] - CONSTANTS[76])*pow(STATES[13], 2.00000);
ALGEBRAIC[49] = ALGEBRAIC[47]+ALGEBRAIC[48];
ALGEBRAIC[50] = CONSTANTS[0]*CONSTANTS[20]*(STATES[0] - CONSTANTS[21])*STATES[14]*STATES[15];
ALGEBRAIC[37] = CONSTANTS[0]*CONSTANTS[14]*(STATES[0] - CONSTANTS[76])*( 0.600000*STATES[8]+ 0.400000*STATES[7])*STATES[9];
ALGEBRAIC[42] = CONSTANTS[0]*CONSTANTS[15]*(STATES[0] - CONSTANTS[77])*pow(STATES[10], 2.00000);
ALGEBRAIC[44] = CONSTANTS[0]*CONSTANTS[16]*(STATES[0] - CONSTANTS[76])*STATES[11]*STATES[12];
ALGEBRAIC[46] = CONSTANTS[0]*CONSTANTS[17]*(STATES[0] - CONSTANTS[76])*STATES[12];
ALGEBRAIC[52] = ( CONSTANTS[0]*CONSTANTS[25])/( (1.00000+pow(CONSTANTS[23]/CONSTANTS[7], 1.20000))*(1.00000+pow(CONSTANTS[24]/CONSTANTS[4], 1.30000))*(1.00000+exp(- ((STATES[0] - CONSTANTS[75])+120.000)/30.0000)));
ALGEBRAIC[54] = exp(( - CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78]));
ALGEBRAIC[60] = 1.00000+ (CONSTANTS[40]/CONSTANTS[39])*(1.00000+exp(( CONSTANTS[30]*STATES[0])/CONSTANTS[78]))+ (CONSTANTS[5]/CONSTANTS[37])*(1.00000+ (CONSTANTS[5]/CONSTANTS[38])*(1.00000+CONSTANTS[5]/CONSTANTS[36]));
ALGEBRAIC[62] = ( (( (CONSTANTS[5]/CONSTANTS[37])*CONSTANTS[5])/CONSTANTS[38])*(1.00000+CONSTANTS[5]/CONSTANTS[36])*exp(( - CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78])))/ALGEBRAIC[60];
ALGEBRAIC[61] = ( (CONSTANTS[40]/CONSTANTS[39])*exp(( CONSTANTS[30]*STATES[0])/CONSTANTS[78]))/ALGEBRAIC[60];
ALGEBRAIC[58] = exp(( CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78]));
ALGEBRAIC[63] = ALGEBRAIC[54]*CONSTANTS[80]*(ALGEBRAIC[62]+ALGEBRAIC[61])+ ALGEBRAIC[61]*ALGEBRAIC[58]*(CONSTANTS[79]+ALGEBRAIC[54]);
ALGEBRAIC[55] = 1.00000+ (STATES[1]/CONSTANTS[32])*(1.00000+exp(( - CONSTANTS[28]*STATES[0])/CONSTANTS[78])+CONSTANTS[4]/CONSTANTS[35])+ (CONSTANTS[4]/CONSTANTS[33])*(1.00000+ (CONSTANTS[4]/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[31]));
ALGEBRAIC[56] = ( (STATES[1]/CONSTANTS[32])*exp(( - CONSTANTS[28]*STATES[0])/CONSTANTS[78]))/ALGEBRAIC[55];
ALGEBRAIC[57] = ( (( (CONSTANTS[4]/CONSTANTS[33])*CONSTANTS[4])/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[31])*exp(( CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78])))/ALGEBRAIC[55];
ALGEBRAIC[59] = ALGEBRAIC[58]*CONSTANTS[79]*(ALGEBRAIC[57]+ALGEBRAIC[56])+ ALGEBRAIC[54]*ALGEBRAIC[56]*(CONSTANTS[80]+ALGEBRAIC[58]);
ALGEBRAIC[64] = ALGEBRAIC[57]*CONSTANTS[79]*(ALGEBRAIC[62]+ALGEBRAIC[61])+ ALGEBRAIC[56]*ALGEBRAIC[62]*(CONSTANTS[79]+ALGEBRAIC[54]);
ALGEBRAIC[65] = ALGEBRAIC[62]*CONSTANTS[80]*(ALGEBRAIC[57]+ALGEBRAIC[56])+ ALGEBRAIC[57]*ALGEBRAIC[61]*(CONSTANTS[80]+ALGEBRAIC[58]);
ALGEBRAIC[66] = ( CONSTANTS[0]*CONSTANTS[27]*( ALGEBRAIC[59]*ALGEBRAIC[61] - ALGEBRAIC[63]*ALGEBRAIC[56]))/(ALGEBRAIC[63]+ALGEBRAIC[59]+ALGEBRAIC[64]+ALGEBRAIC[65]);
ALGEBRAIC[53] = CONSTANTS[0]*CONSTANTS[26]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[51] = CONSTANTS[0]*CONSTANTS[22]*(STATES[0] - CONSTANTS[75]);
RATES[0] = - (ALGEBRAIC[14]+ALGEBRAIC[31]+ALGEBRAIC[49]+ALGEBRAIC[50]+ALGEBRAIC[37]+ALGEBRAIC[42]+ALGEBRAIC[44]+ALGEBRAIC[46]+ALGEBRAIC[52]+ALGEBRAIC[66]+ALGEBRAIC[53]+ALGEBRAIC[51])/CONSTANTS[0];
ALGEBRAIC[69] = CONSTANTS[61]*STATES[1]*(1.00000 - STATES[25]) - CONSTANTS[66]*STATES[25];
RATES[25] = ALGEBRAIC[69];
ALGEBRAIC[67] = CONSTANTS[41]*STATES[17]*(STATES[20] - STATES[1]);
ALGEBRAIC[68] = (STATES[1] - STATES[22])/CONSTANTS[50];
RATES[1] = ( ALGEBRAIC[67]*CONSTANTS[83])/CONSTANTS[82] - (((ALGEBRAIC[14]+ALGEBRAIC[31]+ALGEBRAIC[53]) - 2.00000*ALGEBRAIC[66])/( 2.00000*CONSTANTS[3]*CONSTANTS[82])+ALGEBRAIC[68]+ CONSTANTS[56]*ALGEBRAIC[69]);
ALGEBRAIC[71] = CONSTANTS[58]*STATES[22]*(1.00000 - STATES[26]) - CONSTANTS[63]*STATES[26];
RATES[26] = ALGEBRAIC[71];
ALGEBRAIC[70] = CONSTANTS[52]/(1.00000+CONSTANTS[53]/STATES[22]);
ALGEBRAIC[72] = (STATES[21] - STATES[20])/CONSTANTS[51];
RATES[21] = ALGEBRAIC[70] - ( ALGEBRAIC[72]*CONSTANTS[83])/CONSTANTS[85];
ALGEBRAIC[73] = CONSTANTS[60]*STATES[22]*(1.00000 - (STATES[27]+STATES[23])) - CONSTANTS[64]*STATES[27];
RATES[27] = ALGEBRAIC[73];
ALGEBRAIC[74] = CONSTANTS[62]*STATES[20]*(1.00000 - STATES[28]) - CONSTANTS[67]*STATES[28];
RATES[28] = ALGEBRAIC[74];
RATES[20] = ALGEBRAIC[72] - (ALGEBRAIC[67]+ CONSTANTS[57]*ALGEBRAIC[74]);
ALGEBRAIC[75] = CONSTANTS[61]*STATES[22]*(1.00000 - STATES[24]) - CONSTANTS[66]*STATES[24];
RATES[24] = ALGEBRAIC[75];
RATES[22] = ( ALGEBRAIC[68]*CONSTANTS[82] - ALGEBRAIC[70]*CONSTANTS[85])/CONSTANTS[84] - ( CONSTANTS[56]*ALGEBRAIC[75]+ CONSTANTS[54]*ALGEBRAIC[71]+ CONSTANTS[55]*ALGEBRAIC[73]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[15] = CONSTANTS[59]*CONSTANTS[68]*(1.00000 - (STATES[27]+STATES[23])) - CONSTANTS[65]*STATES[23];
ALGEBRAIC[1] = 1.00000/(1.00000+exp((STATES[0]+35.0000)/7.30000));
ALGEBRAIC[17] = 44.3000+ 257.100*exp(- pow((STATES[0]+32.5000)/13.9000, 2.00000));
ALGEBRAIC[2] = CONSTANTS[11]/(CONSTANTS[11]+STATES[1]);
ALGEBRAIC[18] = ALGEBRAIC[2]/CONSTANTS[10];
ALGEBRAIC[3] = 1.00000/(1.00000+exp(- (STATES[0]+26.3000)/6.00000));
ALGEBRAIC[19] = 1.00000/( 1.06800*exp((STATES[0]+26.3000)/30.0000)+ 1.06800*exp(- (STATES[0]+26.3000)/30.0000));
ALGEBRAIC[4] = 1.00000/(1.00000+exp((STATES[0]+61.7000)/5.60000));
ALGEBRAIC[20] = 1.00000/( 0.0153000*exp(- (STATES[0]+61.7000)/83.3000)+ 0.0150000*exp((STATES[0]+61.7000)/15.3800));
ALGEBRAIC[5] = 1.00000/(1.00000+exp(- (STATES[0]+23.2000)/10.6000));
ALGEBRAIC[21] = 0.846554/( 0.00420000*exp(STATES[0]/17.0000)+ 0.000150000*exp(- STATES[0]/21.6000));
ALGEBRAIC[22] = 0.846554/( 0.0372000*exp(STATES[0]/15.9000)+ 0.000960000*exp(- STATES[0]/22.5000));
ALGEBRAIC[6] = 1.00000/(1.00000+exp((STATES[0]+28.6000)/17.1000));
ALGEBRAIC[23] = 1.00000/( 0.100000*exp(- STATES[0]/54.6450)+ 0.656000*exp(STATES[0]/106.157));
ALGEBRAIC[8] = 1.00000/(1.00000+exp((STATES[0]+49.0000)/13.0000));
ALGEBRAIC[25] = 6.06000+39.1020/( 0.570000*exp( - 0.0800000*(STATES[0]+44.0000))+ 0.0650000*exp( 0.100000*(STATES[0]+45.9300)));
ALGEBRAIC[9] = 1.00000/(1.00000+exp(- (STATES[0] - 19.3000)/15.0000));
ALGEBRAIC[26] = 2.75352+14.4052/( 1.03700*exp( 0.0900000*(STATES[0]+30.6100))+ 0.369000*exp( - 0.120000*(STATES[0]+23.8400)));
ALGEBRAIC[10] = 1.00000/(1.00000+exp((STATES[0] - CONSTANTS[19])/13.5000));
ALGEBRAIC[27] = 0.716653/(exp(- (STATES[0]+386.900)/45.3020)+exp((STATES[0] - 73.0800)/19.2310));
ALGEBRAIC[13] = CONSTANTS[42] - (CONSTANTS[42] - CONSTANTS[43])/(1.00000+pow(CONSTANTS[44]/STATES[20], CONSTANTS[45]));
ALGEBRAIC[30] = CONSTANTS[46]/ALGEBRAIC[13];
ALGEBRAIC[36] = CONSTANTS[47]*ALGEBRAIC[13];
ALGEBRAIC[7] = 0.0140000/(1.00000+exp(- (STATES[0] - 40.0000)/9.00000));
ALGEBRAIC[24] = 0.00100000*exp(- STATES[0]/45.0000);
ALGEBRAIC[33] = ALGEBRAIC[7]/(ALGEBRAIC[7]+ALGEBRAIC[24]);
ALGEBRAIC[39] = 1.00000/(ALGEBRAIC[7]+ALGEBRAIC[24]);
ALGEBRAIC[11] = 1.00000/(1.00000+exp(- (STATES[0]+57.0000)/5.00000));
ALGEBRAIC[28] = 1.00000/( 0.150000*exp(- STATES[0]/11.0000)+ 0.200000*exp(- STATES[0]/700.000));
ALGEBRAIC[34] = 1.00000/( 16.0000*exp(STATES[0]/8.00000)+ 15.0000*exp(STATES[0]/50.0000));
ALGEBRAIC[40] = 1.00000/(ALGEBRAIC[28]+ALGEBRAIC[34]);
ALGEBRAIC[12] = 1.00000/( 3100.00*exp(STATES[0]/13.0000)+ 700.000*exp(STATES[0]/70.0000));
ALGEBRAIC[29] = 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[35] = ALGEBRAIC[12]/(ALGEBRAIC[12]+ALGEBRAIC[29]);
ALGEBRAIC[41] = 6.65000/(ALGEBRAIC[12]+ALGEBRAIC[29]);
ALGEBRAIC[0] = 1.00000/(1.00000+exp(- (STATES[0]+13.5000)/6.00000));
ALGEBRAIC[16] = (STATES[0]==- 35.0000 ? - 35.0000 : STATES[0]==0.00000 ? 1.00000e-05 : STATES[0]);
ALGEBRAIC[32] = ( - 0.0283900*(ALGEBRAIC[16]+35.0000))/(exp(- (ALGEBRAIC[16]+35.0000)/2.50000) - 1.00000) - ( 0.0849000*ALGEBRAIC[16])/(exp(- ALGEBRAIC[16]/4.80000) - 1.00000);
ALGEBRAIC[38] = (STATES[0]==5.00000 ? 5.00001 : STATES[0]);
ALGEBRAIC[43] = ( 0.0114300*(ALGEBRAIC[38] - 5.00000))/(exp((ALGEBRAIC[38] - 5.00000)/2.50000) - 1.00000);
ALGEBRAIC[45] = 1.00000/(ALGEBRAIC[32]+ALGEBRAIC[43]);
ALGEBRAIC[31] = CONSTANTS[0]*CONSTANTS[12]*(STATES[0] - CONSTANTS[13])*STATES[5]*STATES[6];
ALGEBRAIC[14] = CONSTANTS[0]*CONSTANTS[9]*(STATES[0] - CONSTANTS[8])*STATES[2]*STATES[3]*STATES[4];
ALGEBRAIC[47] = CONSTANTS[0]*0.383300*CONSTANTS[18]*(STATES[0] - CONSTANTS[75])*pow(STATES[13], 2.00000);
ALGEBRAIC[48] = CONSTANTS[0]*0.616700*CONSTANTS[18]*(STATES[0] - CONSTANTS[76])*pow(STATES[13], 2.00000);
ALGEBRAIC[49] = ALGEBRAIC[47]+ALGEBRAIC[48];
ALGEBRAIC[50] = CONSTANTS[0]*CONSTANTS[20]*(STATES[0] - CONSTANTS[21])*STATES[14]*STATES[15];
ALGEBRAIC[37] = CONSTANTS[0]*CONSTANTS[14]*(STATES[0] - CONSTANTS[76])*( 0.600000*STATES[8]+ 0.400000*STATES[7])*STATES[9];
ALGEBRAIC[42] = CONSTANTS[0]*CONSTANTS[15]*(STATES[0] - CONSTANTS[77])*pow(STATES[10], 2.00000);
ALGEBRAIC[44] = CONSTANTS[0]*CONSTANTS[16]*(STATES[0] - CONSTANTS[76])*STATES[11]*STATES[12];
ALGEBRAIC[46] = CONSTANTS[0]*CONSTANTS[17]*(STATES[0] - CONSTANTS[76])*STATES[12];
ALGEBRAIC[52] = ( CONSTANTS[0]*CONSTANTS[25])/( (1.00000+pow(CONSTANTS[23]/CONSTANTS[7], 1.20000))*(1.00000+pow(CONSTANTS[24]/CONSTANTS[4], 1.30000))*(1.00000+exp(- ((STATES[0] - CONSTANTS[75])+120.000)/30.0000)));
ALGEBRAIC[54] = exp(( - CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78]));
ALGEBRAIC[60] = 1.00000+ (CONSTANTS[40]/CONSTANTS[39])*(1.00000+exp(( CONSTANTS[30]*STATES[0])/CONSTANTS[78]))+ (CONSTANTS[5]/CONSTANTS[37])*(1.00000+ (CONSTANTS[5]/CONSTANTS[38])*(1.00000+CONSTANTS[5]/CONSTANTS[36]));
ALGEBRAIC[62] = ( (( (CONSTANTS[5]/CONSTANTS[37])*CONSTANTS[5])/CONSTANTS[38])*(1.00000+CONSTANTS[5]/CONSTANTS[36])*exp(( - CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78])))/ALGEBRAIC[60];
ALGEBRAIC[61] = ( (CONSTANTS[40]/CONSTANTS[39])*exp(( CONSTANTS[30]*STATES[0])/CONSTANTS[78]))/ALGEBRAIC[60];
ALGEBRAIC[58] = exp(( CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78]));
ALGEBRAIC[63] = ALGEBRAIC[54]*CONSTANTS[80]*(ALGEBRAIC[62]+ALGEBRAIC[61])+ ALGEBRAIC[61]*ALGEBRAIC[58]*(CONSTANTS[79]+ALGEBRAIC[54]);
ALGEBRAIC[55] = 1.00000+ (STATES[1]/CONSTANTS[32])*(1.00000+exp(( - CONSTANTS[28]*STATES[0])/CONSTANTS[78])+CONSTANTS[4]/CONSTANTS[35])+ (CONSTANTS[4]/CONSTANTS[33])*(1.00000+ (CONSTANTS[4]/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[31]));
ALGEBRAIC[56] = ( (STATES[1]/CONSTANTS[32])*exp(( - CONSTANTS[28]*STATES[0])/CONSTANTS[78]))/ALGEBRAIC[55];
ALGEBRAIC[57] = ( (( (CONSTANTS[4]/CONSTANTS[33])*CONSTANTS[4])/CONSTANTS[34])*(1.00000+CONSTANTS[4]/CONSTANTS[31])*exp(( CONSTANTS[29]*STATES[0])/( 2.00000*CONSTANTS[78])))/ALGEBRAIC[55];
ALGEBRAIC[59] = ALGEBRAIC[58]*CONSTANTS[79]*(ALGEBRAIC[57]+ALGEBRAIC[56])+ ALGEBRAIC[54]*ALGEBRAIC[56]*(CONSTANTS[80]+ALGEBRAIC[58]);
ALGEBRAIC[64] = ALGEBRAIC[57]*CONSTANTS[79]*(ALGEBRAIC[62]+ALGEBRAIC[61])+ ALGEBRAIC[56]*ALGEBRAIC[62]*(CONSTANTS[79]+ALGEBRAIC[54]);
ALGEBRAIC[65] = ALGEBRAIC[62]*CONSTANTS[80]*(ALGEBRAIC[57]+ALGEBRAIC[56])+ ALGEBRAIC[57]*ALGEBRAIC[61]*(CONSTANTS[80]+ALGEBRAIC[58]);
ALGEBRAIC[66] = ( CONSTANTS[0]*CONSTANTS[27]*( ALGEBRAIC[59]*ALGEBRAIC[61] - ALGEBRAIC[63]*ALGEBRAIC[56]))/(ALGEBRAIC[63]+ALGEBRAIC[59]+ALGEBRAIC[64]+ALGEBRAIC[65]);
ALGEBRAIC[53] = CONSTANTS[0]*CONSTANTS[26]*(STATES[0] - CONSTANTS[8]);
ALGEBRAIC[51] = CONSTANTS[0]*CONSTANTS[22]*(STATES[0] - CONSTANTS[75]);
ALGEBRAIC[69] = CONSTANTS[61]*STATES[1]*(1.00000 - STATES[25]) - CONSTANTS[66]*STATES[25];
ALGEBRAIC[67] = CONSTANTS[41]*STATES[17]*(STATES[20] - STATES[1]);
ALGEBRAIC[68] = (STATES[1] - STATES[22])/CONSTANTS[50];
ALGEBRAIC[71] = CONSTANTS[58]*STATES[22]*(1.00000 - STATES[26]) - CONSTANTS[63]*STATES[26];
ALGEBRAIC[70] = CONSTANTS[52]/(1.00000+CONSTANTS[53]/STATES[22]);
ALGEBRAIC[72] = (STATES[21] - STATES[20])/CONSTANTS[51];
ALGEBRAIC[73] = CONSTANTS[60]*STATES[22]*(1.00000 - (STATES[27]+STATES[23])) - CONSTANTS[64]*STATES[27];
ALGEBRAIC[74] = CONSTANTS[62]*STATES[20]*(1.00000 - STATES[28]) - CONSTANTS[67]*STATES[28];
ALGEBRAIC[75] = CONSTANTS[61]*STATES[22]*(1.00000 - STATES[24]) - CONSTANTS[66]*STATES[24];
}