/* There are a total of 108 entries in the algebraic variable array. There are a total of 55 entries in each of the rate and state variable arrays. There are a total of 130 entries in the constant variable array. */ /* * VOI is time in component environment (second). * CONSTANTS[0] is R in component model_parameters (joule_per_kilomole_kelvin). * CONSTANTS[1] is T in component model_parameters (kelvin). * CONSTANTS[2] is F in component model_parameters (coulomb_per_mole). * CONSTANTS[3] is Na_e in component model_parameters (millimolar). * CONSTANTS[4] is Ca_e in component model_parameters (millimolar). * CONSTANTS[5] is K_e in component model_parameters (millimolar). * CONSTANTS[6] is ATP_i in component model_parameters (millimolar). * CONSTANTS[127] is Vt in component model_parameters (cm3). * CONSTANTS[84] is Vd in component model_parameters (cm3). * CONSTANTS[86] is Vmyo in component model_parameters (cm3). * CONSTANTS[80] is Sms in component model_parameters (cm2). * CONSTANTS[90] is Smt in component model_parameters (cm2). * CONSTANTS[129] is Cms in component model_parameters (microF). * CONSTANTS[128] is Cmt in component model_parameters (microF). * CONSTANTS[89] is Rst in component model_parameters (ohm). * CONSTANTS[87] is VSRup in component model_parameters (cm3). * CONSTANTS[85] is VSRrel in component model_parameters (cm3). * CONSTANTS[83] is Vc in component model_parameters (cm3). * CONSTANTS[88] is pt in component model_parameters (dimensionless). * CONSTANTS[7] is Rot in component model_parameters (ohm_cm). * CONSTANTS[8] is rt in component model_parameters (cm). * CONSTANTS[9] is Lt in component model_parameters (cm). * CONSTANTS[10] is rc in component model_parameters (cm). * CONSTANTS[11] is Lc in component model_parameters (cm). * CONSTANTS[12] is ptcm in component model_parameters (per_cm2). * ALGEBRAIC[14] is i_circ in component i_circ (microA). * STATES[0] is Vm_s in component Vm_s (millivolt). * STATES[1] is Vm_t in component Vm_t (millivolt). * CONSTANTS[91] is g_Na_s in component membrane_permeabilities (milliS). * CONSTANTS[92] is g_Naps_s in component membrane_permeabilities (milliS). * CONSTANTS[93] is q_Kr_s in component membrane_permeabilities (milliS). * CONSTANTS[94] is q_Ks_s in component membrane_permeabilities (milliS). * CONSTANTS[95] is g_Kp_s in component membrane_permeabilities (milliS). * CONSTANTS[96] is g_Kto_s in component membrane_permeabilities (milliS). * CONSTANTS[97] is g_K1_s in component membrane_permeabilities (milliS). * CONSTANTS[98] is g_KNa_s in component membrane_permeabilities (milliS). * CONSTANTS[99] is g_KATP_s in component membrane_permeabilities (milliS). * CONSTANTS[100] is g_Nab_s in component membrane_permeabilities (milliS). * CONSTANTS[101] is g_Cab_s in component membrane_permeabilities (milliS). * CONSTANTS[102] is i_pCa_max_s in component membrane_permeabilities (microA). * CONSTANTS[103] is i_NaK_max_s in component membrane_permeabilities (microA). * CONSTANTS[104] is i_NaCa_max_s in component membrane_permeabilities (microA). * CONSTANTS[105] is P_CaL_s in component membrane_permeabilities (litre_per_second). * CONSTANTS[106] is P_KL_s in component membrane_permeabilities (litre_per_second). * CONSTANTS[107] is P_nsNa_s in component membrane_permeabilities (litre_per_second). * CONSTANTS[108] is P_nsK_s in component membrane_permeabilities (litre_per_second). * CONSTANTS[109] is g_Na_t in component membrane_permeabilities (milliS). * CONSTANTS[110] is g_Naps_t in component membrane_permeabilities (milliS). * CONSTANTS[111] is q_Kr_t in component membrane_permeabilities (milliS). * CONSTANTS[112] is q_Ks_t in component membrane_permeabilities (milliS). * CONSTANTS[113] is g_Kp_t in component membrane_permeabilities (milliS). * CONSTANTS[114] is g_Kto_t in component membrane_permeabilities (milliS). * CONSTANTS[115] is g_K1_t in component membrane_permeabilities (milliS). * CONSTANTS[116] is g_KNa_t in component membrane_permeabilities (milliS). * CONSTANTS[117] is g_KATP_t in component membrane_permeabilities (milliS). * CONSTANTS[118] is g_Nab_t in component membrane_permeabilities (milliS). * CONSTANTS[119] is g_Cab_t in component membrane_permeabilities (milliS). * CONSTANTS[120] is i_pCa_max_t in component membrane_permeabilities (microA). * CONSTANTS[121] is i_NaK_max_t in component membrane_permeabilities (microA). * CONSTANTS[122] is i_NaCa_max_t in component membrane_permeabilities (microA). * CONSTANTS[123] is P_CaL_t in component membrane_permeabilities (litre_per_second). * CONSTANTS[124] is P_KL_t in component membrane_permeabilities (litre_per_second). * CONSTANTS[125] is P_nsNa_t in component membrane_permeabilities (litre_per_second). * CONSTANTS[126] is P_nsK_t in component membrane_permeabilities (litre_per_second). * CONSTANTS[13] is g_Na in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[14] is fNat in component membrane_permeabilities (dimensionless). * CONSTANTS[15] is g_Naps in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[16] is fNapst in component membrane_permeabilities (dimensionless). * CONSTANTS[17] is q_Kr in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[18] is fKrt in component membrane_permeabilities (dimensionless). * CONSTANTS[19] is q_Ks in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[20] is fKst in component membrane_permeabilities (dimensionless). * CONSTANTS[21] is g_Kp in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[22] is fKpt in component membrane_permeabilities (dimensionless). * CONSTANTS[23] is g_Kto in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[24] is fKtot in component membrane_permeabilities (dimensionless). * CONSTANTS[25] is g_K1 in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[26] is fK1t in component membrane_permeabilities (dimensionless). * CONSTANTS[27] is g_KNa in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[28] is fKNat in component membrane_permeabilities (dimensionless). * CONSTANTS[29] is g_KATP in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[30] is fKATPt in component membrane_permeabilities (dimensionless). * CONSTANTS[31] is g_Nab in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[32] is fNabt in component membrane_permeabilities (dimensionless). * CONSTANTS[33] is g_Cab in component membrane_permeabilities (milliS_per_cm2). * CONSTANTS[34] is fCabt in component membrane_permeabilities (dimensionless). * CONSTANTS[35] is P_CaL in component membrane_permeabilities (litre_per_second_cm2). * CONSTANTS[36] is fCaLt in component membrane_permeabilities (dimensionless). * CONSTANTS[37] is P_KL in component membrane_permeabilities (litre_per_second_cm2). * CONSTANTS[38] is P_nsNa in component membrane_permeabilities (litre_per_second_cm2). * CONSTANTS[39] is fnsNat in component membrane_permeabilities (dimensionless). * CONSTANTS[40] is P_nsK in component membrane_permeabilities (litre_per_second_cm2). * CONSTANTS[41] is fnsKt in component membrane_permeabilities (dimensionless). * CONSTANTS[42] is i_NaCa_max in component membrane_permeabilities (microA_per_cm2). * CONSTANTS[43] is fNaCat in component membrane_permeabilities (dimensionless). * CONSTANTS[44] is i_NaK_max in component membrane_permeabilities (microA_per_cm2). * CONSTANTS[45] is fNaKt in component membrane_permeabilities (dimensionless). * CONSTANTS[46] is i_pCa_max in component membrane_permeabilities (microA_per_cm2). * CONSTANTS[47] is fpCat in component membrane_permeabilities (dimensionless). * ALGEBRAIC[29] is i_Kext in component i_Kext (microA). * CONSTANTS[48] is stim_Period in component i_Kext (second). * ALGEBRAIC[35] is i_Na_s in component i_Na_s (microA). * ALGEBRAIC[32] is E_Na_s in component i_Na_s (millivolt). * STATES[2] is Na_i in component ion_concentrations (millimolar). * STATES[3] is K_i in component ion_concentrations (millimolar). * STATES[4] is m in component i_Na_s_m_gate (dimensionless). * STATES[5] is h in component i_Na_s_h_gate (dimensionless). * ALGEBRAIC[0] is m_infinity in component i_Na_s_m_gate (dimensionless). * ALGEBRAIC[15] is tau_m in component i_Na_s_m_gate (second). * ALGEBRAIC[1] is h_infinity in component i_Na_s_h_gate (dimensionless). * ALGEBRAIC[16] is tau_h in component i_Na_s_h_gate (second). * ALGEBRAIC[39] is i_Na_t in component i_Na_t (microA). * ALGEBRAIC[38] is E_Na_t in component i_Na_t (millivolt). * STATES[6] is Na_t in component ion_concentrations (millimolar). * STATES[7] is K_t in component ion_concentrations (millimolar). * STATES[8] is m in component i_Na_t_m_gate (dimensionless). * STATES[9] is h in component i_Na_t_h_gate (dimensionless). * ALGEBRAIC[2] is m_infinity in component i_Na_t_m_gate (dimensionless). * ALGEBRAIC[17] is tau_m in component i_Na_t_m_gate (second). * ALGEBRAIC[3] is h_infinity in component i_Na_t_h_gate (dimensionless). * ALGEBRAIC[18] is tau_h in component i_Na_t_h_gate (second). * ALGEBRAIC[40] is i_Naps_s in component i_Naps_s (microA). * ALGEBRAIC[41] is i_Naps_t in component i_Naps_t (microA). * ALGEBRAIC[42] is i_CaL_s in component i_CaL_s (microA). * ALGEBRAIC[44] is i_KL_s in component i_CaL_s (microA). * ALGEBRAIC[4] is alfas in component i_CaL_s (per_second). * ALGEBRAIC[19] is betas in component i_CaL_s (per_second). * ALGEBRAIC[30] is gama in component i_CaL_s (per_second). * CONSTANTS[49] is omega in component i_CaL_s (per_second). * CONSTANTS[50] is a in component i_CaL_s (dimensionless). * CONSTANTS[51] is b in component i_CaL_s (dimensionless). * CONSTANTS[52] is f in component i_CaL_s (per_second). * CONSTANTS[53] is g in component i_CaL_s (per_second). * CONSTANTS[54] is f2 in component i_CaL_s (per_second). * CONSTANTS[55] is g2 in component i_CaL_s (per_second). * ALGEBRAIC[33] is alfa2s in component i_CaL_s (per_second). * ALGEBRAIC[36] is beta2s in component i_CaL_s (per_second). * STATES[10] is Cst in component i_CaL_s (dimensionless). * STATES[11] is C1 in component i_CaL_s (dimensionless). * STATES[12] is C2 in component i_CaL_s (dimensionless). * STATES[13] is C3 in component i_CaL_s (dimensionless). * STATES[14] is C4 in component i_CaL_s (dimensionless). * STATES[15] is Co in component i_CaL_s (dimensionless). * STATES[16] is Ccast in component i_CaL_s (dimensionless). * STATES[17] is Cca1 in component i_CaL_s (dimensionless). * STATES[18] is Cca2 in component i_CaL_s (dimensionless). * STATES[19] is Cca3 in component i_CaL_s (dimensionless). * STATES[20] is Cca4 in component i_CaL_s (dimensionless). * STATES[21] is Ccao in component i_CaL_s (dimensionless). * STATES[22] is Ca_ss in component ion_concentrations (millimolar). * ALGEBRAIC[43] is i_CaL_t in component i_CaL_t (microA). * STATES[23] is y in component i_CaL_s_y_gate (dimensionless). * ALGEBRAIC[5] is y_infinity in component i_CaL_s_y_gate (dimensionless). * ALGEBRAIC[20] is tau_y in component i_CaL_s_y_gate (second). * ALGEBRAIC[45] is i_KL_t in component i_CaL_t (microA). * ALGEBRAIC[21] is alfat in component i_CaL_t (per_second). * ALGEBRAIC[31] is betat in component i_CaL_t (per_second). * ALGEBRAIC[6] is gama in component i_CaL_t (per_second). * CONSTANTS[56] is omega in component i_CaL_t (per_second). * CONSTANTS[57] is a in component i_CaL_t (dimensionless). * CONSTANTS[58] is b in component i_CaL_t (dimensionless). * CONSTANTS[59] is f in component i_CaL_t (per_second). * CONSTANTS[60] is g in component i_CaL_t (per_second). * CONSTANTS[61] is f2 in component i_CaL_t (per_second). * CONSTANTS[62] is g2 in component i_CaL_t (per_second). * ALGEBRAIC[34] is alfa2t in component i_CaL_t (per_second). * ALGEBRAIC[37] is beta2t in component i_CaL_t (per_second). * STATES[24] is TCst in component i_CaL_t (dimensionless). * STATES[25] is TC1 in component i_CaL_t (dimensionless). * STATES[26] is TC2 in component i_CaL_t (dimensionless). * STATES[27] is TC3 in component i_CaL_t (dimensionless). * STATES[28] is TC4 in component i_CaL_t (dimensionless). * STATES[29] is TCo in component i_CaL_t (dimensionless). * STATES[30] is TCcast in component i_CaL_t (dimensionless). * STATES[31] is TCca1 in component i_CaL_t (dimensionless). * STATES[32] is TCca2 in component i_CaL_t (dimensionless). * STATES[33] is TCca3 in component i_CaL_t (dimensionless). * STATES[34] is TCca4 in component i_CaL_t (dimensionless). * STATES[35] is TCcao in component i_CaL_t (dimensionless). * STATES[36] is Ca_t in component ion_concentrations (millimolar). * STATES[37] is y in component i_CaL_t_y_gate (dimensionless). * ALGEBRAIC[7] is y_infinity in component i_CaL_t_y_gate (dimensionless). * ALGEBRAIC[22] is tau_y in component i_CaL_t_y_gate (second). * ALGEBRAIC[48] is i_Kr_s in component i_Kr_s (microA). * ALGEBRAIC[46] is E_Kr_s in component i_Kr_s (millivolt). * STATES[38] is xr in component i_Kr_s_xr_gate (dimensionless). * ALGEBRAIC[47] is xri in component i_Kr_s_xri_gate (dimensionless). * ALGEBRAIC[8] is xr_infinity in component i_Kr_s_xr_gate (dimensionless). * ALGEBRAIC[23] is tau_xr in component i_Kr_s_xr_gate (second). * ALGEBRAIC[51] is i_Kr_t in component i_Kr_t (microA). * ALGEBRAIC[49] is E_Kr_t in component i_Kr_t (millivolt). * STATES[39] is xr in component i_Kr_t_xr_gate (dimensionless). * ALGEBRAIC[50] is xri in component i_Kr_t_xri_gate (dimensionless). * ALGEBRAIC[9] is xr_infinity in component i_Kr_t_xr_gate (dimensionless). * ALGEBRAIC[24] is tau_xr in component i_Kr_t_xr_gate (second). * ALGEBRAIC[53] is i_Ks_s in component i_Ks_s (microA). * ALGEBRAIC[52] is E_Ks_s in component i_Ks_s (millivolt). * CONSTANTS[63] is PRNaK in component i_Ks_s (dimensionless). * STATES[40] is Ca_i in component ion_concentrations (millimolar). * STATES[41] is xs in component i_Ks_s_xs_gate (dimensionless). * ALGEBRAIC[10] is xs_infinity in component i_Ks_s_xs_gate (dimensionless). * ALGEBRAIC[25] is tau_xs in component i_Ks_s_xs_gate (second). * ALGEBRAIC[55] is i_Ks_t in component i_Ks_t (microA). * ALGEBRAIC[54] is E_Ks_t in component i_Ks_t (millivolt). * CONSTANTS[64] is PRNaK in component i_Ks_t (dimensionless). * STATES[42] is xs in component i_Ks_t_xs_gate (dimensionless). * ALGEBRAIC[11] is xs_infinity in component i_Ks_t_xs_gate (dimensionless). * ALGEBRAIC[26] is tau_xs in component i_Ks_t_xs_gate (second). * ALGEBRAIC[59] is i_K1_s in component i_K1_s (microA). * ALGEBRAIC[56] is E_K1_s in component i_K1_s (millivolt). * ALGEBRAIC[57] is aK1s in component i_K1_s (dimensionless). * ALGEBRAIC[58] is bK1s in component i_K1_s (dimensionless). * ALGEBRAIC[63] is i_K1_t in component i_K1_t (microA). * ALGEBRAIC[60] is E_K1_t in component i_K1_t (millivolt). * ALGEBRAIC[61] is aK1t in component i_K1_t (dimensionless). * ALGEBRAIC[62] is bK1t in component i_K1_t (dimensionless). * ALGEBRAIC[66] is i_Kp_s in component i_Kp_s (microA). * ALGEBRAIC[65] is kps in component i_Kp_s (dimensionless). * ALGEBRAIC[64] is E_Kp_s in component i_Kp_s (millivolt). * ALGEBRAIC[69] is i_Kp_t in component i_Kp_t (microA). * ALGEBRAIC[68] is kpt in component i_Kp_t (dimensionless). * ALGEBRAIC[67] is E_Kp_t in component i_Kp_t (millivolt). * ALGEBRAIC[71] is i_Kto_s in component i_Kto_s (microA). * CONSTANTS[65] is fr2 in component i_Kto_s (dimensionless). * ALGEBRAIC[70] is E_Kto_s in component i_Kto_s (millivolt). * STATES[43] is r2 in component i_Kto_s_r2_gate (dimensionless). * STATES[44] is r3 in component i_Kto_s_r3_gate (dimensionless). * CONSTANTS[66] is tau_r2 in component i_Kto_s_r2_gate (second). * ALGEBRAIC[12] is r2_infinity in component i_Kto_s_r2_gate (dimensionless). * CONSTANTS[67] is tau_r3 in component i_Kto_s_r3_gate (second). * ALGEBRAIC[27] is r3_infinity in component i_Kto_s_r3_gate (dimensionless). * ALGEBRAIC[73] is i_Kto_t in component i_Kto_t (microA). * CONSTANTS[68] is fr2 in component i_Kto_t (dimensionless). * ALGEBRAIC[72] is E_Kto_t in component i_Kto_t (millivolt). * STATES[45] is r2 in component i_Kto_t_r2_gate (dimensionless). * STATES[46] is r3 in component i_Kto_t_r3_gate (dimensionless). * CONSTANTS[69] is tau_r2 in component i_Kto_t_r2_gate (second). * ALGEBRAIC[13] is r2_infinity in component i_Kto_t_r2_gate (dimensionless). * ALGEBRAIC[28] is r3_infinity in component i_Kto_t_r3_gate (dimensionless). * CONSTANTS[70] is tau_r3 in component i_Kto_t_r3_gate (second). * ALGEBRAIC[75] is i_KNa_s in component i_KNa_s (microA). * ALGEBRAIC[74] is E_KNa_s in component i_KNa_s (millivolt). * ALGEBRAIC[77] is i_KNa_t in component i_KNa_t (microA). * ALGEBRAIC[76] is E_KNa_t in component i_KNa_t (millivolt). * ALGEBRAIC[78] is i_nsNa_s in component i_nsNa_s (microA). * ALGEBRAIC[79] is i_nsNa_t in component i_nsNa_t (microA). * ALGEBRAIC[80] is i_nsK_s in component i_nsK_s (microA). * ALGEBRAIC[81] is i_nsK_t in component i_nsK_t (microA). * ALGEBRAIC[82] is i_Nab_s in component i_Nab_s (microA). * ALGEBRAIC[83] is i_Nab_t in component i_Nab_t (microA). * ALGEBRAIC[85] is i_Cab_s in component i_Cab_s (microA). * ALGEBRAIC[84] is E_Ca_s in component i_Cab_s (millivolt). * ALGEBRAIC[87] is i_Cab_t in component i_Cab_t (microA). * ALGEBRAIC[86] is E_Ca_t in component i_Cab_t (millivolt). * ALGEBRAIC[88] is i_NaCa_s in component i_NaCa_s (microA). * ALGEBRAIC[89] is i_NaCa_t in component i_NaCa_t (microA). * ALGEBRAIC[90] is i_NaK_s in component i_NaK_s (microA). * ALGEBRAIC[91] is i_NaK_t in component i_NaK_t (microA). * ALGEBRAIC[93] is i_pCa_s in component i_pCa_s (microA). * ALGEBRAIC[94] is i_pCa_t in component i_pCa_t (microA). * ALGEBRAIC[98] is i_KATP_s in component i_KATP_s (microA). * ALGEBRAIC[95] is E_KATP_s in component i_KATP_s (millivolt). * ALGEBRAIC[102] is i_KATP_t in component i_KATP_t (microA). * ALGEBRAIC[100] is E_KATP_t in component i_KATP_t (millivolt). * ALGEBRAIC[92] is JteNa in component t_tubular_ion_fluxes (millimolar_per_second). * ALGEBRAIC[96] is JteCa in component t_tubular_ion_fluxes (millimolar_per_second). * ALGEBRAIC[104] is JteK in component t_tubular_ion_fluxes (millimolar_per_second). * CONSTANTS[71] is tau_Na in component t_tubular_ion_fluxes (second). * CONSTANTS[72] is tau_Ca in component t_tubular_ion_fluxes (second). * CONSTANTS[73] is tau_K in component t_tubular_ion_fluxes (second). * ALGEBRAIC[97] is JCaSRup in component JCaSRup (millimolar_per_second). * ALGEBRAIC[99] is JCaSRleak in component JCaSRleak (millimolar_per_second). * STATES[47] is CaSRup in component CaSRup (millimolar). * ALGEBRAIC[101] is Jtr in component Jtr (millimolar_per_second). * CONSTANTS[74] is tau_tr in component Jtr (second). * STATES[48] is CaSRrel in component CaSRrel (millimolar). * ALGEBRAIC[103] is JCaSRrel in component JCaSRrel (millimolar_per_second). * CONSTANTS[81] is kap in component JCaSRrel (per_millimolar4_per_second). * CONSTANTS[75] is kam in component JCaSRrel (per_second). * CONSTANTS[82] is kbp in component JCaSRrel (per_millimolar3_per_second). * CONSTANTS[76] is kbm in component JCaSRrel (per_second). * CONSTANTS[77] is kcp in component JCaSRrel (per_second). * CONSTANTS[78] is kcm in component JCaSRrel (per_second). * STATES[49] is F1 in component JCaSRrel (dimensionless). * STATES[50] is F2 in component JCaSRrel (dimensionless). * STATES[51] is F3 in component JCaSRrel (dimensionless). * STATES[52] is F4 in component JCaSRrel (dimensionless). * ALGEBRAIC[105] is JCad in component JCad (millimolar_per_second). * CONSTANTS[79] is tau_d in component JCad (second). * STATES[53] is BTRH in component ion_concentrations (millimolar). * STATES[54] is BTRL in component ion_concentrations (millimolar). * ALGEBRAIC[106] is dBTRH in component ion_concentrations (millimolar_per_second). * ALGEBRAIC[107] is dBTRL in component ion_concentrations (millimolar_per_second). * RATES[4] is d/dt m in component i_Na_s_m_gate (dimensionless). * RATES[5] is d/dt h in component i_Na_s_h_gate (dimensionless). * RATES[8] is d/dt m in component i_Na_t_m_gate (dimensionless). * RATES[9] is d/dt h in component i_Na_t_h_gate (dimensionless). * RATES[10] is d/dt Cst in component i_CaL_s (dimensionless). * RATES[11] is d/dt C1 in component i_CaL_s (dimensionless). * RATES[12] is d/dt C2 in component i_CaL_s (dimensionless). * RATES[13] is d/dt C3 in component i_CaL_s (dimensionless). * RATES[14] is d/dt C4 in component i_CaL_s (dimensionless). * RATES[15] is d/dt Co in component i_CaL_s (dimensionless). * RATES[16] is d/dt Ccast in component i_CaL_s (dimensionless). * RATES[17] is d/dt Cca1 in component i_CaL_s (dimensionless). * RATES[18] is d/dt Cca2 in component i_CaL_s (dimensionless). * RATES[19] is d/dt Cca3 in component i_CaL_s (dimensionless). * RATES[20] is d/dt Cca4 in component i_CaL_s (dimensionless). * RATES[21] is d/dt Ccao in component i_CaL_s (dimensionless). * RATES[23] is d/dt y in component i_CaL_s_y_gate (dimensionless). * RATES[24] is d/dt TCst in component i_CaL_t (dimensionless). * RATES[25] is d/dt TC1 in component i_CaL_t (dimensionless). * RATES[26] is d/dt TC2 in component i_CaL_t (dimensionless). * RATES[27] is d/dt TC3 in component i_CaL_t (dimensionless). * RATES[28] is d/dt TC4 in component i_CaL_t (dimensionless). * RATES[29] is d/dt TCo in component i_CaL_t (dimensionless). * RATES[30] is d/dt TCcast in component i_CaL_t (dimensionless). * RATES[31] is d/dt TCca1 in component i_CaL_t (dimensionless). * RATES[32] is d/dt TCca2 in component i_CaL_t (dimensionless). * RATES[33] is d/dt TCca3 in component i_CaL_t (dimensionless). * RATES[34] is d/dt TCca4 in component i_CaL_t (dimensionless). * RATES[35] is d/dt TCcao in component i_CaL_t (dimensionless). * RATES[37] is d/dt y in component i_CaL_t_y_gate (dimensionless). * RATES[38] is d/dt xr in component i_Kr_s_xr_gate (dimensionless). * RATES[39] is d/dt xr in component i_Kr_t_xr_gate (dimensionless). * RATES[41] is d/dt xs in component i_Ks_s_xs_gate (dimensionless). * RATES[42] is d/dt xs in component i_Ks_t_xs_gate (dimensionless). * RATES[43] is d/dt r2 in component i_Kto_s_r2_gate (dimensionless). * RATES[44] is d/dt r3 in component i_Kto_s_r3_gate (dimensionless). * RATES[45] is d/dt r2 in component i_Kto_t_r2_gate (dimensionless). * RATES[46] is d/dt r3 in component i_Kto_t_r3_gate (dimensionless). * RATES[49] is d/dt F1 in component JCaSRrel (dimensionless). * RATES[50] is d/dt F2 in component JCaSRrel (dimensionless). * RATES[51] is d/dt F3 in component JCaSRrel (dimensionless). * RATES[52] is d/dt F4 in component JCaSRrel (dimensionless). * RATES[48] is d/dt CaSRrel in component CaSRrel (millimolar). * RATES[47] is d/dt CaSRup in component CaSRup (millimolar). * RATES[22] is d/dt Ca_ss in component ion_concentrations (millimolar). * RATES[40] is d/dt Ca_i in component ion_concentrations (millimolar). * RATES[53] is d/dt BTRH in component ion_concentrations (millimolar). * RATES[54] is d/dt BTRL in component ion_concentrations (millimolar). * RATES[2] is d/dt Na_i in component ion_concentrations (millimolar). * RATES[3] is d/dt K_i in component ion_concentrations (millimolar). * RATES[6] is d/dt Na_t in component ion_concentrations (millimolar). * RATES[36] is d/dt Ca_t in component ion_concentrations (millimolar). * RATES[7] is d/dt K_t in component ion_concentrations (millimolar). * RATES[0] is d/dt Vm_s in component Vm_s (millivolt). * RATES[1] is d/dt Vm_t in component Vm_t (millivolt). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 8310; CONSTANTS[1] = 310; CONSTANTS[2] = 96500; CONSTANTS[3] = 140; CONSTANTS[4] = 1.8; CONSTANTS[5] = 5.4; CONSTANTS[6] = 6.8; CONSTANTS[7] = 83.33; CONSTANTS[8] = 148.15e-7; CONSTANTS[9] = 5.93e-4; CONSTANTS[10] = 5.93e-4; CONSTANTS[11] = 0.013; CONSTANTS[12] = 21028875; STATES[0] = -8.5207812e1; STATES[1] = -8.5208034e1; CONSTANTS[13] = 30; CONSTANTS[14] = 0.64; CONSTANTS[15] = 0.0053; CONSTANTS[16] = 0.526; CONSTANTS[17] = 0.8; CONSTANTS[18] = 0.526; CONSTANTS[19] = 0.1975; CONSTANTS[20] = 0.526; CONSTANTS[21] = 0.006; CONSTANTS[22] = 0.526; CONSTANTS[23] = 0; CONSTANTS[24] = 0.526; CONSTANTS[25] = 0.75; CONSTANTS[26] = 0.8; CONSTANTS[27] = 0.12848; CONSTANTS[28] = 0.526; CONSTANTS[29] = 1; CONSTANTS[30] = 0.526; CONSTANTS[31] = 0.00075; CONSTANTS[32] = 0.526; CONSTANTS[33] = 0.0021; CONSTANTS[34] = 0.526; CONSTANTS[35] = 48e-4; CONSTANTS[36] = 0.64; CONSTANTS[37] = 48e-7; CONSTANTS[38] = 1.75e-7; CONSTANTS[39] = 0.526; CONSTANTS[40] = 1.75e-7; CONSTANTS[41] = 0.526; CONSTANTS[42] = 0.00025; CONSTANTS[43] = 0.526; CONSTANTS[44] = 1.5; CONSTANTS[45] = 0.526; CONSTANTS[46] = 1.15; CONSTANTS[47] = 0.2; CONSTANTS[48] = 1; STATES[2] = 1.1120279e1; STATES[3] = 1.3678926e2; STATES[4] = 3.4130493e-3; STATES[5] = 8.2699973e-1; STATES[6] = 1.3993239e2; STATES[7] = 5.4140321e0; STATES[8] = 3.4129472e-3; STATES[9] = 8.2700551e-1; CONSTANTS[49] = 10; CONSTANTS[50] = 2; CONSTANTS[51] = 2; CONSTANTS[52] = 300; CONSTANTS[53] = 2000; CONSTANTS[54] = 5; CONSTANTS[55] = 7000; STATES[10] = 9.9797984e-1; STATES[11] = 7.5722514e-5; STATES[12] = 2.1545646e-9; STATES[13] = 2.7463178e-14; STATES[14] = -5.7569102e-17; STATES[15] = -1.3496934e-17; STATES[16] = 1.9155065e-3; STATES[17] = 5.8135383e-7; STATES[18] = 6.615513e-11; STATES[19] = 3.3426555e-15; STATES[20] = 4.6640361e-19; STATES[21] = 3.7986638e-22; STATES[22] = 9.5977033e-5; STATES[23] = 7.8134852e-1; CONSTANTS[56] = 10; CONSTANTS[57] = 2; CONSTANTS[58] = 2; CONSTANTS[59] = 300; CONSTANTS[60] = 2000; CONSTANTS[61] = 5; CONSTANTS[62] = 7000; STATES[24] = 9.9798934e-1; STATES[25] = 7.5720257e-5; STATES[26] = 2.1544156e-9; STATES[27] = 2.7459475e-14; STATES[28] = -5.7367189e-17; STATES[29] = -1.3449454e-17; STATES[30] = 1.9155245e-3; STATES[31] = 5.8133643e-7; STATES[32] = 6.6150547e-11; STATES[33] = 3.3423167e-15; STATES[34] = 4.6165955e-19; STATES[35] = 3.75926e-22; STATES[36] = 1.8476402e0; STATES[37] = 7.813477e-1; STATES[38] = 2.0469344e-4; STATES[39] = 2.0469344e-4; CONSTANTS[63] = 0.01833; STATES[40] = 8.8787034e-5; STATES[41] = 6.1359896e-3; CONSTANTS[64] = 0.01833; STATES[42] = 6.1359896e-3; CONSTANTS[65] = 1; STATES[43] = 9.625701e-1; STATES[44] = 9.6250049e-1; CONSTANTS[66] = 0.02; CONSTANTS[67] = 0.08; CONSTANTS[68] = 1; STATES[45] = 9.625701e-1; STATES[46] = 9.6250049e-1; CONSTANTS[69] = 0.02; CONSTANTS[70] = 0.08; CONSTANTS[71] = 0.2; CONSTANTS[72] = 0.24; CONSTANTS[73] = 0.2; STATES[47] = 1.0200296e0; CONSTANTS[74] = 0.03448; STATES[48] = 1.0326252e0; CONSTANTS[75] = 576; CONSTANTS[76] = 1930; CONSTANTS[77] = 18; CONSTANTS[78] = 0.8; STATES[49] = 8.053082e-1; STATES[50] = 2.6639195e-4; STATES[51] = 4.9455459e-10; STATES[52] = 1.9442578e-1; CONSTANTS[79] = 0.003125; STATES[53] = 9.6700747e-1; STATES[54] = 8.1740868e-2; CONSTANTS[80] = 2.00000*3.14159*pow(CONSTANTS[10], 2.00000)+ 2.00000*3.14159*CONSTANTS[10]*CONSTANTS[11]; CONSTANTS[81] = 3.00000*1.21500e+13; CONSTANTS[82] = 4.05000e+09; CONSTANTS[83] = 3.14159*pow(CONSTANTS[10], 2.00000)*CONSTANTS[11]; CONSTANTS[84] = CONSTANTS[83]*0.680000*6.00000e-05; CONSTANTS[85] = CONSTANTS[83]*0.00316000; CONSTANTS[86] = CONSTANTS[83]*0.680000; CONSTANTS[87] = CONSTANTS[83]*0.0550000; CONSTANTS[88] = CONSTANTS[12]*2.00000*3.14159*CONSTANTS[10]*CONSTANTS[11]; CONSTANTS[89] = (( CONSTANTS[7]*CONSTANTS[9])/2.00000)/( 3.14159*pow(CONSTANTS[8], 2.00000)*CONSTANTS[88]); CONSTANTS[90] = 2.00000*3.14159*CONSTANTS[8]*CONSTANTS[9]*CONSTANTS[88]; CONSTANTS[91] = CONSTANTS[13]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[14]); CONSTANTS[92] = CONSTANTS[15]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[16]); CONSTANTS[93] = CONSTANTS[17]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[18]); CONSTANTS[94] = CONSTANTS[19]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[20]); CONSTANTS[95] = CONSTANTS[21]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[22]); CONSTANTS[96] = CONSTANTS[23]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[24]); CONSTANTS[97] = CONSTANTS[25]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[26]); CONSTANTS[98] = CONSTANTS[27]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[28]); CONSTANTS[99] = CONSTANTS[29]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[30]); CONSTANTS[100] = CONSTANTS[31]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[32]); CONSTANTS[101] = CONSTANTS[33]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[34]); CONSTANTS[102] = CONSTANTS[46]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[47]); CONSTANTS[103] = CONSTANTS[44]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[45]); CONSTANTS[104] = CONSTANTS[42]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[43]); CONSTANTS[105] = CONSTANTS[35]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[36]); CONSTANTS[106] = CONSTANTS[37]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[36]); CONSTANTS[107] = CONSTANTS[38]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[39]); CONSTANTS[108] = CONSTANTS[40]*(CONSTANTS[80]+CONSTANTS[90])*(1.00000 - CONSTANTS[41]); CONSTANTS[109] = CONSTANTS[13]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[14]; CONSTANTS[110] = CONSTANTS[15]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[16]; CONSTANTS[111] = CONSTANTS[17]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[18]; CONSTANTS[112] = CONSTANTS[19]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[20]; CONSTANTS[113] = CONSTANTS[21]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[22]; CONSTANTS[114] = CONSTANTS[23]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[24]; CONSTANTS[115] = CONSTANTS[25]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[26]; CONSTANTS[116] = CONSTANTS[27]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[28]; CONSTANTS[117] = CONSTANTS[29]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[30]; CONSTANTS[118] = CONSTANTS[31]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[32]; CONSTANTS[119] = CONSTANTS[33]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[34]; CONSTANTS[120] = CONSTANTS[46]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[47]; CONSTANTS[121] = CONSTANTS[44]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[45]; CONSTANTS[122] = CONSTANTS[42]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[43]; CONSTANTS[123] = CONSTANTS[35]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[36]; CONSTANTS[124] = CONSTANTS[37]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[36]; CONSTANTS[125] = CONSTANTS[38]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[39]; CONSTANTS[126] = CONSTANTS[40]*(CONSTANTS[80]+CONSTANTS[90])*CONSTANTS[41]; CONSTANTS[127] = 3.14159*pow(CONSTANTS[8], 2.00000)*CONSTANTS[9]*CONSTANTS[88]; CONSTANTS[128] = CONSTANTS[90]*1.00000; CONSTANTS[129] = CONSTANTS[80]*1.00000; } void computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { RATES[15] = CONSTANTS[52]*STATES[14] - CONSTANTS[53]*STATES[15]; RATES[21] = CONSTANTS[54]*STATES[20] - CONSTANTS[55]*STATES[21]; RATES[29] = CONSTANTS[59]*STATES[28] - CONSTANTS[60]*STATES[29]; RATES[35] = CONSTANTS[61]*STATES[34] - CONSTANTS[62]*STATES[35]; RATES[49] = STATES[50]*CONSTANTS[75] - STATES[49]*CONSTANTS[81]*pow(STATES[22], 4.00000); RATES[50] = ( STATES[49]*CONSTANTS[81]*pow(STATES[22], 4.00000)+ STATES[51]*CONSTANTS[76]+ STATES[52]*CONSTANTS[78]) - STATES[50]*(CONSTANTS[75]+ CONSTANTS[82]*pow(STATES[22], 3.00000)+CONSTANTS[77]); RATES[51] = STATES[50]*CONSTANTS[82]*pow(STATES[22], 3.00000) - STATES[51]*CONSTANTS[76]; RATES[52] = STATES[50]*CONSTANTS[77] - STATES[52]*CONSTANTS[78]; ALGEBRAIC[12] = 1.00600/(1.00000+exp((STATES[0]+36.6900+10.0000)/12.4300)); RATES[43] = (ALGEBRAIC[12] - STATES[43])/CONSTANTS[66]; ALGEBRAIC[13] = 1.00600/(1.00000+exp((STATES[1]+36.6900+10.0000)/12.4300)); RATES[45] = (ALGEBRAIC[13] - STATES[45])/CONSTANTS[69]; ALGEBRAIC[0] = 1.00000/(1.00000+exp(((STATES[0]+52.2000) - 9.00000)/- 7.40000)); ALGEBRAIC[15] = (0.00100000/( 101.600*exp( STATES[0]*0.113500)+ 0.0226800*exp( - 0.0717000*STATES[0]))+0.000100000)/4.83000; RATES[4] = (ALGEBRAIC[0] - STATES[4])/ALGEBRAIC[15]; ALGEBRAIC[1] = 1.00000/(1.00000+exp(((STATES[0]+85.6000) - 9.00000)/5.50000)); ALGEBRAIC[16] = (0.00100000/( 1.13810e-06*exp( - 0.101700*STATES[0])+ 6.53700*exp( STATES[0]*0.0801600))+0.000500000)/4.83000; RATES[5] = (ALGEBRAIC[1] - STATES[5])/ALGEBRAIC[16]; ALGEBRAIC[2] = 1.00000/(1.00000+exp(((STATES[1]+52.2000) - 9.00000)/- 7.40000)); ALGEBRAIC[17] = (0.00100000/( 101.600*exp( STATES[1]*0.113500)+ 0.0226800*exp( - 0.0717000*STATES[1]))+0.000100000)/4.83000; RATES[8] = (ALGEBRAIC[2] - STATES[8])/ALGEBRAIC[17]; ALGEBRAIC[3] = 1.00000/(1.00000+exp(((STATES[1]+85.6000) - 9.00000)/5.50000)); ALGEBRAIC[18] = (0.00100000/( 1.13810e-06*exp( - 0.101700*STATES[1])+ 6.53700*exp( STATES[1]*0.0801600))+0.000500000)/4.83000; RATES[9] = (ALGEBRAIC[3] - STATES[9])/ALGEBRAIC[18]; ALGEBRAIC[5] = 1.00000/(1.00000+exp((STATES[0]+35.0000)/6.00000)); ALGEBRAIC[20] = 0.00100000/(0.0200000+ 0.0197000*exp(- pow( (STATES[0]+10.0000)*0.0337000, 2.00000)))+0.550000/(1.00000+exp( ((STATES[0]+40.0000)/9.50000)*4.00000)); RATES[23] = (ALGEBRAIC[5] - STATES[23])/ALGEBRAIC[20]; ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[1]+35.0000)/6.00000)); ALGEBRAIC[22] = 0.00100000/(0.0200000+ 0.0197000*exp(- pow( (STATES[1]+10.0000)*0.0337000, 2.00000)))+0.550000/(1.00000+exp( ((STATES[1]+40.0000)/9.50000)*4.00000)); RATES[37] = (ALGEBRAIC[7] - STATES[37])/ALGEBRAIC[22]; ALGEBRAIC[8] = 1.00000/(1.00000+exp(- (STATES[0]+21.5000)/7.50000)); ALGEBRAIC[23] = 0.00100000/(( 0.00138000*(STATES[0]+14.2000))/(1.00000 - exp( - 0.123000*(STATES[0]+14.2000)))+( 0.000610000*(STATES[0]+38.9000))/(exp( 0.145000*(STATES[0]+38.9000)) - 1.00000)); RATES[38] = (ALGEBRAIC[8] - STATES[38])/ALGEBRAIC[23]; ALGEBRAIC[9] = 1.00000/(1.00000+exp(- (STATES[1]+21.5000)/7.50000)); ALGEBRAIC[24] = 0.00100000/(( 0.00138000*(STATES[1]+14.2000))/(1.00000 - exp( - 0.123000*(STATES[1]+14.2000)))+( 0.000610000*(STATES[1]+38.9000))/(exp( 0.145000*(STATES[1]+38.9000)) - 1.00000)); RATES[39] = (ALGEBRAIC[9] - STATES[39])/ALGEBRAIC[24]; ALGEBRAIC[10] = 1.00000/(1.00000+exp(- (STATES[0] - 1.50000)/16.7000)); ALGEBRAIC[25] = 0.00100000/(( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)))+( 0.000131000*(STATES[0]+30.0000))/(exp( 0.0687000*(STATES[0]+30.0000)) - 1.00000)); RATES[41] = (ALGEBRAIC[10] - STATES[41])/ALGEBRAIC[25]; ALGEBRAIC[11] = 1.00000/(1.00000+exp(- (STATES[1] - 1.50000)/16.7000)); ALGEBRAIC[26] = 0.00100000/(( 7.19000e-05*(STATES[1]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[1]+30.0000)))+( 0.000131000*(STATES[1]+30.0000))/(exp( 0.0687000*(STATES[1]+30.0000)) - 1.00000)); RATES[42] = (ALGEBRAIC[11] - STATES[42])/ALGEBRAIC[26]; ALGEBRAIC[27] = ALGEBRAIC[12]; RATES[44] = (ALGEBRAIC[27] - STATES[44])/CONSTANTS[67]; ALGEBRAIC[28] = ALGEBRAIC[13]; RATES[46] = (ALGEBRAIC[28] - STATES[46])/CONSTANTS[70]; ALGEBRAIC[4] = 400.000*exp((STATES[0]+12.0000)/10.0000); ALGEBRAIC[19] = 50.0000*exp(- (STATES[0]+12.0000)/13.0000); ALGEBRAIC[30] = 187.500*STATES[22]; RATES[10] = ( ALGEBRAIC[19]*STATES[11]+ CONSTANTS[49]*STATES[16]) - ( 4.00000*ALGEBRAIC[4]+ALGEBRAIC[30])*STATES[10]; RATES[11] = ( 4.00000*ALGEBRAIC[4]*STATES[10]+ 2.00000*ALGEBRAIC[19]*STATES[12]+ (CONSTANTS[49]/CONSTANTS[51])*STATES[17]) - (ALGEBRAIC[19]+ 3.00000*ALGEBRAIC[4]+ ALGEBRAIC[30]*CONSTANTS[50])*STATES[11]; RATES[12] = ( 3.00000*ALGEBRAIC[4]*STATES[11]+ 3.00000*ALGEBRAIC[19]*STATES[13]+ (CONSTANTS[49]/pow(CONSTANTS[51], 2.00000))*STATES[18]) - ( 2.00000*ALGEBRAIC[19]+ 2.00000*ALGEBRAIC[4]+ ALGEBRAIC[30]*pow(CONSTANTS[50], 2.00000))*STATES[12]; RATES[13] = ( 2.00000*ALGEBRAIC[4]*STATES[12]+ 4.00000*ALGEBRAIC[19]*STATES[14]+ (CONSTANTS[49]/pow(CONSTANTS[51], 3.00000))*STATES[19]) - ( 3.00000*ALGEBRAIC[19]+ALGEBRAIC[4]+ ALGEBRAIC[30]*pow(CONSTANTS[50], 3.00000))*STATES[13]; RATES[14] = ( ALGEBRAIC[4]*STATES[13]+ CONSTANTS[53]*STATES[15]+ (CONSTANTS[49]/pow(CONSTANTS[51], 4.00000))*STATES[20]) - ( 4.00000*ALGEBRAIC[19]+CONSTANTS[52]+ ALGEBRAIC[30]*pow(CONSTANTS[50], 4.00000))*STATES[14]; ALGEBRAIC[21] = 400.000*exp((STATES[1]+12.0000)/10.0000); ALGEBRAIC[31] = 50.0000*exp(- (STATES[1]+12.0000)/13.0000); ALGEBRAIC[6] = 187.500*STATES[22]; RATES[24] = ( ALGEBRAIC[31]*STATES[25]+ CONSTANTS[56]*STATES[30]) - ( 4.00000*ALGEBRAIC[21]+ALGEBRAIC[6])*STATES[24]; RATES[25] = ( 4.00000*ALGEBRAIC[21]*STATES[24]+ 2.00000*ALGEBRAIC[31]*STATES[26]+ (CONSTANTS[56]/CONSTANTS[58])*STATES[31]) - (ALGEBRAIC[31]+ 3.00000*ALGEBRAIC[21]+ ALGEBRAIC[6]*CONSTANTS[57])*STATES[25]; RATES[26] = ( 3.00000*ALGEBRAIC[21]*STATES[25]+ 3.00000*ALGEBRAIC[31]*STATES[27]+ (CONSTANTS[56]/pow(CONSTANTS[58], 2.00000))*STATES[32]) - ( 2.00000*ALGEBRAIC[31]+ 2.00000*ALGEBRAIC[21]+ ALGEBRAIC[6]*pow(CONSTANTS[57], 2.00000))*STATES[26]; RATES[27] = ( 2.00000*ALGEBRAIC[21]*STATES[26]+ 4.00000*ALGEBRAIC[31]*STATES[28]+ (CONSTANTS[56]/pow(CONSTANTS[58], 3.00000))*STATES[33]) - ( 3.00000*ALGEBRAIC[31]+ALGEBRAIC[21]+ ALGEBRAIC[6]*pow(CONSTANTS[57], 3.00000))*STATES[27]; RATES[28] = ( ALGEBRAIC[21]*STATES[27]+ CONSTANTS[60]*STATES[29]+ (CONSTANTS[56]/pow(CONSTANTS[58], 4.00000))*STATES[34]) - ( 4.00000*ALGEBRAIC[31]+CONSTANTS[59]+ ALGEBRAIC[6]*pow(CONSTANTS[57], 4.00000))*STATES[28]; ALGEBRAIC[33] = ALGEBRAIC[4]*CONSTANTS[50]; ALGEBRAIC[36] = ALGEBRAIC[19]/CONSTANTS[51]; RATES[16] = ( ALGEBRAIC[36]*STATES[17]+ ALGEBRAIC[30]*STATES[10]) - ( 4.00000*ALGEBRAIC[33]+CONSTANTS[49])*STATES[16]; RATES[17] = ( 4.00000*ALGEBRAIC[33]*STATES[16]+ 2.00000*ALGEBRAIC[36]*STATES[18]+ ALGEBRAIC[30]*CONSTANTS[50]*STATES[11]) - (ALGEBRAIC[36]+ 3.00000*ALGEBRAIC[33]+CONSTANTS[49]/CONSTANTS[51])*STATES[17]; RATES[18] = ( 3.00000*ALGEBRAIC[33]*STATES[17]+ 3.00000*ALGEBRAIC[36]*STATES[19]+ ALGEBRAIC[30]*pow(CONSTANTS[50], 2.00000)*STATES[12]) - ( 2.00000*ALGEBRAIC[36]+ 2.00000*ALGEBRAIC[33]+CONSTANTS[49]/pow(CONSTANTS[51], 2.00000))*STATES[18]; RATES[19] = ( 2.00000*ALGEBRAIC[33]*STATES[18]+ 4.00000*ALGEBRAIC[36]*STATES[20]+ ALGEBRAIC[30]*pow(CONSTANTS[50], 3.00000)*STATES[13]) - ( 3.00000*ALGEBRAIC[36]+ALGEBRAIC[33]+CONSTANTS[49]/pow(CONSTANTS[51], 3.00000))*STATES[19]; RATES[20] = ( ALGEBRAIC[33]*STATES[19]+ CONSTANTS[55]*STATES[21]+ ALGEBRAIC[30]*pow(CONSTANTS[50], 4.00000)*STATES[14]) - ( 4.00000*ALGEBRAIC[36]+CONSTANTS[54]+CONSTANTS[49]/pow(CONSTANTS[51], 4.00000))*STATES[20]; ALGEBRAIC[34] = ALGEBRAIC[21]*CONSTANTS[57]; ALGEBRAIC[37] = ALGEBRAIC[31]/CONSTANTS[58]; RATES[30] = ( ALGEBRAIC[37]*STATES[31]+ ALGEBRAIC[6]*STATES[24]) - ( 4.00000*ALGEBRAIC[34]+CONSTANTS[56])*STATES[30]; RATES[31] = ( 4.00000*ALGEBRAIC[34]*STATES[30]+ 2.00000*ALGEBRAIC[37]*STATES[32]+ ALGEBRAIC[6]*CONSTANTS[57]*STATES[25]) - (ALGEBRAIC[37]+ 3.00000*ALGEBRAIC[34]+CONSTANTS[56]/CONSTANTS[58])*STATES[31]; RATES[32] = ( 3.00000*ALGEBRAIC[34]*STATES[31]+ 3.00000*ALGEBRAIC[37]*STATES[33]+ ALGEBRAIC[6]*pow(CONSTANTS[57], 2.00000)*STATES[26]) - ( 2.00000*ALGEBRAIC[37]+ 2.00000*ALGEBRAIC[34]+CONSTANTS[56]/pow(CONSTANTS[58], 2.00000))*STATES[32]; RATES[33] = ( 2.00000*ALGEBRAIC[34]*STATES[32]+ 4.00000*ALGEBRAIC[37]*STATES[34]+ ALGEBRAIC[6]*pow(CONSTANTS[57], 3.00000)*STATES[27]) - ( 3.00000*ALGEBRAIC[37]+ALGEBRAIC[34]+CONSTANTS[56]/pow(CONSTANTS[58], 3.00000))*STATES[33]; RATES[34] = ( ALGEBRAIC[34]*STATES[33]+ CONSTANTS[62]*STATES[35]+ ALGEBRAIC[6]*pow(CONSTANTS[57], 4.00000)*STATES[28]) - ( 4.00000*ALGEBRAIC[37]+CONSTANTS[61]+CONSTANTS[56]/pow(CONSTANTS[58], 4.00000))*STATES[34]; ALGEBRAIC[35] = CONSTANTS[91]*pow(STATES[4], 3.00000)*STATES[5]*(STATES[0] - ( log((CONSTANTS[3]+ 0.120000*CONSTANTS[5])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[39] = CONSTANTS[109]*pow(STATES[8], 3.00000)*STATES[9]*(STATES[1] - ( log((STATES[6]+ 0.120000*STATES[7])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[40] = (CONSTANTS[92]/(1.00000+exp((- 54.0000 - STATES[0])/8.00000)))*(STATES[0] - ( log((CONSTANTS[3]+ 0.120000*CONSTANTS[5])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[41] = (CONSTANTS[110]/(1.00000+exp((- 54.0000 - STATES[1])/8.00000)))*(STATES[1] - ( log((STATES[6]+ 0.120000*STATES[7])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[78] = (( (( 1.00000*CONSTANTS[107]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[2]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*CONSTANTS[3]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[79] = (( (( 1.00000*CONSTANTS[125]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[2]*exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*STATES[6]))/(exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[32] = ( log(CONSTANTS[3]/STATES[2])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[82] = CONSTANTS[100]*(STATES[0] - ALGEBRAIC[32]); ALGEBRAIC[38] = ( log(STATES[6]/STATES[2])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[83] = CONSTANTS[118]*(STATES[1] - ALGEBRAIC[38]); ALGEBRAIC[88] = ( CONSTANTS[104]*exp(( - 0.850000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*CONSTANTS[4] - pow(CONSTANTS[3], 3.00000)*STATES[40]))/(1.00000+ 0.000100000*exp(( - 0.850000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*CONSTANTS[4]+ pow(CONSTANTS[3], 3.00000)*STATES[40])); ALGEBRAIC[89] = ( CONSTANTS[122]*exp(( - 0.850000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*STATES[36] - pow(STATES[6], 3.00000)*STATES[40]))/(1.00000+ 0.000100000*exp(( - 0.850000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*STATES[36]+ pow(STATES[6], 3.00000)*STATES[40])); ALGEBRAIC[90] = ( (( (CONSTANTS[103]/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ (( 0.0365000*1.00000)/7.00000)*(exp(CONSTANTS[3]/67.3000) - 1.00000)*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*1.00000)/(1.00000+pow(10.0000/STATES[2], 1.50000)))*CONSTANTS[5])/(CONSTANTS[5]+1.50000); ALGEBRAIC[91] = ( (( (CONSTANTS[121]/(1.00000+ 0.124500*exp(( - 0.100000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ (( 0.0365000*1.00000)/7.00000)*(exp(STATES[6]/67.3000) - 1.00000)*exp(( - STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*1.00000)/(1.00000+pow(10.0000/STATES[2], 1.50000)))*STATES[7])/(STATES[7]+1.50000); RATES[2] = - (ALGEBRAIC[35]+ALGEBRAIC[39]+ALGEBRAIC[40]+ALGEBRAIC[41]+ALGEBRAIC[78]+ALGEBRAIC[79]+ALGEBRAIC[82]+ALGEBRAIC[83]+ 3.00000*ALGEBRAIC[88]+ 3.00000*ALGEBRAIC[89]+ 3.00000*ALGEBRAIC[90]+ 3.00000*ALGEBRAIC[91])/( CONSTANTS[2]*CONSTANTS[86]); ALGEBRAIC[92] = (( 1.00000*CONSTANTS[127])/CONSTANTS[71])*(STATES[6] - CONSTANTS[3]); RATES[6] = ((ALGEBRAIC[39]+ALGEBRAIC[41]+ALGEBRAIC[79]+ALGEBRAIC[83]+ 3.00000*ALGEBRAIC[89]+ 3.00000*ALGEBRAIC[91])/CONSTANTS[2] - 1.00000*ALGEBRAIC[92])/CONSTANTS[127]; ALGEBRAIC[43] = ( (( 1.00000*CONSTANTS[123]*4.00000*(STATES[29]+STATES[35])*STATES[37]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.00100000*exp(( 2.00000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.341000*STATES[36]))/(exp(( 2.00000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[86] = ( log(STATES[36]/STATES[40])*CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]); ALGEBRAIC[87] = CONSTANTS[119]*(STATES[1] - ALGEBRAIC[86]); ALGEBRAIC[94] = ( CONSTANTS[120]*STATES[40])/(0.000500000+STATES[40]); ALGEBRAIC[96] = (( 1.00000*CONSTANTS[127])/CONSTANTS[72])*(STATES[36] - CONSTANTS[4]); RATES[36] = (( - 2.00000*ALGEBRAIC[89]+ALGEBRAIC[43]+ALGEBRAIC[87]+ALGEBRAIC[94])/( 2.00000*CONSTANTS[2]) - 1.00000*ALGEBRAIC[96])/CONSTANTS[127]; ALGEBRAIC[14] = 1000.00*(STATES[1]/CONSTANTS[89] - STATES[0]/CONSTANTS[89]); ALGEBRAIC[29] = (VOI - floor(VOI/CONSTANTS[48])*CONSTANTS[48]>=0.00000&&VOI - floor(VOI/CONSTANTS[48])*CONSTANTS[48]<=0.00100000 ? 45.0000*(CONSTANTS[80]+CONSTANTS[90]) : 0.00000); ALGEBRAIC[42] = ( (( 1.00000*CONSTANTS[105]*4.00000*(STATES[15]+STATES[21])*STATES[23]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.341000*CONSTANTS[4]))/(exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[44] = ( (( (( 1.00000*CONSTANTS[106])/(1.00000 - (ALGEBRAIC[42]+ALGEBRAIC[43])/( 0.458000*(CONSTANTS[80]+CONSTANTS[90]))))*(STATES[15]+STATES[21])*STATES[23]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( STATES[3]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - CONSTANTS[5]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[46] = ( log(CONSTANTS[5]/STATES[3])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[47] = 1.00000/(1.00000+exp((STATES[0]+9.00000)/22.4000)); ALGEBRAIC[48] = CONSTANTS[93]*0.0261400* pow((CONSTANTS[5]/5.40000), 1.0 / 2)*STATES[38]*ALGEBRAIC[47]*(STATES[0] - ALGEBRAIC[46]); ALGEBRAIC[52] = ( log((CONSTANTS[5]+ CONSTANTS[63]*CONSTANTS[3])/(STATES[3]+ CONSTANTS[63]*STATES[2]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[53] = CONSTANTS[94]*(0.230800+0.769200/(1.00000+exp((- arbitrary_log( 1.00000*STATES[40], 10) - 4.20000)/0.600000)))*pow(STATES[41], 2.00000)*(STATES[0] - ALGEBRAIC[52]); ALGEBRAIC[56] = ALGEBRAIC[46]; ALGEBRAIC[57] = 1020.00/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[56]) - 59.2150))); ALGEBRAIC[58] = ( 1000.00*( 0.491240*exp( 0.0803200*((STATES[0] - ALGEBRAIC[56])+5.47600))+exp( 0.0617500*((STATES[0] - ALGEBRAIC[56]) - 594.310))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[56])+4.75300))); ALGEBRAIC[59] = (( CONSTANTS[97]* pow((CONSTANTS[5]/5.40000), 1.0 / 2)*ALGEBRAIC[57])/(ALGEBRAIC[57]+ALGEBRAIC[58]))*(STATES[0] - ALGEBRAIC[56]); ALGEBRAIC[65] = 1.00000/(1.00000+exp((20.0000 - STATES[0])/5.00000)); ALGEBRAIC[64] = ALGEBRAIC[46]; ALGEBRAIC[66] = CONSTANTS[95]*ALGEBRAIC[65]*(STATES[0] - ALGEBRAIC[64]); ALGEBRAIC[70] = ALGEBRAIC[46]; ALGEBRAIC[71] = (CONSTANTS[96]/( (1.00000+exp((STATES[0]+57.5300)/- 5.86300))*(1.00000+exp((STATES[0] - 45.8000)/25.8700))))*( STATES[43]*CONSTANTS[65]+ STATES[44]*(1.00000 - CONSTANTS[65]))*(STATES[0] - ALGEBRAIC[70]); ALGEBRAIC[74] = ALGEBRAIC[46]; ALGEBRAIC[75] = (( CONSTANTS[98]*(0.800000 - 0.650000/(1.00000+exp((STATES[0]+125.000)/15.0000)))*0.850000)/(1.00000+pow(66.0000/STATES[2], 2.80000)))*(STATES[0] - ALGEBRAIC[74]); ALGEBRAIC[80] = (( (( 1.00000*CONSTANTS[108]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[3]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*CONSTANTS[5]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[84] = ( log(CONSTANTS[4]/STATES[40])*CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]); ALGEBRAIC[85] = CONSTANTS[101]*(STATES[0] - ALGEBRAIC[84]); ALGEBRAIC[93] = ( CONSTANTS[102]*STATES[40])/(0.000500000+STATES[40]); ALGEBRAIC[95] = ALGEBRAIC[46]; ALGEBRAIC[98] = (( CONSTANTS[99]*1.00000)/(1.00000+pow(CONSTANTS[6]/0.114000, 2.00000)))*pow(CONSTANTS[5]/4.00000, 0.240000)*(STATES[0] - ALGEBRAIC[95]); RATES[0] = ( 1000.00*(((((((((((((((((((ALGEBRAIC[29]+ALGEBRAIC[14]) - ALGEBRAIC[35]) - ALGEBRAIC[40]) - ALGEBRAIC[42]) - ALGEBRAIC[44]) - ALGEBRAIC[53]) - ALGEBRAIC[48]) - ALGEBRAIC[59]) - ALGEBRAIC[66]) - ALGEBRAIC[75]) - ALGEBRAIC[98]) - ALGEBRAIC[78]) - ALGEBRAIC[80]) - ALGEBRAIC[82]) - ALGEBRAIC[85]) - ALGEBRAIC[88]) - ALGEBRAIC[90]) - ALGEBRAIC[93]) - ALGEBRAIC[71]))/CONSTANTS[129]; ALGEBRAIC[97] = ( 1.00000*1000.00*CONSTANTS[86]*0.00180000*pow(STATES[40], 2.00000))/(pow(STATES[40], 2.00000)+pow(0.000500000, 2.00000)); ALGEBRAIC[99] = 1250.00*CONSTANTS[86]*5.80000e-05*(STATES[47] - STATES[40]); ALGEBRAIC[101] = (( 1.00000*CONSTANTS[85])/CONSTANTS[74])*(STATES[47] - STATES[48]); RATES[47] = ((ALGEBRAIC[97] - ALGEBRAIC[99]) - ALGEBRAIC[101])/( 1.00000*CONSTANTS[87]); ALGEBRAIC[103] = 1800.00*CONSTANTS[85]*(STATES[50]+STATES[51])*(STATES[48] - STATES[22]); RATES[48] = ( (1.00000/(1.00000+( 15.0000*0.800000)/pow(0.800000+STATES[48], 2.00000)))*(ALGEBRAIC[101] - ALGEBRAIC[103]))/CONSTANTS[85]; ALGEBRAIC[45] = ( (( (( 1.00000*CONSTANTS[124])/(1.00000 - (ALGEBRAIC[42]+ALGEBRAIC[43])/( 0.458000*(CONSTANTS[80]+CONSTANTS[90]))))*(STATES[29]+STATES[35])*STATES[37]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( STATES[3]*exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - STATES[7]))/(exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[49] = ( log(STATES[7]/STATES[3])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[50] = 1.00000/(1.00000+exp((STATES[1]+9.00000)/22.4000)); ALGEBRAIC[51] = CONSTANTS[111]*0.0261400* pow((STATES[7]/5.40000), 1.0 / 2)*STATES[39]*ALGEBRAIC[50]*(STATES[1] - ALGEBRAIC[49]); ALGEBRAIC[54] = ( log((STATES[7]+ CONSTANTS[64]*STATES[6])/(STATES[3]+ CONSTANTS[64]*STATES[2]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[55] = CONSTANTS[112]*(0.230800+0.769200/(1.00000+exp((- arbitrary_log( 1.00000*STATES[40], 10) - 4.20000)/0.600000)))*pow(STATES[42], 2.00000)*(STATES[1] - ALGEBRAIC[54]); ALGEBRAIC[60] = ALGEBRAIC[49]; ALGEBRAIC[61] = 1020.00/(1.00000+exp( 0.238500*((STATES[1] - ALGEBRAIC[60]) - 59.2150))); ALGEBRAIC[62] = ( 1000.00*( 0.491240*exp( 0.0803200*((STATES[1] - ALGEBRAIC[60])+5.47600))+exp( 0.0617500*((STATES[1] - ALGEBRAIC[60]) - 594.310))))/(1.00000+exp( - 0.514300*((STATES[1] - ALGEBRAIC[60])+4.75300))); ALGEBRAIC[63] = (( CONSTANTS[115]* pow((STATES[7]/5.40000), 1.0 / 2)*ALGEBRAIC[61])/(ALGEBRAIC[61]+ALGEBRAIC[62]))*(STATES[1] - ALGEBRAIC[60]); ALGEBRAIC[68] = 1.00000/(1.00000+exp((20.0000 - STATES[1])/5.00000)); ALGEBRAIC[67] = ALGEBRAIC[49]; ALGEBRAIC[69] = CONSTANTS[113]*ALGEBRAIC[68]*(STATES[1] - ALGEBRAIC[67]); ALGEBRAIC[72] = ALGEBRAIC[49]; ALGEBRAIC[73] = (CONSTANTS[114]/( (1.00000+exp((STATES[1]+57.5300)/- 5.86300))*(1.00000+exp((STATES[1] - 45.8000)/25.8700))))*( STATES[45]*CONSTANTS[68]+ STATES[46]*(1.00000 - CONSTANTS[68]))*(STATES[1] - ALGEBRAIC[72]); ALGEBRAIC[76] = ALGEBRAIC[49]; ALGEBRAIC[77] = (( CONSTANTS[116]*(0.800000 - 0.650000/(1.00000+exp((STATES[1]+125.000)/15.0000)))*0.850000)/(1.00000+pow(66.0000/STATES[2], 2.80000)))*(STATES[1] - ALGEBRAIC[76]); ALGEBRAIC[81] = (( (( 1.00000*CONSTANTS[126]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[3]*exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*STATES[7]))/(exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[100] = ALGEBRAIC[49]; ALGEBRAIC[102] = (( CONSTANTS[117]*1.00000)/(1.00000+pow(CONSTANTS[6]/0.114000, 2.00000)))*pow(STATES[7]/4.00000, 0.240000)*(STATES[1] - ALGEBRAIC[100]); RATES[3] = - (((( - 0.00000*ALGEBRAIC[29]+ALGEBRAIC[53]+ALGEBRAIC[55]+ALGEBRAIC[48]+ALGEBRAIC[51]+ALGEBRAIC[59]+ALGEBRAIC[63]+ALGEBRAIC[66]+ALGEBRAIC[69]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[75]+ALGEBRAIC[77]+ALGEBRAIC[80]+ALGEBRAIC[81]) - 2.00000*ALGEBRAIC[90]) - 2.00000*ALGEBRAIC[91])+ALGEBRAIC[98]+ALGEBRAIC[102]+ALGEBRAIC[71]+ALGEBRAIC[73])/( CONSTANTS[2]*CONSTANTS[86]); RATES[1] = ( 1000.00*((((((((((((((((((- ALGEBRAIC[14] - ALGEBRAIC[39]) - ALGEBRAIC[41]) - ALGEBRAIC[43]) - ALGEBRAIC[45]) - ALGEBRAIC[55]) - ALGEBRAIC[51]) - ALGEBRAIC[63]) - ALGEBRAIC[69]) - ALGEBRAIC[77]) - ALGEBRAIC[102]) - ALGEBRAIC[79]) - ALGEBRAIC[81]) - ALGEBRAIC[83]) - ALGEBRAIC[87]) - ALGEBRAIC[89]) - ALGEBRAIC[91]) - ALGEBRAIC[94]) - ALGEBRAIC[73]))/CONSTANTS[128]; ALGEBRAIC[105] = ( 1.00000*(STATES[22] - STATES[40])*CONSTANTS[86])/CONSTANTS[79]; RATES[22] = (1.00000/(1.00000+( 0.0500000*0.00238000)/pow(0.00238000+STATES[22], 2.00000)))*((- (ALGEBRAIC[42]+ALGEBRAIC[43])/( 2.00000*CONSTANTS[2]*CONSTANTS[84])+( 1.00000*ALGEBRAIC[103])/CONSTANTS[84]) - ( 1.00000*ALGEBRAIC[105])/CONSTANTS[84]); ALGEBRAIC[104] = (( 1.00000*CONSTANTS[127])/CONSTANTS[73])*(STATES[7] - CONSTANTS[5]); RATES[7] = ((((ALGEBRAIC[55]+ALGEBRAIC[51]+ALGEBRAIC[63]+ALGEBRAIC[69]+ALGEBRAIC[77]+ALGEBRAIC[45]+ALGEBRAIC[81]) - 2.00000*ALGEBRAIC[91])+ALGEBRAIC[102]+ALGEBRAIC[73])/CONSTANTS[2] - 1.00000*ALGEBRAIC[104])/CONSTANTS[127]; ALGEBRAIC[106] = 20000.0*STATES[40]*(1.00000 - STATES[53]) - 0.0700000*STATES[53]; RATES[53] = ALGEBRAIC[106]; ALGEBRAIC[107] = 40000.0*STATES[40]*(1.00000 - STATES[54]) - STATES[54]*40.0000; RATES[40] = (1.00000/(1.00000+( 0.0500000*0.00238000)/pow(0.00238000+STATES[40], 2.00000)))*(((( - 2.00000*(ALGEBRAIC[88]+ALGEBRAIC[89])+ALGEBRAIC[85]+ALGEBRAIC[87]+ALGEBRAIC[93]+ALGEBRAIC[94])/( - 2.00000*CONSTANTS[2]*CONSTANTS[86])+( 1.00000*((ALGEBRAIC[105]+ALGEBRAIC[99]) - ALGEBRAIC[97]))/CONSTANTS[86]) - ALGEBRAIC[106]*0.140000) - ALGEBRAIC[107]*0.0700000); RATES[54] = ALGEBRAIC[107]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[12] = 1.00600/(1.00000+exp((STATES[0]+36.6900+10.0000)/12.4300)); ALGEBRAIC[13] = 1.00600/(1.00000+exp((STATES[1]+36.6900+10.0000)/12.4300)); ALGEBRAIC[0] = 1.00000/(1.00000+exp(((STATES[0]+52.2000) - 9.00000)/- 7.40000)); ALGEBRAIC[15] = (0.00100000/( 101.600*exp( STATES[0]*0.113500)+ 0.0226800*exp( - 0.0717000*STATES[0]))+0.000100000)/4.83000; ALGEBRAIC[1] = 1.00000/(1.00000+exp(((STATES[0]+85.6000) - 9.00000)/5.50000)); ALGEBRAIC[16] = (0.00100000/( 1.13810e-06*exp( - 0.101700*STATES[0])+ 6.53700*exp( STATES[0]*0.0801600))+0.000500000)/4.83000; ALGEBRAIC[2] = 1.00000/(1.00000+exp(((STATES[1]+52.2000) - 9.00000)/- 7.40000)); ALGEBRAIC[17] = (0.00100000/( 101.600*exp( STATES[1]*0.113500)+ 0.0226800*exp( - 0.0717000*STATES[1]))+0.000100000)/4.83000; ALGEBRAIC[3] = 1.00000/(1.00000+exp(((STATES[1]+85.6000) - 9.00000)/5.50000)); ALGEBRAIC[18] = (0.00100000/( 1.13810e-06*exp( - 0.101700*STATES[1])+ 6.53700*exp( STATES[1]*0.0801600))+0.000500000)/4.83000; ALGEBRAIC[5] = 1.00000/(1.00000+exp((STATES[0]+35.0000)/6.00000)); ALGEBRAIC[20] = 0.00100000/(0.0200000+ 0.0197000*exp(- pow( (STATES[0]+10.0000)*0.0337000, 2.00000)))+0.550000/(1.00000+exp( ((STATES[0]+40.0000)/9.50000)*4.00000)); ALGEBRAIC[7] = 1.00000/(1.00000+exp((STATES[1]+35.0000)/6.00000)); ALGEBRAIC[22] = 0.00100000/(0.0200000+ 0.0197000*exp(- pow( (STATES[1]+10.0000)*0.0337000, 2.00000)))+0.550000/(1.00000+exp( ((STATES[1]+40.0000)/9.50000)*4.00000)); ALGEBRAIC[8] = 1.00000/(1.00000+exp(- (STATES[0]+21.5000)/7.50000)); ALGEBRAIC[23] = 0.00100000/(( 0.00138000*(STATES[0]+14.2000))/(1.00000 - exp( - 0.123000*(STATES[0]+14.2000)))+( 0.000610000*(STATES[0]+38.9000))/(exp( 0.145000*(STATES[0]+38.9000)) - 1.00000)); ALGEBRAIC[9] = 1.00000/(1.00000+exp(- (STATES[1]+21.5000)/7.50000)); ALGEBRAIC[24] = 0.00100000/(( 0.00138000*(STATES[1]+14.2000))/(1.00000 - exp( - 0.123000*(STATES[1]+14.2000)))+( 0.000610000*(STATES[1]+38.9000))/(exp( 0.145000*(STATES[1]+38.9000)) - 1.00000)); ALGEBRAIC[10] = 1.00000/(1.00000+exp(- (STATES[0] - 1.50000)/16.7000)); ALGEBRAIC[25] = 0.00100000/(( 7.19000e-05*(STATES[0]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[0]+30.0000)))+( 0.000131000*(STATES[0]+30.0000))/(exp( 0.0687000*(STATES[0]+30.0000)) - 1.00000)); ALGEBRAIC[11] = 1.00000/(1.00000+exp(- (STATES[1] - 1.50000)/16.7000)); ALGEBRAIC[26] = 0.00100000/(( 7.19000e-05*(STATES[1]+30.0000))/(1.00000 - exp( - 0.148000*(STATES[1]+30.0000)))+( 0.000131000*(STATES[1]+30.0000))/(exp( 0.0687000*(STATES[1]+30.0000)) - 1.00000)); ALGEBRAIC[27] = ALGEBRAIC[12]; ALGEBRAIC[28] = ALGEBRAIC[13]; ALGEBRAIC[4] = 400.000*exp((STATES[0]+12.0000)/10.0000); ALGEBRAIC[19] = 50.0000*exp(- (STATES[0]+12.0000)/13.0000); ALGEBRAIC[30] = 187.500*STATES[22]; ALGEBRAIC[21] = 400.000*exp((STATES[1]+12.0000)/10.0000); ALGEBRAIC[31] = 50.0000*exp(- (STATES[1]+12.0000)/13.0000); ALGEBRAIC[6] = 187.500*STATES[22]; ALGEBRAIC[33] = ALGEBRAIC[4]*CONSTANTS[50]; ALGEBRAIC[36] = ALGEBRAIC[19]/CONSTANTS[51]; ALGEBRAIC[34] = ALGEBRAIC[21]*CONSTANTS[57]; ALGEBRAIC[37] = ALGEBRAIC[31]/CONSTANTS[58]; ALGEBRAIC[35] = CONSTANTS[91]*pow(STATES[4], 3.00000)*STATES[5]*(STATES[0] - ( log((CONSTANTS[3]+ 0.120000*CONSTANTS[5])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[39] = CONSTANTS[109]*pow(STATES[8], 3.00000)*STATES[9]*(STATES[1] - ( log((STATES[6]+ 0.120000*STATES[7])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[40] = (CONSTANTS[92]/(1.00000+exp((- 54.0000 - STATES[0])/8.00000)))*(STATES[0] - ( log((CONSTANTS[3]+ 0.120000*CONSTANTS[5])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[41] = (CONSTANTS[110]/(1.00000+exp((- 54.0000 - STATES[1])/8.00000)))*(STATES[1] - ( log((STATES[6]+ 0.120000*STATES[7])/(STATES[2]+ 0.120000*STATES[3]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]); ALGEBRAIC[78] = (( (( 1.00000*CONSTANTS[107]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[2]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*CONSTANTS[3]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[79] = (( (( 1.00000*CONSTANTS[125]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[2]*exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*STATES[6]))/(exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[32] = ( log(CONSTANTS[3]/STATES[2])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[82] = CONSTANTS[100]*(STATES[0] - ALGEBRAIC[32]); ALGEBRAIC[38] = ( log(STATES[6]/STATES[2])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[83] = CONSTANTS[118]*(STATES[1] - ALGEBRAIC[38]); ALGEBRAIC[88] = ( CONSTANTS[104]*exp(( - 0.850000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*CONSTANTS[4] - pow(CONSTANTS[3], 3.00000)*STATES[40]))/(1.00000+ 0.000100000*exp(( - 0.850000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*CONSTANTS[4]+ pow(CONSTANTS[3], 3.00000)*STATES[40])); ALGEBRAIC[89] = ( CONSTANTS[122]*exp(( - 0.850000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*STATES[36] - pow(STATES[6], 3.00000)*STATES[40]))/(1.00000+ 0.000100000*exp(( - 0.850000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*( exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))*pow(STATES[2], 3.00000)*STATES[36]+ pow(STATES[6], 3.00000)*STATES[40])); ALGEBRAIC[90] = ( (( (CONSTANTS[103]/(1.00000+ 0.124500*exp(( - 0.100000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ (( 0.0365000*1.00000)/7.00000)*(exp(CONSTANTS[3]/67.3000) - 1.00000)*exp(( - STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*1.00000)/(1.00000+pow(10.0000/STATES[2], 1.50000)))*CONSTANTS[5])/(CONSTANTS[5]+1.50000); ALGEBRAIC[91] = ( (( (CONSTANTS[121]/(1.00000+ 0.124500*exp(( - 0.100000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))+ (( 0.0365000*1.00000)/7.00000)*(exp(STATES[6]/67.3000) - 1.00000)*exp(( - STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1]))))*1.00000)/(1.00000+pow(10.0000/STATES[2], 1.50000)))*STATES[7])/(STATES[7]+1.50000); ALGEBRAIC[92] = (( 1.00000*CONSTANTS[127])/CONSTANTS[71])*(STATES[6] - CONSTANTS[3]); ALGEBRAIC[43] = ( (( 1.00000*CONSTANTS[123]*4.00000*(STATES[29]+STATES[35])*STATES[37]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.00100000*exp(( 2.00000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.341000*STATES[36]))/(exp(( 2.00000*STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[86] = ( log(STATES[36]/STATES[40])*CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]); ALGEBRAIC[87] = CONSTANTS[119]*(STATES[1] - ALGEBRAIC[86]); ALGEBRAIC[94] = ( CONSTANTS[120]*STATES[40])/(0.000500000+STATES[40]); ALGEBRAIC[96] = (( 1.00000*CONSTANTS[127])/CONSTANTS[72])*(STATES[36] - CONSTANTS[4]); ALGEBRAIC[14] = 1000.00*(STATES[1]/CONSTANTS[89] - STATES[0]/CONSTANTS[89]); ALGEBRAIC[29] = (VOI - floor(VOI/CONSTANTS[48])*CONSTANTS[48]>=0.00000&&VOI - floor(VOI/CONSTANTS[48])*CONSTANTS[48]<=0.00100000 ? 45.0000*(CONSTANTS[80]+CONSTANTS[90]) : 0.00000); ALGEBRAIC[42] = ( (( 1.00000*CONSTANTS[105]*4.00000*(STATES[15]+STATES[21])*STATES[23]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.00100000*exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.341000*CONSTANTS[4]))/(exp(( 2.00000*STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[44] = ( (( (( 1.00000*CONSTANTS[106])/(1.00000 - (ALGEBRAIC[42]+ALGEBRAIC[43])/( 0.458000*(CONSTANTS[80]+CONSTANTS[90]))))*(STATES[15]+STATES[21])*STATES[23]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( STATES[3]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - CONSTANTS[5]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[46] = ( log(CONSTANTS[5]/STATES[3])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[47] = 1.00000/(1.00000+exp((STATES[0]+9.00000)/22.4000)); ALGEBRAIC[48] = CONSTANTS[93]*0.0261400* pow((CONSTANTS[5]/5.40000), 1.0 / 2)*STATES[38]*ALGEBRAIC[47]*(STATES[0] - ALGEBRAIC[46]); ALGEBRAIC[52] = ( log((CONSTANTS[5]+ CONSTANTS[63]*CONSTANTS[3])/(STATES[3]+ CONSTANTS[63]*STATES[2]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[53] = CONSTANTS[94]*(0.230800+0.769200/(1.00000+exp((- arbitrary_log( 1.00000*STATES[40], 10) - 4.20000)/0.600000)))*pow(STATES[41], 2.00000)*(STATES[0] - ALGEBRAIC[52]); ALGEBRAIC[56] = ALGEBRAIC[46]; ALGEBRAIC[57] = 1020.00/(1.00000+exp( 0.238500*((STATES[0] - ALGEBRAIC[56]) - 59.2150))); ALGEBRAIC[58] = ( 1000.00*( 0.491240*exp( 0.0803200*((STATES[0] - ALGEBRAIC[56])+5.47600))+exp( 0.0617500*((STATES[0] - ALGEBRAIC[56]) - 594.310))))/(1.00000+exp( - 0.514300*((STATES[0] - ALGEBRAIC[56])+4.75300))); ALGEBRAIC[59] = (( CONSTANTS[97]* pow((CONSTANTS[5]/5.40000), 1.0 / 2)*ALGEBRAIC[57])/(ALGEBRAIC[57]+ALGEBRAIC[58]))*(STATES[0] - ALGEBRAIC[56]); ALGEBRAIC[65] = 1.00000/(1.00000+exp((20.0000 - STATES[0])/5.00000)); ALGEBRAIC[64] = ALGEBRAIC[46]; ALGEBRAIC[66] = CONSTANTS[95]*ALGEBRAIC[65]*(STATES[0] - ALGEBRAIC[64]); ALGEBRAIC[70] = ALGEBRAIC[46]; ALGEBRAIC[71] = (CONSTANTS[96]/( (1.00000+exp((STATES[0]+57.5300)/- 5.86300))*(1.00000+exp((STATES[0] - 45.8000)/25.8700))))*( STATES[43]*CONSTANTS[65]+ STATES[44]*(1.00000 - CONSTANTS[65]))*(STATES[0] - ALGEBRAIC[70]); ALGEBRAIC[74] = ALGEBRAIC[46]; ALGEBRAIC[75] = (( CONSTANTS[98]*(0.800000 - 0.650000/(1.00000+exp((STATES[0]+125.000)/15.0000)))*0.850000)/(1.00000+pow(66.0000/STATES[2], 2.80000)))*(STATES[0] - ALGEBRAIC[74]); ALGEBRAIC[80] = (( (( 1.00000*CONSTANTS[108]*STATES[0]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[3]*exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*CONSTANTS[5]))/(exp(( STATES[0]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[84] = ( log(CONSTANTS[4]/STATES[40])*CONSTANTS[0]*CONSTANTS[1])/( 2.00000*CONSTANTS[2]); ALGEBRAIC[85] = CONSTANTS[101]*(STATES[0] - ALGEBRAIC[84]); ALGEBRAIC[93] = ( CONSTANTS[102]*STATES[40])/(0.000500000+STATES[40]); ALGEBRAIC[95] = ALGEBRAIC[46]; ALGEBRAIC[98] = (( CONSTANTS[99]*1.00000)/(1.00000+pow(CONSTANTS[6]/0.114000, 2.00000)))*pow(CONSTANTS[5]/4.00000, 0.240000)*(STATES[0] - ALGEBRAIC[95]); ALGEBRAIC[97] = ( 1.00000*1000.00*CONSTANTS[86]*0.00180000*pow(STATES[40], 2.00000))/(pow(STATES[40], 2.00000)+pow(0.000500000, 2.00000)); ALGEBRAIC[99] = 1250.00*CONSTANTS[86]*5.80000e-05*(STATES[47] - STATES[40]); ALGEBRAIC[101] = (( 1.00000*CONSTANTS[85])/CONSTANTS[74])*(STATES[47] - STATES[48]); ALGEBRAIC[103] = 1800.00*CONSTANTS[85]*(STATES[50]+STATES[51])*(STATES[48] - STATES[22]); ALGEBRAIC[45] = ( (( (( 1.00000*CONSTANTS[124])/(1.00000 - (ALGEBRAIC[42]+ALGEBRAIC[43])/( 0.458000*(CONSTANTS[80]+CONSTANTS[90]))))*(STATES[29]+STATES[35])*STATES[37]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( STATES[3]*exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - STATES[7]))/(exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000); ALGEBRAIC[49] = ( log(STATES[7]/STATES[3])*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[50] = 1.00000/(1.00000+exp((STATES[1]+9.00000)/22.4000)); ALGEBRAIC[51] = CONSTANTS[111]*0.0261400* pow((STATES[7]/5.40000), 1.0 / 2)*STATES[39]*ALGEBRAIC[50]*(STATES[1] - ALGEBRAIC[49]); ALGEBRAIC[54] = ( log((STATES[7]+ CONSTANTS[64]*STATES[6])/(STATES[3]+ CONSTANTS[64]*STATES[2]))*CONSTANTS[0]*CONSTANTS[1])/CONSTANTS[2]; ALGEBRAIC[55] = CONSTANTS[112]*(0.230800+0.769200/(1.00000+exp((- arbitrary_log( 1.00000*STATES[40], 10) - 4.20000)/0.600000)))*pow(STATES[42], 2.00000)*(STATES[1] - ALGEBRAIC[54]); ALGEBRAIC[60] = ALGEBRAIC[49]; ALGEBRAIC[61] = 1020.00/(1.00000+exp( 0.238500*((STATES[1] - ALGEBRAIC[60]) - 59.2150))); ALGEBRAIC[62] = ( 1000.00*( 0.491240*exp( 0.0803200*((STATES[1] - ALGEBRAIC[60])+5.47600))+exp( 0.0617500*((STATES[1] - ALGEBRAIC[60]) - 594.310))))/(1.00000+exp( - 0.514300*((STATES[1] - ALGEBRAIC[60])+4.75300))); ALGEBRAIC[63] = (( CONSTANTS[115]* pow((STATES[7]/5.40000), 1.0 / 2)*ALGEBRAIC[61])/(ALGEBRAIC[61]+ALGEBRAIC[62]))*(STATES[1] - ALGEBRAIC[60]); ALGEBRAIC[68] = 1.00000/(1.00000+exp((20.0000 - STATES[1])/5.00000)); ALGEBRAIC[67] = ALGEBRAIC[49]; ALGEBRAIC[69] = CONSTANTS[113]*ALGEBRAIC[68]*(STATES[1] - ALGEBRAIC[67]); ALGEBRAIC[72] = ALGEBRAIC[49]; ALGEBRAIC[73] = (CONSTANTS[114]/( (1.00000+exp((STATES[1]+57.5300)/- 5.86300))*(1.00000+exp((STATES[1] - 45.8000)/25.8700))))*( STATES[45]*CONSTANTS[68]+ STATES[46]*(1.00000 - CONSTANTS[68]))*(STATES[1] - ALGEBRAIC[72]); ALGEBRAIC[76] = ALGEBRAIC[49]; ALGEBRAIC[77] = (( CONSTANTS[116]*(0.800000 - 0.650000/(1.00000+exp((STATES[1]+125.000)/15.0000)))*0.850000)/(1.00000+pow(66.0000/STATES[2], 2.80000)))*(STATES[1] - ALGEBRAIC[76]); ALGEBRAIC[81] = (( (( 1.00000*CONSTANTS[126]*STATES[1]*pow(CONSTANTS[2], 2.00000))/( CONSTANTS[0]*CONSTANTS[1]))*( 0.750000*STATES[3]*exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 0.750000*STATES[7]))/(exp(( STATES[1]*CONSTANTS[2])/( CONSTANTS[0]*CONSTANTS[1])) - 1.00000))/(1.00000+pow(0.00250000/STATES[40], 3.00000)); ALGEBRAIC[100] = ALGEBRAIC[49]; ALGEBRAIC[102] = (( CONSTANTS[117]*1.00000)/(1.00000+pow(CONSTANTS[6]/0.114000, 2.00000)))*pow(STATES[7]/4.00000, 0.240000)*(STATES[1] - ALGEBRAIC[100]); ALGEBRAIC[105] = ( 1.00000*(STATES[22] - STATES[40])*CONSTANTS[86])/CONSTANTS[79]; ALGEBRAIC[104] = (( 1.00000*CONSTANTS[127])/CONSTANTS[73])*(STATES[7] - CONSTANTS[5]); ALGEBRAIC[106] = 20000.0*STATES[40]*(1.00000 - STATES[53]) - 0.0700000*STATES[53]; ALGEBRAIC[107] = 40000.0*STATES[40]*(1.00000 - STATES[54]) - STATES[54]*40.0000; }