Personal tools
You are here: Home Model Repository
 
Views

Emergent Properties of Networks of Biological Signaling Pathways

Generated code

Programming Language

Download this code
  • The required array sizes are:
    • VARIABLES[9]
    • CONSTANTS[20]
    • BOUND[1]
    • RATES[9]
  • The following variables are defined:
    • Variable t(BOUND[0])
    • Variable kf1(CONSTANTS[0])
    • Variable kb1(CONSTANTS[1])
    • Variable kf2(CONSTANTS[2])
    • Variable kb2(CONSTANTS[3])
    • Variable kf3(CONSTANTS[4])
    • Variable kb3(CONSTANTS[5])
    • Variable kf4(CONSTANTS[6])
    • Variable kb4(CONSTANTS[7])
    • Variable kf5(CONSTANTS[8])
    • Variable kb5(CONSTANTS[9])
    • Variable kf6(CONSTANTS[10])
    • Variable kb6(CONSTANTS[11])
    • Variable kf7(CONSTANTS[12])
    • Variable DAG(RATES[0] and VARIABLES[0])
    • Variable Ca(RATES[1] and VARIABLES[1])
    • Variable DAGPKC(RATES[2] and VARIABLES[2])
    • Variable AADAGPKC(RATES[3] and VARIABLES[3])
    • Variable DAGCaPKC(RATES[4] and VARIABLES[4])
    • Variable CaPKC(RATES[5] and VARIABLES[5])
    • Variable PKC_a(RATES[6] and VARIABLES[6])
    • Variable kb7(CONSTANTS[13])
    • Variable kf8(CONSTANTS[14])
    • Variable kb8(CONSTANTS[15])
    • Variable kf9(CONSTANTS[16])
    • Variable kb9(CONSTANTS[17])
    • Variable kf10(CONSTANTS[18])
    • Variable kb10(CONSTANTS[19])
    • Variable PKC_i(RATES[7] and VARIABLES[7])
    • Variable AA(RATES[8] and VARIABLES[8])
  • No equation needed for Newton-Raphson evaluation.
void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 50;
CONSTANTS[2] = 2e-10;
CONSTANTS[3] = 0.1;
CONSTANTS[4] = 1.2705;
CONSTANTS[5] = 3.5026;
CONSTANTS[6] = 2e-09;
CONSTANTS[7] = 0.1;
CONSTANTS[8] = 1;
CONSTANTS[9] = 0.1;
CONSTANTS[10] = 2;
CONSTANTS[11] = 0.2;
CONSTANTS[12] = 1e-06;
CONSTANTS[13] = 0.5;
CONSTANTS[14] = 1.3333e-08;
CONSTANTS[15] = 8.6348;
CONSTANTS[16] = 1e-09;
CONSTANTS[17] = 0.1;
CONSTANTS[18] = 3e-08;
CONSTANTS[19] = 2;
}

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

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

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