Emergent Properties of Networks of Biological Signaling Pathways

Generated code

• The required array sizes are:
  • VARIABLES[10]
  • CONSTANTS[16]
  • BOUND[1]
  • RATES[10]
• The following variables are defined:
  • Variable t(BOUND[0])
  • Variable km1(CONSTANTS[0])
  • Variable Vmax1(CONSTANTS[1])
  • Variable kf6(CONSTANTS[2])
  • Variable kf7(CONSTANTS[3])
  • Variable kb7(CONSTANTS[4])
  • Variable kf8(CONSTANTS[5])
  • Variable kb8(CONSTANTS[6])
  • Variable kf9(CONSTANTS[7])
  • Variable kb9(CONSTANTS[8])
  • Variable p2(CONSTANTS[9])
  • Variable p3(CONSTANTS[10])
  • Variable p4(CONSTANTS[11])
  • Variable p5(CONSTANTS[12])
  • Variable Ca(RATES[0] and VARIABLES[0])
  • Variable Ca_ext(RATES[1] and VARIABLES[1])
  • Variable Extracell_Leak(CONSTANTS[13])
  • Variable Cap_channel(RATES[2] and VARIABLES[2])
  • Variable Cap_channel_star(RATES[3] and VARIABLES[3])
  • Variable Ca2Ca_transp(RATES[4] and VARIABLES[4])
  • Variable Ca_stores(RATES[5] and VARIABLES[5])
  • Variable Stores_Leak(CONSTANTS[14])
  • Variable Ca_transp(RATES[6] and VARIABLES[6])
  • Variable IP3R(RATES[7] and VARIABLES[7])
  • Variable IP3R_star(RATES[8] and VARIABLES[8])
  • Variable IP3(RATES[9] and VARIABLES[9])
  • Variable Ca_pump(CONSTANTS[15])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 0.2;
CONSTANTS[1] = 72;
CONSTANTS[2] = 25;
CONSTANTS[3] = 1e-08;
CONSTANTS[4] = 144;
CONSTANTS[5] = 1.2e-11;
CONSTANTS[6] = 1;
CONSTANTS[7] = 1e-20;
CONSTANTS[8] = 1;
CONSTANTS[9] = 0.004;
CONSTANTS[10] = 8;
CONSTANTS[11] = 0.005;
CONSTANTS[12] = 19.2;
CONSTANTS[13] = 0.00083333;
CONSTANTS[14] = 0.001;
CONSTANTS[15] = 0.00083333;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[0] = 0.08;
VARIABLES[1] = 4000;
VARIABLES[2] = 0.00083333;
VARIABLES[3] = 0;
VARIABLES[4] = 0.08;
VARIABLES[5] = 6.3328;
VARIABLES[6] = 0.0083838;
VARIABLES[7] = 0.016643;
VARIABLES[8] = 0;
VARIABLES[9] = 0;
}

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

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