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 56 entries in the algebraic variable array.
   There are a total of 14 entries in each of the rate and state variable arrays.
   There are a total of 43 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (millisecond).
 * STATES[0] is Vm_n in component Vm_n (millivolt).
 * ALGEBRAIC[0] is i_app in component Vm_n (nanoA_per_cm2).
 * CONSTANTS[0] is Cm in component Vm_n (microF_per_cm2).
 * CONSTANTS[1] is t0 in component Vm_n (millisecond).
 * CONSTANTS[2] is t1 in component Vm_n (millisecond).
 * ALGEBRAIC[36] is i_NaT in component i_NaT (nanoA_per_cm2).
 * ALGEBRAIC[38] is i_NaP in component i_NaP (nanoA_per_cm2).
 * ALGEBRAIC[40] is i_leakNa in component i_leakNa (nanoA_per_cm2).
 * ALGEBRAIC[25] is i_KDR in component i_KDR (nanoA_per_cm2).
 * ALGEBRAIC[26] is i_KA in component i_KA (nanoA_per_cm2).
 * ALGEBRAIC[27] is i_leakK in component i_leakK (nanoA_per_cm2).
 * ALGEBRAIC[8] is i_leakf in component i_leakf (nanoA_per_cm2).
 * ALGEBRAIC[34] is i_NaKATPase_n in component i_NaKATPase_n (nanoA_per_cm2).
 * ALGEBRAIC[37] is J_NaT in component i_NaT (picomole_per_cm2_millisecond).
 * CONSTANTS[3] is gNaT in component i_NaT (microS_per_cm2).
 * CONSTANTS[4] is F in component model_parameters (coulomb_per_mole).
 * ALGEBRAIC[35] is ENa_n in component electric_potentials (millivolt).
 * STATES[1] is m in component i_NaT_m_gate (dimensionless).
 * STATES[2] is h in component i_NaT_h_gate (dimensionless).
 * ALGEBRAIC[1] is alpha_m in component i_NaT_m_gate (per_millisecond).
 * ALGEBRAIC[9] is beta_m in component i_NaT_m_gate (per_millisecond).
 * ALGEBRAIC[2] is alpha_h in component i_NaT_h_gate (per_millisecond).
 * ALGEBRAIC[10] is beta_h in component i_NaT_h_gate (per_millisecond).
 * ALGEBRAIC[39] is J_NaP in component i_NaP (picomole_per_cm2_millisecond).
 * CONSTANTS[5] is gNaP in component i_NaP (microS_per_cm2).
 * STATES[3] is m in component i_NaP_m_gate (dimensionless).
 * STATES[4] is h in component i_NaP_h_gate (dimensionless).
 * ALGEBRAIC[3] is alpha_m in component i_NaP_m_gate (per_millisecond).
 * ALGEBRAIC[11] is beta_m in component i_NaP_m_gate (per_millisecond).
 * CONSTANTS[6] is tau_activation in component i_NaP_m_gate (millisecond).
 * ALGEBRAIC[4] is alpha_h in component i_NaP_h_gate (per_millisecond).
 * ALGEBRAIC[12] is beta_h in component i_NaP_h_gate (per_millisecond).
 * ALGEBRAIC[22] is J_KDR in component i_KDR (picomole_per_cm2_millisecond).
 * CONSTANTS[7] is gKDR in component i_KDR (microS_per_cm2).
 * ALGEBRAIC[21] is EK_n in component electric_potentials (millivolt).
 * STATES[5] is n in component i_KDR_n_gate (dimensionless).
 * ALGEBRAIC[5] is alpha_n in component i_KDR_n_gate (per_millisecond).
 * ALGEBRAIC[13] is beta_n in component i_KDR_n_gate (per_millisecond).
 * ALGEBRAIC[23] is J_KA in component i_KA (picomole_per_cm2_millisecond).
 * CONSTANTS[8] is gKA in component i_KA (microS_per_cm2).
 * STATES[6] is m in component i_KA_m_gate (dimensionless).
 * STATES[7] is h in component i_KA_h_gate (dimensionless).
 * ALGEBRAIC[6] is alpha_m in component i_KA_m_gate (per_millisecond).
 * ALGEBRAIC[14] is beta_m in component i_KA_m_gate (per_millisecond).
 * ALGEBRAIC[7] is alpha_h in component i_KA_h_gate (per_millisecond).
 * ALGEBRAIC[15] is beta_h in component i_KA_h_gate (per_millisecond).
 * ALGEBRAIC[33] is J_NaKATPase_n in component i_NaKATPase_n (picomole_per_cm2_millisecond).
 * CONSTANTS[9] is I_NaKATPase_n_max in component i_NaKATPase_n (picomole_per_cm2_millisecond).
 * CONSTANTS[10] is KmNa in component model_parameters (millimolar).
 * CONSTANTS[11] is KmK in component model_parameters (millimolar).
 * ALGEBRAIC[32] is Nan in component ion_concentrations (millimolar).
 * ALGEBRAIC[18] is Ko in component ion_concentrations (millimolar).
 * ALGEBRAIC[41] is J_leakNa in component i_leakNa (picomole_per_cm2_millisecond).
 * CONSTANTS[12] is gleakNa in component i_leakNa (microS_per_cm2).
 * ALGEBRAIC[24] is J_leakK in component i_leakK (picomole_per_cm2_millisecond).
 * CONSTANTS[13] is gleakK in component i_leakK (microS_per_cm2).
 * CONSTANTS[14] is gleakf in component i_leakf (microS_per_cm2).
 * CONSTANTS[15] is Ef_n in component electric_potentials (millivolt).
 * ALGEBRAIC[50] is Vm_g in component Vm_g (millivolt).
 * ALGEBRAIC[31] is ENa_g in component electric_potentials (millivolt).
 * ALGEBRAIC[19] is EK_g in component electric_potentials (millivolt).
 * ALGEBRAIC[49] is ECl_g in component electric_potentials (millivolt).
 * ALGEBRAIC[45] is ENBC_g in component electric_potentials (millivolt).
 * CONSTANTS[16] is gNa in component J_Na (microS_per_cm2).
 * CONSTANTS[17] is gK in component J_K (microS_per_cm2).
 * CONSTANTS[18] is gCl in component model_parameters (microS_per_cm2).
 * CONSTANTS[19] is gNBC in component J_NBC (microS_per_cm2).
 * ALGEBRAIC[29] is J_NaKATPase_g in component J_NaKATPase_g (picomole_per_cm2_millisecond).
 * ALGEBRAIC[51] is J_Na in component J_Na (picomole_per_cm2_millisecond).
 * ALGEBRAIC[52] is J_K in component J_K (picomole_per_cm2_millisecond).
 * CONSTANTS[20] is I_NaKATPase_g_max in component J_NaKATPase_g (picomole_per_cm2_millisecond).
 * ALGEBRAIC[28] is Nag in component ion_concentrations (millimolar).
 * ALGEBRAIC[53] is J_NBC in component J_NBC (picomole_per_cm2_millisecond).
 * ALGEBRAIC[48] is J_NKCC1 in component J_NKCC1 (picomole_per_cm2_millisecond).
 * CONSTANTS[21] is gNKCC1 in component J_NKCC1 (microS_per_cm2).
 * CONSTANTS[40] is Psi in component model_parameters (millivolt).
 * ALGEBRAIC[30] is Nao in component ion_concentrations (millimolar).
 * CONSTANTS[22] is P_Ko in component ion_concentrations (millimolar).
 * ALGEBRAIC[17] is Kg in component ion_concentrations (millimolar).
 * ALGEBRAIC[47] is Clo in component ion_concentrations (millimolar).
 * ALGEBRAIC[46] is Clg in component ion_concentrations (millimolar).
 * STATES[8] is N_Nag in component N_Nag (nanomole_per_cm2).
 * ALGEBRAIC[55] is dN_Nag_dt in component N_Nag (nanomole_per_cm2_millisecond).
 * STATES[9] is N_Kg in component N_Kg (nanomole_per_cm2).
 * ALGEBRAIC[54] is dN_Kg_dt in component N_Kg (nanomole_per_cm2_millisecond).
 * STATES[10] is wg in component wg (micrometre).
 * CONSTANTS[23] is Lp in component model_parameters (cm_per_second_millimolar).
 * CONSTANTS[24] is Xg in component model_parameters (millimolar_micrometre).
 * ALGEBRAIC[43] is HCO3o in component ion_concentrations (millimolar).
 * ALGEBRAIC[44] is HCO3g in component ion_concentrations (millimolar).
 * ALGEBRAIC[16] is wo in component wo (micrometre).
 * CONSTANTS[25] is P_wo in component model_parameters (micrometre).
 * CONSTANTS[26] is P_wg in component model_parameters (micrometre).
 * STATES[11] is N_Nao in component N_Nao (nanomole_per_cm2).
 * STATES[12] is N_Ko in component N_Ko (nanomole_per_cm2).
 * STATES[13] is N_HCO3o in component N_HCO3o (nanomole_per_cm2).
 * ALGEBRAIC[20] is Kn in component ion_concentrations (millimolar).
 * ALGEBRAIC[42] is Cln in component ion_concentrations (millimolar).
 * CONSTANTS[27] is P_Vm_g in component ion_concentrations (millivolt).
 * CONSTANTS[39] is P_Clo in component ion_concentrations (millimolar).
 * CONSTANTS[28] is P_Cln in component ion_concentrations (millimolar).
 * CONSTANTS[41] is P_Clg in component ion_concentrations (millimolar).
 * CONSTANTS[29] is P_Nan in component ion_concentrations (millimolar).
 * CONSTANTS[30] is P_Nao in component ion_concentrations (millimolar).
 * CONSTANTS[31] is P_Nag in component ion_concentrations (millimolar).
 * CONSTANTS[32] is P_Kn in component ion_concentrations (millimolar).
 * CONSTANTS[33] is P_Kg in component ion_concentrations (millimolar).
 * CONSTANTS[42] is P_HCO3g in component ion_concentrations (millimolar).
 * CONSTANTS[34] is P_HCO3o in component ion_concentrations (millimolar).
 * CONSTANTS[35] is P_wn in component ion_concentrations (micrometre).
 * CONSTANTS[36] is rho in component model_parameters (dimensionless).
 * CONSTANTS[37] is R in component model_parameters (joule_per_mole_kelvin).
 * CONSTANTS[38] is T in component model_parameters (kelvin).
 * RATES[0] is d/dt Vm_n in component Vm_n (millivolt).
 * RATES[1] is d/dt m in component i_NaT_m_gate (dimensionless).
 * RATES[2] is d/dt h in component i_NaT_h_gate (dimensionless).
 * RATES[3] is d/dt m in component i_NaP_m_gate (dimensionless).
 * RATES[4] is d/dt h in component i_NaP_h_gate (dimensionless).
 * RATES[5] is d/dt n in component i_KDR_n_gate (dimensionless).
 * RATES[6] is d/dt m in component i_KA_m_gate (dimensionless).
 * RATES[7] is d/dt h in component i_KA_h_gate (dimensionless).
 * RATES[8] is d/dt N_Nag in component N_Nag (nanomole_per_cm2).
 * RATES[9] is d/dt N_Kg in component N_Kg (nanomole_per_cm2).
 * RATES[10] is d/dt wg in component wg (micrometre).
 * RATES[11] is d/dt N_Nao in component N_Nao (nanomole_per_cm2).
 * RATES[12] is d/dt N_Ko in component N_Ko (nanomole_per_cm2).
 * RATES[13] is d/dt N_HCO3o in component N_HCO3o (nanomole_per_cm2).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -70;
