A kinetic model of the thiazide-sensitive Na-Cl cotransporter

Generated code

• The required array sizes are:
  • VARIABLES[12]
  • CONSTANTS[38]
  • BOUND[1]
  • RATES[12]
• The following variables are defined:
  • Variable ENaD_(RATES[0] and VARIABLES[0])
  • Variable ED(RATES[1] and VARIABLES[1])
  • Variable time(BOUND[0])
  • Variable E(RATES[2] and VARIABLES[2])
  • Variable ED_(RATES[3] and VARIABLES[3])
  • Variable ENaD(RATES[4] and VARIABLES[4])
  • Variable ENaCl(RATES[5] and VARIABLES[5])
  • Variable ENaCl_(RATES[6] and VARIABLES[6])
  • Variable k1(CONSTANTS[0])
  • Variable k2(CONSTANTS[1])
  • Variable k3(CONSTANTS[2])
  • Variable k4(CONSTANTS[3])
  • Variable k5(CONSTANTS[4])
  • Variable k6(CONSTANTS[5])
  • Variable k7(CONSTANTS[6])
  • Variable k8(CONSTANTS[7])
  • Variable k9(CONSTANTS[8])
  • Variable k10(CONSTANTS[9])
  • Variable k11(CONSTANTS[10])
  • Variable k12(CONSTANTS[11])
  • Variable k13(CONSTANTS[12])
  • Variable k14(CONSTANTS[13])
  • Variable k15(CONSTANTS[14])
  • Variable k16(CONSTANTS[15])
  • Variable k17(CONSTANTS[16])
  • Variable k18(CONSTANTS[17])
  • Variable k19(CONSTANTS[18])
  • Variable k20(CONSTANTS[19])
  • Variable k21(CONSTANTS[20])
  • Variable k22(CONSTANTS[21])
  • Variable k23(CONSTANTS[22])
  • Variable k24(CONSTANTS[23])
  • Variable k25(CONSTANTS[24])
  • Variable k26(CONSTANTS[25])
  • Variable k27(CONSTANTS[26])
  • Variable k28(CONSTANTS[27])
  • Variable k29(CONSTANTS[28])
  • Variable ENa(RATES[7] and VARIABLES[7])
  • Variable ENa_(RATES[8] and VARIABLES[8])
  • Variable E_(RATES[9] and VARIABLES[9])
  • Variable ECl(RATES[10] and VARIABLES[10])
  • Variable ECl_(RATES[11] and VARIABLES[11])
  • Variable k30(CONSTANTS[29])
  • Variable k31(CONSTANTS[30])
  • Variable k32(CONSTANTS[31])
  • Variable Na(CONSTANTS[32])
  • Variable Cl(CONSTANTS[33])
  • Variable D(CONSTANTS[34])
  • Variable Cl_(CONSTANTS[35])
  • Variable Na_(CONSTANTS[36])
  • Variable D_(CONSTANTS[37])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 1e+08;
CONSTANTS[1] = 418300;
CONSTANTS[2] = 1e+08;
CONSTANTS[3] = 4.928e+06;
CONSTANTS[4] = 1e+08;
CONSTANTS[5] = 1.065e+06;
CONSTANTS[6] = 1e+08;
CONSTANTS[7] = 89400;
CONSTANTS[8] = 1e+08;
CONSTANTS[9] = 418300;
CONSTANTS[10] = 1e+08;
CONSTANTS[11] = 4.982e+06;
CONSTANTS[12] = 1e+08;
CONSTANTS[13] = 1.065e+06;
CONSTANTS[14] = 1e+08;
CONSTANTS[15] = 89400;
CONSTANTS[16] = 458700;
CONSTANTS[17] = 100000;
CONSTANTS[18] = 1000;
CONSTANTS[19] = 218;
CONSTANTS[20] = 1e+08;
CONSTANTS[21] = 418300;
CONSTANTS[22] = 1e+08;
CONSTANTS[23] = 31.92;
CONSTANTS[24] = 1e+08;
CONSTANTS[25] = 31.92;
CONSTANTS[26] = 1e+08;
CONSTANTS[27] = 138900;
CONSTANTS[28] = 1e+08;
CONSTANTS[29] = 0.3514;
CONSTANTS[30] = 1e+08;
CONSTANTS[31] = 0.1166;
CONSTANTS[32] = 50;
CONSTANTS[33] = 96;
CONSTANTS[34] = 1;
CONSTANTS[35] = 40;
CONSTANTS[36] = 10;
CONSTANTS[37] = 1;
}

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

