/* There are a total of 73 entries in the algebraic variable array. There are a total of 29 entries in each of the rate and state variable arrays. There are a total of 55 entries in the constant variable array. */ /* * VOI is time in component environment (second). * STATES[0] is V in component membrane (millivolt). * CONSTANTS[0] is R in component membrane (joule_per_kilomole_kelvin). * CONSTANTS[1] is T in component membrane (kelvin). * CONSTANTS[2] is F in component membrane (coulomb_per_mole). * CONSTANTS[3] is Cm in component membrane (microF). * ALGEBRAIC[64] is i_Na in component sodium_current (nanoA). * ALGEBRAIC[20] is i_Ca_T in component T_type_Ca_channel (nanoA). * ALGEBRAIC[65] is i_Ca_L in component L_type_Ca_channel (nanoA). * ALGEBRAIC[66] is i_K in component delayed_rectifying_potassium_current (nanoA). * ALGEBRAIC[37] is i_f in component hyperpolarisation_activated_current (nanoA). * ALGEBRAIC[70] is i_B in component linear_background_current (nanoA). * ALGEBRAIC[42] is i_NaK in component sodium_potassium_pump (nanoA). * ALGEBRAIC[44] is i_NaCa in component sodium_calcium_pump (nanoA). * ALGEBRAIC[43] is i_Ca_P in component calcium_pump_current (nanoA). * ALGEBRAIC[72] is i_K_ACh in component muscarinic_potassium_current (nanoA). * CONSTANTS[4] is P_Na in component sodium_current (mul_per_second). * ALGEBRAIC[61] is E_Na in component reversal_potentials (millivolt). * STATES[1] is Na_c in component cleft_space_equations (millimolar). * CONSTANTS[45] is F_ACh_Na in component sodium_current (dimensionless). * CONSTANTS[5] is ACh in component cAMP_balance (millimolar). * STATES[2] is m in component sodium_current_m_gate (dimensionless). * STATES[3] is h1 in component sodium_current_h_gate (dimensionless). * STATES[4] is h2 in component sodium_current_h_gate (dimensionless). * ALGEBRAIC[2] is m_infinity in component sodium_current_m_gate (dimensionless). * ALGEBRAIC[3] is tau_m in component sodium_current_m_gate (second). * ALGEBRAIC[0] is alpha_m in component sodium_current_m_gate (per_second). * ALGEBRAIC[1] is beta_m in component sodium_current_m_gate (per_second). * ALGEBRAIC[6] is h1_infinity in component sodium_current_h_gate (dimensionless). * ALGEBRAIC[8] is h2_infinity in component sodium_current_h_gate (dimensionless). * ALGEBRAIC[7] is tau_h1 in component sodium_current_h_gate (second). * ALGEBRAIC[9] is tau_h2 in component sodium_current_h_gate (second). * ALGEBRAIC[4] is alpha_h1 in component sodium_current_h_gate (per_second). * ALGEBRAIC[5] is beta_h1 in component sodium_current_h_gate (per_second). * ALGEBRAIC[11] is g_Ca_L in component L_type_Ca_channel (microS). * CONSTANTS[6] is g_Ca_L_cont in component L_type_Ca_channel (microS). * ALGEBRAIC[10] is F_cAMP_CaL in component L_type_Ca_channel (dimensionless). * STATES[5] is cAMP in component cAMP_balance (millimolar). * CONSTANTS[7] is E_Ca_L in component L_type_Ca_channel (millivolt). * STATES[6] is d_L in component L_type_Ca_channel_d_gate (dimensionless). * ALGEBRAIC[15] is d_L_infinity in component L_type_Ca_channel_d_gate (dimensionless). * STATES[7] is f_L in component L_type_Ca_channel_f_gate (dimensionless). * ALGEBRAIC[12] is alpha_d_L in component L_type_Ca_channel_d_gate (per_second). * ALGEBRAIC[13] is beta_d_L in component L_type_Ca_channel_d_gate (per_second). * ALGEBRAIC[14] is tau_d_L in component L_type_Ca_channel_d_gate (second). * ALGEBRAIC[16] is alpha_f_L in component L_type_Ca_channel_f_gate (per_second). * ALGEBRAIC[17] is beta_f_L in component L_type_Ca_channel_f_gate (per_second). * ALGEBRAIC[19] is f_L_infinity in component L_type_Ca_channel_f_gate (dimensionless). * ALGEBRAIC[18] is tau_f_L in component L_type_Ca_channel_f_gate (second). * CONSTANTS[8] is g_Ca_T in component T_type_Ca_channel (microS). * CONSTANTS[9] is E_Ca_T in component T_type_Ca_channel (millivolt). * STATES[8] is d_T in component T_type_Ca_channel_d_gate (dimensionless). * STATES[9] is f_T in component T_type_Ca_channel_f_gate (dimensionless). * ALGEBRAIC[21] is alpha_d_T in component T_type_Ca_channel_d_gate (per_second). * ALGEBRAIC[22] is beta_d_T in component T_type_Ca_channel_d_gate (per_second). * ALGEBRAIC[24] is d_T_infinity in component T_type_Ca_channel_d_gate (dimensionless). * ALGEBRAIC[23] is tau_d_T in component T_type_Ca_channel_d_gate (second). * ALGEBRAIC[25] is alpha_f_T in component T_type_Ca_channel_f_gate (per_second). * ALGEBRAIC[26] is beta_f_T in component T_type_Ca_channel_f_gate (per_second). * ALGEBRAIC[28] is f_T_infinity in component T_type_Ca_channel_f_gate (dimensionless). * ALGEBRAIC[27] is tau_f_T in component T_type_Ca_channel_f_gate (second). * ALGEBRAIC[30] is g_K in component delayed_rectifying_potassium_current (microS). * ALGEBRAIC[29] is F_cAMP_K in component delayed_rectifying_potassium_current (dimensionless). * ALGEBRAIC[62] is E_K in component reversal_potentials (millivolt). * CONSTANTS[10] is K_b in component cleft_space_equations (millimolar). * STATES[10] is P_a in component delayed_rectifying_potassium_current_P_a_gate (dimensionless). * STATES[11] is P_i in component delayed_rectifying_potassium_current_P_i_gate (dimensionless). * ALGEBRAIC[32] is tau_P_a in component delayed_rectifying_potassium_current_P_a_gate (second). * ALGEBRAIC[31] is P_a_infinity in component delayed_rectifying_potassium_current_P_a_gate (dimensionless). * ALGEBRAIC[33] is alpha_P_i in component delayed_rectifying_potassium_current_P_i_gate (per_second). * ALGEBRAIC[34] is beta_P_i in component delayed_rectifying_potassium_current_P_i_gate (per_second). * ALGEBRAIC[67] is i_B_Na in component linear_background_current (nanoA). * ALGEBRAIC[68] is i_B_Ca in component linear_background_current (nanoA). * ALGEBRAIC[69] is i_B_K in component linear_background_current (nanoA). * CONSTANTS[11] is g_B_Na in component linear_background_current (microS). * CONSTANTS[12] is g_B_Ca in component linear_background_current (microS). * CONSTANTS[13] is g_B_K in component linear_background_current (microS). * ALGEBRAIC[63] is E_Ca in component reversal_potentials (millivolt). * CONSTANTS[46] is F_ACh_bNa in component linear_background_current (dimensionless). * ALGEBRAIC[35] is i_f_Na in component hyperpolarisation_activated_current (nanoA). * ALGEBRAIC[36] is i_f_K in component hyperpolarisation_activated_current (nanoA). * CONSTANTS[14] is g_f_Na in component hyperpolarisation_activated_current (microS). * CONSTANTS[15] is g_f_K in component hyperpolarisation_activated_current (microS). * STATES[12] is y in component hyperpolarisation_activated_current_y_gate (dimensionless). * ALGEBRAIC[39] is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless). * ALGEBRAIC[38] is V_half in component hyperpolarisation_activated_current_y_gate (millivolt). * ALGEBRAIC[40] is tau_y in component hyperpolarisation_activated_current_y_gate (second). * CONSTANTS[16] is K_m_Na in component sodium_potassium_pump (millimolar). * CONSTANTS[17] is K_m_K in component sodium_potassium_pump (millimolar). * CONSTANTS[18] is i_NaK_max in component sodium_potassium_pump (nanoA). * ALGEBRAIC[41] is F_cAMP_NaK in component sodium_potassium_pump (dimensionless). * STATES[13] is Na_i in component intracellular_concentrations_and_buffer_equations (millimolar). * STATES[14] is K_c in component cleft_space_equations (millimolar). * STATES[15] is Ca_i in component intracellular_concentrations_and_buffer_equations (millimolar). * CONSTANTS[19] is i_Ca_P_max in component calcium_pump_current (nanoA). * CONSTANTS[20] is K_NaCa in component sodium_calcium_pump (nanoA). * CONSTANTS[21] is d_NaCa in component sodium_calcium_pump (dimensionless). * CONSTANTS[22] is gamma in component sodium_calcium_pump (dimensionless). * STATES[16] is Ca_c in component cleft_space_equations (millimolar). * ALGEBRAIC[71] is I_K_ACh in component muscarinic_potassium_current (nanoA). * CONSTANTS[52] is g_K_ACh in component muscarinic_potassium_current (microS). * CONSTANTS[23] is g_K_ACh_base in component muscarinic_potassium_current (microS). * CONSTANTS[47] is P_M2_KACh in component muscarinic_potassium_current (dimensionless). * STATES[17] is a in component muscarinic_potassium_current (dimensionless). * ALGEBRAIC[45] is alpha_a in component muscarinic_potassium_current (per_second). * CONSTANTS[48] is beta_a in component muscarinic_potassium_current (per_second). * CONSTANTS[24] is f_Vagal in component cAMP_balance (per_second). * STATES[18] is K_i in component intracellular_concentrations_and_buffer_equations (millimolar). * STATES[19] is Ca_Calmod in component intracellular_concentrations_and_buffer_equations (dimensionless). * STATES[20] is Ca_Trop in component intracellular_concentrations_and_buffer_equations (dimensionless). * STATES[21] is Ca_Mg_Trop in component intracellular_concentrations_and_buffer_equations (dimensionless). * STATES[22] is Mg_Mg_Trop in component intracellular_concentrations_and_buffer_equations (dimensionless). * ALGEBRAIC[46] is phi_C in component intracellular_concentrations_and_buffer_equations (per_second). * ALGEBRAIC[47] is phi_TC in component intracellular_concentrations_and_buffer_equations (per_second). * ALGEBRAIC[48] is phi_TMgC in component intracellular_concentrations_and_buffer_equations (per_second). * ALGEBRAIC[49] is phi_TMgM in component intracellular_concentrations_and_buffer_equations (per_second). * ALGEBRAIC[53] is phi_B in component intracellular_concentrations_and_buffer_equations (millimolar_per_second). * CONSTANTS[25] is Mg_i in component intracellular_concentrations_and_buffer_equations (millimolar). * ALGEBRAIC[50] is F_C in component intracellular_concentrations_and_buffer_equations (millimolar_per_second). * ALGEBRAIC[51] is F_TC in component intracellular_concentrations_and_buffer_equations (millimolar_per_second). * ALGEBRAIC[52] is F_TMgC in component intracellular_concentrations_and_buffer_equations (millimolar_per_second). * CONSTANTS[26] is Vol in component cleft_space_equations (mm_cubed). * CONSTANTS[49] is V_i in component intracellular_concentrations_and_buffer_equations (mm_cubed). * ALGEBRAIC[60] is i_up in component SR_Ca_uptake_and_release (nanoA). * ALGEBRAIC[58] is i_rel in component SR_Ca_uptake_and_release (nanoA). * CONSTANTS[27] is Na_b in component cleft_space_equations (millimolar). * CONSTANTS[28] is Ca_b in component cleft_space_equations (millimolar). * CONSTANTS[50] is V_c in component cleft_space_equations (mm_cubed). * CONSTANTS[29] is tau_p in component cleft_space_equations (second). * STATES[23] is Ca_up in component SR_Ca_uptake_and_release (millimolar). * CONSTANTS[30] is alpha_up in component SR_Ca_uptake_and_release (nanoA). * CONSTANTS[31] is beta_up in component SR_Ca_uptake_and_release (nanoA). * STATES[24] is Ca_rel in component SR_Ca_uptake_and_release (millimolar). * CONSTANTS[32] is alpha_rel in component SR_Ca_uptake_and_release (nanoA_per_millimolar). * ALGEBRAIC[59] is i_tr in component SR_Ca_uptake_and_release (nanoA). * CONSTANTS[51] is K1 in component SR_Ca_uptake_and_release (dimensionless). * ALGEBRAIC[54] is K2 in component SR_Ca_uptake_and_release (millimolar). * CONSTANTS[33] is k_cyca in component SR_Ca_uptake_and_release (millimolar). * CONSTANTS[34] is k_xcs in component SR_Ca_uptake_and_release (dimensionless). * CONSTANTS[35] is k_SRCa in component SR_Ca_uptake_and_release (millimolar). * CONSTANTS[36] is k_rel in component SR_Ca_uptake_and_release (millimolar). * ALGEBRAIC[55] is r_act in component SR_Ca_uptake_and_release (per_second). * ALGEBRAIC[56] is r_inact in component SR_Ca_uptake_and_release (per_second). * STATES[25] is Ca_Calse in component SR_Ca_uptake_and_release (dimensionless). * ALGEBRAIC[57] is phi_Calse in component SR_Ca_uptake_and_release (per_second). * STATES[26] is F1 in component SR_Ca_uptake_and_release (dimensionless). * STATES[27] is F2 in component SR_Ca_uptake_and_release (dimensionless). * STATES[28] is F3 in component SR_Ca_uptake_and_release (dimensionless). * CONSTANTS[53] is V_up in component SR_Ca_uptake_and_release (mm_cubed). * CONSTANTS[54] is V_rel in component SR_Ca_uptake_and_release (mm_cubed). * CONSTANTS[37] is cGMP in component cAMP_balance (millimolar). * CONSTANTS[38] is Iso in component cAMP_balance (millimolar). * CONSTANTS[39] is Km_Iso in component cAMP_balance (millimolar). * CONSTANTS[40] is Km_ACh in component cAMP_balance (millimolar). * CONSTANTS[41] is K_PDE in component cAMP_balance (dimensionless). * CONSTANTS[42] is K_ADC in component cAMP_balance (millimolar_per_second). * CONSTANTS[43] is V_PDE in component cAMP_balance (per_second). * CONSTANTS[44] is P_M2_ADC in component cAMP_balance (dimensionless). * RATES[0] is d/dt V in component membrane (millivolt). * RATES[2] is d/dt m in component sodium_current_m_gate (dimensionless). * RATES[3] is d/dt h1 in component sodium_current_h_gate (dimensionless). * RATES[4] is d/dt h2 in component sodium_current_h_gate (dimensionless). * RATES[6] is d/dt d_L in component L_type_Ca_channel_d_gate (dimensionless). * RATES[7] is d/dt f_L in component L_type_Ca_channel_f_gate (dimensionless). * RATES[8] is d/dt d_T in component T_type_Ca_channel_d_gate (dimensionless). * RATES[9] is d/dt f_T in component T_type_Ca_channel_f_gate (dimensionless). * RATES[10] is d/dt P_a in component delayed_rectifying_potassium_current_P_a_gate (dimensionless). * RATES[11] is d/dt P_i in component delayed_rectifying_potassium_current_P_i_gate (dimensionless). * RATES[12] is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless). * RATES[17] is d/dt a in component muscarinic_potassium_current (dimensionless). * RATES[19] is d/dt Ca_Calmod in component intracellular_concentrations_and_buffer_equations (dimensionless). * RATES[20] is d/dt Ca_Trop in component intracellular_concentrations_and_buffer_equations (dimensionless). * RATES[21] is d/dt Ca_Mg_Trop in component intracellular_concentrations_and_buffer_equations (dimensionless). * RATES[22] is d/dt Mg_Mg_Trop in component intracellular_concentrations_and_buffer_equations (dimensionless). * RATES[13] is d/dt Na_i in component intracellular_concentrations_and_buffer_equations (millimolar). * RATES[18] is d/dt K_i in component intracellular_concentrations_and_buffer_equations (millimolar). * RATES[15] is d/dt Ca_i in component intracellular_concentrations_and_buffer_equations (millimolar). * RATES[1] is d/dt Na_c in component cleft_space_equations (millimolar). * RATES[14] is d/dt K_c in component cleft_space_equations (millimolar). * RATES[16] is d/dt Ca_c in component cleft_space_equations (millimolar). * RATES[25] is d/dt Ca_Calse in component SR_Ca_uptake_and_release (dimensionless). * RATES[26] is d/dt F1 in component SR_Ca_uptake_and_release (dimensionless). * RATES[27] is d/dt F2 in component SR_Ca_uptake_and_release (dimensionless). * RATES[28] is d/dt F3 in component SR_Ca_uptake_and_release (dimensionless). * RATES[23] is d/dt Ca_up in component SR_Ca_uptake_and_release (millimolar). * RATES[24] is d/dt Ca_rel in component SR_Ca_uptake_and_release (millimolar). * RATES[5] is d/dt cAMP in component cAMP_balance (millimolar). * There are a total of 1 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = -49.54105; CONSTANTS[0] = 8314.472; CONSTANTS[1] = 310; CONSTANTS[2] = 96485.3415; CONSTANTS[3] = 5.5e-5; CONSTANTS[4] = 0.00344; STATES[1] = 139.9988; CONSTANTS[5] = 0; STATES[2] = 0.250113; STATES[3] = 0.001386897; STATES[4] = 0.002065463; CONSTANTS[6] = 0.02115; STATES[5] = 3e-3; CONSTANTS[7] = 46.4; STATES[6] = 0.002572773; STATES[7] = 0.98651; CONSTANTS[8] = 0.02521; CONSTANTS[9] = 45; STATES[8] = 0.02012114; STATES[9] = 0.1945111; CONSTANTS[10] = 5.4; STATES[10] = 0.02302278; STATES[11] = 0.3777728; CONSTANTS[11] = 0.00016; CONSTANTS[12] = 0.0000364; CONSTANTS[13] = 0.0000694; CONSTANTS[14] = 0.0067478; CONSTANTS[15] = 0.0128821; STATES[12] = 0.09227776; CONSTANTS[16] = 5.46; CONSTANTS[17] = 0.621; CONSTANTS[18] = 0.2192; STATES[13] = 9.701621; STATES[14] = 5.389014; STATES[15] = 3.787018e-4; CONSTANTS[19] = 0.02869; CONSTANTS[20] = 0.00001248; CONSTANTS[21] = 0.0001; CONSTANTS[22] = 0.5; STATES[16] = 2.00474; CONSTANTS[23] = 7.833e-3; STATES[17] = 0; CONSTANTS[24] = 200; STATES[18] = 1.407347e2; STATES[19] = 0.1411678; STATES[20] = 0.07331396; STATES[21] = 0.7618549; STATES[22] = 0.2097049; CONSTANTS[25] = 2.5; CONSTANTS[26] = 3.497e-6; CONSTANTS[27] = 140; CONSTANTS[28] = 2; CONSTANTS[29] = 0.01; STATES[23] = 16.95311; CONSTANTS[30] = 0.08; CONSTANTS[31] = 0.072; STATES[24] = 16.85024; CONSTANTS[32] = 0.5; CONSTANTS[33] = 0.00005; CONSTANTS[34] = 0.9; CONSTANTS[35] = 22; CONSTANTS[36] = 0.004; STATES[25] = 0.9528726; STATES[26] = 0.1133251; STATES[27] = 0.0007594214; STATES[28] = 0.8859153; CONSTANTS[37] = 2e-3; CONSTANTS[38] = 0; CONSTANTS[39] = 0.14e-3; CONSTANTS[40] = 0.14e-3; CONSTANTS[41] = 6; CONSTANTS[42] = 8e-3; CONSTANTS[43] = 20; CONSTANTS[44] = 0.02; CONSTANTS[45] = 1.00000 - CONSTANTS[5]/(CONSTANTS[5]+0.00100000); CONSTANTS[46] = 1.00000 - CONSTANTS[5]/(CONSTANTS[5]+0.500000); CONSTANTS[47] = (CONSTANTS[24]<100.000&&CONSTANTS[24]>25.0000 ? 1.02600/(1.00000+exp((CONSTANTS[24]+11.0500)/- 7.50950)) - 0.990000 : 0.000600000); CONSTANTS[48] = 12.3200/(1.00000+0.00420000/CONSTANTS[5]); CONSTANTS[49] = 0.465000*CONSTANTS[26]; CONSTANTS[50] = 0.136000*CONSTANTS[26]; CONSTANTS[51] = ( CONSTANTS[33]*CONSTANTS[34])/CONSTANTS[35]; CONSTANTS[52] = CONSTANTS[47]*CONSTANTS[23]; CONSTANTS[53] = 0.0116600*CONSTANTS[49]; CONSTANTS[54] = 0.00129600*CONSTANTS[49]; RATES[0] = 0.1001; RATES[2] = 0.1001; RATES[3] = 0.1001; RATES[4] = 0.1001; RATES[6] = 0.1001; RATES[7] = 0.1001; RATES[8] = 0.1001; RATES[9] = 0.1001; RATES[10] = 0.1001; RATES[11] = 0.1001; RATES[12] = 0.1001; RATES[17] = 0.1001; RATES[19] = 0.1001; RATES[20] = 0.1001; RATES[21] = 0.1001; RATES[22] = 0.1001; RATES[13] = 0.1001; RATES[18] = 0.1001; RATES[15] = 0.1001; RATES[1] = 0.1001; RATES[14] = 0.1001; RATES[16] = 0.1001; RATES[25] = 0.1001; RATES[26] = 0.1001; RATES[27] = 0.1001; RATES[28] = 0.1001; RATES[23] = 0.1001; RATES[24] = 0.1001; RATES[5] = 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[64]+ALGEBRAIC[20]+ALGEBRAIC[65]+ALGEBRAIC[66]+ALGEBRAIC[37]+ALGEBRAIC[70]+ALGEBRAIC[42]+ALGEBRAIC[44]+ALGEBRAIC[43]+ALGEBRAIC[72])/CONSTANTS[3]; resid[1] = RATES[2] - (ALGEBRAIC[2] - STATES[2])/ALGEBRAIC[3]; resid[2] = RATES[3] - (ALGEBRAIC[6] - STATES[3])/ALGEBRAIC[7]; resid[3] = RATES[4] - (ALGEBRAIC[8] - STATES[4])/ALGEBRAIC[9]; resid[4] = RATES[6] - (ALGEBRAIC[15] - STATES[6])/ALGEBRAIC[14]; resid[5] = RATES[7] - (ALGEBRAIC[19] - STATES[7])/ALGEBRAIC[18]; resid[6] = RATES[8] - (ALGEBRAIC[24] - STATES[8])/ALGEBRAIC[23]; resid[7] = RATES[9] - (ALGEBRAIC[28] - STATES[9])/ALGEBRAIC[27]; resid[8] = RATES[10] - (ALGEBRAIC[31] - STATES[10])/ALGEBRAIC[32]; resid[9] = RATES[11] - ALGEBRAIC[33]*(1.00000 - STATES[11]) - ALGEBRAIC[34]*STATES[11]; resid[10] = RATES[12] - (ALGEBRAIC[39] - STATES[12])/ALGEBRAIC[40]; resid[11] = RATES[17] - CONSTANTS[48]*(1.00000 - STATES[17]) - ALGEBRAIC[45]*STATES[17]; resid[12] = RATES[19] - ALGEBRAIC[46]; resid[13] = RATES[20] - ALGEBRAIC[47]; resid[14] = RATES[21] - ALGEBRAIC[48]; resid[15] = RATES[22] - ALGEBRAIC[49]; resid[16] = RATES[13] - - ( 3.00000*ALGEBRAIC[42]+ 3.00000*ALGEBRAIC[44]+ALGEBRAIC[67]+ALGEBRAIC[35]+ALGEBRAIC[64])/( CONSTANTS[2]*CONSTANTS[49]); resid[17] = RATES[18] - ( 2.00000*ALGEBRAIC[42] - (ALGEBRAIC[66]+ALGEBRAIC[36]+ALGEBRAIC[69]))/( CONSTANTS[2]*CONSTANTS[49]); resid[18] = RATES[15] - (( 2.00000*ALGEBRAIC[44]+ALGEBRAIC[58]) - (ALGEBRAIC[65]+ALGEBRAIC[20]+ALGEBRAIC[43]+ALGEBRAIC[68]+ALGEBRAIC[60]))/( 2.00000*CONSTANTS[49]*CONSTANTS[2]) - ALGEBRAIC[53]; resid[19] = RATES[1] - (CONSTANTS[27] - STATES[1])/CONSTANTS[29]+(ALGEBRAIC[64]+ 3.00000*ALGEBRAIC[44]+ 3.00000*ALGEBRAIC[42]+ALGEBRAIC[67]+ALGEBRAIC[35])/( CONSTANTS[2]*CONSTANTS[50]); resid[20] = RATES[14] - (CONSTANTS[10] - STATES[14])/CONSTANTS[29]+( - 2.00000*ALGEBRAIC[42]+ALGEBRAIC[66]+ALGEBRAIC[69]+ALGEBRAIC[36])/( CONSTANTS[2]*CONSTANTS[50]); resid[21] = RATES[16] - (CONSTANTS[28] - STATES[16])/CONSTANTS[29]+( - 2.00000*ALGEBRAIC[44]+ALGEBRAIC[65]+ALGEBRAIC[20]+ALGEBRAIC[43]+ALGEBRAIC[68])/( 2.00000*CONSTANTS[2]*CONSTANTS[50]); resid[22] = RATES[25] - ALGEBRAIC[57]; resid[23] = RATES[26] - 0.960000*STATES[28] - ALGEBRAIC[55]*STATES[26]; resid[24] = RATES[27] - ALGEBRAIC[55]*STATES[26] - ALGEBRAIC[56]*STATES[27]; resid[25] = RATES[28] - ALGEBRAIC[56]*STATES[27] - 0.960000*STATES[28]; resid[26] = RATES[23] - (ALGEBRAIC[60] - ALGEBRAIC[59])/( 2.00000*CONSTANTS[53]*CONSTANTS[2]); resid[27] = RATES[24] - (ALGEBRAIC[59] - ALGEBRAIC[58])/( 2.00000*CONSTANTS[54]*CONSTANTS[2]) - 11.4800*ALGEBRAIC[57]; resid[28] = RATES[5] - 1.00000*( CONSTANTS[42]*((1.00000+CONSTANTS[38]/(CONSTANTS[38]+CONSTANTS[39])) - ( CONSTANTS[44]*CONSTANTS[5])/( CONSTANTS[44]*CONSTANTS[5]+CONSTANTS[40])) - ( CONSTANTS[43]*CONSTANTS[37]*STATES[5])/(STATES[5]+ CONSTANTS[41]*CONSTANTS[37])); } 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] = ( - 824.000*(STATES[0]+51.9000))/(exp((STATES[0]+51.9000)/- 8.90000) - 1.00000); ALGEBRAIC[1] = 32960.0*exp((STATES[0]+51.9000)/- 8.90000); ALGEBRAIC[2] = ALGEBRAIC[0]/(ALGEBRAIC[0]+ALGEBRAIC[1]); ALGEBRAIC[3] = 1.00000/(ALGEBRAIC[0]+ALGEBRAIC[1])+1.50000e-05; ALGEBRAIC[4] = 165.000*exp((STATES[0]+101.300)/- 12.6000); ALGEBRAIC[5] = 12360.0/( 320.000*exp((STATES[0]+101.300)/- 12.6000)+1.00000); ALGEBRAIC[6] = ALGEBRAIC[4]/(ALGEBRAIC[4]+ALGEBRAIC[5]); ALGEBRAIC[7] = 1.00000/(ALGEBRAIC[4]+ALGEBRAIC[5]); ALGEBRAIC[8] = ALGEBRAIC[6]; ALGEBRAIC[9] = 20.0000*ALGEBRAIC[7]; ALGEBRAIC[12] = ( - 28.3900*(STATES[0]+35.0000))/(exp((STATES[0]+35.0000)/- 2.50000) - 1.00000)+( - 84.9000*STATES[0])/(exp( - 0.208000*STATES[0]) - 1.00000); ALGEBRAIC[13] = ( 11.4300*(STATES[0] - 5.00000))/(exp( 0.400000*(STATES[0] - 5.00000)) - 1.00000); ALGEBRAIC[14] = 1.00000/(ALGEBRAIC[12]+ALGEBRAIC[13]); ALGEBRAIC[15] = 1.00000/(1.00000+exp((STATES[0]+14.1000)/- 6.00000)); ALGEBRAIC[16] = ( 3.75000*(STATES[0]+28.0000))/(exp((STATES[0]+28.0000)/4.00000) - 1.00000); ALGEBRAIC[17] = 30.0000/(1.00000+exp((STATES[0]+28.0000)/- 4.00000)); ALGEBRAIC[18] = 1.00000/(ALGEBRAIC[16]+ALGEBRAIC[17]); ALGEBRAIC[19] = 1.00000/(1.00000+exp((STATES[0]+30.0000)/5.00000)); ALGEBRAIC[20] = CONSTANTS[8]*STATES[8]*STATES[9]*(STATES[0] - CONSTANTS[9]); ALGEBRAIC[21] = 1068.00*exp((STATES[0]+26.3000)/30.0000); ALGEBRAIC[22] = 1068.00*exp((STATES[0]+26.3000)/- 30.0000); ALGEBRAIC[23] = 1.00000/(ALGEBRAIC[21]+ALGEBRAIC[22]); ALGEBRAIC[24] = 1.00000/(1.00000+exp((STATES[0]+26.3000)/- 6.00000)); ALGEBRAIC[25] = 15.3000*exp((STATES[0]+61.7000)/- 83.3000); ALGEBRAIC[26] = 15.0000*exp((STATES[0]+61.7000)/15.3800); ALGEBRAIC[27] = 1.00000/(ALGEBRAIC[25]+ALGEBRAIC[26]); ALGEBRAIC[28] = 1.00000/(1.00000+exp((STATES[0]+61.7000)/5.60000)); ALGEBRAIC[31] = 1.00000/(1.00000+exp((STATES[0]+5.10000)/- 7.40000)); ALGEBRAIC[32] = 1.00000/( 17.0000*exp( 0.0398000*STATES[0])+ 2.11000*exp( - 0.0510000*STATES[0])); ALGEBRAIC[33] = 100.000*exp( - 0.0183000*STATES[0]); ALGEBRAIC[34] = 656.000*exp( 0.00942000*STATES[0]); ALGEBRAIC[35] = CONSTANTS[14]*pow(STATES[12], 2.00000)*(STATES[0] - 75.0000); ALGEBRAIC[36] = CONSTANTS[15]*pow(STATES[12], 2.00000)*(STATES[0]+85.0000); ALGEBRAIC[37] = ALGEBRAIC[35]+ALGEBRAIC[36]; ALGEBRAIC[38] = 20.5000/(1.00000+exp((STATES[5] - 0.00340000)/- 0.000500000)) - 78.5600; ALGEBRAIC[39] = 1.00000/(1.00000+exp((STATES[0] - ALGEBRAIC[38])/9.00000)); ALGEBRAIC[40] = 1.00000/( 1.64830*exp((STATES[0]+54.0600)/- 24.3300)+14.0106/(0.700000+exp((STATES[0]+60.0000)/- 5.50000))); ALGEBRAIC[41] = 1.60000/(1.00000+exp((STATES[5] - 0.00375000)/- 0.000150000))+0.990000; ALGEBRAIC[42] = (CONDVAR[0]<0.00000 ? ( (2.00000 - ALGEBRAIC[41])*CONSTANTS[18]*pow(STATES[13]/(CONSTANTS[16]+STATES[13]), 3.00000)*pow(STATES[14]/(CONSTANTS[17]+STATES[14]), 2.00000)*1.60000)/(1.50000+exp((STATES[0]+60.0000)/- 40.0000)) : ( ALGEBRAIC[41]*CONSTANTS[18]*pow(STATES[13]/(CONSTANTS[16]+STATES[13]), 3.00000)*pow(STATES[14]/(CONSTANTS[17]+STATES[14]), 2.00000)*1.60000)/(1.50000+exp((STATES[0]+60.0000)/- 40.0000))); ALGEBRAIC[43] = ( CONSTANTS[19]*STATES[15])/(STATES[15]+0.000400000); ALGEBRAIC[44] = ( CONSTANTS[20]*( pow(STATES[13], 3.00000)*STATES[16]*exp( 0.0374300*STATES[0]*CONSTANTS[22]) - pow(STATES[1], 3.00000)*STATES[15]*exp( 0.0374300*STATES[0]*(CONSTANTS[22] - 1.00000))))/(1.00000+ CONSTANTS[21]*( STATES[15]*pow(STATES[1], 3.00000)+ STATES[16]*pow(STATES[13], 3.00000))); ALGEBRAIC[45] = 17.0000*exp( 0.0133000*(STATES[0]+40.0000)); ALGEBRAIC[46] = 129000.*STATES[15]*(1.00000 - STATES[19]) - 307.000*STATES[19]; ALGEBRAIC[47] = 50500.0*STATES[15]*(1.00000 - STATES[20]) - 252.000*STATES[20]; ALGEBRAIC[48] = 129000.*STATES[15]*(1.00000 - (STATES[21]+STATES[22])) - 4.25000*STATES[21]; ALGEBRAIC[49] = 1290.00*CONSTANTS[25]*(1.00000 - (STATES[21]+STATES[22])) - 429.000*STATES[22]; ALGEBRAIC[50] = 0.0900000*ALGEBRAIC[46]; ALGEBRAIC[51] = 0.0310000*ALGEBRAIC[47]; ALGEBRAIC[52] = 0.0620000*ALGEBRAIC[48]; ALGEBRAIC[53] = ALGEBRAIC[50]+ALGEBRAIC[51]+ALGEBRAIC[52]; ALGEBRAIC[55] = 240.000*exp( (STATES[0] - 40.0000)*0.0800000)+ 240.000*pow(STATES[15]/(STATES[15]+CONSTANTS[36]), 4.00000); ALGEBRAIC[56] = 40.0000+ 240.000*pow(STATES[15]/(STATES[15]+CONSTANTS[36]), 4.00000); ALGEBRAIC[57] = 770.000*STATES[24]*(1.00000 - STATES[25]) - 641.000*STATES[25]; ALGEBRAIC[58] = CONSTANTS[32]*pow(STATES[27]/(STATES[27]+0.250000), 2.00000)*STATES[24]; ALGEBRAIC[59] = ( (STATES[23] - STATES[24])*2.00000*CONSTANTS[2]*CONSTANTS[53])/0.0641800; ALGEBRAIC[54] = STATES[15]+ STATES[23]*CONSTANTS[51]+ CONSTANTS[33]*CONSTANTS[34]+CONSTANTS[33]; ALGEBRAIC[60] = ( CONSTANTS[30]*STATES[15] - CONSTANTS[31]*STATES[23]*CONSTANTS[51])/ALGEBRAIC[54]; ALGEBRAIC[61] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(STATES[1]/STATES[13]); ALGEBRAIC[64] = ( (( CONSTANTS[45]*CONSTANTS[4]*pow(STATES[2], 3.00000)*STATES[3]*STATES[4]*STATES[1]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*(exp(( (STATES[0] - ALGEBRAIC[61])*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[10] = 0.400000*(1.00000+( 4.50000*STATES[5])/(STATES[5]+0.00650000))+0.0315700; ALGEBRAIC[11] = CONSTANTS[6]*ALGEBRAIC[10]; ALGEBRAIC[65] = ALGEBRAIC[11]*( STATES[7]*STATES[6]+ 0.0950000*ALGEBRAIC[15])*(STATES[0] - CONSTANTS[7]); ALGEBRAIC[29] = 0.620000*(1.00000+( 2.61290*STATES[5])/(STATES[5]+0.00900000)) - 0.0250000; ALGEBRAIC[30] = ALGEBRAIC[29]*0.00693000*pow(CONSTANTS[10]/1.00000, 0.590000); ALGEBRAIC[62] = (( CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(STATES[14]/STATES[18]); ALGEBRAIC[66] = ALGEBRAIC[30]*STATES[10]*STATES[11]*(STATES[0] - ALGEBRAIC[62]); ALGEBRAIC[67] = CONSTANTS[46]*CONSTANTS[11]*(STATES[0] - ALGEBRAIC[61]); ALGEBRAIC[63] = (( 0.500000*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2])*log(STATES[16]/STATES[15]); ALGEBRAIC[68] = CONSTANTS[12]*(STATES[0] - ALGEBRAIC[63]); ALGEBRAIC[69] = CONSTANTS[13]*(STATES[0] - ALGEBRAIC[62]); ALGEBRAIC[70] = ALGEBRAIC[67]+ALGEBRAIC[68]+ALGEBRAIC[69]; ALGEBRAIC[71] = CONSTANTS[52]*(STATES[0] - ALGEBRAIC[62]); ALGEBRAIC[72] = STATES[17]*ALGEBRAIC[71]; } 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; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = STATES[15] - 0.000150000; }