Comparison of Simulated and Measured Calcium Sparks in Intact Skeletal Muscle Fibers of the Frog

Generated code

• The required array sizes are:
  • VARIABLES[6]
  • CONSTANTS[8]
  • BOUND[1]
  • RATES[6]
• The following variables are defined:
  • Variable PrFluo(RATES[0] and VARIABLES[0])
  • Variable CaFluo(RATES[1] and VARIABLES[1])
  • Variable CaPrFluo(RATES[2] and VARIABLES[2])
  • Variable time(BOUND[0])
  • Variable Ca(RATES[3] and VARIABLES[3])
  • Variable k1(CONSTANTS[0])
  • Variable k1_(CONSTANTS[1])
  • Variable k2(CONSTANTS[2])
  • Variable k2_(CONSTANTS[3])
  • Variable k3(CONSTANTS[4])
  • Variable k3_(CONSTANTS[5])
  • Variable k4(CONSTANTS[6])
  • Variable k4_(CONSTANTS[7])
  • Variable Fluo(RATES[4] and VARIABLES[4])
  • Variable Pr(RATES[5] and VARIABLES[5])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 2.676e+08;
CONSTANTS[1] = 137;
CONSTANTS[2] = 1.72e+07;
CONSTANTS[3] = 32.9;
CONSTANTS[4] = 1.149e+07;
CONSTANTS[5] = 4216;
CONSTANTS[6] = 1.149e+07;
CONSTANTS[7] = 15777;
}

void SetupComputedConstants(double* CONSTANTS, double* VARIABLES)
{
VARIABLES[0] = 88.12;
VARIABLES[1] = 0;
VARIABLES[2] = 0;
VARIABLES[3] = 0.05;
VARIABLES[4] = 11.88;
VARIABLES[5] = 3000;
}

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

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