void ComputeRates(double* BOUND, double* RATES, double* CONSTANTS, double* VARIABLES)
{
RATES[2] = (((CONSTANTS[23]*VARIABLES[1])+(CONSTANTS[17]*VARIABLES[9])+(CONSTANTS[1]*VARIABLES[7])+(CONSTANTS[3]*VARIABLES[10]))-((CONSTANTS[22]*CONSTANTS[34]*VARIABLES[2])+(CONSTANTS[16]*VARIABLES[2])+(CONSTANTS[0]*CONSTANTS[32]*VARIABLES[2])+(CONSTANTS[2]*CONSTANTS[33]*VARIABLES[2])));
RATES[9] = (((CONSTANTS[29]*VARIABLES[3])+(CONSTANTS[16]*VARIABLES[2])+(CONSTANTS[9]*VARIABLES[8])+(CONSTANTS[11]*VARIABLES[11]))-((CONSTANTS[28]*CONSTANTS[37]*VARIABLES[9])+(CONSTANTS[17]*VARIABLES[9])+(CONSTANTS[8]*CONSTANTS[36]*VARIABLES[9])+(CONSTANTS[10]*CONSTANTS[35]*VARIABLES[9])));
RATES[1] = (((CONSTANTS[22]*VARIABLES[2]*CONSTANTS[34])+(CONSTANTS[21]*VARIABLES[4]))-((CONSTANTS[23]*VARIABLES[1])+(CONSTANTS[20]*CONSTANTS[32]*VARIABLES[1])));
RATES[3] = (((CONSTANTS[28]*VARIABLES[9]*CONSTANTS[37])+((CONSTANTS[26]) * 0.00100000*VARIABLES[0]))-((CONSTANTS[29]*VARIABLES[3])+(CONSTANTS[26]*CONSTANTS[36]*VARIABLES[3])));
RATES[4] = (((CONSTANTS[20]*CONSTANTS[32]*VARIABLES[1])+(CONSTANTS[24]*CONSTANTS[34]*VARIABLES[7]))-((CONSTANTS[21]*VARIABLES[4])+(CONSTANTS[25]*VARIABLES[4])));
RATES[0] = (((CONSTANTS[26]*CONSTANTS[36]*VARIABLES[3])+(CONSTANTS[30]*CONSTANTS[37]*VARIABLES[8]))-((CONSTANTS[27]*VARIABLES[0])+(CONSTANTS[31]*VARIABLES[0])));
RATES[7] = (((CONSTANTS[0]*CONSTANTS[32]*VARIABLES[2])+(CONSTANTS[25]*VARIABLES[4])+(CONSTANTS[5]*VARIABLES[5]))-((CONSTANTS[1]*VARIABLES[7])+(CONSTANTS[4]*CONSTANTS[33]*VARIABLES[7])+(CONSTANTS[24]*CONSTANTS[34]*VARIABLES[7])));
RATES[8] = (((CONSTANTS[8]*CONSTANTS[36]*VARIABLES[9])+(CONSTANTS[31]*VARIABLES[0])+(CONSTANTS[13]*VARIABLES[6]))-((CONSTANTS[9]*VARIABLES[8])+(CONSTANTS[12]*CONSTANTS[35]*VARIABLES[8])+(CONSTANTS[30]*CONSTANTS[37]*VARIABLES[8])));
RATES[10] = (((CONSTANTS[2]*CONSTANTS[33]*VARIABLES[2])+(CONSTANTS[7]*VARIABLES[5]))-((CONSTANTS[6]*CONSTANTS[32]*VARIABLES[10])+(CONSTANTS[3]*VARIABLES[10])));
RATES[11] = (((CONSTANTS[10]*CONSTANTS[35]*VARIABLES[9])+(CONSTANTS[15]*VARIABLES[5]))-((CONSTANTS[14]*CONSTANTS[36]*VARIABLES[11])+(CONSTANTS[11]*VARIABLES[11])));
RATES[5] = (((CONSTANTS[4]*CONSTANTS[33]*VARIABLES[7])+(CONSTANTS[6]*CONSTANTS[32]*VARIABLES[10])+(CONSTANTS[19]*VARIABLES[6]))-((CONSTANTS[5]*VARIABLES[5])+(CONSTANTS[7]*VARIABLES[5])+(CONSTANTS[18]*VARIABLES[5])));
RATES[6] = (((CONSTANTS[12]*CONSTANTS[35]*VARIABLES[8])+(CONSTANTS[14]*CONSTANTS[36]*VARIABLES[11])+(CONSTANTS[18]*VARIABLES[5]))-((CONSTANTS[13]*VARIABLES[6])+(CONSTANTS[15]*VARIABLES[6])+(CONSTANTS[19]*VARIABLES[6])));
}

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