Emergent Properties of Networks of Biological Signaling Pathways

Generated code

• The required array sizes are:
  • VARIABLES[11]
  • CONSTANTS[15]
  • BOUND[1]
  • RATES[11]
• The following variables are defined:
  • Variable kf8(CONSTANTS[0])
  • Variable kf9(CONSTANTS[1])
  • Variable km6(CONSTANTS[2])
  • Variable Vmax6(CONSTANTS[3])
  • Variable km7(CONSTANTS[4])
  • Variable Vmax7(CONSTANTS[5])
  • Variable PIP2(RATES[0] and VARIABLES[0])
  • Variable DAG(RATES[1] and VARIABLES[1])
  • Variable IP3(RATES[2] and VARIABLES[2])
  • Variable Inositol(RATES[3] and VARIABLES[3])
  • Variable PC(RATES[4] and VARIABLES[4])
  • Variable kf1(CONSTANTS[6])
  • Variable kb1(CONSTANTS[7])
  • Variable kf2(CONSTANTS[8])
  • Variable kb2(CONSTANTS[9])
  • Variable kf3(CONSTANTS[10])
  • Variable kb3(CONSTANTS[11])
  • Variable kf4(CONSTANTS[12])
  • Variable kb4(CONSTANTS[13])
  • Variable kf5(CONSTANTS[14])
  • Variable PLC(RATES[5] and VARIABLES[5])
  • Variable Gq(RATES[6] and VARIABLES[6])
  • Variable Ca(RATES[7] and VARIABLES[7])
  • Variable GqPLC(RATES[8] and VARIABLES[8])
  • Variable CaGqPLC(RATES[9] and VARIABLES[9])
  • Variable CaPLC(RATES[10] and VARIABLES[10])
  • Variable t(BOUND[0])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 0.15;
CONSTANTS[1] = 2.5;
CONSTANTS[2] = 5;
CONSTANTS[3] = 48;
CONSTANTS[4] = 19.8417;
CONSTANTS[5] = 10;
CONSTANTS[6] = 5e-06;
CONSTANTS[7] = 1;
CONSTANTS[8] = 4.2e-06;
CONSTANTS[9] = 1;
CONSTANTS[10] = 5e-05;
CONSTANTS[11] = 1;
CONSTANTS[12] = 4.2e-06;
CONSTANTS[13] = 1;
CONSTANTS[14] = 0.0133;
}

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

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

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