/* There are a total of 71 entries in the algebraic variable array. There are a total of 41 entries in each of the rate and state variable arrays. There are a total of 73 entries in the constant variable array. */ /* * VOI is time in component environment (millisecond). * STATES[0] is V in component membrane (millivolt). * CONSTANTS[0] is Cm in component membrane (microF_per_cm2). * CONSTANTS[1] is Vmyo in component membrane (microlitre). * CONSTANTS[2] is VJSR in component membrane (microlitre). * CONSTANTS[3] is VNSR in component membrane (microlitre). * CONSTANTS[4] is Vss in component membrane (microlitre). * CONSTANTS[5] is Acap in component membrane (cm2). * CONSTANTS[6] is Ko in component membrane (micromolar). * CONSTANTS[7] is Nao in component membrane (micromolar). * CONSTANTS[8] is Cao in component membrane (micromolar). * CONSTANTS[9] is R in component membrane (joule_per_mole_kelvin). * CONSTANTS[10] is T in component membrane (kelvin). * CONSTANTS[11] is F in component membrane (coulomb_per_millimole). * ALGEBRAIC[0] is i_stim in component membrane (picoA_per_picoF). * ALGEBRAIC[11] is i_CaL in component L_type_calcium_current (picoA_per_picoF). * ALGEBRAIC[17] is i_pCa in component calcium_pump_current (picoA_per_picoF). * ALGEBRAIC[18] is i_NaCa in component sodium_calcium_exchange_current (picoA_per_picoF). * ALGEBRAIC[66] is i_Cab in component calcium_background_current (picoA_per_picoF). * ALGEBRAIC[67] is i_Na in component fast_sodium_current (picoA_per_picoF). * ALGEBRAIC[36] is i_Nab in component sodium_background_current (picoA_per_picoF). * ALGEBRAIC[69] is i_NaK in component sodium_potassium_pump_current (picoA_per_picoF). * ALGEBRAIC[68] is i_Kto_f in component fast_transient_outward_potassium_current (picoA_per_picoF). * ALGEBRAIC[42] is i_Kto_s in component slow_transient_outward_potassium_current (picoA_per_picoF). * ALGEBRAIC[47] is i_K1 in component time_independent_potassium_current (picoA_per_picoF). * ALGEBRAIC[48] is i_Ks in component slow_delayed_rectifier_potassium_current (picoA_per_picoF). * ALGEBRAIC[51] is i_Kur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (picoA_per_picoF). * ALGEBRAIC[54] is i_Kss in component non_inactivating_steady_state_potassium_current (picoA_per_picoF). * ALGEBRAIC[70] is i_ClCa in component calcium_activated_chloride_current (picoA_per_picoF). * ALGEBRAIC[56] is i_Kr in component rapid_delayed_rectifier_potassium_current (picoA_per_picoF). * CONSTANTS[12] is stim_start in component membrane (millisecond). * CONSTANTS[13] is stim_end in component membrane (millisecond). * CONSTANTS[14] is stim_period in component membrane (millisecond). * CONSTANTS[15] is stim_duration in component membrane (millisecond). * CONSTANTS[16] is stim_amplitude in component membrane (picoA_per_picoF). * STATES[1] is Cai in component calcium_concentration (micromolar). * STATES[2] is Cass in component calcium_concentration (micromolar). * STATES[3] is CaJSR in component calcium_concentration (micromolar). * STATES[4] is CaNSR in component calcium_concentration (micromolar). * ALGEBRAIC[1] is Bi in component calcium_concentration (dimensionless). * ALGEBRAIC[2] is Bss in component calcium_concentration (dimensionless). * ALGEBRAIC[3] is BJSR in component calcium_concentration (dimensionless). * CONSTANTS[17] is CMDN_tot in component calcium_concentration (micromolar). * CONSTANTS[18] is CSQN_tot in component calcium_concentration (micromolar). * CONSTANTS[19] is Km_CMDN in component calcium_concentration (micromolar). * CONSTANTS[20] is Km_CSQN in component calcium_concentration (micromolar). * ALGEBRAIC[7] is J_leak in component calcium_fluxes (micromolar_per_millisecond). * ALGEBRAIC[4] is J_rel in component calcium_fluxes (micromolar_per_millisecond). * ALGEBRAIC[8] is J_up in component calcium_fluxes (micromolar_per_millisecond). * ALGEBRAIC[5] is J_tr in component calcium_fluxes (micromolar_per_millisecond). * ALGEBRAIC[9] is J_trpn in component calcium_fluxes (micromolar_per_millisecond). * ALGEBRAIC[6] is J_xfer in component calcium_fluxes (micromolar_per_millisecond). * CONSTANTS[21] is k_plus_htrpn in component calcium_fluxes (per_micromolar_millisecond). * CONSTANTS[22] is k_minus_htrpn in component calcium_fluxes (per_millisecond). * CONSTANTS[23] is k_plus_ltrpn in component calcium_fluxes (per_micromolar_millisecond). * CONSTANTS[24] is k_minus_ltrpn in component calcium_fluxes (per_millisecond). * STATES[5] is P_RyR in component calcium_fluxes (dimensionless). * CONSTANTS[25] is v1 in component calcium_fluxes (per_millisecond). * CONSTANTS[26] is tau_tr in component calcium_fluxes (millisecond). * CONSTANTS[27] is v2 in component calcium_fluxes (per_millisecond). * CONSTANTS[28] is tau_xfer in component calcium_fluxes (millisecond). * CONSTANTS[29] is v3 in component calcium_fluxes (micromolar_per_millisecond). * CONSTANTS[30] is Km_up in component calcium_fluxes (micromolar). * CONSTANTS[31] is LTRPN_tot in component calcium_buffering (micromolar). * CONSTANTS[32] is HTRPN_tot in component calcium_buffering (micromolar). * STATES[6] is LTRPN_Ca in component calcium_buffering (micromolar). * STATES[7] is HTRPN_Ca in component calcium_buffering (micromolar). * CONSTANTS[33] is i_CaL_max in component L_type_calcium_current (picoA_per_picoF). * STATES[8] is P_O1 in component ryanodine_receptors (dimensionless). * STATES[9] is P_O2 in component ryanodine_receptors (dimensionless). * ALGEBRAIC[10] is P_C1 in component ryanodine_receptors (dimensionless). * STATES[10] is P_C2 in component ryanodine_receptors (dimensionless). * CONSTANTS[34] is k_plus_a in component ryanodine_receptors (micromolar4_per_millisecond). * CONSTANTS[35] is k_minus_a in component ryanodine_receptors (per_millisecond). * CONSTANTS[36] is k_plus_b in component ryanodine_receptors (micromolar3_per_millisecond). * CONSTANTS[37] is k_minus_b in component ryanodine_receptors (per_millisecond). * CONSTANTS[38] is k_plus_c in component ryanodine_receptors (per_millisecond). * CONSTANTS[39] is k_minus_c in component ryanodine_receptors (per_millisecond). * CONSTANTS[40] is m in component ryanodine_receptors (dimensionless). * CONSTANTS[41] is n in component ryanodine_receptors (dimensionless). * CONSTANTS[42] is E_CaL in component L_type_calcium_current (millivolt). * CONSTANTS[43] is g_CaL in component L_type_calcium_current (milliS_per_microF). * STATES[11] is O in component L_type_calcium_current (dimensionless). * ALGEBRAIC[12] is C1 in component L_type_calcium_current (dimensionless). * STATES[12] is C2 in component L_type_calcium_current (dimensionless). * STATES[13] is C3 in component L_type_calcium_current (dimensionless). * STATES[14] is C4 in component L_type_calcium_current (dimensionless). * STATES[15] is I1 in component L_type_calcium_current (dimensionless). * STATES[16] is I2 in component L_type_calcium_current (dimensionless). * STATES[17] is I3 in component L_type_calcium_current (dimensionless). * ALGEBRAIC[13] is alpha in component L_type_calcium_current (per_millisecond). * ALGEBRAIC[14] is beta in component L_type_calcium_current (per_millisecond). * ALGEBRAIC[15] is gamma in component L_type_calcium_current (per_millisecond). * ALGEBRAIC[16] is Kpcf in component L_type_calcium_current (per_millisecond). * CONSTANTS[44] is Kpcb in component L_type_calcium_current (per_millisecond). * CONSTANTS[45] is Kpc_max in component L_type_calcium_current (per_millisecond). * CONSTANTS[46] is Kpc_half in component L_type_calcium_current (micromolar). * CONSTANTS[47] is i_pCa_max in component calcium_pump_current (picoA_per_picoF). * CONSTANTS[48] is Km_pCa in component calcium_pump_current (micromolar). * CONSTANTS[49] is k_NaCa in component sodium_calcium_exchange_current (picoA_per_picoF). * CONSTANTS[50] is K_mNa in component sodium_calcium_exchange_current (micromolar). * CONSTANTS[51] is K_mCa in component sodium_calcium_exchange_current (micromolar). * CONSTANTS[52] is k_sat in component sodium_calcium_exchange_current (dimensionless). * CONSTANTS[53] is eta in component sodium_calcium_exchange_current (dimensionless). * STATES[18] is Nai in component sodium_concentration (micromolar). * CONSTANTS[54] is g_Cab in component calcium_background_current (milliS_per_microF). * ALGEBRAIC[19] is E_CaN in component calcium_background_current (millivolt). * ALGEBRAIC[20] is E_Na in component fast_sodium_current (millivolt). * CONSTANTS[55] is g_Na in component fast_sodium_current (milliS_per_microF). * STATES[19] is O_Na in component fast_sodium_current (dimensionless). * STATES[20] is C_Na1 in component fast_sodium_current (dimensionless). * STATES[21] is C_Na2 in component fast_sodium_current (dimensionless). * ALGEBRAIC[21] is C_Na3 in component fast_sodium_current (dimensionless). * STATES[22] is I1_Na in component fast_sodium_current (dimensionless). * STATES[23] is I2_Na in component fast_sodium_current (dimensionless). * STATES[24] is IF_Na in component fast_sodium_current (dimensionless). * STATES[25] is IC_Na2 in component fast_sodium_current (dimensionless). * STATES[26] is IC_Na3 in component fast_sodium_current (dimensionless). * ALGEBRAIC[22] is alpha_Na11 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[25] is beta_Na11 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[23] is alpha_Na12 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[26] is beta_Na12 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[24] is alpha_Na13 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[27] is beta_Na13 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[28] is alpha_Na3 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[29] is beta_Na3 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[30] is alpha_Na2 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[31] is beta_Na2 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[32] is alpha_Na4 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[33] is beta_Na4 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[34] is alpha_Na5 in component fast_sodium_current (per_millisecond). * ALGEBRAIC[35] is beta_Na5 in component fast_sodium_current (per_millisecond). * STATES[27] is Ki in component potassium_concentration (micromolar). * CONSTANTS[56] is g_Nab in component sodium_background_current (milliS_per_microF). * ALGEBRAIC[37] is E_K in component fast_transient_outward_potassium_current (millivolt). * CONSTANTS[57] is g_Kto_f in component fast_transient_outward_potassium_current (milliS_per_microF). * STATES[28] is ato_f in component fast_transient_outward_potassium_current (dimensionless). * STATES[29] is ito_f in component fast_transient_outward_potassium_current (dimensionless). * ALGEBRAIC[38] is alpha_a in component fast_transient_outward_potassium_current (per_millisecond). * ALGEBRAIC[39] is beta_a in component fast_transient_outward_potassium_current (per_millisecond). * ALGEBRAIC[40] is alpha_i in component fast_transient_outward_potassium_current (per_millisecond). * ALGEBRAIC[41] is beta_i in component fast_transient_outward_potassium_current (per_millisecond). * ALGEBRAIC[43] is ass in component slow_transient_outward_potassium_current (dimensionless). * ALGEBRAIC[44] is iss in component slow_transient_outward_potassium_current (dimensionless). * CONSTANTS[58] is g_Kto_s in component slow_transient_outward_potassium_current (milliS_per_microF). * STATES[30] is ato_s in component slow_transient_outward_potassium_current (dimensionless). * STATES[31] is ito_s in component slow_transient_outward_potassium_current (dimensionless). * ALGEBRAIC[45] is tau_ta_s in component slow_transient_outward_potassium_current (millisecond). * ALGEBRAIC[46] is tau_ti_s in component slow_transient_outward_potassium_current (millisecond). * CONSTANTS[59] is g_Ks in component slow_delayed_rectifier_potassium_current (milliS_per_microF). * STATES[32] is nKs in component slow_delayed_rectifier_potassium_current (dimensionless). * ALGEBRAIC[49] is alpha_n in component slow_delayed_rectifier_potassium_current (per_millisecond). * ALGEBRAIC[50] is beta_n in component slow_delayed_rectifier_potassium_current (per_millisecond). * CONSTANTS[60] is g_Kur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (milliS_per_microF). * STATES[33] is aur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless). * STATES[34] is iur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless). * ALGEBRAIC[52] is tau_aur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (millisecond). * ALGEBRAIC[53] is tau_iur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (millisecond). * CONSTANTS[61] is g_Kss in component non_inactivating_steady_state_potassium_current (milliS_per_microF). * STATES[35] is aKss in component non_inactivating_steady_state_potassium_current (dimensionless). * STATES[36] is iKss in component non_inactivating_steady_state_potassium_current (dimensionless). * ALGEBRAIC[55] is tau_Kss in component non_inactivating_steady_state_potassium_current (millisecond). * CONSTANTS[62] is g_Kr in component rapid_delayed_rectifier_potassium_current (milliS_per_microF). * STATES[37] is O_K in component rapid_delayed_rectifier_potassium_current (dimensionless). * STATES[38] is C_K1 in component rapid_delayed_rectifier_potassium_current (dimensionless). * STATES[39] is C_K2 in component rapid_delayed_rectifier_potassium_current (dimensionless). * ALGEBRAIC[57] is C_K0 in component rapid_delayed_rectifier_potassium_current (dimensionless). * STATES[40] is I_K in component rapid_delayed_rectifier_potassium_current (dimensionless). * ALGEBRAIC[58] is alpha_a0 in component rapid_delayed_rectifier_potassium_current (per_millisecond). * ALGEBRAIC[59] is beta_a0 in component rapid_delayed_rectifier_potassium_current (per_millisecond). * CONSTANTS[63] is kb in component rapid_delayed_rectifier_potassium_current (per_millisecond). * CONSTANTS[64] is kf in component rapid_delayed_rectifier_potassium_current (per_millisecond). * ALGEBRAIC[60] is alpha_a1 in component rapid_delayed_rectifier_potassium_current (per_millisecond). * ALGEBRAIC[61] is beta_a1 in component rapid_delayed_rectifier_potassium_current (per_millisecond). * ALGEBRAIC[62] is alpha_i in component rapid_delayed_rectifier_potassium_current (per_millisecond). * ALGEBRAIC[63] is beta_i in component rapid_delayed_rectifier_potassium_current (per_millisecond). * CONSTANTS[65] is i_NaK_max in component sodium_potassium_pump_current (picoA_per_picoF). * CONSTANTS[66] is Km_Nai in component sodium_potassium_pump_current (micromolar). * CONSTANTS[67] is Km_Ko in component sodium_potassium_pump_current (micromolar). * ALGEBRAIC[64] is f_NaK in component sodium_potassium_pump_current (dimensionless). * CONSTANTS[71] is sigma in component sodium_potassium_pump_current (dimensionless). * CONSTANTS[68] is g_ClCa in component calcium_activated_chloride_current (milliS_per_microF). * ALGEBRAIC[65] is O_ClCa in component calcium_activated_chloride_current (dimensionless). * CONSTANTS[69] is E_Cl in component calcium_activated_chloride_current (millivolt). * CONSTANTS[70] is Km_Cl in component calcium_activated_chloride_current (micromolar). * RATES[0] is d/dt V in component membrane (millivolt). * RATES[1] is d/dt Cai in component calcium_concentration (micromolar). * RATES[2] is d/dt Cass in component calcium_concentration (micromolar). * RATES[3] is d/dt CaJSR in component calcium_concentration (micromolar). * RATES[4] is d/dt CaNSR in component calcium_concentration (micromolar). * RATES[5] is d/dt P_RyR in component calcium_fluxes (dimensionless). * RATES[6] is d/dt LTRPN_Ca in component calcium_buffering (micromolar). * RATES[7] is d/dt HTRPN_Ca in component calcium_buffering (micromolar). * RATES[8] is d/dt P_O1 in component ryanodine_receptors (dimensionless). * RATES[9] is d/dt P_O2 in component ryanodine_receptors (dimensionless). * RATES[10] is d/dt P_C2 in component ryanodine_receptors (dimensionless). * RATES[11] is d/dt O in component L_type_calcium_current (dimensionless). * RATES[12] is d/dt C2 in component L_type_calcium_current (dimensionless). * RATES[13] is d/dt C3 in component L_type_calcium_current (dimensionless). * RATES[14] is d/dt C4 in component L_type_calcium_current (dimensionless). * RATES[15] is d/dt I1 in component L_type_calcium_current (dimensionless). * RATES[16] is d/dt I2 in component L_type_calcium_current (dimensionless). * RATES[17] is d/dt I3 in component L_type_calcium_current (dimensionless). * RATES[18] is d/dt Nai in component sodium_concentration (micromolar). * RATES[21] is d/dt C_Na2 in component fast_sodium_current (dimensionless). * RATES[20] is d/dt C_Na1 in component fast_sodium_current (dimensionless). * RATES[19] is d/dt O_Na in component fast_sodium_current (dimensionless). * RATES[24] is d/dt IF_Na in component fast_sodium_current (dimensionless). * RATES[22] is d/dt I1_Na in component fast_sodium_current (dimensionless). * RATES[23] is d/dt I2_Na in component fast_sodium_current (dimensionless). * RATES[25] is d/dt IC_Na2 in component fast_sodium_current (dimensionless). * RATES[26] is d/dt IC_Na3 in component fast_sodium_current (dimensionless). * RATES[27] is d/dt Ki in component potassium_concentration (micromolar). * RATES[28] is d/dt ato_f in component fast_transient_outward_potassium_current (dimensionless). * RATES[29] is d/dt ito_f in component fast_transient_outward_potassium_current (dimensionless). * RATES[30] is d/dt ato_s in component slow_transient_outward_potassium_current (dimensionless). * RATES[31] is d/dt ito_s in component slow_transient_outward_potassium_current (dimensionless). * RATES[32] is d/dt nKs in component slow_delayed_rectifier_potassium_current (dimensionless). * RATES[33] is d/dt aur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless). * RATES[34] is d/dt iur in component ultra_rapidly_activating_delayed_rectifier_potassium_current (dimensionless). * RATES[35] is d/dt aKss in component non_inactivating_steady_state_potassium_current (dimensionless). * RATES[36] is d/dt iKss in component non_inactivating_steady_state_potassium_current (dimensionless). * RATES[39] is d/dt C_K2 in component rapid_delayed_rectifier_potassium_current (dimensionless). * RATES[38] is d/dt C_K1 in component rapid_delayed_rectifier_potassium_current (dimensionless). * RATES[37] is d/dt O_K in component rapid_delayed_rectifier_potassium_current (dimensionless). * RATES[40] is d/dt I_K in component rapid_delayed_rectifier_potassium_current (dimensionless). * There are a total of 3 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = -82.4202; CONSTANTS[0] = 1; CONSTANTS[1] = 25.84e-6; CONSTANTS[2] = 0.12e-6; CONSTANTS[3] = 2.098e-6; CONSTANTS[4] = 1.485e-9; CONSTANTS[5] = 1.534e-4; CONSTANTS[6] = 5400; CONSTANTS[7] = 140000; CONSTANTS[8] = 1800; CONSTANTS[9] = 8.314; CONSTANTS[10] = 298; CONSTANTS[11] = 96.5; CONSTANTS[12] = 20; CONSTANTS[13] = 100000; CONSTANTS[14] = 71.43; CONSTANTS[15] = 0.5; CONSTANTS[16] = -80; STATES[1] = 0.115001; STATES[2] = 0.115001; STATES[3] = 1299.5; STATES[4] = 1299.5; CONSTANTS[17] = 50; CONSTANTS[18] = 15000; CONSTANTS[19] = 0.238; CONSTANTS[20] = 800; CONSTANTS[21] = 0.00237; CONSTANTS[22] = 3.2e-5; CONSTANTS[23] = 0.0327; CONSTANTS[24] = 0.0196; STATES[5] = 0; CONSTANTS[25] = 4.5; CONSTANTS[26] = 20; CONSTANTS[27] = 1.74e-5; CONSTANTS[28] = 8; CONSTANTS[29] = 0.45; CONSTANTS[30] = 0.5; CONSTANTS[31] = 70; CONSTANTS[32] = 140; STATES[6] = 11.2684; STATES[7] = 125.29; CONSTANTS[33] = 7; STATES[8] = 0.149102e-4; STATES[9] = 0.951726e-10; STATES[10] = 0.16774e-3; CONSTANTS[34] = 0.006075; CONSTANTS[35] = 0.07125; CONSTANTS[36] = 0.00405; CONSTANTS[37] = 0.965; CONSTANTS[38] = 0.009; CONSTANTS[39] = 0.0008; CONSTANTS[40] = 3; CONSTANTS[41] = 4; CONSTANTS[42] = 63; CONSTANTS[43] = 0.1729; STATES[11] = 0.930308e-18; STATES[12] = 0.124216e-3; STATES[13] = 0.578679e-8; STATES[14] = 0.119816e-12; STATES[15] = 0.497923e-18; STATES[16] = 0.345847e-13; STATES[17] = 0.185106e-13; CONSTANTS[44] = 0.0005; CONSTANTS[45] = 0.23324; CONSTANTS[46] = 20; CONSTANTS[47] = 1; CONSTANTS[48] = 0.5; CONSTANTS[49] = 292.8; CONSTANTS[50] = 87500; CONSTANTS[51] = 1380; CONSTANTS[52] = 0.1; CONSTANTS[53] = 0.35; STATES[18] = 14237.1; CONSTANTS[54] = 0.000367; CONSTANTS[55] = 13; STATES[19] = 0.713483e-6; STATES[20] = 0.279132e-3; STATES[21] = 0.020752; STATES[22] = 0.673345e-6; STATES[23] = 0.155787e-8; STATES[24] = 0.153176e-3; STATES[25] = 0.0113879; STATES[26] = 0.34278; STATES[27] = 143720; CONSTANTS[56] = 0.0026; CONSTANTS[57] = 0.0798; STATES[28] = 0.265563e-2; STATES[29] = 0.999977; CONSTANTS[58] = 0.0629; STATES[30] = 0.417069e-3; STATES[31] = 0.998543; CONSTANTS[59] = 0.00575; STATES[32] = 0.262753e-3; CONSTANTS[60] = 0.0975; STATES[33] = 0.417069e-3; STATES[34] = 0.998543; CONSTANTS[61] = 0.0324; STATES[35] = 0.417069e-3; STATES[36] = 1; CONSTANTS[62] = 0.078; STATES[37] = 0.175298e-3; STATES[38] = 0.992513e-3; STATES[39] = 0.641229e-3; STATES[40] = 0.319129e-4; CONSTANTS[63] = 0.036778; CONSTANTS[64] = 0.023761; CONSTANTS[65] = 0.88; CONSTANTS[66] = 21000; CONSTANTS[67] = 1500; CONSTANTS[68] = 10; CONSTANTS[69] = -40; CONSTANTS[70] = 10; CONSTANTS[71] = (1.00000/7.00000)*(exp(CONSTANTS[7]/67300.0) - 1.00000); CONSTANTS[72] = 0.00000; RATES[0] = 0.1001; RATES[1] = 0.1001; RATES[2] = 0.1001; RATES[3] = 0.1001; RATES[4] = 0.1001; RATES[5] = 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[13] = 0.1001; RATES[14] = 0.1001; RATES[15] = 0.1001; RATES[16] = 0.1001; RATES[17] = 0.1001; RATES[18] = 0.1001; RATES[21] = 0.1001; RATES[20] = 0.1001; RATES[19] = 0.1001; RATES[24] = 0.1001; RATES[22] = 0.1001; RATES[23] = 0.1001; RATES[25] = 0.1001; RATES[26] = 0.1001; RATES[27] = 0.1001; RATES[28] = 0.1001; RATES[29] = 0.1001; RATES[30] = 0.1001; RATES[31] = 0.1001; RATES[32] = 0.1001; RATES[33] = 0.1001; RATES[34] = 0.1001; RATES[35] = 0.1001; RATES[39] = 0.1001; RATES[38] = 0.1001; RATES[37] = 0.1001; RATES[40] = 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[11]+ALGEBRAIC[17]+ALGEBRAIC[18]+ALGEBRAIC[66]+ALGEBRAIC[67]+ALGEBRAIC[36]+ALGEBRAIC[69]+ALGEBRAIC[68]+ALGEBRAIC[42]+ALGEBRAIC[47]+ALGEBRAIC[48]+ALGEBRAIC[51]+ALGEBRAIC[54]+ALGEBRAIC[56]+ALGEBRAIC[70]+ALGEBRAIC[0]); resid[1] = RATES[1] - ALGEBRAIC[1]*((ALGEBRAIC[7]+ALGEBRAIC[6]) - (ALGEBRAIC[8]+ALGEBRAIC[9]+( ((ALGEBRAIC[66]+ALGEBRAIC[17]) - 2.00000*ALGEBRAIC[18])*CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[1]*CONSTANTS[11]))); resid[2] = RATES[2] - ALGEBRAIC[2]*(( ALGEBRAIC[4]*CONSTANTS[2])/CONSTANTS[4] - (( ALGEBRAIC[6]*CONSTANTS[1])/CONSTANTS[4]+( ALGEBRAIC[11]*CONSTANTS[5]*CONSTANTS[0])/( 2.00000*CONSTANTS[4]*CONSTANTS[11]))); resid[3] = RATES[3] - ALGEBRAIC[3]*(ALGEBRAIC[5] - ALGEBRAIC[4]); resid[4] = RATES[4] - ( (ALGEBRAIC[8] - ALGEBRAIC[7])*CONSTANTS[1])/CONSTANTS[3] - ( ALGEBRAIC[5]*CONSTANTS[2])/CONSTANTS[3]; resid[5] = RATES[5] - - 0.0400000*STATES[5] - (( 0.100000*ALGEBRAIC[11])/CONSTANTS[33])*exp(- pow(STATES[0] - 5.00000, 2.00000)/648.000); resid[6] = RATES[6] - CONSTANTS[23]*STATES[1]*(CONSTANTS[31] - STATES[6]) - CONSTANTS[24]*STATES[6]; resid[7] = RATES[7] - CONSTANTS[21]*STATES[1]*(CONSTANTS[32] - STATES[7]) - CONSTANTS[22]*STATES[7]; resid[8] = RATES[8] - ( CONSTANTS[34]*pow(STATES[2], CONSTANTS[41])*ALGEBRAIC[10]+ CONSTANTS[37]*STATES[9]+ CONSTANTS[39]*STATES[10]) - ( CONSTANTS[35]*STATES[8]+ CONSTANTS[36]*pow(STATES[2], CONSTANTS[40])*STATES[8]+ CONSTANTS[38]*STATES[8]); resid[9] = RATES[9] - CONSTANTS[36]*pow(STATES[2], CONSTANTS[40])*STATES[8] - CONSTANTS[37]*STATES[9]; resid[10] = RATES[10] - CONSTANTS[38]*STATES[8] - CONSTANTS[39]*STATES[10]; resid[11] = RATES[11] - ( ALGEBRAIC[13]*STATES[14]+ CONSTANTS[44]*STATES[15]+ 0.00100000*( ALGEBRAIC[13]*STATES[16] - ALGEBRAIC[16]*STATES[11])) - ( 4.00000*ALGEBRAIC[14]*STATES[11]+ ALGEBRAIC[15]*STATES[11]); resid[12] = RATES[12] - ( 4.00000*ALGEBRAIC[13]*ALGEBRAIC[12]+ 2.00000*ALGEBRAIC[14]*STATES[13]) - ( ALGEBRAIC[14]*STATES[12]+ 3.00000*ALGEBRAIC[13]*STATES[12]); resid[13] = RATES[13] - ( 3.00000*ALGEBRAIC[13]*STATES[12]+ 3.00000*ALGEBRAIC[14]*STATES[14]) - ( 2.00000*ALGEBRAIC[14]*STATES[13]+ 2.00000*ALGEBRAIC[13]*STATES[13]); resid[14] = RATES[14] - ( 2.00000*ALGEBRAIC[13]*STATES[13]+ 4.00000*ALGEBRAIC[14]*STATES[11]+ 0.0100000*( 4.00000*CONSTANTS[44]*ALGEBRAIC[14]*STATES[15] - ALGEBRAIC[13]*ALGEBRAIC[15]*STATES[14])+ 0.00200000*( 4.00000*ALGEBRAIC[14]*STATES[16] - ALGEBRAIC[16]*STATES[14])+ 4.00000*ALGEBRAIC[14]*CONSTANTS[44]*STATES[17]) - ( 3.00000*ALGEBRAIC[14]*STATES[14]+ ALGEBRAIC[13]*STATES[14]+ 1.00000*ALGEBRAIC[15]*ALGEBRAIC[16]*STATES[14]); resid[15] = RATES[15] - ( ALGEBRAIC[15]*STATES[11]+ 0.00100000*( ALGEBRAIC[13]*STATES[17] - ALGEBRAIC[16]*STATES[15])+ 0.0100000*( ALGEBRAIC[13]*ALGEBRAIC[15]*STATES[14] - 4.00000*ALGEBRAIC[14]*ALGEBRAIC[16]*STATES[15])) - CONSTANTS[44]*STATES[15]; resid[16] = RATES[16] - ( 0.00100000*( ALGEBRAIC[16]*STATES[11] - ALGEBRAIC[13]*STATES[16])+ CONSTANTS[44]*STATES[17]+ 0.00200000*( ALGEBRAIC[16]*STATES[14] - 4.00000*ALGEBRAIC[14]*STATES[16])) - ALGEBRAIC[15]*STATES[16]; resid[17] = RATES[17] - ( 0.00100000*( ALGEBRAIC[16]*STATES[15] - ALGEBRAIC[13]*STATES[17])+ ALGEBRAIC[15]*STATES[16]+ 1.00000*ALGEBRAIC[15]*ALGEBRAIC[16]*STATES[14]) - ( 4.00000*ALGEBRAIC[14]*CONSTANTS[44]*STATES[17]+ CONSTANTS[44]*STATES[17]); resid[18] = RATES[18] - ( - (ALGEBRAIC[67]+ALGEBRAIC[36]+ 3.00000*ALGEBRAIC[69]+ 3.00000*ALGEBRAIC[18])*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[11]); resid[19] = RATES[21] - ( ALGEBRAIC[22]*ALGEBRAIC[21]+ ALGEBRAIC[26]*STATES[20]+ ALGEBRAIC[28]*STATES[25]) - ( ALGEBRAIC[25]*STATES[21]+ ALGEBRAIC[23]*STATES[21]+ ALGEBRAIC[29]*STATES[21]); resid[20] = RATES[20] - ( ALGEBRAIC[23]*STATES[21]+ ALGEBRAIC[27]*STATES[19]+ ALGEBRAIC[28]*STATES[24]) - ( ALGEBRAIC[26]*STATES[20]+ ALGEBRAIC[24]*STATES[20]+ ALGEBRAIC[29]*STATES[20]); resid[21] = RATES[19] - ( ALGEBRAIC[24]*STATES[20]+ ALGEBRAIC[31]*STATES[24]) - ( ALGEBRAIC[27]*STATES[19]+ ALGEBRAIC[30]*STATES[19]); resid[22] = RATES[24] - ( ALGEBRAIC[30]*STATES[19]+ ALGEBRAIC[29]*STATES[20]+ ALGEBRAIC[33]*STATES[22]+ ALGEBRAIC[23]*STATES[25]) - ( ALGEBRAIC[31]*STATES[24]+ ALGEBRAIC[28]*STATES[24]+ ALGEBRAIC[32]*STATES[24]+ ALGEBRAIC[26]*STATES[24]); resid[23] = RATES[22] - ( ALGEBRAIC[32]*STATES[24]+ ALGEBRAIC[35]*STATES[23]) - ( ALGEBRAIC[33]*STATES[22]+ ALGEBRAIC[34]*STATES[22]); resid[24] = RATES[23] - ALGEBRAIC[34]*STATES[22] - ALGEBRAIC[35]*STATES[23]; resid[25] = RATES[25] - ( ALGEBRAIC[22]*STATES[26]+ ALGEBRAIC[26]*STATES[24]+ ALGEBRAIC[29]*STATES[21]) - ( ALGEBRAIC[25]*STATES[25]+ ALGEBRAIC[23]*STATES[25]+ ALGEBRAIC[28]*STATES[25]); resid[26] = RATES[26] - ( ALGEBRAIC[25]*STATES[25]+ ALGEBRAIC[29]*ALGEBRAIC[21]) - ( ALGEBRAIC[22]*STATES[26]+ ALGEBRAIC[28]*STATES[26]); resid[27] = RATES[27] - ( - ((ALGEBRAIC[68]+ALGEBRAIC[42]+ALGEBRAIC[47]+ALGEBRAIC[48]+ALGEBRAIC[54]+ALGEBRAIC[51]+ALGEBRAIC[56]) - 2.00000*ALGEBRAIC[69])*CONSTANTS[5]*CONSTANTS[0])/( CONSTANTS[1]*CONSTANTS[11]); resid[28] = RATES[28] - ALGEBRAIC[38]*(1.00000 - STATES[28]) - ALGEBRAIC[39]*STATES[28]; resid[29] = RATES[29] - ALGEBRAIC[40]*(1.00000 - STATES[29]) - ALGEBRAIC[41]*STATES[29]; resid[30] = RATES[30] - (ALGEBRAIC[43] - STATES[30])/ALGEBRAIC[45]; resid[31] = RATES[31] - (ALGEBRAIC[44] - STATES[31])/ALGEBRAIC[46]; resid[32] = RATES[32] - ALGEBRAIC[49]*(1.00000 - STATES[32]) - ALGEBRAIC[50]*STATES[32]; resid[33] = RATES[33] - (ALGEBRAIC[43] - STATES[33])/ALGEBRAIC[52]; resid[34] = RATES[34] - (ALGEBRAIC[44] - STATES[34])/ALGEBRAIC[53]; resid[35] = RATES[35] - (ALGEBRAIC[43] - STATES[35])/ALGEBRAIC[55]; resid[36] = RATES[39] - ( CONSTANTS[64]*STATES[38]+ ALGEBRAIC[61]*STATES[37]) - ( CONSTANTS[63]*STATES[39]+ ALGEBRAIC[60]*STATES[39]); resid[37] = RATES[38] - ( ALGEBRAIC[58]*ALGEBRAIC[57]+ CONSTANTS[63]*STATES[39]) - ( ALGEBRAIC[59]*STATES[38]+ CONSTANTS[64]*STATES[38]); resid[38] = RATES[37] - ( ALGEBRAIC[60]*STATES[39]+ ALGEBRAIC[63]*STATES[40]) - ( ALGEBRAIC[61]*STATES[37]+ ALGEBRAIC[62]*STATES[37]); resid[39] = RATES[40] - ALGEBRAIC[62]*STATES[37] - ALGEBRAIC[63]*STATES[40]; resid[40] = RATES[36] - CONSTANTS[72]; } 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[16] : 0.00000); ALGEBRAIC[1] = pow(1.00000+( CONSTANTS[17]*CONSTANTS[19])/pow(CONSTANTS[19]+STATES[1], 2.00000), - 1.00000); ALGEBRAIC[2] = pow(1.00000+( CONSTANTS[17]*CONSTANTS[19])/pow(CONSTANTS[19]+STATES[2], 2.00000), - 1.00000); ALGEBRAIC[3] = pow(1.00000+( CONSTANTS[18]*CONSTANTS[20])/pow(CONSTANTS[20]+STATES[3], 2.00000), - 1.00000); ALGEBRAIC[4] = CONSTANTS[25]*(STATES[8]+STATES[9])*(STATES[3] - STATES[2])*STATES[5]; ALGEBRAIC[5] = (STATES[4] - STATES[3])/CONSTANTS[26]; ALGEBRAIC[6] = (STATES[2] - STATES[1])/CONSTANTS[28]; ALGEBRAIC[7] = CONSTANTS[27]*(STATES[4] - STATES[1]); ALGEBRAIC[8] = ( CONSTANTS[29]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[30], 2.00000)+pow(STATES[1], 2.00000)); ALGEBRAIC[9] = ( CONSTANTS[21]*STATES[1]*(CONSTANTS[32] - STATES[7])+ CONSTANTS[23]*STATES[1]*(CONSTANTS[31] - STATES[6])) - ( CONSTANTS[22]*STATES[7]+ CONSTANTS[24]*STATES[6]); ALGEBRAIC[10] = 1.00000 - (STATES[10]+STATES[8]+STATES[9]); ALGEBRAIC[11] = CONSTANTS[43]*STATES[11]*(STATES[0] - CONSTANTS[42]); ALGEBRAIC[12] = 1.00000 - (STATES[11]+STATES[12]+STATES[13]+STATES[14]+STATES[15]+STATES[16]+STATES[17]); ALGEBRAIC[13] = ( 0.400000*exp((STATES[0]+12.0000)/10.0000)*((1.00000+ 0.700000*exp(- pow(STATES[0]+40.0000, 2.00000)/10.0000)) - 0.750000*exp(- pow(STATES[0]+20.0000, 2.00000)/400.000)))/(1.00000+ 0.120000*exp((STATES[0]+12.0000)/10.0000)); ALGEBRAIC[14] = 0.0500000*exp(- (STATES[0]+12.0000)/13.0000); ALGEBRAIC[15] = ( CONSTANTS[45]*STATES[2])/(CONSTANTS[46]+STATES[2]); ALGEBRAIC[16] = 13.0000*(1.00000 - exp(- pow(STATES[0]+14.5000, 2.00000)/100.000)); ALGEBRAIC[17] = ( CONSTANTS[47]*pow(STATES[1], 2.00000))/(pow(CONSTANTS[48], 2.00000)+pow(STATES[1], 2.00000)); ALGEBRAIC[18] = (( (( (( CONSTANTS[49]*1.00000)/(pow(CONSTANTS[50], 3.00000)+pow(CONSTANTS[7], 3.00000)))*1.00000)/(CONSTANTS[51]+CONSTANTS[8]))*1.00000)/(1.00000+ CONSTANTS[52]*exp(( (CONSTANTS[53] - 1.00000)*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))))*( exp(( CONSTANTS[53]*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))*pow(STATES[18], 3.00000)*CONSTANTS[8] - exp(( (CONSTANTS[53] - 1.00000)*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))*pow(CONSTANTS[7], 3.00000)*STATES[1]); ALGEBRAIC[21] = 1.00000 - (STATES[19]+STATES[20]+STATES[21]+STATES[24]+STATES[22]+STATES[23]+STATES[25]+STATES[26]); ALGEBRAIC[22] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/17.0000)+ 0.200000*exp(- (STATES[0]+2.50000)/150.000)); ALGEBRAIC[23] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/15.0000)+ 0.230000*exp(- (STATES[0]+2.50000)/150.000)); ALGEBRAIC[24] = 3.80200/( 0.102700*exp(- (STATES[0]+2.50000)/12.0000)+ 0.250000*exp(- (STATES[0]+2.50000)/150.000)); ALGEBRAIC[25] = 0.191700*exp(- (STATES[0]+2.50000)/20.3000); ALGEBRAIC[26] = 0.200000*exp(- (STATES[0] - 2.50000)/20.3000); ALGEBRAIC[27] = 0.220000*exp(- (STATES[0] - 7.50000)/20.3000); ALGEBRAIC[28] = 7.00000e-07*exp(- (STATES[0]+7.00000)/7.70000); ALGEBRAIC[29] = 0.00840000+ 2.00000e-05*(STATES[0]+7.00000); ALGEBRAIC[30] = 1.00000/( 0.188495*exp(- (STATES[0]+7.00000)/16.6000)+0.393956); ALGEBRAIC[31] = ( ALGEBRAIC[24]*ALGEBRAIC[30]*ALGEBRAIC[28])/( ALGEBRAIC[27]*ALGEBRAIC[29]); ALGEBRAIC[32] = ALGEBRAIC[30]/1000.00; ALGEBRAIC[33] = ALGEBRAIC[28]; ALGEBRAIC[34] = ALGEBRAIC[30]/95000.0; ALGEBRAIC[35] = ALGEBRAIC[28]/50.0000; ALGEBRAIC[20] = (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(( 0.900000*CONSTANTS[7]+ 0.100000*CONSTANTS[6])/( 0.900000*STATES[18]+ 0.100000*STATES[27])); ALGEBRAIC[36] = CONSTANTS[56]*(STATES[0] - ALGEBRAIC[20]); ALGEBRAIC[38] = 0.180640*exp( 0.0357700*(STATES[0]+30.0000)); ALGEBRAIC[39] = 0.395600*exp( - 0.0623700*(STATES[0]+30.0000)); ALGEBRAIC[40] = ( 0.000152000*exp(- (STATES[0]+13.5000)/7.00000))/( 0.00670830*exp(- (STATES[0]+33.5000)/7.00000)+1.00000); ALGEBRAIC[41] = ( 0.000950000*exp((STATES[0]+33.5000)/7.00000))/( 0.0513350*exp((STATES[0]+33.5000)/7.00000)+1.00000); ALGEBRAIC[37] = (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(CONSTANTS[6]/STATES[27]); ALGEBRAIC[42] = CONSTANTS[58]*STATES[30]*STATES[31]*(STATES[0] - ALGEBRAIC[37]); ALGEBRAIC[43] = 1.00000/(1.00000+exp(- (STATES[0]+22.5000)/7.70000)); ALGEBRAIC[44] = 1.00000/(1.00000+exp((STATES[0]+45.2000)/5.70000)); ALGEBRAIC[45] = 0.493000*exp( - 0.0629000*STATES[0])+2.05800; ALGEBRAIC[46] = 270.000+1050.00/(1.00000+exp((STATES[0]+45.2000)/5.70000)); ALGEBRAIC[47] = ( (( 0.293800*CONSTANTS[6])/(CONSTANTS[6]+210.000))*(STATES[0] - ALGEBRAIC[37]))/(1.00000+exp( 0.0896000*(STATES[0] - ALGEBRAIC[37]))); ALGEBRAIC[48] = CONSTANTS[59]*pow(STATES[32], 2.00000)*(STATES[0] - ALGEBRAIC[37]); ALGEBRAIC[49] = ( 4.81333e-06*(STATES[0]+26.5000))/(1.00000 - exp( - 0.128000*(STATES[0]+26.5000))); ALGEBRAIC[50] = 9.53333e-05*exp( - 0.0380000*(STATES[0]+26.5000)); ALGEBRAIC[51] = CONSTANTS[60]*STATES[33]*STATES[34]*(STATES[0] - ALGEBRAIC[37]); ALGEBRAIC[52] = 0.493000*exp( - 0.0629000*STATES[0])+2.05800; ALGEBRAIC[53] = 1200.00 - 170.000/(1.00000+exp((STATES[0]+45.2000)/5.70000)); ALGEBRAIC[54] = CONSTANTS[61]*STATES[35]*STATES[36]*(STATES[0] - ALGEBRAIC[37]); ALGEBRAIC[55] = 39.3000*exp( - 0.0862000*STATES[0])+13.1700; ALGEBRAIC[56] = CONSTANTS[62]*STATES[37]*(STATES[0] - (( CONSTANTS[9]*CONSTANTS[10])/CONSTANTS[11])*log(( 0.980000*CONSTANTS[6]+ 0.0200000*CONSTANTS[7])/( 0.980000*STATES[27]+ 0.0200000*STATES[18]))); ALGEBRAIC[57] = 1.00000 - (STATES[38]+STATES[39]+STATES[37]+STATES[40]); ALGEBRAIC[58] = 0.0223480*exp( 0.0117600*STATES[0]); ALGEBRAIC[59] = 0.0470020*exp( - 0.0631000*STATES[0]); ALGEBRAIC[60] = 0.0137330*exp( 0.0381980*STATES[0]); ALGEBRAIC[61] = 6.89000e-05*exp( - 0.0417800*STATES[0]); ALGEBRAIC[62] = 0.0908210*exp( 0.0233910*(STATES[0]+5.00000)); ALGEBRAIC[63] = 0.00649700*exp( - 0.0326800*(STATES[0]+5.00000)); ALGEBRAIC[19] = (( CONSTANTS[9]*CONSTANTS[10])/( 2.00000*CONSTANTS[11]))*log(CONSTANTS[8]/STATES[1]); ALGEBRAIC[66] = CONSTANTS[54]*(STATES[0] - ALGEBRAIC[19]); ALGEBRAIC[67] = CONSTANTS[55]*STATES[19]*(STATES[0] - ALGEBRAIC[20]); ALGEBRAIC[68] = CONSTANTS[57]*pow(STATES[28], 3.00000)*STATES[29]*(STATES[0] - ALGEBRAIC[37]); ALGEBRAIC[64] = 1.00000/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))+ 0.0365000*CONSTANTS[71]*exp(( - STATES[0]*CONSTANTS[11])/( CONSTANTS[9]*CONSTANTS[10]))); ALGEBRAIC[69] = ( (( CONSTANTS[65]*ALGEBRAIC[64]*1.00000)/(1.00000+pow(CONSTANTS[66]/STATES[18], 1.50000)))*CONSTANTS[6])/(CONSTANTS[6]+CONSTANTS[67]); ALGEBRAIC[65] = 0.200000/(1.00000+exp(- (STATES[0] - 46.7000)/7.80000)); ALGEBRAIC[70] = (( CONSTANTS[68]*ALGEBRAIC[65]*STATES[1])/(STATES[1]+CONSTANTS[70]))*(STATES[0] - CONSTANTS[69]); } 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; SI[36] = 1.0; SI[37] = 1.0; SI[38] = 1.0; SI[39] = 1.0; SI[40] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = VOI - CONSTANTS[12]; CONDVAR[1] = VOI - CONSTANTS[13]; CONDVAR[2] = ((VOI - CONSTANTS[12]) - floor((VOI - CONSTANTS[12])/CONSTANTS[14])*CONSTANTS[14]) - CONSTANTS[15]; }