/*
   There are a total of 3 entries in the algebraic variable array.
   There are a total of 1 entries in each of the rate and state variable arrays.
   There are a total of 14 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (minute).
 * STATES[0] is Pha in component Pha (dimensionless).
 * ALGEBRAIC[1] is V_1 in component V_1 (dimensionless).
 * ALGEBRAIC[0] is Z in component Z (dimensionless).
 * ALGEBRAIC[2] is K_1 in component K_1 (dimensionless).
 * CONSTANTS[0] is V_M2 in component Pha (per_minute).
 * CONSTANTS[1] is alpha in component Pha (dimensionless).
 * CONSTANTS[2] is Glc in component Pha (micromolar).
 * CONSTANTS[3] is K_2 in component Pha (dimensionless).
 * CONSTANTS[4] is K_a1 in component Pha (micromolar).
 * CONSTANTS[5] is K_a2 in component Pha (micromolar).
 * CONSTANTS[6] is V_M1 in component V_1 (per_minute).
 * CONSTANTS[7] is gamma in component V_1 (dimensionless).
 * CONSTANTS[8] is K_a5 in component V_1 (micromolar).
 * CONSTANTS[9] is K_11 in component K_1 (dimensionless).
 * CONSTANTS[10] is K_a6 in component K_1 (micromolar).
 * CONSTANTS[11] is A in component Z (micromolar).
 * CONSTANTS[12] is B in component Z (micromolar).
 * CONSTANTS[13] is v in component Z (per_minute).
 * RATES[0] is d/dt Pha in component Pha (dimensionless).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 0.098;
CONSTANTS[0] = 0.6;
CONSTANTS[1] = 9;
CONSTANTS[2] = 10000;
CONSTANTS[3] = 0.2;
CONSTANTS[4] = 10000;
CONSTANTS[5] = 10000;
CONSTANTS[6] = 1.5;
CONSTANTS[7] = 9;
CONSTANTS[8] = 0.5;
CONSTANTS[9] = 0.1;
CONSTANTS[10] = 0.5;
CONSTANTS[11] = 0.3;
CONSTANTS[12] = 0.25;
CONSTANTS[13] = 1;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = CONSTANTS[11]+ CONSTANTS[12]* sin( 2.00000* 3.14159265358979*CONSTANTS[13]*VOI);
ALGEBRAIC[1] =  CONSTANTS[6]*(1.00000+ CONSTANTS[7]*(pow(ALGEBRAIC[0], 4.00000)/pow(CONSTANTS[8], 4.00000)+pow(ALGEBRAIC[0], 4.00000)));
ALGEBRAIC[2] = CONSTANTS[9]/(1.00000+pow(ALGEBRAIC[0], 4.00000)/pow(CONSTANTS[10], 4.00000));
RATES[0] =  ALGEBRAIC[1]*ALGEBRAIC[0]*((1.00000 - STATES[0])/(( ALGEBRAIC[2]*ALGEBRAIC[0]+1.00000) - STATES[0])) - ( CONSTANTS[0]*(1.00000+( CONSTANTS[1]*CONSTANTS[2])/(CONSTANTS[4]+CONSTANTS[2]))*STATES[0])/(CONSTANTS[3]/(1.00000+CONSTANTS[2]/CONSTANTS[5])+STATES[0]);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] = CONSTANTS[11]+ CONSTANTS[12]* sin( 2.00000* 3.14159265358979*CONSTANTS[13]*VOI);
ALGEBRAIC[1] =  CONSTANTS[6]*(1.00000+ CONSTANTS[7]*(pow(ALGEBRAIC[0], 4.00000)/pow(CONSTANTS[8], 4.00000)+pow(ALGEBRAIC[0], 4.00000)));
ALGEBRAIC[2] = CONSTANTS[9]/(1.00000+pow(ALGEBRAIC[0], 4.00000)/pow(CONSTANTS[10], 4.00000));
}