CONSTANTS[0] = 1;
CONSTANTS[1] = 100;
CONSTANTS[2] = 700;
CONSTANTS[3] = 5000;
CONSTANTS[4] = 9.649e4;
STATES[1] = 0.005;
STATES[2] = 0.9961;
CONSTANTS[5] = 150;
STATES[3] = 0.0129;
STATES[4] = 0.9718;
CONSTANTS[6] = 6;
CONSTANTS[7] = 5000;
STATES[5] = 0.0012;
CONSTANTS[8] = 1000;
STATES[6] = 0.1193;
STATES[7] = 0.1205;
CONSTANTS[9] = 0.0289;
CONSTANTS[10] = 10;
CONSTANTS[11] = 1.5;
CONSTANTS[12] = 20;
CONSTANTS[13] = 66.06;
CONSTANTS[14] = 10;
CONSTANTS[15] = -70;
CONSTANTS[16] = 100.0;
CONSTANTS[17] = 1696.0;
CONSTANTS[18] = 50;
CONSTANTS[19] = 80;
CONSTANTS[20] = 0.1151;
CONSTANTS[21] = 2;
CONSTANTS[22] = 3;
STATES[8] = 0.75;
STATES[9] = 5;
STATES[10] = 0.05;
CONSTANTS[23] = 2e-8;
CONSTANTS[24] = 8.45;
CONSTANTS[25] = 0.025;
CONSTANTS[26] = 0.05;
STATES[11] = 3.65;
STATES[12] = 0.075;
STATES[13] = 0.375;
CONSTANTS[27] = -85;
CONSTANTS[28] = 50;
CONSTANTS[29] = 10;
CONSTANTS[30] = 146;
CONSTANTS[31] = 15;
CONSTANTS[32] = 130;
CONSTANTS[33] = 100;
CONSTANTS[34] = 15;
CONSTANTS[35] = 0.05;
CONSTANTS[36] = 0.5975;
CONSTANTS[37] = 8.315;
CONSTANTS[38] = 300;
CONSTANTS[39] = CONSTANTS[30]+CONSTANTS[22] - CONSTANTS[34];
CONSTANTS[40] =  1000.00*CONSTANTS[37]*CONSTANTS[38]/CONSTANTS[4];
CONSTANTS[41] =  CONSTANTS[39]*exp(CONSTANTS[27]/CONSTANTS[40]);
CONSTANTS[42] =  CONSTANTS[34]* pow( CONSTANTS[30]/CONSTANTS[31]*exp(CONSTANTS[27]/CONSTANTS[40]), 1.0 / 2);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] =  0.320000*(- STATES[0] - 51.9000)/(exp(- ( 0.250000*STATES[0]+12.9750)) - 1.00000);
ALGEBRAIC[9] =  0.280000*(STATES[0]+24.8900)/(exp( 0.200000*STATES[0]+4.97800) - 1.00000);
RATES[1] =  ALGEBRAIC[1]*(1.00000 - STATES[1]) -  ALGEBRAIC[9]*STATES[1];
ALGEBRAIC[2] =  0.128000*exp(- ( 0.0560000*STATES[0]+2.94000));
ALGEBRAIC[10] = 4.00000/(exp(- ( 0.200000*STATES[0]+6.00000))+1.00000);
RATES[2] =  ALGEBRAIC[2]*(1.00000 - STATES[2]) -  ALGEBRAIC[10]*STATES[2];
ALGEBRAIC[3] =  1.00000/CONSTANTS[6]*1.00000/(exp(- ( 0.143000*STATES[0]+5.67000))+1.00000);
ALGEBRAIC[11] =  1.00000/CONSTANTS[6]*exp(- ( 0.143000*STATES[0]+5.67000))/(exp(- ( 0.143000*STATES[0]+5.67000))+1.00000);
RATES[3] =  ALGEBRAIC[3]*(1.00000 - STATES[3]) -  ALGEBRAIC[11]*STATES[3];
ALGEBRAIC[4] =  5.12000e-08*exp(- ( 0.0560000*STATES[0]+2.94000));
ALGEBRAIC[12] = 1.60000e-06/(exp(- ( 0.200000*STATES[0]+8.00000))+1.00000);
RATES[4] =  ALGEBRAIC[4]*(1.00000 - STATES[4]) -  ALGEBRAIC[12]*STATES[4];
ALGEBRAIC[5] =  0.0160000*(- STATES[0] - 34.9000)/(exp(- ( 0.200000*STATES[0]+6.98000)) - 1.00000);
ALGEBRAIC[13] =  0.250000*exp(- ( 0.0250000*STATES[0]+1.25000));
RATES[5] =  ALGEBRAIC[5]*(1.00000 - STATES[5]) -  ALGEBRAIC[13]*STATES[5];
ALGEBRAIC[6] =  0.0200000*(- STATES[0] - 56.9000)/(exp(- ( 0.100000*STATES[0]+5.69000)) - 1.00000);
ALGEBRAIC[14] =  0.0175000*(STATES[0]+29.9000)/(exp( 0.100000*STATES[0]+2.99000) - 1.00000);
RATES[6] =  ALGEBRAIC[6]*(1.00000 - STATES[6]) -  ALGEBRAIC[14]*STATES[6];
ALGEBRAIC[7] =  0.0160000*exp(- ( 0.0560000*STATES[0]+4.61000));
ALGEBRAIC[15] = 0.500000/(exp(- ( 0.200000*STATES[0]+11.9800))+1.00000);
RATES[7] =  ALGEBRAIC[7]*(1.00000 - STATES[7]) -  ALGEBRAIC[15]*STATES[7];
ALGEBRAIC[0] = (VOI>=CONSTANTS[1]&&VOI<=CONSTANTS[2] ? 3500.00 : 0.00000);
ALGEBRAIC[32] = ( CONSTANTS[35]*CONSTANTS[29]+ CONSTANTS[26]*CONSTANTS[31]+ CONSTANTS[25]*CONSTANTS[30] - (STATES[11]+STATES[8]))/CONSTANTS[35];
ALGEBRAIC[16] = CONSTANTS[26]+CONSTANTS[25] - STATES[10];
ALGEBRAIC[30] = STATES[11]/ALGEBRAIC[16];
ALGEBRAIC[35] =  CONSTANTS[40]*log(ALGEBRAIC[30]/ALGEBRAIC[32]);
ALGEBRAIC[36] =  CONSTANTS[3]*pow(STATES[1], 3.00000)*STATES[2]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[38] =  CONSTANTS[5]*pow(STATES[3], 2.00000)*STATES[4]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[40] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[18] = STATES[12]/ALGEBRAIC[16];
ALGEBRAIC[20] = ( CONSTANTS[35]*CONSTANTS[32]+ CONSTANTS[26]*CONSTANTS[33]+ CONSTANTS[25]*CONSTANTS[22] - (STATES[12]+STATES[9]))/CONSTANTS[35];
ALGEBRAIC[21] =  CONSTANTS[40]*log(ALGEBRAIC[18]/ALGEBRAIC[20]);
ALGEBRAIC[25] =  CONSTANTS[7]*pow(STATES[5], 2.00000)*(STATES[0] - ALGEBRAIC[21]);
ALGEBRAIC[26] =  CONSTANTS[8]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[21]);
ALGEBRAIC[27] =  CONSTANTS[13]*(STATES[0] - ALGEBRAIC[21]);
ALGEBRAIC[8] =  CONSTANTS[14]*(STATES[0] - CONSTANTS[15]);
ALGEBRAIC[33] =   CONSTANTS[9]*pow(ALGEBRAIC[32], 1.50000)/(pow(ALGEBRAIC[32], 1.50000)+pow(CONSTANTS[10], 1.50000))*ALGEBRAIC[18]/(ALGEBRAIC[18]+CONSTANTS[11]);
ALGEBRAIC[34] =  ALGEBRAIC[33]*CONSTANTS[4];
RATES[0] =  0.00100000*(- (ALGEBRAIC[36]+ALGEBRAIC[38]+ALGEBRAIC[40]+ALGEBRAIC[25]+ALGEBRAIC[26]+ALGEBRAIC[27]+ALGEBRAIC[8]+ALGEBRAIC[34])+ALGEBRAIC[0])/CONSTANTS[0];
ALGEBRAIC[28] = STATES[8]/STATES[10];
ALGEBRAIC[17] = STATES[9]/STATES[10];
ALGEBRAIC[44] = ( CONSTANTS[26]*CONSTANTS[42]+ CONSTANTS[25]*CONSTANTS[34] - STATES[13])/STATES[10];
ALGEBRAIC[46] = ALGEBRAIC[28]+ALGEBRAIC[17] - (ALGEBRAIC[44]+ CONSTANTS[36]*CONSTANTS[24]/STATES[10]);
ALGEBRAIC[42] = CONSTANTS[28]+ALGEBRAIC[32] - CONSTANTS[29]+(ALGEBRAIC[20] - CONSTANTS[32]);
ALGEBRAIC[47] = ( CONSTANTS[26]*CONSTANTS[41]+ CONSTANTS[25]*CONSTANTS[39]+ CONSTANTS[35]*CONSTANTS[28] - ( STATES[10]*ALGEBRAIC[46]+ CONSTANTS[35]*ALGEBRAIC[42]))/ALGEBRAIC[16];
ALGEBRAIC[43] = STATES[13]/ALGEBRAIC[16];
RATES[10] =  10.0000*CONSTANTS[23]*(ALGEBRAIC[28]+ALGEBRAIC[17]+ALGEBRAIC[46]+ALGEBRAIC[44]+CONSTANTS[24]/STATES[10] - (ALGEBRAIC[30]+ALGEBRAIC[18]+ALGEBRAIC[47]+ALGEBRAIC[43]));
ALGEBRAIC[29] =   CONSTANTS[20]*pow(ALGEBRAIC[28], 1.50000)/(pow(ALGEBRAIC[28], 1.50000)+pow(CONSTANTS[10], 1.50000))*ALGEBRAIC[18]/(ALGEBRAIC[18]+CONSTANTS[11]);
ALGEBRAIC[31] =  CONSTANTS[40]*log(ALGEBRAIC[30]/ALGEBRAIC[28]);
ALGEBRAIC[19] =  CONSTANTS[40]*log(ALGEBRAIC[18]/ALGEBRAIC[17]);
ALGEBRAIC[49] =  - 1.00000*CONSTANTS[40]*log(ALGEBRAIC[47]/ALGEBRAIC[46]);
ALGEBRAIC[45] =  - CONSTANTS[40]*log( ALGEBRAIC[30]/ALGEBRAIC[28]*pow(ALGEBRAIC[43]/ALGEBRAIC[44], 2.00000));
ALGEBRAIC[50] = ( CONSTANTS[16]*ALGEBRAIC[31]+ CONSTANTS[17]*ALGEBRAIC[19]+ CONSTANTS[18]*ALGEBRAIC[49]+ CONSTANTS[19]*ALGEBRAIC[45])/(CONSTANTS[16]+CONSTANTS[17]+CONSTANTS[18]+CONSTANTS[19]) -  CONSTANTS[4]*ALGEBRAIC[29]/(CONSTANTS[16]+CONSTANTS[17]+CONSTANTS[18]+CONSTANTS[19]);
ALGEBRAIC[52] =  CONSTANTS[17]*(ALGEBRAIC[50] - ALGEBRAIC[19])/CONSTANTS[4];
ALGEBRAIC[48] = (ALGEBRAIC[18]>CONSTANTS[22] ?  pow(ALGEBRAIC[18] - CONSTANTS[22], 10.0000)/(pow(ALGEBRAIC[18] - CONSTANTS[22], 10.0000)+pow(0.0300000, 10.0000))*CONSTANTS[40]*CONSTANTS[21]*log(  ALGEBRAIC[18]/ALGEBRAIC[17]*ALGEBRAIC[30]/ALGEBRAIC[28]*pow(ALGEBRAIC[47]/ALGEBRAIC[46], 2.00000))/CONSTANTS[4] : 0.00000);
ALGEBRAIC[54] =  0.0100000*(- ALGEBRAIC[52]+ 2.00000*ALGEBRAIC[29]+ALGEBRAIC[48]);
RATES[9] = ALGEBRAIC[54];
ALGEBRAIC[37] =  CONSTANTS[3]*pow(STATES[1], 3.00000)*STATES[2]*(STATES[0] - ALGEBRAIC[35])/CONSTANTS[4];
ALGEBRAIC[39] =  CONSTANTS[5]*pow(STATES[3], 2.00000)*STATES[4]*(STATES[0] - ALGEBRAIC[35])/CONSTANTS[4];
ALGEBRAIC[41] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[35])/CONSTANTS[4];
RATES[12] =  -0.0100000*(ALGEBRAIC[37]+ALGEBRAIC[39]+ALGEBRAIC[41]+ 3.00000*ALGEBRAIC[33]) - ALGEBRAIC[54];
ALGEBRAIC[53] =  CONSTANTS[19]*(ALGEBRAIC[50] - ALGEBRAIC[45])/CONSTANTS[4];
RATES[13] =  - 0.0100000*2.00000*ALGEBRAIC[53];
ALGEBRAIC[51] =  CONSTANTS[16]*(ALGEBRAIC[50] - ALGEBRAIC[31])/CONSTANTS[4];
ALGEBRAIC[55] =  0.0100000*(- ALGEBRAIC[51] -  3.00000*ALGEBRAIC[29]+ALGEBRAIC[48]+ALGEBRAIC[53]);
RATES[8] = ALGEBRAIC[55];
RATES[11] =  0.0100000*( 3.00000*ALGEBRAIC[33]+ALGEBRAIC[37]+ALGEBRAIC[39]+ALGEBRAIC[41] -  100.000*ALGEBRAIC[55]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] =  0.320000*(- STATES[0] - 51.9000)/(exp(- ( 0.250000*STATES[0]+12.9750)) - 1.00000);
ALGEBRAIC[9] =  0.280000*(STATES[0]+24.8900)/(exp( 0.200000*STATES[0]+4.97800) - 1.00000);
ALGEBRAIC[2] =  0.128000*exp(- ( 0.0560000*STATES[0]+2.94000));
ALGEBRAIC[10] = 4.00000/(exp(- ( 0.200000*STATES[0]+6.00000))+1.00000);
ALGEBRAIC[3] =  1.00000/CONSTANTS[6]*1.00000/(exp(- ( 0.143000*STATES[0]+5.67000))+1.00000);
ALGEBRAIC[11] =  1.00000/CONSTANTS[6]*exp(- ( 0.143000*STATES[0]+5.67000))/(exp(- ( 0.143000*STATES[0]+5.67000))+1.00000);
ALGEBRAIC[4] =  5.12000e-08*exp(- ( 0.0560000*STATES[0]+2.94000));
ALGEBRAIC[12] = 1.60000e-06/(exp(- ( 0.200000*STATES[0]+8.00000))+1.00000);
ALGEBRAIC[5] =  0.0160000*(- STATES[0] - 34.9000)/(exp(- ( 0.200000*STATES[0]+6.98000)) - 1.00000);
ALGEBRAIC[13] =  0.250000*exp(- ( 0.0250000*STATES[0]+1.25000));
ALGEBRAIC[6] =  0.0200000*(- STATES[0] - 56.9000)/(exp(- ( 0.100000*STATES[0]+5.69000)) - 1.00000);
ALGEBRAIC[14] =  0.0175000*(STATES[0]+29.9000)/(exp( 0.100000*STATES[0]+2.99000) - 1.00000);
ALGEBRAIC[7] =  0.0160000*exp(- ( 0.0560000*STATES[0]+4.61000));
ALGEBRAIC[15] = 0.500000/(exp(- ( 0.200000*STATES[0]+11.9800))+1.00000);
ALGEBRAIC[0] = (VOI>=CONSTANTS[1]&&VOI<=CONSTANTS[2] ? 3500.00 : 0.00000);
ALGEBRAIC[32] = ( CONSTANTS[35]*CONSTANTS[29]+ CONSTANTS[26]*CONSTANTS[31]+ CONSTANTS[25]*CONSTANTS[30] - (STATES[11]+STATES[8]))/CONSTANTS[35];
ALGEBRAIC[16] = CONSTANTS[26]+CONSTANTS[25] - STATES[10];
ALGEBRAIC[30] = STATES[11]/ALGEBRAIC[16];
ALGEBRAIC[35] =  CONSTANTS[40]*log(ALGEBRAIC[30]/ALGEBRAIC[32]);
ALGEBRAIC[36] =  CONSTANTS[3]*pow(STATES[1], 3.00000)*STATES[2]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[38] =  CONSTANTS[5]*pow(STATES[3], 2.00000)*STATES[4]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[40] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[35]);
ALGEBRAIC[18] = STATES[12]/ALGEBRAIC[16];
ALGEBRAIC[20] = ( CONSTANTS[35]*CONSTANTS[32]+ CONSTANTS[26]*CONSTANTS[33]+ CONSTANTS[25]*CONSTANTS[22] - (STATES[12]+STATES[9]))/CONSTANTS[35];
ALGEBRAIC[21] =  CONSTANTS[40]*log(ALGEBRAIC[18]/ALGEBRAIC[20]);
ALGEBRAIC[25] =  CONSTANTS[7]*pow(STATES[5], 2.00000)*(STATES[0] - ALGEBRAIC[21]);
ALGEBRAIC[26] =  CONSTANTS[8]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[21]);
ALGEBRAIC[27] =  CONSTANTS[13]*(STATES[0] - ALGEBRAIC[21]);
ALGEBRAIC[8] =  CONSTANTS[14]*(STATES[0] - CONSTANTS[15]);
ALGEBRAIC[33] =   CONSTANTS[9]*pow(ALGEBRAIC[32], 1.50000)/(pow(ALGEBRAIC[32], 1.50000)+pow(CONSTANTS[10], 1.50000))*ALGEBRAIC[18]/(ALGEBRAIC[18]+CONSTANTS[11]);
ALGEBRAIC[34] =  ALGEBRAIC[33]*CONSTANTS[4];
ALGEBRAIC[28] = STATES[8]/STATES[10];
ALGEBRAIC[17] = STATES[9]/STATES[10];
ALGEBRAIC[44] = ( CONSTANTS[26]*CONSTANTS[42]+ CONSTANTS[25]*CONSTANTS[34] - STATES[13])/STATES[10];
ALGEBRAIC[46] = ALGEBRAIC[28]+ALGEBRAIC[17] - (ALGEBRAIC[44]+ CONSTANTS[36]*CONSTANTS[24]/STATES[10]);
ALGEBRAIC[42] = CONSTANTS[28]+ALGEBRAIC[32] - CONSTANTS[29]+(ALGEBRAIC[20] - CONSTANTS[32]);
ALGEBRAIC[47] = ( CONSTANTS[26]*CONSTANTS[41]+ CONSTANTS[25]*CONSTANTS[39]+ CONSTANTS[35]*CONSTANTS[28] - ( STATES[10]*ALGEBRAIC[46]+ CONSTANTS[35]*ALGEBRAIC[42]))/ALGEBRAIC[16];
ALGEBRAIC[43] = STATES[13]/ALGEBRAIC[16];
ALGEBRAIC[29] =   CONSTANTS[20]*pow(ALGEBRAIC[28], 1.50000)/(pow(ALGEBRAIC[28], 1.50000)+pow(CONSTANTS[10], 1.50000))*ALGEBRAIC[18]/(ALGEBRAIC[18]+CONSTANTS[11]);
ALGEBRAIC[31] =  CONSTANTS[40]*log(ALGEBRAIC[30]/ALGEBRAIC[28]);
ALGEBRAIC[19] =  CONSTANTS[40]*log(ALGEBRAIC[18]/ALGEBRAIC[17]);
ALGEBRAIC[49] =  - 1.00000*CONSTANTS[40]*log(ALGEBRAIC[47]/ALGEBRAIC[46]);
ALGEBRAIC[45] =  - CONSTANTS[40]*log( ALGEBRAIC[30]/ALGEBRAIC[28]*pow(ALGEBRAIC[43]/ALGEBRAIC[44], 2.00000));
ALGEBRAIC[50] = ( CONSTANTS[16]*ALGEBRAIC[31]+ CONSTANTS[17]*ALGEBRAIC[19]+ CONSTANTS[18]*ALGEBRAIC[49]+ CONSTANTS[19]*ALGEBRAIC[45])/(CONSTANTS[16]+CONSTANTS[17]+CONSTANTS[18]+CONSTANTS[19]) -  CONSTANTS[4]*ALGEBRAIC[29]/(CONSTANTS[16]+CONSTANTS[17]+CONSTANTS[18]+CONSTANTS[19]);
ALGEBRAIC[52] =  CONSTANTS[17]*(ALGEBRAIC[50] - ALGEBRAIC[19])/CONSTANTS[4];
ALGEBRAIC[48] = (ALGEBRAIC[18]>CONSTANTS[22] ?  pow(ALGEBRAIC[18] - CONSTANTS[22], 10.0000)/(pow(ALGEBRAIC[18] - CONSTANTS[22], 10.0000)+pow(0.0300000, 10.0000))*CONSTANTS[40]*CONSTANTS[21]*log(  ALGEBRAIC[18]/ALGEBRAIC[17]*ALGEBRAIC[30]/ALGEBRAIC[28]*pow(ALGEBRAIC[47]/ALGEBRAIC[46], 2.00000))/CONSTANTS[4] : 0.00000);
ALGEBRAIC[54] =  0.0100000*(- ALGEBRAIC[52]+ 2.00000*ALGEBRAIC[29]+ALGEBRAIC[48]);
ALGEBRAIC[37] =  CONSTANTS[3]*pow(STATES[1], 3.00000)*STATES[2]*(STATES[0] - ALGEBRAIC[35])/CONSTANTS[4];
ALGEBRAIC[39] =  CONSTANTS[5]*pow(STATES[3], 2.00000)*STATES[4]*(STATES[0] - ALGEBRAIC[35])/CONSTANTS[4];
ALGEBRAIC[41] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[35])/CONSTANTS[4];
ALGEBRAIC[53] =  CONSTANTS[19]*(ALGEBRAIC[50] - ALGEBRAIC[45])/CONSTANTS[4];
ALGEBRAIC[51] =  CONSTANTS[16]*(ALGEBRAIC[50] - ALGEBRAIC[31])/CONSTANTS[4];
ALGEBRAIC[55] =  0.0100000*(- ALGEBRAIC[51] -  3.00000*ALGEBRAIC[29]+ALGEBRAIC[48]+ALGEBRAIC[53]);
ALGEBRAIC[22] =  CONSTANTS[7]*pow(STATES[5], 2.00000)*(STATES[0] - ALGEBRAIC[21])/CONSTANTS[4];
ALGEBRAIC[23] =  CONSTANTS[8]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[21])/CONSTANTS[4];
ALGEBRAIC[24] =  CONSTANTS[13]*(STATES[0] - ALGEBRAIC[21])/CONSTANTS[4];
}