/*
   There are a total of 0 entries in the algebraic variable array.
   There are a total of 0 entries in each of the rate and state variable arrays.
   There are a total of 21 entries in the constant variable array.
 */
/*
 * CONSTANTS[16] is v_cystathionine in component v_cystathionine (flux).
 * CONSTANTS[0] is Cys in component Cys (micromolar).
 * CONSTANTS[1] is CGS in component CGS (micromolar).
 * CONSTANTS[2] is Pi in component Pi (micromolar).
 * CONSTANTS[3] is Phser in component Phser (micromolar).
 * CONSTANTS[15] is Km_CGS_app_Cys in component v_cystathionine (micromolar).
 * CONSTANTS[4] is Km_CGS_Cys in component v_cystathionine (micromolar).
 * CONSTANTS[5] is kcat_CGS in component v_cystathionine (first_order_rate_constant).
 * CONSTANTS[14] is kcat_CGS_app_Cys in component v_cystathionine (first_order_rate_constant).
 * CONSTANTS[6] is Km_CGS_Phser in component v_cystathionine (micromolar).
 * CONSTANTS[7] is Ki_CGS_Pi in component v_cystathionine (micromolar).
 * CONSTANTS[19] is v_Thr in component v_Thr (flux).
 * CONSTANTS[8] is TS in component TS (micromolar).
 * CONSTANTS[9] is AdoMet in component AdoMet (micromolar).
 * CONSTANTS[18] is Km_TS in component v_Thr (micromolar).
 * CONSTANTS[10] is kcat_TS_noAdoMet in component v_Thr (first_order_rate_constant).
 * CONSTANTS[11] is kcat_TS_AdoMet in component v_Thr (first_order_rate_constant).
 * CONSTANTS[17] is kcat_TS in component v_Thr (first_order_rate_constant).
 * CONSTANTS[12] is K1K2 in component v_Thr (micromolar2).
 * CONSTANTS[13] is Ki_TS_Pi in component v_Thr (micromolar).
 * CONSTANTS[20] is J_Phser in component J_Phser (flux).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 15.0;
CONSTANTS[1] = 0.7;
CONSTANTS[2] = 10000.0;
CONSTANTS[3] = 500.0;
CONSTANTS[4] = 460.0;
CONSTANTS[5] = 30.0;
CONSTANTS[6] = 2500.0;
CONSTANTS[7] = 2500.0;
CONSTANTS[8] = 5.0;
CONSTANTS[9] = 20.0;
CONSTANTS[10] = 0.42;
CONSTANTS[11] = 3.5;
CONSTANTS[12] = 73.0;
CONSTANTS[13] = 1000.0;
CONSTANTS[14] = CONSTANTS[5]/(1.00000+ (CONSTANTS[6]/CONSTANTS[3])*(1.00000+CONSTANTS[2]/CONSTANTS[7]));
CONSTANTS[15] = CONSTANTS[4]/(1.00000+ (CONSTANTS[6]/CONSTANTS[3])*(1.00000+CONSTANTS[2]/CONSTANTS[7]));
CONSTANTS[16] = ( CONSTANTS[14]*CONSTANTS[1]*CONSTANTS[0])/(CONSTANTS[15]+CONSTANTS[0]);
CONSTANTS[17] = (CONSTANTS[10]+ CONSTANTS[11]*(pow(CONSTANTS[9], 2.00000)/CONSTANTS[12]))/(1.00000+pow(CONSTANTS[9], 2.00000)/CONSTANTS[12]);
CONSTANTS[18] =  (( 250.000*((1.00000+CONSTANTS[9]/0.500000)/(1.00000+CONSTANTS[9]/1.10000)))/(1.00000+pow(CONSTANTS[9], 2.00000)/140.000))*(1.00000+CONSTANTS[2]/CONSTANTS[13]);
CONSTANTS[19] = ( CONSTANTS[8]*CONSTANTS[17]*CONSTANTS[3])/(CONSTANTS[18]+CONSTANTS[3]);
CONSTANTS[20] = CONSTANTS[16]+CONSTANTS[19];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}