/*
   There are a total of 2 entries in the algebraic variable array.
   There are a total of 5 entries in each of the rate and state variable arrays.
   There are a total of 26 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * STATES[0] is R in component R (dimensionless).
 * CONSTANTS[19] is phi_1 in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[20] is phi_2 in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[21] is phi_2b in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[0] is k_1b in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[1] is l_2b in component reaction_rate_constants (first_order_rate_constant).
 * STATES[1] is I_1 in component I_1 (dimensionless).
 * STATES[2] is O in component O (dimensionless).
 * CONSTANTS[2] is p in component reaction_rate_constants (micromolar).
 * CONSTANTS[22] is phi_3 in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[23] is phi_4 in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[24] is phi_4b in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[3] is k_3b in component reaction_rate_constants (first_order_rate_constant).
 * STATES[3] is A in component A (dimensionless).
 * ALGEBRAIC[0] is S in component S (dimensionless).
 * STATES[4] is I_2 in component I_2 (dimensionless).
 * CONSTANTS[25] is phi_5 in component reaction_rate_constants (first_order_rate_constant).
 * ALGEBRAIC[1] is open_probability in component open_probability (dimensionless).
 * CONSTANTS[4] is k_1a in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[5] is k_2a in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[6] is k_2b in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[7] is k_3a in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[8] is k_4a in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[9] is k_4b in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[10] is l_2a in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[11] is l_4a in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[12] is l_4b in component reaction_rate_constants (second_order_rate_constant).
 * CONSTANTS[13] is l_6a in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[14] is l_6b in component reaction_rate_constants (first_order_rate_constant).
 * CONSTANTS[15] is L_1 in component reaction_rate_constants (micromolar).
 * CONSTANTS[16] is L_3 in component reaction_rate_constants (micromolar).
 * CONSTANTS[17] is L_5 in component reaction_rate_constants (micromolar).
 * CONSTANTS[18] is c in component reaction_rate_constants (micromolar).
 * RATES[0] is d/dt R in component R (dimensionless).
 * RATES[2] is d/dt O in component O (dimensionless).
 * RATES[1] is d/dt I_1 in component I_1 (dimensionless).
 * RATES[4] is d/dt I_2 in component I_2 (dimensionless).
 * RATES[3] is d/dt A in component A (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 1;
CONSTANTS[0] = 0.04;
CONSTANTS[1] = 0.8;
STATES[1] = 0;
STATES[2] = 0;
CONSTANTS[2] = 10;
CONSTANTS[3] = 29.8;
STATES[3] = 0;
STATES[4] = 0;
CONSTANTS[4] = 0.64;
CONSTANTS[5] = 37.4;
CONSTANTS[6] = 1.4;
CONSTANTS[7] = 0.11;
CONSTANTS[8] = 4;
CONSTANTS[9] = 0.54;
CONSTANTS[10] = 1.7;
CONSTANTS[11] = 1.7;
CONSTANTS[12] = 2.5;
CONSTANTS[13] = 4707;
CONSTANTS[14] = 11.4;
CONSTANTS[15] = 0.12;
CONSTANTS[16] = 0.025;
CONSTANTS[17] = 54.7;
CONSTANTS[18] = 1;
CONSTANTS[19] = ( ( CONSTANTS[4]*CONSTANTS[15]+CONSTANTS[10])*CONSTANTS[18])/(CONSTANTS[15]+ CONSTANTS[18]*(1.00000+CONSTANTS[15]/CONSTANTS[16]));
CONSTANTS[20] = ( CONSTANTS[5]*CONSTANTS[16]+ CONSTANTS[11]*CONSTANTS[18])/(CONSTANTS[16]+ CONSTANTS[18]*(1.00000+CONSTANTS[16]/CONSTANTS[15]));
CONSTANTS[21] = (CONSTANTS[6]+ CONSTANTS[12]*CONSTANTS[18])/(1.00000+CONSTANTS[18]/CONSTANTS[17]);
CONSTANTS[22] = ( CONSTANTS[7]*CONSTANTS[17])/(CONSTANTS[18]+CONSTANTS[17]);
CONSTANTS[23] = ( ( CONSTANTS[8]*CONSTANTS[17]+CONSTANTS[13])*CONSTANTS[18])/(CONSTANTS[18]+CONSTANTS[17]);
CONSTANTS[24] = ( CONSTANTS[15]*(CONSTANTS[9]+CONSTANTS[14]))/(CONSTANTS[18]+CONSTANTS[15]);
CONSTANTS[25] = ( ( CONSTANTS[4]*CONSTANTS[15]+CONSTANTS[10])*CONSTANTS[18])/(CONSTANTS[18]+CONSTANTS[15]);
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[0] = ( CONSTANTS[21]*STATES[2]+ (CONSTANTS[0]+CONSTANTS[1])*STATES[1]) - ( CONSTANTS[20]*CONSTANTS[2]*STATES[0]+ CONSTANTS[19]*STATES[0]);
RATES[1] =  CONSTANTS[19]*STATES[0] -  (CONSTANTS[0]+CONSTANTS[1])*STATES[1];
RATES[4] =  CONSTANTS[25]*STATES[3] -  (CONSTANTS[0]+CONSTANTS[1])*STATES[4];
RATES[3] = ( CONSTANTS[23]*STATES[2]+ (CONSTANTS[0]+CONSTANTS[1])*STATES[4]) - ( CONSTANTS[24]*STATES[3]+ CONSTANTS[25]*STATES[3]);
ALGEBRAIC[0] = 1.00000 - (STATES[0]+STATES[2]+STATES[3]+STATES[1]+STATES[4]);
RATES[2] = ( CONSTANTS[20]*CONSTANTS[2]*STATES[0]+ CONSTANTS[24]*STATES[3]+ CONSTANTS[3]*ALGEBRAIC[0]) -  (CONSTANTS[21]+CONSTANTS[23]+ 1.00000*CONSTANTS[22])*STATES[2];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = 1.00000 - (STATES[0]+STATES[2]+STATES[3]+STATES[1]+STATES[4]);
ALGEBRAIC[1] = pow( 0.100000*STATES[2]+ 0.900000*STATES[3], 4.00000);
}