/*
   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[35] is i_NaT in component i_NaT (nanoA_per_cm2).
 * ALGEBRAIC[37] is i_NaP in component i_NaP (nanoA_per_cm2).
 * ALGEBRAIC[44] is i_leakNa in component i_leakNa (nanoA_per_cm2).
 * ALGEBRAIC[39] is i_KDR in component i_KDR (nanoA_per_cm2).
 * ALGEBRAIC[41] is i_KA in component i_KA (nanoA_per_cm2).
 * ALGEBRAIC[46] is i_leakK in component i_leakK (nanoA_per_cm2).
 * ALGEBRAIC[15] is i_leakf in component i_leakf (nanoA_per_cm2).
 * ALGEBRAIC[43] is i_NaKATPase_n in component i_NaKATPase_n (nanoA_per_cm2).
 * ALGEBRAIC[36] 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[28] 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[2] is beta_m in component i_NaT_m_gate (per_millisecond).
 * ALGEBRAIC[3] is alpha_h in component i_NaT_h_gate (per_millisecond).
 * ALGEBRAIC[4] is beta_h in component i_NaT_h_gate (per_millisecond).
 * ALGEBRAIC[38] 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[5] is alpha_m in component i_NaP_m_gate (per_millisecond).
 * ALGEBRAIC[6] is beta_m in component i_NaP_m_gate (per_millisecond).
 * CONSTANTS[6] is tau_activation in component i_NaP_m_gate (millisecond).
 * ALGEBRAIC[7] is alpha_h in component i_NaP_h_gate (per_millisecond).
 * ALGEBRAIC[8] is beta_h in component i_NaP_h_gate (per_millisecond).
 * ALGEBRAIC[40] is J_KDR in component i_KDR (picomole_per_cm2_millisecond).
 * CONSTANTS[7] is gKDR in component i_KDR (microS_per_cm2).
 * ALGEBRAIC[29] is EK_n in component electric_potentials (millivolt).
 * STATES[5] is n in component i_KDR_n_gate (dimensionless).
 * ALGEBRAIC[9] is alpha_n in component i_KDR_n_gate (per_millisecond).
 * ALGEBRAIC[10] is beta_n in component i_KDR_n_gate (per_millisecond).
 * ALGEBRAIC[42] 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[11] is alpha_m in component i_KA_m_gate (per_millisecond).
 * ALGEBRAIC[12] is beta_m in component i_KA_m_gate (per_millisecond).
 * ALGEBRAIC[13] is alpha_h in component i_KA_h_gate (per_millisecond).
 * ALGEBRAIC[14] is beta_h in component i_KA_h_gate (per_millisecond).
 * ALGEBRAIC[26] 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[22] is Nan in component ion_concentrations (millimolar).
 * ALGEBRAIC[18] is Ko in component ion_concentrations (millimolar).
 * ALGEBRAIC[45] is J_leakNa in component i_leakNa (picomole_per_cm2_millisecond).
 * CONSTANTS[12] is gleakNa in component i_leakNa (microS_per_cm2).
 * ALGEBRAIC[47] 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[30] is ENa_g in component electric_potentials (millivolt).
 * ALGEBRAIC[31] is EK_g in component electric_potentials (millivolt).
 * ALGEBRAIC[49] is ECl_g in component electric_potentials (millivolt).
 * ALGEBRAIC[32] 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[27] 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[20] 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[21] 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[34] is Clo in component ion_concentrations (millimolar).
 * ALGEBRAIC[33] is Clg in component ion_concentrations (millimolar).
 * STATES[8] is N_Nag in component N_Nag (nanomole_per_cm2).
 * ALGEBRAIC[54] 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[55] 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[24] is HCO3o in component ion_concentrations (millimolar).
 * ALGEBRAIC[25] 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[19] is Kn in component ion_concentrations (millimolar).
 * ALGEBRAIC[23] 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).
 * There are a total of 3 condition variables.
 */
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);
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[13] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - ( 0.00100000*(- (ALGEBRAIC[35]+ALGEBRAIC[37]+ALGEBRAIC[44]+ALGEBRAIC[39]+ALGEBRAIC[41]+ALGEBRAIC[46]+ALGEBRAIC[15]+ALGEBRAIC[43])+ALGEBRAIC[0]))/CONSTANTS[0];
resid[1] = RATES[1] -  ALGEBRAIC[1]*(1.00000 - STATES[1]) -  ALGEBRAIC[2]*STATES[1];
resid[2] = RATES[2] -  ALGEBRAIC[3]*(1.00000 - STATES[2]) -  ALGEBRAIC[4]*STATES[2];
resid[3] = RATES[3] -  ALGEBRAIC[5]*(1.00000 - STATES[3]) -  ALGEBRAIC[6]*STATES[3];
resid[4] = RATES[4] -  ALGEBRAIC[7]*(1.00000 - STATES[4]) -  ALGEBRAIC[8]*STATES[4];
resid[5] = RATES[5] -  ALGEBRAIC[9]*(1.00000 - STATES[5]) -  ALGEBRAIC[10]*STATES[5];
resid[6] = RATES[6] -  ALGEBRAIC[11]*(1.00000 - STATES[6]) -  ALGEBRAIC[12]*STATES[6];
resid[7] = RATES[7] -  ALGEBRAIC[13]*(1.00000 - STATES[7]) -  ALGEBRAIC[14]*STATES[7];
resid[8] = RATES[8] - ALGEBRAIC[54];
resid[9] = RATES[9] - ALGEBRAIC[55];
resid[10] = RATES[10] -  10.0000*CONSTANTS[23]*((ALGEBRAIC[20]+ALGEBRAIC[17]+ALGEBRAIC[33]+ALGEBRAIC[25]+CONSTANTS[24]/STATES[10]) - (ALGEBRAIC[21]+ALGEBRAIC[18]+ALGEBRAIC[34]+ALGEBRAIC[24]));
resid[11] = RATES[11] -  0.0100000*(( 3.00000*ALGEBRAIC[26]+ALGEBRAIC[36]+ALGEBRAIC[38]+ALGEBRAIC[45]) -  100.000*ALGEBRAIC[54]);
resid[12] = RATES[12] -  -0.0100000*(ALGEBRAIC[36]+ALGEBRAIC[38]+ALGEBRAIC[45]+ 3.00000*ALGEBRAIC[26]) - ALGEBRAIC[55];
resid[13] = RATES[13] -  - 0.0100000*2.00000*ALGEBRAIC[53];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[40] = ( CONSTANTS[7]*pow(STATES[5], 2.00000)*(STATES[0] - ALGEBRAIC[29]))/CONSTANTS[4];
ALGEBRAIC[42] = ( CONSTANTS[8]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[29]))/CONSTANTS[4];
ALGEBRAIC[47] = ( CONSTANTS[13]*(STATES[0] - ALGEBRAIC[29]))/CONSTANTS[4];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<=0.00000 ? 3500.00 : 0.00000);
ALGEBRAIC[1] = ( 0.320000*(- STATES[0] - 51.9000))/(exp(- ( 0.250000*STATES[0]+12.9750)) - 1.00000);
ALGEBRAIC[2] = ( 0.280000*(STATES[0]+24.8900))/(exp( 0.200000*STATES[0]+4.97800) - 1.00000);
ALGEBRAIC[3] =  0.128000*exp(- ( 0.0560000*STATES[0]+2.94000));
ALGEBRAIC[4] = 4.00000/(exp(- ( 0.200000*STATES[0]+6.00000))+1.00000);
ALGEBRAIC[5] = ( (1.00000/CONSTANTS[6])*1.00000)/(exp(- ( 0.143000*STATES[0]+5.67000))+1.00000);
ALGEBRAIC[6] = ( (1.00000/CONSTANTS[6])*exp(- ( 0.143000*STATES[0]+5.67000)))/(exp(- ( 0.143000*STATES[0]+5.67000))+1.00000);
ALGEBRAIC[7] =  5.12000e-08*exp(- ( 0.0560000*STATES[0]+2.94000));
ALGEBRAIC[8] = 1.60000e-06/(exp(- ( 0.200000*STATES[0]+8.00000))+1.00000);
ALGEBRAIC[9] = ( 0.0160000*(- STATES[0] - 34.9000))/(exp(- ( 0.200000*STATES[0]+6.98000)) - 1.00000);
ALGEBRAIC[10] =  0.250000*exp(- ( 0.0250000*STATES[0]+1.25000));
ALGEBRAIC[11] = ( 0.0200000*(- STATES[0] - 56.9000))/(exp(- ( 0.100000*STATES[0]+5.69000)) - 1.00000);
ALGEBRAIC[12] = ( 0.0175000*(STATES[0]+29.9000))/(exp( 0.100000*STATES[0]+2.99000) - 1.00000);
ALGEBRAIC[13] =  0.0160000*exp(- ( 0.0560000*STATES[0]+4.61000));
ALGEBRAIC[14] = 0.500000/(exp(- ( 0.200000*STATES[0]+11.9800))+1.00000);
ALGEBRAIC[15] =  CONSTANTS[14]*(STATES[0] - CONSTANTS[15]);
ALGEBRAIC[17] = STATES[9]/STATES[10];
ALGEBRAIC[16] = (CONSTANTS[26]+CONSTANTS[25]) - STATES[10];
ALGEBRAIC[18] = STATES[12]/ALGEBRAIC[16];
ALGEBRAIC[20] = STATES[8]/STATES[10];
ALGEBRAIC[21] = STATES[11]/ALGEBRAIC[16];
ALGEBRAIC[24] = STATES[13]/ALGEBRAIC[16];
ALGEBRAIC[25] = (( CONSTANTS[26]*CONSTANTS[42]+ CONSTANTS[25]*CONSTANTS[34]) - STATES[13])/STATES[10];
ALGEBRAIC[22] = (( CONSTANTS[35]*CONSTANTS[29]+ CONSTANTS[26]*CONSTANTS[31]+ CONSTANTS[25]*CONSTANTS[30]) - (STATES[11]+STATES[8]))/CONSTANTS[35];
ALGEBRAIC[26] = ( (( CONSTANTS[9]*pow(ALGEBRAIC[22], 1.50000))/(pow(ALGEBRAIC[22], 1.50000)+pow(CONSTANTS[10], 1.50000)))*ALGEBRAIC[18])/(ALGEBRAIC[18]+CONSTANTS[11]);
ALGEBRAIC[33] = (ALGEBRAIC[20]+ALGEBRAIC[17]) - (ALGEBRAIC[25]+( CONSTANTS[36]*CONSTANTS[24])/STATES[10]);
ALGEBRAIC[19] = (( CONSTANTS[35]*CONSTANTS[32]+ CONSTANTS[26]*CONSTANTS[33]+ CONSTANTS[25]*CONSTANTS[22]) - (STATES[12]+STATES[9]))/CONSTANTS[35];
ALGEBRAIC[23] = ((CONSTANTS[28]+ALGEBRAIC[22]) - CONSTANTS[29])+(ALGEBRAIC[19] - CONSTANTS[32]);
ALGEBRAIC[34] = (( CONSTANTS[26]*CONSTANTS[41]+ CONSTANTS[25]*CONSTANTS[39]+ CONSTANTS[35]*CONSTANTS[28]) - ( STATES[10]*ALGEBRAIC[33]+ CONSTANTS[35]*ALGEBRAIC[23]))/ALGEBRAIC[16];
ALGEBRAIC[28] =  CONSTANTS[40]*log(ALGEBRAIC[21]/ALGEBRAIC[22]);
ALGEBRAIC[35] =  CONSTANTS[3]*pow(STATES[1], 3.00000)*STATES[2]*(STATES[0] - ALGEBRAIC[28]);
ALGEBRAIC[36] = ( CONSTANTS[3]*pow(STATES[1], 3.00000)*STATES[2]*(STATES[0] - ALGEBRAIC[28]))/CONSTANTS[4];
ALGEBRAIC[37] =  CONSTANTS[5]*pow(STATES[3], 2.00000)*STATES[4]*(STATES[0] - ALGEBRAIC[28]);
ALGEBRAIC[38] = ( CONSTANTS[5]*pow(STATES[3], 2.00000)*STATES[4]*(STATES[0] - ALGEBRAIC[28]))/CONSTANTS[4];
ALGEBRAIC[29] =  CONSTANTS[40]*log(ALGEBRAIC[18]/ALGEBRAIC[19]);
ALGEBRAIC[39] =  CONSTANTS[7]*pow(STATES[5], 2.00000)*(STATES[0] - ALGEBRAIC[29]);
ALGEBRAIC[41] =  CONSTANTS[8]*pow(STATES[6], 2.00000)*STATES[7]*(STATES[0] - ALGEBRAIC[29]);
ALGEBRAIC[43] =  ALGEBRAIC[26]*CONSTANTS[4];
ALGEBRAIC[44] =  CONSTANTS[12]*(STATES[0] - ALGEBRAIC[28]);
ALGEBRAIC[45] = ( CONSTANTS[12]*(STATES[0] - ALGEBRAIC[28]))/CONSTANTS[4];
ALGEBRAIC[46] =  CONSTANTS[13]*(STATES[0] - ALGEBRAIC[29]);
ALGEBRAIC[30] =  CONSTANTS[40]*log(ALGEBRAIC[21]/ALGEBRAIC[20]);
ALGEBRAIC[31] =  CONSTANTS[40]*log(ALGEBRAIC[18]/ALGEBRAIC[17]);
ALGEBRAIC[49] =  - 1.00000*CONSTANTS[40]*log(ALGEBRAIC[34]/ALGEBRAIC[33]);
ALGEBRAIC[32] =  - CONSTANTS[40]*log( (ALGEBRAIC[21]/ALGEBRAIC[20])*pow(ALGEBRAIC[24]/ALGEBRAIC[25], 2.00000));
ALGEBRAIC[27] = ( (( CONSTANTS[20]*pow(ALGEBRAIC[20], 1.50000))/(pow(ALGEBRAIC[20], 1.50000)+pow(CONSTANTS[10], 1.50000)))*ALGEBRAIC[18])/(ALGEBRAIC[18]+CONSTANTS[11]);
ALGEBRAIC[50] = ( CONSTANTS[16]*ALGEBRAIC[30]+ CONSTANTS[17]*ALGEBRAIC[31]+ CONSTANTS[18]*ALGEBRAIC[49]+ CONSTANTS[19]*ALGEBRAIC[32])/(CONSTANTS[16]+CONSTANTS[17]+CONSTANTS[18]+CONSTANTS[19]) - ( CONSTANTS[4]*ALGEBRAIC[27])/(CONSTANTS[16]+CONSTANTS[17]+CONSTANTS[18]+CONSTANTS[19]);
ALGEBRAIC[53] = ( CONSTANTS[19]*(ALGEBRAIC[50] - ALGEBRAIC[32]))/CONSTANTS[4];
ALGEBRAIC[51] = ( CONSTANTS[16]*(ALGEBRAIC[50] - ALGEBRAIC[30]))/CONSTANTS[4];
ALGEBRAIC[48] = (CONDVAR[2]>0.00000 ? ( (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[21])/ALGEBRAIC[20])*pow(ALGEBRAIC[34]/ALGEBRAIC[33], 2.00000)))/CONSTANTS[4] : 0.00000);
ALGEBRAIC[54] =  0.0100000*((- ALGEBRAIC[51] -  3.00000*ALGEBRAIC[27])+ALGEBRAIC[48]+ALGEBRAIC[53]);
ALGEBRAIC[52] = ( CONSTANTS[17]*(ALGEBRAIC[50] - ALGEBRAIC[31]))/CONSTANTS[4];
ALGEBRAIC[55] =  0.0100000*(- ALGEBRAIC[52]+ 2.00000*ALGEBRAIC[27]+ALGEBRAIC[48]);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[1];
CONDVAR[1] = VOI - CONSTANTS[2];
CONDVAR[2] = ALGEBRAIC[18] - CONSTANTS[22];
}