/*
   There are a total of 34 entries in the algebraic variable array.
   There are a total of 9 entries in each of the rate and state variable arrays.
   There are a total of 54 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (millisecond).
 * CONSTANTS[0] is isotonic in component parameters (dimensionless).
 * CONSTANTS[1] is alpha_1 in component parameters (per_um).
 * CONSTANTS[2] is beta_1 in component parameters (mN).
 * CONSTANTS[3] is alpha_2 in component parameters (per_um).
 * CONSTANTS[4] is beta_2 in component parameters (mN).
 * CONSTANTS[5] is alpha_3 in component parameters (per_um).
 * CONSTANTS[6] is beta_3 in component parameters (mN).
 * CONSTANTS[7] is lambda in component parameters (mN).
 * CONSTANTS[8] is A_half in component parameters (dimensionless).
 * CONSTANTS[9] is mu in component parameters (dimensionless).
 * CONSTANTS[10] is chi in component parameters (dimensionless).
 * CONSTANTS[11] is chi_0 in component parameters (dimensionless).
 * CONSTANTS[12] is m_0 in component parameters (dimensionless).
 * CONSTANTS[13] is v_max in component parameters (um_per_msec).
 * CONSTANTS[14] is a in component parameters (dimensionless).
 * CONSTANTS[15] is d_h in component parameters (dimensionless).
 * CONSTANTS[16] is alpha_P in component parameters (dimensionless).
 * ALGEBRAIC[7] is l in component length (um).
 * ALGEBRAIC[25] is F_muscle in component force (mN).
 * ALGEBRAIC[23] is flag in component isotonic (dimensionless).
 * CONSTANTS[17] is F_afterload in component isotonic (mN).
 * ALGEBRAIC[29] is isotonic_mode in component isotonic (dimensionless).
 * CONSTANTS[18] is l_0 in component isotonic (um).
 * CONSTANTS[19] is S_0 in component parameters_izakov_et_al_1991 (um).
 * ALGEBRAIC[0] is q_v in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[20] is q_1 in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[21] is q_2 in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[22] is q_3 in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[23] is q_4 in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[24] is v_star in component parameters_izakov_et_al_1991 (um_per_msec).
 * CONSTANTS[50] is v_1 in component parameters_izakov_et_al_1991 (um_per_msec).
 * CONSTANTS[25] is alpha_G in component parameters_izakov_et_al_1991 (dimensionless).
 * CONSTANTS[26] is a_on in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[27] is a_off in component parameters_izakov_et_al_1991 (per_millisecond).
 * CONSTANTS[28] is k_A in component parameters_izakov_et_al_1991 (dimensionless).
 * STATES[0] is v in component CE_velocity (um_per_msec).
 * CONSTANTS[29] is alpha_Q in component parameters_izakov_et_al_1991 (dimensionless).
 * CONSTANTS[30] is beta_Q in component parameters_izakov_et_al_1991 (dimensionless).
 * ALGEBRAIC[24] is F_CE in component force (mN).
 * ALGEBRAIC[3] is F_XSE in component force (mN).
 * ALGEBRAIC[1] is F_SE in component force (mN).
 * ALGEBRAIC[2] is F_PE in component force (mN).
 * STATES[1] is N in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[10] is k_P_vis in component CE_velocity (mN_msec_per_um).
 * ALGEBRAIC[12] is k_S_vis in component PE_velocity (mN_msec_per_um).
 * STATES[2] is w in component PE_velocity (um_per_msec).
 * STATES[3] is l_1 in component length (um).
 * STATES[4] is l_2 in component length (um).
 * STATES[5] is l_3 in component length (um).
 * ALGEBRAIC[17] is p_v in component average_crossbridge_force (dimensionless).
 * ALGEBRAIC[28] is K_chi in component crossbridge_kinetics (per_millisecond).
 * ALGEBRAIC[4] is M_A in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[5] is n_1 in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[6] is L_oz in component crossbridge_kinetics (dimensionless).
 * ALGEBRAIC[26] is k_p_v in component crossbridge_kinetics (per_millisecond).
 * ALGEBRAIC[27] is k_m_v in component crossbridge_kinetics (per_millisecond).
 * STATES[6] is A in component calcium_handling (dimensionless).
 * ALGEBRAIC[14] is G_star in component average_crossbridge_force (dimensionless).
 * ALGEBRAIC[13] is P_star in component average_crossbridge_force (dimensionless).
 * CONSTANTS[31] is v_0 in component crossbridge_kinetics (um_per_msec).
 * CONSTANTS[32] is q_star in component crossbridge_kinetics (per_millisecond).
 * ALGEBRAIC[8] is dl_1_dt in component length (um_per_msec).
 * ALGEBRAIC[32] is dl_2_dt in component length (um_per_msec).
 * ALGEBRAIC[33] is dl_3_dt in component length (um_per_msec).
 * ALGEBRAIC[31] is phi_chi_2 in component CE_velocity (um_per_msec).
 * ALGEBRAIC[30] is phi_chi in component CE_velocity (um_per_msec2).
 * ALGEBRAIC[18] is p_prime_v in component average_crossbridge_force (msec_per_um).
 * CONSTANTS[33] is alpha_P_lengthening in component CE_velocity (per_um).
 * CONSTANTS[34] is beta_P_lengthening in component CE_velocity (mN_msec_per_um).
 * CONSTANTS[35] is alpha_P_shortening in component CE_velocity (per_um).
 * CONSTANTS[36] is beta_P_shortening in component CE_velocity (mN_msec_per_um).
 * ALGEBRAIC[9] is alp_p in component CE_velocity (per_um).
 * CONSTANTS[37] is alpha_S_lengthening in component PE_velocity (per_um).
 * CONSTANTS[38] is beta_S_lengthening in component PE_velocity (mN_msec_per_um).
 * CONSTANTS[39] is alpha_S_shortening in component PE_velocity (per_um).
 * CONSTANTS[40] is beta_S_shortening in component PE_velocity (mN_msec_per_um).
 * ALGEBRAIC[11] is alp_s in component PE_velocity (per_um).
 * CONSTANTS[53] is gamma in component average_crossbridge_force (dimensionless).
 * CONSTANTS[51] is case_1 in component average_crossbridge_force (msec_per_um).
 * ALGEBRAIC[15] is case_2 in component average_crossbridge_force (msec_per_um).
 * CONSTANTS[52] is case_3 in component average_crossbridge_force (msec_per_um).
 * ALGEBRAIC[16] is case_4 in component average_crossbridge_force (msec_per_um).
 * ALGEBRAIC[21] is dA_dt in component calcium_handling (per_millisecond).
 * ALGEBRAIC[19] is N_A in component calcium_handling (dimensionless).
 * ALGEBRAIC[20] is pi_N_A in component calcium_handling (dimensionless).
 * STATES[7] is B in component calcium_handling (dimensionless).
 * ALGEBRAIC[22] is dB_dt in component calcium_handling (per_millisecond).
 * STATES[8] is Ca_C in component calcium_handling (dimensionless).
 * CONSTANTS[41] is A_tot in component calcium_handling (dimensionless).
 * CONSTANTS[42] is B_tot in component calcium_handling (dimensionless).
 * CONSTANTS[43] is b_on in component calcium_handling (per_millisecond).
 * CONSTANTS[44] is b_off in component calcium_handling (per_millisecond).
 * CONSTANTS[45] is a_c in component calcium_handling (per_millisecond2).
 * CONSTANTS[46] is r_Ca in component calcium_handling (per_millisecond).
 * CONSTANTS[47] is q_Ca in component calcium_handling (dimensionless).
 * CONSTANTS[48] is t_d in component calcium_handling (millisecond).
 * CONSTANTS[49] is Ca_m in component calcium_handling (dimensionless).
 * RATES[1] is d/dt N in component crossbridge_kinetics (dimensionless).
 * RATES[3] is d/dt l_1 in component length (um).
 * RATES[4] is d/dt l_2 in component length (um).
 * RATES[5] is d/dt l_3 in component length (um).
 * RATES[0] is d/dt v in component CE_velocity (um_per_msec).
 * RATES[2] is d/dt w in component PE_velocity (um_per_msec).
 * RATES[6] is d/dt A in component calcium_handling (dimensionless).
 * RATES[7] is d/dt B in component calcium_handling (dimensionless).
 * RATES[8] is d/dt Ca_C in component calcium_handling (dimensionless).
 * There are a total of 24 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 0;
CONSTANTS[1] = 19;
CONSTANTS[2] = 0.29;
CONSTANTS[3] = 14.6;
CONSTANTS[4] = 0.000924;
CONSTANTS[5] = 48;
CONSTANTS[6] = 0.01;
CONSTANTS[7] = 96;
CONSTANTS[8] = 0.6;
CONSTANTS[9] = 3;
CONSTANTS[10] = 0.705;
CONSTANTS[11] = 3;
CONSTANTS[12] = 0.9;
CONSTANTS[13] = 0.0056;
CONSTANTS[14] = 0.25;
CONSTANTS[15] = 0.5;
CONSTANTS[16] = 4;
CONSTANTS[17] = 2;
CONSTANTS[18] = 0.527;
CONSTANTS[19] = 1.14;
CONSTANTS[20] = 0.0173;
CONSTANTS[21] = 0.259;
CONSTANTS[22] = 0.0173;
CONSTANTS[23] = 0.015;
CONSTANTS[24] = 0.0053035675;
CONSTANTS[25] = 1;
CONSTANTS[26] = 2.3;
CONSTANTS[27] = 0.29;
CONSTANTS[28] = 2.8;
STATES[0] = 0;
CONSTANTS[29] = 10;
CONSTANTS[30] = 5000;
STATES[1] = 0.0001;
STATES[2] = 0;
STATES[3] = 0.437;
STATES[4] = 0.439;
STATES[5] = 0.089;
STATES[6] = 0.01;
CONSTANTS[31] = 0.0056;
CONSTANTS[32] = 1;
CONSTANTS[33] = 16;
CONSTANTS[34] = 1.5;
CONSTANTS[35] = 16;
CONSTANTS[36] = 1.5;
CONSTANTS[37] = 39;
CONSTANTS[38] = 8;
CONSTANTS[39] = 46;
CONSTANTS[40] = 6;
STATES[7] = 0;
STATES[8] = 0;
CONSTANTS[41] = 1;
CONSTANTS[42] = 0.4;
CONSTANTS[43] = 2.6;
CONSTANTS[44] = 0.182;
CONSTANTS[45] = 0.0052;
CONSTANTS[46] = 0.65;
CONSTANTS[47] = 50;
CONSTANTS[48] = 33;
CONSTANTS[49] = 0.03;
CONSTANTS[50] = CONSTANTS[13]/10.0000;
CONSTANTS[51] = ( CONSTANTS[14]*(0.400000+ 0.400000*CONSTANTS[14]))/( CONSTANTS[13]*pow( (CONSTANTS[14]+1.00000)*0.400000, 2.00000));
CONSTANTS[52] = ( 0.400000*CONSTANTS[14]+1.00000)/( CONSTANTS[14]*CONSTANTS[13]);
CONSTANTS[53] = ( CONSTANTS[14]*CONSTANTS[15]*pow(CONSTANTS[50]/CONSTANTS[13], 2.00000))/( 3.00000*CONSTANTS[14]*CONSTANTS[15] - ( (CONSTANTS[14]+1.00000)*CONSTANTS[50])/CONSTANTS[13]);
RATES[1] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[0] = 0.1001;
RATES[2] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[1] - ALGEBRAIC[28];
resid[1] = RATES[3] - ALGEBRAIC[8];
resid[2] = RATES[4] - ALGEBRAIC[32];
resid[3] = RATES[5] - ALGEBRAIC[33];
resid[4] = RATES[0] - (ALGEBRAIC[12]==0.00000 ? ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*(ALGEBRAIC[31] - STATES[0]) - ( CONSTANTS[7]*ALGEBRAIC[28]*ALGEBRAIC[17]+ ALGEBRAIC[9]*ALGEBRAIC[10]*pow(STATES[0], 2.00000)))/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[18]+ALGEBRAIC[10]) : ALGEBRAIC[30]);
resid[5] = RATES[2] - (ALGEBRAIC[29]==1.00000 ? (( ALGEBRAIC[12]*(ALGEBRAIC[30] -  ALGEBRAIC[11]*pow(STATES[2] - STATES[0], 2.00000)) -  CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*(STATES[2] - STATES[0])) -  CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])*STATES[2])/ALGEBRAIC[12] : (ALGEBRAIC[30] -  ALGEBRAIC[11]*pow(STATES[2] - STATES[0], 2.00000)) - ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*(STATES[2] - STATES[0])+ ( CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5]))*STATES[2])/ALGEBRAIC[12]);
resid[6] = RATES[6] - ALGEBRAIC[21];
resid[7] = RATES[7] - ALGEBRAIC[22];
resid[8] = RATES[8] - (CONDVAR[23]<0.00000 ?  4.00000*CONSTANTS[45]*CONSTANTS[49]*VOI*(1.00000 - exp( - CONSTANTS[45]*pow(VOI, 2.00000)))*exp( - CONSTANTS[45]*pow(VOI, 2.00000)) : (- ALGEBRAIC[21] - ALGEBRAIC[22]) -  CONSTANTS[46]*exp( - CONSTANTS[47]*STATES[8])*STATES[8]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[1] =  CONSTANTS[2]*(exp( CONSTANTS[1]*(STATES[4] - STATES[3])) - 1.00000);
ALGEBRAIC[2] =  CONSTANTS[4]*(exp( CONSTANTS[3]*STATES[4]) - 1.00000);
ALGEBRAIC[24] =  CONSTANTS[7]*ALGEBRAIC[17]*STATES[1];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[8] = STATES[0];
ALGEBRAIC[9] = (CONDVAR[8]<=0.00000 ? CONSTANTS[33] : CONSTANTS[35]);
ALGEBRAIC[10] = (CONDVAR[9]<=0.00000 ?  CONSTANTS[34]*exp( CONSTANTS[33]*STATES[3]) :  CONSTANTS[36]*exp( CONSTANTS[35]*STATES[3]));
ALGEBRAIC[11] = (CONDVAR[10]<=0.00000 ? CONSTANTS[37] : CONSTANTS[39]);
ALGEBRAIC[12] = (CONDVAR[11]<=0.00000 ?  CONSTANTS[38]*exp( CONSTANTS[37]*(STATES[4] - STATES[3])) :  CONSTANTS[40]*exp( CONSTANTS[39]*(STATES[4] - STATES[3])));
ALGEBRAIC[13] = (CONDVAR[12]<=0.00000 ? ( CONSTANTS[14]*(1.00000+STATES[0]/CONSTANTS[13]))/(CONSTANTS[14] - STATES[0]/CONSTANTS[13]) : (1.00000+CONSTANTS[15]) - ( pow(CONSTANTS[15], 2.00000)*CONSTANTS[14])/( (( CONSTANTS[14]*CONSTANTS[15])/CONSTANTS[53])*pow(STATES[0]/CONSTANTS[13], 2.00000)+( (CONSTANTS[14]+1.00000)*STATES[0])/CONSTANTS[13]+ CONSTANTS[14]*CONSTANTS[15]));
ALGEBRAIC[14] = (CONDVAR[13]<=0.00000&&CONDVAR[14]<=0.00000 ? 1.00000+( 0.600000*STATES[0])/CONSTANTS[13] : CONDVAR[15]<0.00000&&CONDVAR[16]<=0.00000 ? ALGEBRAIC[13]/(( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13])+1.00000) : ( ALGEBRAIC[13]*exp( - CONSTANTS[25]*pow((STATES[0] - CONSTANTS[50])/CONSTANTS[13], CONSTANTS[16])))/(( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13])+1.00000));
ALGEBRAIC[17] = ALGEBRAIC[13]/ALGEBRAIC[14];
ALGEBRAIC[15] = ( CONSTANTS[14]*1.00000*(1.00000+ 0.400000*CONSTANTS[14]+( 1.20000*STATES[0])/CONSTANTS[13]+ 0.600000*pow(STATES[0]/CONSTANTS[13], 2.00000)))/( CONSTANTS[13]*pow( (CONSTANTS[14] - STATES[0]/CONSTANTS[13])*(1.00000+( 0.600000*STATES[0])/CONSTANTS[13]), 2.00000));
ALGEBRAIC[16] =  (1.00000/CONSTANTS[13])*exp( - CONSTANTS[25]*pow(STATES[0]/CONSTANTS[13] - CONSTANTS[50]/CONSTANTS[13], CONSTANTS[16]))*(( 0.400000*CONSTANTS[14]+1.00000)/CONSTANTS[14]+ CONSTANTS[25]*CONSTANTS[16]*(1.00000+( ( 0.400000*CONSTANTS[14]+1.00000)*STATES[0])/( CONSTANTS[14]*CONSTANTS[13]))*pow(STATES[0]/CONSTANTS[13] - CONSTANTS[50]/CONSTANTS[13], CONSTANTS[16] - 1.00000));
ALGEBRAIC[18] = (CONDVAR[17]<=0.00000 ? CONSTANTS[51] : CONDVAR[18]<0.00000&&CONDVAR[19]<=0.00000 ? ALGEBRAIC[15] : CONDVAR[20]<0.00000&&CONDVAR[21]<=0.00000 ? CONSTANTS[52] : ALGEBRAIC[16]);
ALGEBRAIC[6] = (STATES[3]+CONSTANTS[19])/(0.460000+CONSTANTS[19]);
ALGEBRAIC[19] = STATES[1]/( ALGEBRAIC[6]*STATES[6]);
ALGEBRAIC[20] = (CONDVAR[22]>=0.00000 ? 1.00000 : pow(0.0200000, ALGEBRAIC[19]));
ALGEBRAIC[21] =  CONSTANTS[26]*(CONSTANTS[41] - STATES[6])*STATES[8] -  CONSTANTS[27]*exp( - CONSTANTS[28]*STATES[6])*ALGEBRAIC[20]*STATES[6];
ALGEBRAIC[22] =  CONSTANTS[43]*(CONSTANTS[42] - STATES[7])*STATES[8] -  CONSTANTS[44]*STATES[7];
ALGEBRAIC[4] = pow(STATES[6], CONSTANTS[9])/(pow(STATES[6], CONSTANTS[9])+pow(CONSTANTS[8], CONSTANTS[9]));
ALGEBRAIC[5] =  0.600000*STATES[3]+0.500000;
ALGEBRAIC[0] = (CONDVAR[4]<=0.00000 ? CONSTANTS[20] - ( CONSTANTS[21]*STATES[0])/CONSTANTS[13] : CONDVAR[5]<=0.00000&&CONDVAR[6]<0.00000 ? ( (CONSTANTS[23] - CONSTANTS[22])*STATES[0])/CONSTANTS[24]+CONSTANTS[22] : CONSTANTS[23]/pow(1.00000+( CONSTANTS[30]*(STATES[0] - CONSTANTS[24]))/CONSTANTS[13], CONSTANTS[29]));
ALGEBRAIC[26] =  CONSTANTS[10]*CONSTANTS[11]*ALGEBRAIC[0]*CONSTANTS[12]*ALGEBRAIC[14];
ALGEBRAIC[27] = (CONDVAR[7]<=0.00000 ?  CONSTANTS[11]*ALGEBRAIC[0]*(1.00000 -  CONSTANTS[10]*CONSTANTS[12]*ALGEBRAIC[14]) :  CONSTANTS[11]*( CONSTANTS[23]*(1.00000 -  CONSTANTS[10]*CONSTANTS[12]*ALGEBRAIC[14])+( CONSTANTS[32]*(STATES[0] - CONSTANTS[24]))/(CONSTANTS[31] - CONSTANTS[24])));
ALGEBRAIC[28] =  ALGEBRAIC[26]*ALGEBRAIC[4]*ALGEBRAIC[5]*ALGEBRAIC[6]*(1.00000 - STATES[1]) -  ALGEBRAIC[27]*STATES[1];
ALGEBRAIC[7] = STATES[4]+STATES[5];
ALGEBRAIC[3] =  CONSTANTS[6]*(exp( CONSTANTS[5]*STATES[5]) - 1.00000);
ALGEBRAIC[25] = ALGEBRAIC[3];
ALGEBRAIC[23] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<0.00000 ? 0.00000 : 1.00000);
ALGEBRAIC[29] = (CONSTANTS[0]==0.00000 ? 0.00000 : CONSTANTS[0]==1.00000&&CONDVAR[2]>=0.00000 ? 1.00000 : CONSTANTS[0]==1.00000&&CONDVAR[3]>=0.00000&&ALGEBRAIC[23]==1.00000 ? 0.00000 : 0.0/0.0);
ALGEBRAIC[30] = (ALGEBRAIC[29]==1.00000 ? - ( CONSTANTS[7]*ALGEBRAIC[28]*ALGEBRAIC[17]+ ALGEBRAIC[9]*ALGEBRAIC[10]*pow(STATES[0], 2.00000)+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])*STATES[2])/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[18]+ALGEBRAIC[10]) : - ( CONSTANTS[7]*ALGEBRAIC[28]*ALGEBRAIC[17]+ ALGEBRAIC[9]*ALGEBRAIC[10]*pow(STATES[0], 2.00000)+ ( CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5]))*STATES[2])/( CONSTANTS[7]*STATES[1]*ALGEBRAIC[18]+ALGEBRAIC[10]));
ALGEBRAIC[31] = (ALGEBRAIC[29]==1.00000 ? ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*STATES[0])/( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])) : ( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))*STATES[0])/( CONSTANTS[1]*CONSTANTS[2]*exp( CONSTANTS[1]*(STATES[4] - STATES[3]))+ CONSTANTS[3]*CONSTANTS[4]*exp( CONSTANTS[3]*STATES[4])+ CONSTANTS[5]*CONSTANTS[6]*exp( CONSTANTS[5]*STATES[5])));
ALGEBRAIC[32] = (ALGEBRAIC[12]==0.00000 ? ALGEBRAIC[31] : STATES[2]);
ALGEBRAIC[33] = (ALGEBRAIC[29]==1.00000 ? 0.00000 : ALGEBRAIC[29]==0.00000&&ALGEBRAIC[12]==0.00000 ? - ALGEBRAIC[31] : - STATES[2]);
}
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;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = ALGEBRAIC[7] - CONSTANTS[18];
CONDVAR[1] = VOI - 150.000;
CONDVAR[2] = ALGEBRAIC[25] - CONSTANTS[17];
CONDVAR[3] = ALGEBRAIC[7] - CONSTANTS[18];
CONDVAR[4] = STATES[0] - 0.00000;
CONDVAR[5] = STATES[0] - CONSTANTS[24];
CONDVAR[6] = 0.00000 - STATES[0];
CONDVAR[7] = STATES[0] - CONSTANTS[24];
CONDVAR[8] = STATES[0] - 0.00000;
CONDVAR[9] = STATES[0] - 0.00000;
CONDVAR[10] = STATES[2] - STATES[0];
CONDVAR[11] = STATES[2] - STATES[0];
CONDVAR[12] = STATES[0] - 0.00000;
CONDVAR[13] = - CONSTANTS[13] - STATES[0];
CONDVAR[14] = STATES[0] - 0.00000;
CONDVAR[15] = 0.00000 - STATES[0];
CONDVAR[16] = STATES[0] - CONSTANTS[50];
CONDVAR[17] = STATES[0] - - CONSTANTS[13];
CONDVAR[18] = - CONSTANTS[13] - STATES[0];
CONDVAR[19] = STATES[0] - 0.00000;
CONDVAR[20] = 0.00000 - STATES[0];
CONDVAR[21] = STATES[0] - CONSTANTS[50];
CONDVAR[22] = ALGEBRAIC[19] - 1.00000;
CONDVAR[23] = VOI - CONSTANTS[48];
}