/*
   There are a total of 0 entries in the algebraic variable array.
   There are a total of 4 entries in each of the rate and state variable arrays.
   There are a total of 16 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * STATES[0] is a in component a (micromolar).
 * CONSTANTS[0] is k1 in component model_parameters (dimensionless).
 * CONSTANTS[1] is k3 in component model_parameters (dimensionless).
 * STATES[1] is b in component b (micromolar).
 * STATES[2] is c in component c (micromolar).
 * CONSTANTS[2] is K4 in component model_parameters (dimensionless).
 * CONSTANTS[3] is k5 in component model_parameters (dimensionless).
 * CONSTANTS[4] is K6 in component model_parameters (dimensionless).
 * CONSTANTS[5] is k7 in component model_parameters (dimensionless).
 * CONSTANTS[6] is k8 in component model_parameters (dimensionless).
 * CONSTANTS[7] is K9 in component model_parameters (dimensionless).
 * CONSTANTS[8] is k10 in component model_parameters (dimensionless).
 * STATES[3] is d in component d (micromolar).
 * CONSTANTS[9] is K11 in component model_parameters (dimensionless).
 * CONSTANTS[10] is k12 in component model_parameters (dimensionless).
 * CONSTANTS[11] is k13 in component model_parameters (dimensionless).
 * CONSTANTS[12] is k14 in component model_parameters (dimensionless).
 * CONSTANTS[13] is K15 in component model_parameters (dimensionless).
 * CONSTANTS[14] is k16 in component model_parameters (dimensionless).
 * CONSTANTS[15] is K17 in component model_parameters (dimensionless).
 * RATES[0] is d/dt a in component a (micromolar).
 * RATES[1] is d/dt b in component b (micromolar).
 * RATES[2] is d/dt c in component c (micromolar).
 * RATES[3] is d/dt d in component d (micromolar).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 0.01;
CONSTANTS[0] = 1.5;
CONSTANTS[1] = 0.00001;
STATES[1] = 0.01;
STATES[2] = 0.01;
CONSTANTS[2] = 0.788;
CONSTANTS[3] = 1.45;
CONSTANTS[4] = 0.18;
CONSTANTS[5] = 5.82;
CONSTANTS[6] = 32.24;
CONSTANTS[7] = 29.09;
CONSTANTS[8] = 5;
STATES[3] = 20;
CONSTANTS[9] = 2.67;
CONSTANTS[10] = 0.7;
CONSTANTS[11] = 0.12;
CONSTANTS[12] = 153;
CONSTANTS[13] = 0.16;
CONSTANTS[14] = 4.85;
CONSTANTS[15] = 0.05;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[0] = (CONSTANTS[0] - ( CONSTANTS[1]*STATES[0]*STATES[1])/(STATES[0]+CONSTANTS[2])) - ( CONSTANTS[3]*STATES[0]*STATES[2])/(STATES[0]+CONSTANTS[4]);
RATES[1] =  CONSTANTS[5]*STATES[0] - ( CONSTANTS[6]*STATES[1])/(STATES[1]+CONSTANTS[7]);
RATES[2] = ((( CONSTANTS[8]*STATES[2]*STATES[1]*STATES[3])/(STATES[3]+CONSTANTS[9])+ CONSTANTS[10]*STATES[1]+ CONSTANTS[11]*STATES[0]) - ( CONSTANTS[12]*STATES[2])/(STATES[2]+CONSTANTS[13])) - ( CONSTANTS[14]*STATES[2])/(STATES[2]+CONSTANTS[15]);
RATES[3] = ( - CONSTANTS[8]*STATES[2]*STATES[1]*STATES[3])/(STATES[3]+CONSTANTS[9])+( CONSTANTS[14]*STATES[2])/(STATES[2]+CONSTANTS[15]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}