Generated Code

The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

/*
   There are a total of 65 entries in the algebraic variable array.
   There are a total of 36 entries in each of the rate and state variable arrays.
   There are a total of 100 entries in the constant variable array.
 */
/*
 * VOI is time in component environment (ms).
 * CONSTANTS[0] is lambda in component environment (dimensionless).
 * CONSTANTS[1] is dlambdadt in component environment (per_ms).
 * STATES[0] is V in component membrane (mV).
 * CONSTANTS[2] is R in component membrane (gas_constant_units).
 * CONSTANTS[3] is T in component membrane (kelvin).
 * CONSTANTS[4] is F in component membrane (faradays_constant_units).
 * ALGEBRAIC[54] is i_Na in component fast_sodium_current (uA_per_mm2).
 * ALGEBRAIC[55] is i_Ca_L_Ca in component L_type_Ca_channel (uA_per_mm2).
 * ALGEBRAIC[64] is i_Ca_L_K in component L_type_Ca_channel (uA_per_mm2).
 * ALGEBRAIC[57] is i_K in component time_dependent_potassium_current (uA_per_mm2).
 * ALGEBRAIC[58] is i_K1 in component time_independent_potassium_current (uA_per_mm2).
 * ALGEBRAIC[29] is i_NaCa in component Na_Ca_exchanger (uA_per_mm2).
 * ALGEBRAIC[28] is i_Kp in component plateau_potassium_current (uA_per_mm2).
 * ALGEBRAIC[30] is i_p_Ca in component sarcolemmal_calcium_pump (uA_per_mm2).
 * ALGEBRAIC[32] is i_Na_b in component sodium_background_current (uA_per_mm2).
 * ALGEBRAIC[34] is i_Ca_b in component calcium_background_current (uA_per_mm2).
 * ALGEBRAIC[36] is i_NaK in component sodium_potassium_pump (uA_per_mm2).
 * ALGEBRAIC[61] is i_ns_Ca in component non_specific_calcium_activated_current (uA_per_mm2).
 * CONSTANTS[5] is Cm in component membrane (uF_per_mm2).
 * ALGEBRAIC[0] is I_stim in component membrane (uA_per_mm2).
 * CONSTANTS[6] is stim_start in component membrane (ms).
 * CONSTANTS[7] is stim_end in component membrane (ms).
 * CONSTANTS[8] is stim_period in component membrane (ms).
 * CONSTANTS[9] is stim_duration in component membrane (ms).
 * CONSTANTS[10] is stim_amplitude in component membrane (uA_per_mm2).
 * ALGEBRAIC[1] is E_Na in component fast_sodium_current (mV).
 * CONSTANTS[11] is g_Na in component fast_sodium_current (mS_per_mm2).
 * STATES[1] is Nai in component ionic_concentrations (mM).
 * CONSTANTS[12] is Nao in component ionic_concentrations (mM).
 * STATES[2] is m in component fast_sodium_current_m_gate (dimensionless).
 * STATES[3] is h in component fast_sodium_current_h_gate (dimensionless).
 * STATES[4] is j in component fast_sodium_current_j_gate (dimensionless).
 * ALGEBRAIC[2] is alpha_m in component fast_sodium_current_m_gate (per_ms).
 * ALGEBRAIC[3] is beta_m in component fast_sodium_current_m_gate (per_ms).
 * ALGEBRAIC[4] is alpha_h in component fast_sodium_current_h_gate (per_ms).
 * ALGEBRAIC[5] is beta_h in component fast_sodium_current_h_gate (per_ms).
 * ALGEBRAIC[6] is alpha_j in component fast_sodium_current_j_gate (per_ms).
 * ALGEBRAIC[7] is beta_j in component fast_sodium_current_j_gate (per_ms).
 * CONSTANTS[13] is P_Ca in component L_type_Ca_channel (mm_per_ms).
 * CONSTANTS[14] is P_K in component L_type_Ca_channel (mm_per_ms).
 * ALGEBRAIC[56] is p_k in component L_type_Ca_channel (mm_per_ms).
 * CONSTANTS[15] is i_Ca_L_Ca_half in component L_type_Ca_channel (uA_per_mm2).
 * ALGEBRAIC[8] is i_Ca_L_Ca_max in component L_type_Ca_channel (uA_per_mm2).
 * STATES[5] is O in component L_type_Ca_channel (dimensionless).
 * STATES[6] is O_Ca in component L_type_Ca_channel (dimensionless).
 * ALGEBRAIC[9] is alpha in component L_type_Ca_channel (per_ms).
 * ALGEBRAIC[10] is beta in component L_type_Ca_channel (per_ms).
 * ALGEBRAIC[13] is gamma in component L_type_Ca_channel (per_ms).
 * ALGEBRAIC[11] is alpha_a in component L_type_Ca_channel (per_ms).
 * ALGEBRAIC[12] is beta_b in component L_type_Ca_channel (per_ms).
 * CONSTANTS[16] is a in component L_type_Ca_channel (dimensionless).
 * CONSTANTS[17] is b in component L_type_Ca_channel (dimensionless).
 * CONSTANTS[18] is g in component L_type_Ca_channel (per_ms).
 * CONSTANTS[19] is f in component L_type_Ca_channel (per_ms).
 * CONSTANTS[20] is g_ in component L_type_Ca_channel (per_ms).
 * CONSTANTS[21] is f_ in component L_type_Ca_channel (per_ms).
 * CONSTANTS[22] is omega in component L_type_Ca_channel (per_ms).
 * STATES[7] is C0 in component L_type_Ca_channel (dimensionless).
 * STATES[8] is C1 in component L_type_Ca_channel (dimensionless).
 * STATES[9] is C2 in component L_type_Ca_channel (dimensionless).
 * STATES[10] is C3 in component L_type_Ca_channel (dimensionless).
 * STATES[11] is C4 in component L_type_Ca_channel (dimensionless).
 * STATES[12] is C_Ca0 in component L_type_Ca_channel (dimensionless).
 * STATES[13] is C_Ca1 in component L_type_Ca_channel (dimensionless).
 * STATES[14] is C_Ca2 in component L_type_Ca_channel (dimensionless).
 * STATES[15] is C_Ca3 in component L_type_Ca_channel (dimensionless).
 * STATES[16] is C_Ca4 in component L_type_Ca_channel (dimensionless).
 * STATES[17] is Ca_SS in component calcium_subsystem (mM).
 * CONSTANTS[23] is Cao in component ionic_concentrations (mM).
 * STATES[18] is Ko in component ionic_concentrations (mM).
 * STATES[19] is Ki in component ionic_concentrations (mM).
 * STATES[20] is y in component L_type_Ca_channel_y_gate (dimensionless).
 * ALGEBRAIC[14] is y_infinity in component L_type_Ca_channel_y_gate (dimensionless).
 * ALGEBRAIC[15] is tau_y in component L_type_Ca_channel_y_gate (ms).
 * ALGEBRAIC[16] is g_K in component time_dependent_potassium_current (mS_per_mm2).
 * CONSTANTS[24] is g_K_max in component time_dependent_potassium_current (mS_per_mm2).
 * ALGEBRAIC[17] is E_K in component time_dependent_potassium_current (mV).
 * CONSTANTS[25] is P_NaK in component time_dependent_potassium_current (dimensionless).
 * STATES[21] is X in component time_dependent_potassium_current_X_gate (dimensionless).
 * ALGEBRAIC[20] is Xi in component time_dependent_potassium_current_Xi_gate (dimensionless).
 * ALGEBRAIC[18] is alpha_X in component time_dependent_potassium_current_X_gate (per_ms).
 * ALGEBRAIC[19] is beta_X in component time_dependent_potassium_current_X_gate (per_ms).
 * ALGEBRAIC[22] is E_K1 in component time_independent_potassium_current (mV).
 * ALGEBRAIC[21] is g_K1 in component time_independent_potassium_current (mS_per_mm2).
 * CONSTANTS[26] is g_K1_max in component time_independent_potassium_current (mS_per_mm2).
 * ALGEBRAIC[25] is K1_infinity in component time_independent_potassium_current_K1_gate (dimensionless).
 * ALGEBRAIC[23] is alpha_K1 in component time_independent_potassium_current_K1_gate (per_ms).
 * ALGEBRAIC[24] is beta_K1 in component time_independent_potassium_current_K1_gate (per_ms).
 * ALGEBRAIC[26] is E_Kp in component plateau_potassium_current (mV).
 * CONSTANTS[27] is g_Kp in component plateau_potassium_current (mS_per_mm2).
 * ALGEBRAIC[27] is Kp in component plateau_potassium_current (dimensionless).
 * CONSTANTS[28] is k_NaCa in component Na_Ca_exchanger (uA_per_mm2).
 * CONSTANTS[29] is K_mNa in component Na_Ca_exchanger (mM).
 * CONSTANTS[30] is K_mCa in component Na_Ca_exchanger (mM).
 * CONSTANTS[31] is k_sat in component Na_Ca_exchanger (dimensionless).
 * CONSTANTS[32] is eta in component Na_Ca_exchanger (dimensionless).
 * STATES[22] is Cai in component calcium_subsystem (mM).
 * CONSTANTS[33] is K_mpCa in component sarcolemmal_calcium_pump (mM).
 * CONSTANTS[34] is I_pCa in component sarcolemmal_calcium_pump (uA_per_mm2).
 * CONSTANTS[35] is g_Nab in component sodium_background_current (mS_per_mm2).
 * ALGEBRAIC[31] is E_NaN in component sodium_background_current (mV).
 * CONSTANTS[36] is g_Cab in component calcium_background_current (mS_per_mm2).
 * ALGEBRAIC[33] is E_CaN in component calcium_background_current (mV).
 * CONSTANTS[37] is I_NaK in component sodium_potassium_pump (uA_per_mm2).
 * ALGEBRAIC[35] is f_NaK in component sodium_potassium_pump (dimensionless).
 * CONSTANTS[38] is K_mNai in component sodium_potassium_pump (mM).
 * CONSTANTS[39] is K_mKo in component sodium_potassium_pump (mM).
 * CONSTANTS[90] is sigma in component sodium_potassium_pump (dimensionless).
 * ALGEBRAIC[59] is i_ns_Na in component non_specific_calcium_activated_current (uA_per_mm2).
 * ALGEBRAIC[60] is i_ns_K in component non_specific_calcium_activated_current (uA_per_mm2).
 * ALGEBRAIC[39] is I_ns_Na in component non_specific_calcium_activated_current (uA_per_mm2).
 * ALGEBRAIC[40] is I_ns_K in component non_specific_calcium_activated_current (uA_per_mm2).
 * CONSTANTS[40] is K_m_ns_Ca in component non_specific_calcium_activated_current (mM).
 * CONSTANTS[41] is P_ns_Ca in component non_specific_calcium_activated_current (mm_per_ms).
 * ALGEBRAIC[37] is EnsCa in component non_specific_calcium_activated_current (mV).
 * ALGEBRAIC[38] is VnsCa in component non_specific_calcium_activated_current (mV).
 * CONSTANTS[42] is Am in component calcium_subsystem (per_mm).
 * CONSTANTS[43] is V_myo in component calcium_subsystem (dimensionless).
 * ALGEBRAIC[41] is RyR_open in component calcium_subsystem (dimensionless).
 * STATES[23] is P_O1 in component calcium_subsystem (dimensionless).
 * STATES[24] is P_O2 in component calcium_subsystem (dimensionless).
 * STATES[25] is P_C1 in component calcium_subsystem (dimensionless).
 * STATES[26] is P_C2 in component calcium_subsystem (dimensionless).
 * CONSTANTS[44] is v1 in component calcium_subsystem (per_ms).
 * CONSTANTS[45] is v2 in component calcium_subsystem (per_ms).
 * CONSTANTS[46] is v3 in component calcium_subsystem (mM_per_ms).
 * CONSTANTS[47] is nCa in component calcium_subsystem (dimensionless).
 * CONSTANTS[48] is mCa in component calcium_subsystem (dimensionless).
 * CONSTANTS[49] is k_a_plus in component calcium_subsystem (per_mM4_per_ms).
 * CONSTANTS[50] is k_a_minus in component calcium_subsystem (per_ms).
 * CONSTANTS[51] is k_b_plus in component calcium_subsystem (per_mM3_per_ms).
 * CONSTANTS[52] is k_b_minus in component calcium_subsystem (per_ms).
 * CONSTANTS[53] is k_c_plus in component calcium_subsystem (per_ms).
 * CONSTANTS[54] is k_c_minus in component calcium_subsystem (per_ms).
 * CONSTANTS[55] is k_htrpn_plus in component calcium_subsystem (per_mM_per_ms).
 * CONSTANTS[56] is k_htrpn_minus in component calcium_subsystem (per_ms).
 * CONSTANTS[57] is k_ltrpn_plus in component calcium_subsystem (per_mM_per_ms).
 * CONSTANTS[58] is k_ltrpn_minus in component calcium_subsystem (per_ms).
 * CONSTANTS[59] is tau_tr in component calcium_subsystem (ms).
 * STATES[27] is Ca_JSR in component calcium_subsystem (mM).
 * STATES[28] is Ca_NSR in component calcium_subsystem (mM).
 * CONSTANTS[96] is V_JSR in component calcium_subsystem (dimensionless).
 * CONSTANTS[94] is V_NSR in component calcium_subsystem (dimensionless).
 * CONSTANTS[91] is V_SS in component calcium_subsystem (dimensionless).
 * CONSTANTS[60] is K_mup in component calcium_subsystem (mM).
 * CONSTANTS[61] is K_mCMDN in component calcium_subsystem (mM).
 * CONSTANTS[62] is K_mCSQN in component calcium_subsystem (mM).
 * CONSTANTS[63] is tau_xfer in component calcium_subsystem (ms).
 * CONSTANTS[64] is HTRPN_tot in component calcium_subsystem (mM).
 * CONSTANTS[65] is LTRPN_tot in component calcium_subsystem (mM).
 * STATES[29] is HTRPNCa in component calcium_subsystem (mM).
 * STATES[30] is LTRPNCa in component calcium_subsystem (mM).
 * CONSTANTS[66] is CSQN_tot in component calcium_subsystem (mM).
 * CONSTANTS[67] is CMDN_tot in component calcium_subsystem (mM).
 * ALGEBRAIC[49] is Bi in component calcium_subsystem (dimensionless).
 * ALGEBRAIC[50] is B_SS in component calcium_subsystem (dimensionless).
 * ALGEBRAIC[51] is B_JSR in component calcium_subsystem (dimensionless).
 * ALGEBRAIC[62] is J_rel in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[42] is J_leak in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[43] is J_up in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[44] is J_tr in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[45] is J_xfer in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[48] is J_trpn in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[46] is J_htrpn in component calcium_subsystem (mM_per_ms).
 * ALGEBRAIC[47] is J_ltrpn in component calcium_subsystem (mM_per_ms).
 * CONSTANTS[68] is Tref in component Myofilaments (N_per_mm2).
 * CONSTANTS[69] is beta0 in component Myofilaments (dimensionless).
 * CONSTANTS[70] is a in component Myofilaments (dimensionless).
 * STATES[31] is Q1 in component Myofilaments (dimensionless).
 * STATES[32] is Q2 in component Myofilaments (dimensionless).
 * STATES[33] is Q3 in component Myofilaments (dimensionless).
 * CONSTANTS[71] is A1 in component Myofilaments (dimensionless).
 * CONSTANTS[72] is A2 in component Myofilaments (dimensionless).
 * CONSTANTS[73] is A3 in component Myofilaments (dimensionless).
 * CONSTANTS[74] is alpha1 in component Myofilaments (dimensionless).
 * CONSTANTS[75] is alpha2 in component Myofilaments (dimensionless).
 * CONSTANTS[76] is alpha3 in component Myofilaments (dimensionless).
 * CONSTANTS[77] is Ca50ref in component Myofilaments (mM).
 * CONSTANTS[78] is zp in component Myofilaments (dimensionless).
 * CONSTANTS[79] is beta1 in component Myofilaments (dimensionless).
 * CONSTANTS[80] is alpha0 in component Myofilaments (per_ms).
 * CONSTANTS[81] is alphar1 in component Myofilaments (per_ms).
 * CONSTANTS[82] is alphar2 in component Myofilaments (per_ms).
 * CONSTANTS[83] is nRel in component Myofilaments (dimensionless).
 * CONSTANTS[84] is Kz in component Myofilaments (dimensionless).
 * CONSTANTS[85] is nHill in component Myofilaments (dimensionless).
 * CONSTANTS[86] is kon in component Myofilaments (per_mM_per_ms).
 * CONSTANTS[87] is koff in component Myofilaments (per_ms).
 * CONSTANTS[88] is gamma_trpn in component Myofilaments (dimensionless).
 * CONSTANTS[89] is TRPN_tot in component Myofilaments (mM).
 * ALGEBRAIC[52] is T0 in component Myofilaments (N_per_mm2).
 * CONSTANTS[92] is T0max in component Myofilaments (N_per_mm2).
 * STATES[34] is z in component Myofilaments (dimensionless).
 * CONSTANTS[99] is z_max in component Myofilaments (dimensionless).
 * ALGEBRAIC[53] is Q in component Myofilaments (dimensionless).
 * STATES[35] is Cab in component Myofilaments (mM).
 * CONSTANTS[97] is Ca50 in component Myofilaments (mM).
 * CONSTANTS[98] is CaTRPN50 in component Myofilaments (mM).
 * CONSTANTS[93] is K_2 in component Myofilaments (dimensionless).
 * CONSTANTS[95] is K_1 in component Myofilaments (dimensionless).
 * ALGEBRAIC[63] is Tension in component Myofilaments (N_per_mm2).
 * RATES[0] is d/dt V in component membrane (mV).
 * RATES[2] is d/dt m in component fast_sodium_current_m_gate (dimensionless).
 * RATES[3] is d/dt h in component fast_sodium_current_h_gate (dimensionless).
 * RATES[4] is d/dt j in component fast_sodium_current_j_gate (dimensionless).
 * RATES[7] is d/dt C0 in component L_type_Ca_channel (dimensionless).
 * RATES[8] is d/dt C1 in component L_type_Ca_channel (dimensionless).
 * RATES[9] is d/dt C2 in component L_type_Ca_channel (dimensionless).
 * RATES[10] is d/dt C3 in component L_type_Ca_channel (dimensionless).
 * RATES[11] is d/dt C4 in component L_type_Ca_channel (dimensionless).
 * RATES[5] is d/dt O in component L_type_Ca_channel (dimensionless).
 * RATES[12] is d/dt C_Ca0 in component L_type_Ca_channel (dimensionless).
 * RATES[13] is d/dt C_Ca1 in component L_type_Ca_channel (dimensionless).
 * RATES[14] is d/dt C_Ca2 in component L_type_Ca_channel (dimensionless).
 * RATES[15] is d/dt C_Ca3 in component L_type_Ca_channel (dimensionless).
 * RATES[16] is d/dt C_Ca4 in component L_type_Ca_channel (dimensionless).
 * RATES[6] is d/dt O_Ca in component L_type_Ca_channel (dimensionless).
 * RATES[20] is d/dt y in component L_type_Ca_channel_y_gate (dimensionless).
 * RATES[21] is d/dt X in component time_dependent_potassium_current_X_gate (dimensionless).
 * RATES[25] is d/dt P_C1 in component calcium_subsystem (dimensionless).
 * RATES[23] is d/dt P_O1 in component calcium_subsystem (dimensionless).
 * RATES[24] is d/dt P_O2 in component calcium_subsystem (dimensionless).
 * RATES[26] is d/dt P_C2 in component calcium_subsystem (dimensionless).
 * RATES[29] is d/dt HTRPNCa in component calcium_subsystem (mM).
 * RATES[30] is d/dt LTRPNCa in component calcium_subsystem (mM).
 * RATES[22] is d/dt Cai in component calcium_subsystem (mM).
 * RATES[17] is d/dt Ca_SS in component calcium_subsystem (mM).
 * RATES[27] is d/dt Ca_JSR in component calcium_subsystem (mM).
 * RATES[28] is d/dt Ca_NSR in component calcium_subsystem (mM).
 * RATES[1] is d/dt Nai in component ionic_concentrations (mM).
 * RATES[19] is d/dt Ki in component ionic_concentrations (mM).
 * RATES[18] is d/dt Ko in component ionic_concentrations (mM).
 * RATES[35] is d/dt Cab in component Myofilaments (mM).
 * RATES[34] is d/dt z in component Myofilaments (dimensionless).
 * RATES[31] is d/dt Q1 in component Myofilaments (dimensionless).
 * RATES[32] is d/dt Q2 in component Myofilaments (dimensionless).
 * RATES[33] is d/dt Q3 in component Myofilaments (dimensionless).
 * There are a total of 8 condition variables.
 */
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1;
CONSTANTS[1] = 0;
STATES[0] = -84.1638;
CONSTANTS[2] = 8.3145e3;
CONSTANTS[3] = 310;
CONSTANTS[4] = 9.6845e4;
CONSTANTS[5] = 0.01;
CONSTANTS[6] = 100;
CONSTANTS[7] = 9000;
CONSTANTS[8] = 750;
CONSTANTS[9] = 1;
CONSTANTS[10] = -100;
CONSTANTS[11] = 0.128;
STATES[1] = 10.2042;
CONSTANTS[12] = 140;
STATES[2] = 0.0328302;
STATES[3] = 0.988354;
STATES[4] = 0.99254;
CONSTANTS[13] = 33.75e-6;
CONSTANTS[14] = 1e-9;
CONSTANTS[15] = -4.58e-3;
STATES[5] = 9.84546e-21;
STATES[6] = 0;
CONSTANTS[16] = 2;
CONSTANTS[17] = 2;
CONSTANTS[18] = 2;
CONSTANTS[19] = 0.3;
CONSTANTS[20] = 0;
CONSTANTS[21] = 0;
CONSTANTS[22] = 0.01;
STATES[7] = 0.997208;
STATES[8] = 6.38897e-5;
STATES[9] = 1.535e-9;
STATES[10] = 1.63909e-14;
STATES[11] = 6.56337e-20;
STATES[12] = 2.72826e-3;
STATES[13] = 6.99215e-7;
STATES[14] = 6.71989e-11;
STATES[15] = 2.87031e-15;
STATES[16] = 4.59752e-20;
STATES[17] = 1.36058e-4;
CONSTANTS[23] = 1.8;
STATES[18] = 5.4;
STATES[19] = 143.727;
STATES[20] = 0.998983;
CONSTANTS[24] = 0.001128;
CONSTANTS[25] = 0.01833;
STATES[21] = 0.000928836;
CONSTANTS[26] = 7.5e-3;
CONSTANTS[27] = 8.28e-5;
CONSTANTS[28] = 50;
CONSTANTS[29] = 87.5;
CONSTANTS[30] = 1.38;
CONSTANTS[31] = 0.1;
CONSTANTS[32] = 0.35;
STATES[22] = 9.94893e-11;
CONSTANTS[33] = 0.5e-3;
CONSTANTS[34] = 1.15e-2;
CONSTANTS[35] = 1.41e-5;
CONSTANTS[36] = 6.032e-5;
CONSTANTS[37] = 0.013;
CONSTANTS[38] = 10;
CONSTANTS[39] = 1.5;
CONSTANTS[40] = 1.2e-3;
CONSTANTS[41] = 1.75e-9;
CONSTANTS[42] = 546.69;
CONSTANTS[43] = 0.92;
STATES[23] = 1.19168e-3;
STATES[24] = 6.30613e-9;
STATES[25] = 0.762527;
STATES[26] = 0.236283;
CONSTANTS[44] = 1.8;
CONSTANTS[45] = 0.58e-4;
CONSTANTS[46] = 1.8e-3;
CONSTANTS[47] = 4;
CONSTANTS[48] = 3;
CONSTANTS[49] = 1.215e10;
CONSTANTS[50] = 0.1425;
CONSTANTS[51] = 4.05e7;
CONSTANTS[52] = 1.93;
CONSTANTS[53] = 0.018;
CONSTANTS[54] = 0.0008;
CONSTANTS[55] = 20;
CONSTANTS[56] = 0.066e-3;
CONSTANTS[57] = 40;
CONSTANTS[58] = 0.04;
CONSTANTS[59] = 34.48;
STATES[27] = 1.17504;
STATES[28] = 1.243891;
CONSTANTS[60] = 0.5e-3;
CONSTANTS[61] = 2.38e-3;
CONSTANTS[62] = 0.8;
CONSTANTS[63] = 3.125;
CONSTANTS[64] = 0.14;
CONSTANTS[65] = 0.07;
STATES[29] = 0.13598;
STATES[30] = 0.00635;
CONSTANTS[66] = 15;
CONSTANTS[67] = 0.05;
CONSTANTS[68] = 56.2;
CONSTANTS[69] = 4.9;
CONSTANTS[70] = 0.35;
STATES[31] = 0;
STATES[32] = 0;
STATES[33] = 0;
CONSTANTS[71] = -29;
CONSTANTS[72] = 138;
CONSTANTS[73] = 129;
CONSTANTS[74] = 0.03;
CONSTANTS[75] = 0.13;
CONSTANTS[76] = 0.625;
CONSTANTS[77] = 1.05e-3;
CONSTANTS[78] = 0.85;
CONSTANTS[79] = -4;
CONSTANTS[80] = 8e-3;
CONSTANTS[81] = 2e-3;
CONSTANTS[82] = 1.75e-3;
CONSTANTS[83] = 3;
CONSTANTS[84] = 0.15;
CONSTANTS[85] = 3;
CONSTANTS[86] = 100;
CONSTANTS[87] = 0.2;
CONSTANTS[88] = 2;
CONSTANTS[89] = 0.07;
STATES[34] = 0;
STATES[35] = 0;
CONSTANTS[90] =  (1.00000/7.00000)*(exp(CONSTANTS[12]/67.3000) - 1.00000);
CONSTANTS[91] =  5.82800e-05*CONSTANTS[43];
CONSTANTS[92] =  CONSTANTS[68]*(1.00000+ CONSTANTS[69]*(CONSTANTS[0] - 1.00000));
CONSTANTS[93] =  1.00000*( (( CONSTANTS[82]*pow(CONSTANTS[78], CONSTANTS[83]))/(pow(CONSTANTS[78], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83])))*(1.00000 - ( CONSTANTS[83]*pow(CONSTANTS[84], CONSTANTS[83]))/(pow(CONSTANTS[78], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83]))));
CONSTANTS[94] =  0.0810000*CONSTANTS[43];
CONSTANTS[95] =  1.00000*(( CONSTANTS[82]*pow(CONSTANTS[78], CONSTANTS[83] - 1.00000)*CONSTANTS[83]*pow(CONSTANTS[84], CONSTANTS[83]))/pow(pow(CONSTANTS[78], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83]), 2.00000));
CONSTANTS[96] =  0.00464000*CONSTANTS[43];
CONSTANTS[97] =  CONSTANTS[77]*(1.00000+ CONSTANTS[79]*(CONSTANTS[0] - 1.00000));
CONSTANTS[98] = ( CONSTANTS[97]*CONSTANTS[89])/(CONSTANTS[97]+ (CONSTANTS[87]/CONSTANTS[86])*(1.00000 - ( (1.00000+ CONSTANTS[69]*(CONSTANTS[0] - 1.00000))*0.500000)/CONSTANTS[88]));
CONSTANTS[99] = ( CONSTANTS[80]*1.00000)/pow(CONSTANTS[98]/CONSTANTS[89], CONSTANTS[85]) - ( CONSTANTS[93]*1.00000)/(CONSTANTS[81]+ CONSTANTS[95]*1.00000+CONSTANTS[80]/pow(CONSTANTS[98]/CONSTANTS[89], CONSTANTS[85]));
RATES[0] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[4] = 0.1001;
RATES[7] = 0.1001;
RATES[8] = 0.1001;
RATES[9] = 0.1001;
RATES[10] = 0.1001;
RATES[11] = 0.1001;
RATES[5] = 0.1001;
RATES[12] = 0.1001;
RATES[13] = 0.1001;
RATES[14] = 0.1001;
RATES[15] = 0.1001;
RATES[16] = 0.1001;
RATES[6] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 0.1001;
RATES[25] = 0.1001;
RATES[23] = 0.1001;
RATES[24] = 0.1001;
RATES[26] = 0.1001;
RATES[29] = 0.1001;
RATES[30] = 0.1001;
RATES[22] = 0.1001;
RATES[17] = 0.1001;
RATES[27] = 0.1001;
RATES[28] = 0.1001;
RATES[1] = 0.1001;
RATES[19] = 0.1001;
RATES[18] = 0.1001;
RATES[35] = 0.1001;
RATES[34] = 0.1001;
RATES[31] = 0.1001;
RATES[32] = 0.1001;
RATES[33] = 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[54]+ALGEBRAIC[55]+ALGEBRAIC[64]+ALGEBRAIC[57]+ALGEBRAIC[29]+ALGEBRAIC[58]+ALGEBRAIC[28]+ALGEBRAIC[30]+ALGEBRAIC[32]+ALGEBRAIC[34]+ALGEBRAIC[36]+ALGEBRAIC[61]+ALGEBRAIC[0])/CONSTANTS[5];
resid[1] = RATES[2] -  ALGEBRAIC[2]*(1.00000 - STATES[2]) -  ALGEBRAIC[3]*STATES[2];
resid[2] = RATES[3] -  ALGEBRAIC[4]*(1.00000 - STATES[3]) -  ALGEBRAIC[5]*STATES[3];
resid[3] = RATES[4] -  ALGEBRAIC[6]*(1.00000 - STATES[4]) -  ALGEBRAIC[7]*STATES[4];
resid[4] = RATES[7] - ( ALGEBRAIC[10]*STATES[8]+ CONSTANTS[22]*STATES[12]) -  ( 4.00000*ALGEBRAIC[9]+ALGEBRAIC[13])*STATES[7];
resid[5] = RATES[8] - ( 4.00000*ALGEBRAIC[9]*STATES[7]+ 2.00000*ALGEBRAIC[10]*STATES[9]+ (CONSTANTS[22]/CONSTANTS[17])*STATES[13]) -  (ALGEBRAIC[10]+ 3.00000*ALGEBRAIC[9]+ ALGEBRAIC[13]*CONSTANTS[16])*STATES[8];
resid[6] = RATES[9] - ( 3.00000*ALGEBRAIC[9]*STATES[8]+ 3.00000*ALGEBRAIC[10]*STATES[10]+ (CONSTANTS[22]/pow(CONSTANTS[17], 2.00000))*STATES[14]) -  ( ALGEBRAIC[10]*2.00000+ 2.00000*ALGEBRAIC[9]+ ALGEBRAIC[13]*pow(CONSTANTS[16], 2.00000))*STATES[9];
resid[7] = RATES[10] - ( 2.00000*ALGEBRAIC[9]*STATES[9]+ 4.00000*ALGEBRAIC[10]*STATES[11]+ (CONSTANTS[22]/pow(CONSTANTS[17], 3.00000))*STATES[15]) -  ( ALGEBRAIC[10]*3.00000+ALGEBRAIC[9]+ ALGEBRAIC[13]*pow(CONSTANTS[16], 3.00000))*STATES[10];
resid[8] = RATES[11] - ( ALGEBRAIC[9]*STATES[10]+ CONSTANTS[18]*STATES[5]+ (CONSTANTS[22]/pow(CONSTANTS[17], 4.00000))*STATES[16]) -  ( ALGEBRAIC[10]*4.00000+CONSTANTS[19]+ ALGEBRAIC[13]*pow(CONSTANTS[16], 4.00000))*STATES[11];
resid[9] = RATES[5] -  CONSTANTS[19]*STATES[11] -  CONSTANTS[18]*STATES[5];
resid[10] = RATES[12] - ( ALGEBRAIC[12]*STATES[13]+ ALGEBRAIC[13]*STATES[12]) -  ( 4.00000*ALGEBRAIC[11]+CONSTANTS[22])*STATES[12];
resid[11] = RATES[13] - ( 4.00000*ALGEBRAIC[11]*STATES[12]+ 2.00000*ALGEBRAIC[12]*STATES[14]+ ALGEBRAIC[13]*CONSTANTS[16]*STATES[8]) -  (ALGEBRAIC[12]+ 3.00000*ALGEBRAIC[11]+CONSTANTS[22]/CONSTANTS[17])*STATES[13];
resid[12] = RATES[14] - ( 3.00000*ALGEBRAIC[11]*STATES[13]+ 3.00000*ALGEBRAIC[12]*STATES[15]+ ALGEBRAIC[13]*pow(CONSTANTS[16], 2.00000)*STATES[9]) -  ( ALGEBRAIC[12]*2.00000+ 2.00000*ALGEBRAIC[11]+CONSTANTS[22]/pow(CONSTANTS[17], 2.00000))*STATES[14];
resid[13] = RATES[15] - ( 2.00000*ALGEBRAIC[11]*STATES[14]+ 4.00000*ALGEBRAIC[12]*STATES[16]+ ALGEBRAIC[13]*pow(CONSTANTS[16], 3.00000)*STATES[10]) -  ( ALGEBRAIC[12]*3.00000+ALGEBRAIC[11]+CONSTANTS[22]/pow(CONSTANTS[17], 3.00000))*STATES[15];
resid[14] = RATES[16] - ( ALGEBRAIC[11]*STATES[15]+ CONSTANTS[20]*STATES[6]+ ALGEBRAIC[13]*pow(CONSTANTS[16], 4.00000)*STATES[11]) -  ( ALGEBRAIC[12]*4.00000+CONSTANTS[21]+CONSTANTS[22]/pow(CONSTANTS[17], 4.00000))*STATES[16];
resid[15] = RATES[6] -  CONSTANTS[21]*STATES[16] -  CONSTANTS[20]*STATES[6];
resid[16] = RATES[20] - (ALGEBRAIC[14] - STATES[20])/ALGEBRAIC[15];
resid[17] = RATES[21] -  ALGEBRAIC[18]*(1.00000 - STATES[21]) -  ALGEBRAIC[19]*STATES[21];
resid[18] = RATES[25] -  - CONSTANTS[49]*pow(STATES[17], CONSTANTS[47])*STATES[25]+ CONSTANTS[50]*STATES[23];
resid[19] = RATES[23] - ( CONSTANTS[49]*pow(STATES[17], CONSTANTS[47])*STATES[25] - ( CONSTANTS[50]*STATES[23]+ CONSTANTS[51]*pow(STATES[17], CONSTANTS[48])*STATES[23]+ CONSTANTS[53]*STATES[23]))+ CONSTANTS[52]*STATES[24]+ CONSTANTS[54]*STATES[26];
resid[20] = RATES[24] -  CONSTANTS[51]*pow(STATES[17], CONSTANTS[48])*STATES[23] -  CONSTANTS[52]*STATES[24];
resid[21] = RATES[26] -  CONSTANTS[53]*STATES[23] -  CONSTANTS[54]*STATES[26];
resid[22] = RATES[29] - ALGEBRAIC[46];
resid[23] = RATES[30] - ALGEBRAIC[47];
resid[24] = RATES[22] -  ALGEBRAIC[49]*((ALGEBRAIC[42]+ALGEBRAIC[45]) - (ALGEBRAIC[43]+ALGEBRAIC[48]+( ((ALGEBRAIC[34] -  2.00000*ALGEBRAIC[29])+ALGEBRAIC[30])*CONSTANTS[42])/( 2.00000*CONSTANTS[43]*CONSTANTS[4])));
resid[25] = RATES[17] -  ALGEBRAIC[50]*((( ALGEBRAIC[62]*CONSTANTS[96])/CONSTANTS[91] - ( ALGEBRAIC[45]*CONSTANTS[43])/CONSTANTS[91]) - ( ALGEBRAIC[55]*CONSTANTS[42])/( 2.00000*CONSTANTS[91]*CONSTANTS[4]));
resid[26] = RATES[27] -  ALGEBRAIC[51]*(ALGEBRAIC[44] - ALGEBRAIC[62]);
resid[27] = RATES[28] - ( (ALGEBRAIC[43] - ALGEBRAIC[42])*CONSTANTS[43])/CONSTANTS[94] - ( ALGEBRAIC[44]*CONSTANTS[96])/CONSTANTS[94];
resid[28] = RATES[1] - ( - (ALGEBRAIC[54]+ALGEBRAIC[32]+ALGEBRAIC[59]+ ALGEBRAIC[29]*3.00000+ ALGEBRAIC[36]*3.00000)*CONSTANTS[42])/( CONSTANTS[43]*CONSTANTS[4]);
resid[29] = RATES[19] - ( - (ALGEBRAIC[64]+ALGEBRAIC[57]+ALGEBRAIC[58]+ALGEBRAIC[28]+ALGEBRAIC[60]+ - ALGEBRAIC[36]*2.00000)*CONSTANTS[42])/( CONSTANTS[43]*CONSTANTS[4]);
resid[30] = RATES[18] - ( (ALGEBRAIC[64]+ALGEBRAIC[57]+ALGEBRAIC[58]+ALGEBRAIC[28]+ALGEBRAIC[60]+ - ALGEBRAIC[36]*2.00000)*CONSTANTS[42])/( CONSTANTS[43]*CONSTANTS[4]);
resid[31] = RATES[35] -  CONSTANTS[86]*STATES[22]*(CONSTANTS[89] - STATES[35]) -  CONSTANTS[87]*(1.00000 - ALGEBRAIC[63]/( CONSTANTS[88]*CONSTANTS[68]))*STATES[35];
resid[32] = RATES[34] -  CONSTANTS[80]*pow(STATES[35]/CONSTANTS[98], CONSTANTS[85])*(1.00000 - STATES[34])+ - STATES[34]*CONSTANTS[81]+( - CONSTANTS[82]*pow(STATES[34], CONSTANTS[83]))/(pow(STATES[34], CONSTANTS[83])+pow(CONSTANTS[84], CONSTANTS[83]));
resid[33] = RATES[31] -  CONSTANTS[71]*CONSTANTS[1] -  1.00000*( CONSTANTS[74]*STATES[31]);
resid[34] = RATES[32] -  CONSTANTS[72]*CONSTANTS[1] -  1.00000*( CONSTANTS[75]*STATES[32]);
resid[35] = RATES[33] -  CONSTANTS[73]*CONSTANTS[1] -  1.00000*( CONSTANTS[76]*STATES[33]);
}
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] = (CONDVAR[0]>=0.00000&&CONDVAR[1]<=0.00000&&CONDVAR[2]<=0.00000 ? CONSTANTS[10] : 0.00000);
ALGEBRAIC[2] = ( 0.320000*(STATES[0]+47.1300))/(1.00000 - exp( - 0.100000*(STATES[0]+47.1300)));
ALGEBRAIC[3] =  0.0800000*exp(- STATES[0]/11.0000);
ALGEBRAIC[4] = (CONDVAR[3]<0.00000 ?  0.135000*exp((80.0000+STATES[0])/- 6.80000) : 0.00000);
ALGEBRAIC[5] = (CONDVAR[4]<0.00000 ?  3.56000*exp( 0.0790000*STATES[0])+ 310000.*exp( 0.350000*STATES[0]) : 1.00000/( 0.130000*(1.00000+exp((STATES[0]+10.6600)/- 11.1000))));
ALGEBRAIC[6] = (CONDVAR[5]<0.00000 ? ( ( - 127140.*exp( 0.244400*STATES[0]) -  3.47400e-05*exp( - 0.0439100*STATES[0]))*(STATES[0]+37.7800))/(1.00000+exp( 0.311000*(STATES[0]+79.2300))) : 0.00000);
ALGEBRAIC[7] = (CONDVAR[6]<0.00000 ? ( 0.121200*exp( - 0.0105200*STATES[0]))/(1.00000+exp( - 0.137800*(STATES[0]+40.1400))) : ( 0.300000*exp( - 2.53500e-07*STATES[0]))/(1.00000+exp( - 0.100000*(STATES[0]+32.0000))));
ALGEBRAIC[9] =  0.400000*exp((STATES[0]+12.0000)/10.0000);
ALGEBRAIC[10] =  0.0500000*exp((STATES[0]+12.0000)/- 13.0000);
ALGEBRAIC[11] =  ALGEBRAIC[9]*CONSTANTS[16];
ALGEBRAIC[12] = ALGEBRAIC[10]/CONSTANTS[17];
ALGEBRAIC[13] = ( 0.187500*STATES[17])/1.00000;
ALGEBRAIC[14] = 1.00000/(1.00000+exp((STATES[0]+55.0000)/7.50000))+0.100000/(1.00000+exp((- STATES[0]+21.0000)/6.00000));
ALGEBRAIC[15] = 20.0000+600.000/(1.00000+exp((STATES[0]+30.0000)/9.50000));
ALGEBRAIC[18] = ( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)));
ALGEBRAIC[19] = ( 0.000131000*(STATES[0]+30.0000))/(- 1.00000+exp( 0.0687000*(STATES[0]+30.0000)));
ALGEBRAIC[22] =  (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log(STATES[18]/STATES[19]);
ALGEBRAIC[26] = ALGEBRAIC[22];
ALGEBRAIC[27] = 1.00000/(1.00000+exp((7.48800 - STATES[0])/5.98000));
ALGEBRAIC[28] =  CONSTANTS[27]*ALGEBRAIC[27]*(STATES[0] - ALGEBRAIC[26]);
ALGEBRAIC[29] =  (( (( (( CONSTANTS[28]*1.00000)/(pow(CONSTANTS[29], 3.00000)+pow(CONSTANTS[12], 3.00000)))*1.00000)/(CONSTANTS[30]+CONSTANTS[23]))*1.00000)/(1.00000+ CONSTANTS[31]*exp(( (CONSTANTS[32] - 1.00000)*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))))*( exp(( CONSTANTS[32]*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))*pow(STATES[1], 3.00000)*CONSTANTS[23] -  exp(( (CONSTANTS[32] - 1.00000)*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))*pow(CONSTANTS[12], 3.00000)*STATES[22]);
ALGEBRAIC[30] = ( CONSTANTS[34]*STATES[22])/(CONSTANTS[33]+STATES[22]);
ALGEBRAIC[1] =  (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log(CONSTANTS[12]/STATES[1]);
ALGEBRAIC[31] = ALGEBRAIC[1];
ALGEBRAIC[32] =  CONSTANTS[35]*(STATES[0] - ALGEBRAIC[31]);
ALGEBRAIC[33] =  (( CONSTANTS[2]*CONSTANTS[3])/( 2.00000*CONSTANTS[4]))*log(CONSTANTS[23]/STATES[22]);
ALGEBRAIC[34] =  CONSTANTS[36]*(STATES[0] - ALGEBRAIC[33]);
ALGEBRAIC[35] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3]))+ 0.0365000*CONSTANTS[90]*exp(( - STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])));
ALGEBRAIC[36] = ( (( CONSTANTS[37]*ALGEBRAIC[35]*1.00000)/(1.00000+pow(CONSTANTS[38]/STATES[1], 1.50000)))*STATES[18])/(STATES[18]+CONSTANTS[39]);
ALGEBRAIC[42] =  CONSTANTS[45]*(STATES[28] - STATES[22]);
ALGEBRAIC[43] = ( CONSTANTS[46]*pow(STATES[22], 2.00000))/(pow(CONSTANTS[60], 2.00000)+pow(STATES[22], 2.00000));
ALGEBRAIC[44] = (STATES[28] - STATES[27])/CONSTANTS[59];
ALGEBRAIC[45] = (STATES[17] - STATES[22])/CONSTANTS[63];
ALGEBRAIC[46] =  CONSTANTS[55]*STATES[22]*(CONSTANTS[64] - STATES[29]) -  CONSTANTS[56]*STATES[29];
ALGEBRAIC[47] =  CONSTANTS[57]*STATES[22]*(CONSTANTS[65] - STATES[30]) -  CONSTANTS[58]*STATES[30];
ALGEBRAIC[48] = ALGEBRAIC[46]+ALGEBRAIC[47];
ALGEBRAIC[49] = 1.00000/(1.00000+( CONSTANTS[67]*CONSTANTS[61])/pow(CONSTANTS[61]+STATES[22], 2.00000));
ALGEBRAIC[50] = 1.00000/(1.00000+( CONSTANTS[67]*CONSTANTS[61])/pow(CONSTANTS[61]+STATES[17], 2.00000));
ALGEBRAIC[51] = 1.00000/(1.00000+( CONSTANTS[66]*CONSTANTS[62])/pow(CONSTANTS[62]+STATES[27], 2.00000));
ALGEBRAIC[54] =  CONSTANTS[11]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*(STATES[0] - ALGEBRAIC[1]);
ALGEBRAIC[8] =  1.00000*(( (( CONSTANTS[13]*4.00000*STATES[0]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) -  0.341000*CONSTANTS[23]))/(exp(( 2.00000*STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000));
ALGEBRAIC[55] =  ALGEBRAIC[8]*STATES[20]*(STATES[5]+STATES[6]);
ALGEBRAIC[16] =  CONSTANTS[24]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[17] =  (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log((STATES[18]+ CONSTANTS[25]*CONSTANTS[12])/(STATES[19]+ CONSTANTS[25]*STATES[1]));
ALGEBRAIC[20] = 1.00000/(1.00000+exp((STATES[0] - 56.2600)/32.1000));
ALGEBRAIC[57] =  ALGEBRAIC[16]*ALGEBRAIC[20]*pow(STATES[21], 2.00000)*(STATES[0] - ALGEBRAIC[17]);
ALGEBRAIC[21] =  CONSTANTS[26]* pow((STATES[18]/5.40000), 1.0 / 2);
ALGEBRAIC[23] = 1.02000/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[22]) - 59.2150)));
ALGEBRAIC[24] = ( 0.491240*(exp( 0.0803200*((STATES[0]+5.47600) - ALGEBRAIC[22]))+exp( 0.0617500*(STATES[0] - (ALGEBRAIC[22]+594.310)))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[22])+4.75300)));
ALGEBRAIC[25] = ALGEBRAIC[23]/(ALGEBRAIC[23]+ALGEBRAIC[24]);
ALGEBRAIC[58] =  ALGEBRAIC[21]*ALGEBRAIC[25]*(STATES[0] - ALGEBRAIC[22]);
ALGEBRAIC[37] =  (( CONSTANTS[2]*CONSTANTS[3])/CONSTANTS[4])*log((STATES[18]+CONSTANTS[12])/(STATES[19]+STATES[1]));
ALGEBRAIC[38] = STATES[0] - ALGEBRAIC[37];
ALGEBRAIC[39] = ( (( 1.00000*CONSTANTS[41]*pow(1.00000, 2.00000)*ALGEBRAIC[38]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.750000*STATES[1]*exp(( ALGEBRAIC[38]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) -  0.750000*CONSTANTS[12]))/(exp(( ALGEBRAIC[38]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[59] = ( ALGEBRAIC[39]*1.00000)/(1.00000+pow(CONSTANTS[40]/STATES[22], 3.00000));
ALGEBRAIC[40] = ( (( 1.00000*CONSTANTS[41]*pow(1.00000, 2.00000)*ALGEBRAIC[38]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( 0.750000*STATES[19]*exp(( ALGEBRAIC[38]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) -  0.750000*STATES[18]))/(exp(( ALGEBRAIC[38]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
ALGEBRAIC[60] = ( ALGEBRAIC[40]*1.00000)/(1.00000+pow(CONSTANTS[40]/STATES[22], 3.00000));
ALGEBRAIC[61] = ALGEBRAIC[59]+ALGEBRAIC[60];
ALGEBRAIC[41] = STATES[23]+STATES[24];
ALGEBRAIC[62] =  CONSTANTS[44]*ALGEBRAIC[41]*(STATES[27] - STATES[17]);
ALGEBRAIC[52] = ( CONSTANTS[92]*STATES[34])/CONSTANTS[99];
ALGEBRAIC[53] = STATES[31]+STATES[32]+STATES[33];
ALGEBRAIC[63] = (CONDVAR[7]<0.00000 ? ( ALGEBRAIC[52]*( CONSTANTS[70]*ALGEBRAIC[53]+1.00000))/(1.00000 - ALGEBRAIC[53]) : ( ALGEBRAIC[52]*(1.00000 -  (CONSTANTS[70]+2.00000)*ALGEBRAIC[53]))/(1.00000+ALGEBRAIC[53]));
ALGEBRAIC[56] = CONSTANTS[14]/(1.00000+ALGEBRAIC[8]/CONSTANTS[15]);
ALGEBRAIC[64] = ( (( 1.00000*ALGEBRAIC[56]*STATES[20]*(STATES[5]+STATES[6])*STATES[0]*pow(CONSTANTS[4], 2.00000))/( CONSTANTS[2]*CONSTANTS[3]))*( STATES[19]*exp(( STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - STATES[18]))/(exp(( STATES[0]*CONSTANTS[4])/( CONSTANTS[2]*CONSTANTS[3])) - 1.00000);
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[16] = 1.0;
SI[17] = 1.0;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[21] = 1.0;
SI[22] = 1.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[25] = 1.0;
SI[26] = 1.0;
SI[27] = 1.0;
SI[28] = 1.0;
SI[29] = 1.0;
SI[30] = 1.0;
SI[31] = 1.0;
SI[32] = 1.0;
SI[33] = 1.0;
SI[34] = 1.0;
SI[35] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
             double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[6];
CONDVAR[1] = VOI - CONSTANTS[7];
CONDVAR[2] = ((VOI - CONSTANTS[6]) -  floor((VOI - CONSTANTS[6])/CONSTANTS[8])*CONSTANTS[8]) - CONSTANTS[9];
CONDVAR[3] = STATES[0] - - 40.0000;
CONDVAR[4] = STATES[0] - - 40.0000;
CONDVAR[5] = STATES[0] - - 40.0000;
CONDVAR[6] = STATES[0] - - 40.0000;
CONDVAR[7] = ALGEBRAIC[53] - 0.00000;
}