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

Generated code

• The required array sizes are:
  • VARIABLES[5]
  • CONSTANTS[4]
  • BOUND[1]
  • RATES[5]
• The following variables are defined:
  • Variable MgParv(RATES[0] and VARIABLES[0])
  • Variable k1(CONSTANTS[0])
  • Variable k1_(CONSTANTS[1])
  • Variable k2(CONSTANTS[2])
  • Variable k2_(CONSTANTS[3])
  • Variable time(BOUND[0])
  • Variable Ca(RATES[1] and VARIABLES[1])
  • Variable Parv(RATES[2] and VARIABLES[2])
  • Variable CaParv(RATES[3] and VARIABLES[3])
  • Variable Mg(RATES[4] and VARIABLES[4])
• No equation needed for Newton-Raphson evaluation.

void SetupFixedConstants(double* CONSTANTS)
{
CONSTANTS[0] = 4.79e+07;
CONSTANTS[1] = 0.574;
CONSTANTS[2] = 39700;
CONSTANTS[3] = 3.45;
}

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

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

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