Personal tools
You are here: Home Model Repository
 
Views

Risks and benefits of structured antiretroviral drug therapy interruptions in HIV-1 infection

Generated code

Programming Language

Download this code
  • The required array sizes are:
    • VARIABLES[4]
    • CONSTANTS[13]
    • BOUND[1]
    • RATES[4]
  • The following variables are defined:
    • Variable time(BOUND[0])
    • Variable T(RATES[0] and VARIABLES[0])
    • Variable s(CONSTANTS[0])
    • Variable dT(CONSTANTS[1])
    • Variable I(RATES[1] and VARIABLES[1])
    • Variable al(CONSTANTS[2])
    • Variable E(RATES[2] and VARIABLES[2])
    • Variable c(CONSTANTS[3])
    • Variable dE(CONSTANTS[4])
    • Variable K(CONSTANTS[5])
    • Variable p(CONSTANTS[6])
    • Variable dI(CONSTANTS[7])
    • Variable Il(RATES[3] and VARIABLES[3])
    • Variable dE_(CONSTANTS[8])
    • Variable K_(CONSTANTS[9])
    • Variable b(CONSTANTS[10])
    • Variable qa(CONSTANTS[11])
    • Variable ql(CONSTANTS[12])
  • No equation needed for Newton-Raphson evaluation.
void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 10;
CONSTANTS[1] = 0.01;
CONSTANTS[2] = 0.01;
CONSTANTS[3] = 0.3;
CONSTANTS[4] = 0.1;
CONSTANTS[5] = 0.1;
CONSTANTS[6] = 0.05;
CONSTANTS[7] = 0.3;
CONSTANTS[8] = 0.25;
CONSTANTS[9] = 5;
CONSTANTS[10] = 0.001;
CONSTANTS[11] = 0.001;
CONSTANTS[12] = 0.001;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[0] = 0;
VARIABLES[1] = 0;
VARIABLES[2] = 0;
VARIABLES[3] = 0;
}

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

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