Generated Code

/*
   There are a total of 0 entries in the algebraic variable array.
   There are a total of 2 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 (minute).
 * STATES[0] is Pha in component Pha (millimolar).
 * CONSTANTS[0] is V_1 in component Pha (per_minute).
 * CONSTANTS[1] is K_1 in component Pha (dimensionless).
 * CONSTANTS[2] is K_2 in component Pha (dimensionless).
 * CONSTANTS[3] is V_M2 in component Pha (per_minute).
 * CONSTANTS[4] is alpha in component Pha (dimensionless).
 * CONSTANTS[5] is GLC in component Pha (millimolar).
 * CONSTANTS[6] is K_a1 in component Pha (millimolar).
 * CONSTANTS[7] is K_a2 in component Pha (millimolar).
 * STATES[1] is GSa in component GSa (millimolar).
 * CONSTANTS[8] is V_M3 in component GSa (per_minute).
 * CONSTANTS[9] is beta in component GSa (dimensionless).
 * CONSTANTS[10] is G6P in component GSa (millimolar).
 * CONSTANTS[11] is K_a3 in component GSa (millimolar).
 * CONSTANTS[12] is K_3 in component GSa (dimensionless).
 * CONSTANTS[13] is K_4 in component GSa (dimensionless).
 * CONSTANTS[14] is K_a4 in component GSa (millimolar).
 * CONSTANTS[15] is V_4 in component GSa (per_minute).
 * RATES[0] is d/dt Pha in component Pha (millimolar).
 * RATES[1] is d/dt GSa in component GSa (millimolar).
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = 0.845;
CONSTANTS[0] = 1.25;
CONSTANTS[1] = 0.1;
CONSTANTS[2] = 0.2;
CONSTANTS[3] = 0.22;
CONSTANTS[4] = 9;
CONSTANTS[5] = 32.5;
CONSTANTS[6] = 10;
CONSTANTS[7] = 10;
STATES[1] = 0.02;
CONSTANTS[8] = 0.05;
CONSTANTS[9] = 9;
CONSTANTS[10] = 0.4;
CONSTANTS[11] = 0.5;
CONSTANTS[12] = 0.4;
CONSTANTS[13] = 0.4;
CONSTANTS[14] = 0.5;
CONSTANTS[15] = 0.2;
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
RATES[0] = ( CONSTANTS[0]*(1.00000 - STATES[0]))/((CONSTANTS[1]+1.00000) - STATES[0]) - ( CONSTANTS[3]*STATES[0]*(1.00000+( CONSTANTS[4]*CONSTANTS[5])/(CONSTANTS[6]+CONSTANTS[5])))/(CONSTANTS[2]/(1.00000+CONSTANTS[5]/CONSTANTS[7])+STATES[0]);
RATES[1] = ( CONSTANTS[8]*(1.00000+( CONSTANTS[9]*CONSTANTS[10])/(CONSTANTS[11]+CONSTANTS[10]))*(CONSTANTS[1]/(CONSTANTS[1]+STATES[0]))*(1.00000 - STATES[1]))/((CONSTANTS[12]/(1.00000+CONSTANTS[10]/CONSTANTS[14])+1.00000) - STATES[1]) - ( CONSTANTS[15]*STATES[1])/(CONSTANTS[13]+ STATES[1]*1.00000);
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}