Generated Code

/*
   There are a total of 5 entries in the algebraic variable array.
   There are a total of 3 entries in each of the rate and state variable arrays.
   There are a total of 14 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (second).
 * STATES[0] is V in component membrane (millivolt).
 * CONSTANTS[0] is Cm in component membrane (picoF).
 * ALGEBRAIC[0] is i_s in component calcium_channel (femtoA).
 * ALGEBRAIC[1] is i_K in component potassium_channel (femtoA).
 * ALGEBRAIC[2] is i_K_ACh in component acetyl_choline_activated_potassium_channel (femtoA).
 * ALGEBRAIC[4] is i_j in component coupling_current (femtoA).
 * CONSTANTS[1] is g_s in component calcium_channel (picoS).
 * CONSTANTS[2] is V_s in component calcium_channel (millivolt).
 * CONSTANTS[3] is V_1 in component calcium_channel (millivolt).
 * CONSTANTS[4] is V_2 in component calcium_channel (millivolt).
 * CONSTANTS[5] is g_K in component potassium_channel (picoS).
 * CONSTANTS[6] is V_K in component potassium_channel (millivolt).
 * STATES[1] is w in component potassium_channel_w_gate (dimensionless).
 * CONSTANTS[7] is lambda_w in component potassium_channel_w_gate (per_second).
 * CONSTANTS[8] is V_3 in component potassium_channel_w_gate (millivolt).
 * CONSTANTS[9] is V_4 in component potassium_channel_w_gate (millivolt).
 * STATES[2] is u in component acetyl_choline_activated_potassium_channel_u_gate (dimensionless).
 * CONSTANTS[13] is alpha in component acetyl_choline_activated_potassium_channel_u_gate (per_second).
 * ALGEBRAIC[3] is beta in component acetyl_choline_activated_potassium_channel_u_gate (per_second).
 * CONSTANTS[10] is ACh in component acetyl_choline_activated_potassium_channel_u_gate (molar).
 * CONSTANTS[11] is g_j in component coupling_current (picoS).
 * CONSTANTS[12] is V_B in component coupling_current (millivolt).
 * RATES[0] is d/dt V in component membrane (millivolt).
 * RATES[1] is d/dt w in component potassium_channel_w_gate (dimensionless).
 * RATES[2] is d/dt u in component acetyl_choline_activated_potassium_channel_u_gate (dimensionless).
 * There are a total of 0 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
STATES[0] = -52.07606;
CONSTANTS[0] = 60;
CONSTANTS[1] = 382.9118;
CONSTANTS[2] = 214.1429;
CONSTANTS[3] = -35.9358;
CONSTANTS[4] = 7.8589;
CONSTANTS[5] = 536.1093;
CONSTANTS[6] = -259.0783;
STATES[1] = 0.0008971;
CONSTANTS[7] = 20.7796;
CONSTANTS[8] = -27.9375;
CONSTANTS[9] = 6.321;
STATES[2] = 0.2344555;
CONSTANTS[10] = 1e-6;
CONSTANTS[11] = 0;
CONSTANTS[12] = -50;
CONSTANTS[13] = 0.0123320/(1.00000+4.20000e-06/CONSTANTS[10]);
RATES[0] = 0.1001;
RATES[1] = 0.1001;
RATES[2] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
                 double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - - (ALGEBRAIC[0]+ALGEBRAIC[1]+ALGEBRAIC[2]+ALGEBRAIC[4])/CONSTANTS[0];
resid[1] = RATES[1] -  CONSTANTS[7]*cosh((STATES[0] - CONSTANTS[8])/( 2.00000*CONSTANTS[9]))*( (1.00000/2.00000)*(1.00000+ tanh((STATES[0] - CONSTANTS[8])/CONSTANTS[9])) - STATES[1]);
resid[2] = RATES[2] -  CONSTANTS[13]*(1.00000 - STATES[2]) -  ALGEBRAIC[3]*STATES[2];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[0] =  (1.00000/2.00000)*CONSTANTS[1]*(1.00000+ tanh((STATES[0] - CONSTANTS[3])/CONSTANTS[4]))*(STATES[0] - CONSTANTS[2]);
ALGEBRAIC[1] =  CONSTANTS[5]*STATES[1]*(STATES[0] - CONSTANTS[6]);
ALGEBRAIC[2] =  1.00000*0.270000*STATES[2]*(STATES[0]+90.0000);
ALGEBRAIC[3] =  0.0100000*exp( 0.0133000*(STATES[0]+40.0000));
ALGEBRAIC[4] =  CONSTANTS[11]*(STATES[0] - CONSTANTS[12]);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
}