Emergent Properties of Networks of Biological Signaling Pathways

Generated code

• The required array sizes are:
  • VARIABLES[8]
  • CONSTANTS[12]
  • BOUND[1]
  • RATES[8]
• The following variables are defined:
  • Variable t(BOUND[0])
  • Variable km2(CONSTANTS[0])
  • Variable km6(CONSTANTS[1])
  • Variable Vmax2(CONSTANTS[2])
  • Variable Vmax6(CONSTANTS[3])
  • Variable PIP2(RATES[0] and VARIABLES[0])
  • Variable DAG(RATES[1] and VARIABLES[1])
  • Variable IP3(RATES[2] and VARIABLES[2])
  • Variable kf1(CONSTANTS[4])
  • Variable kb1(CONSTANTS[5])
  • Variable Ca(RATES[3] and VARIABLES[3])
  • Variable EGF_EGFR(CONSTANTS[6])
  • Variable kf3(CONSTANTS[7])
  • Variable kf5(CONSTANTS[8])
  • Variable kb5(CONSTANTS[9])
  • Variable km4(CONSTANTS[10])
  • Variable Vmax4(CONSTANTS[11])
  • Variable PLCg(RATES[4] and VARIABLES[4])
  • Variable PLCg_star(RATES[5] and VARIABLES[5])
  • Variable CaPLCg(RATES[6] and VARIABLES[6])
  • Variable CaPLCg_star(RATES[7] and VARIABLES[7])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 97;
CONSTANTS[1] = 19.7917;
CONSTANTS[2] = 14;
CONSTANTS[3] = 57;
CONSTANTS[4] = 0.0003;
CONSTANTS[5] = 10;
CONSTANTS[6] = 0.1;
CONSTANTS[7] = 0.05;
CONSTANTS[8] = 2e-05;
CONSTANTS[9] = 10;
CONSTANTS[10] = 0.333333;
CONSTANTS[11] = 0.2;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[0] = 10;
VARIABLES[1] = 0;
VARIABLES[2] = 0;
VARIABLES[3] = 0.1;
VARIABLES[4] = 0.82;
VARIABLES[5] = 0;
VARIABLES[6] = 0;
VARIABLES[7] = 0;
}

void ComputeRates(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
RATES[0] = ((((-VARIABLES[6])*CONSTANTS[2]*VARIABLES[0])/(CONSTANTS[0]+VARIABLES[0]))-((VARIABLES[7]*CONSTANTS[3]*VARIABLES[0])/(CONSTANTS[1]+VARIABLES[0])));
RATES[1] = (((VARIABLES[6]*CONSTANTS[2]*VARIABLES[0])/(CONSTANTS[0]+VARIABLES[0]))+((VARIABLES[7]*CONSTANTS[3]*VARIABLES[0])/(CONSTANTS[1]+VARIABLES[0])));
RATES[2] = (((VARIABLES[6]*CONSTANTS[2]*VARIABLES[0])/(CONSTANTS[0]+VARIABLES[0]))+((VARIABLES[7]*CONSTANTS[3]*VARIABLES[0])/(CONSTANTS[1]+VARIABLES[0])));
RATES[4] = (((-VARIABLES[3])*VARIABLES[4]*CONSTANTS[4])+(VARIABLES[6]*CONSTANTS[5]));
RATES[6] = (((-(((-VARIABLES[3])*VARIABLES[4]*CONSTANTS[4])+(VARIABLES[6]*CONSTANTS[5])))+(CONSTANTS[7]*VARIABLES[7]))-((CONSTANTS[11]*CONSTANTS[6]*VARIABLES[6])/(CONSTANTS[10]+VARIABLES[6])));
RATES[7] = (((-((CONSTANTS[7]*VARIABLES[7])-((CONSTANTS[11]*CONSTANTS[6]*VARIABLES[6])/(CONSTANTS[10]+VARIABLES[6]))))+(VARIABLES[3]*VARIABLES[5]*CONSTANTS[8]))-(CONSTANTS[9]*VARIABLES[7]));
RATES[5] = (((-VARIABLES[3])*VARIABLES[5]*CONSTANTS[8])+(CONSTANTS[9]*VARIABLES[7]));
RATES[3] = ((((-VARIABLES[3])*VARIABLES[4]*CONSTANTS[4])+(VARIABLES[6]*CONSTANTS[5])+(CONSTANTS[9]*VARIABLES[7]))-(VARIABLES[3]*VARIABLES[5]*CONSTANTS[8]));
}

void ComputeVariables(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
#ifndef VARIABLES_FOR_RATES_ONLY
#endif
}