/* There are a total of 233 entries in the algebraic variable array. There are a total of 44 entries in each of the rate and state variable arrays. There are a total of 186 entries in the constant variable array. */ /* * CONSTANTS[0] is kappa_R_TRPNCa in component BG_parameters (fmol_per_sec). * CONSTANTS[1] is kappa_R_BU in component BG_parameters (fmol_per_sec). * CONSTANTS[2] is kappa_R_UW in component BG_parameters (fmol_per_sec). * CONSTANTS[3] is kappa_R_WS in component BG_parameters (fmol_per_sec). * CONSTANTS[4] is kappa_R_SU in component BG_parameters (fmol_per_sec). * CONSTANTS[5] is kappa_LCC_Ca1 in component BG_parameters (fmol_per_sec). * CONSTANTS[6] is kappa_LCC_Ca2 in component BG_parameters (fmol_per_sec). * CONSTANTS[7] is kappa_LCC_K1 in component BG_parameters (fmol_per_sec). * CONSTANTS[8] is kappa_LCC_K2 in component BG_parameters (fmol_per_sec). * CONSTANTS[9] is kappa_d000 in component BG_parameters (fmol_per_sec). * CONSTANTS[10] is kappa_d010 in component BG_parameters (fmol_per_sec). * CONSTANTS[11] is kappa_d020 in component BG_parameters (fmol_per_sec). * CONSTANTS[12] is kappa_d001 in component BG_parameters (fmol_per_sec). * CONSTANTS[13] is kappa_d011 in component BG_parameters (fmol_per_sec). * CONSTANTS[14] is kappa_d021 in component BG_parameters (fmol_per_sec). * CONSTANTS[15] is kappa_f1_000 in component BG_parameters (fmol_per_sec). * CONSTANTS[16] is kappa_f1_100 in component BG_parameters (fmol_per_sec). * CONSTANTS[17] is kappa_f1_001 in component BG_parameters (fmol_per_sec). * CONSTANTS[18] is kappa_f1_101 in component BG_parameters (fmol_per_sec). * CONSTANTS[19] is kappa_f2_000 in component BG_parameters (fmol_per_sec). * CONSTANTS[20] is kappa_f2_100 in component BG_parameters (fmol_per_sec). * CONSTANTS[21] is kappa_f2_001 in component BG_parameters (fmol_per_sec). * CONSTANTS[22] is kappa_f2_101 in component BG_parameters (fmol_per_sec). * CONSTANTS[23] is kappa_f3_010 in component BG_parameters (fmol_per_sec). * CONSTANTS[24] is kappa_f3_110 in component BG_parameters (fmol_per_sec). * CONSTANTS[25] is kappa_f3_011 in component BG_parameters (fmol_per_sec). * CONSTANTS[26] is kappa_f3_111 in component BG_parameters (fmol_per_sec). * CONSTANTS[27] is kappa_fCa000 in component BG_parameters (fmol_per_sec). * CONSTANTS[28] is kappa_fCa100 in component BG_parameters (fmol_per_sec). * CONSTANTS[29] is kappa_fCa010 in component BG_parameters (fmol_per_sec). * CONSTANTS[30] is kappa_fCa110 in component BG_parameters (fmol_per_sec). * CONSTANTS[31] is kappa_fCa020 in component BG_parameters (fmol_per_sec). * CONSTANTS[32] is kappa_fCa120 in component BG_parameters (fmol_per_sec). * CONSTANTS[33] is kappa_PLBph1 in component BG_parameters (fmol_per_sec). * CONSTANTS[34] is kappa_PLBph2 in component BG_parameters (fmol_per_sec). * CONSTANTS[35] is kappa_PLBd1 in component BG_parameters (fmol_per_sec). * CONSTANTS[36] is kappa_PLBd2 in component BG_parameters (fmol_per_sec). * CONSTANTS[37] is kappa_Inh in component BG_parameters (fmol_per_sec). * CONSTANTS[38] is kappa_RyR in component BG_parameters (fmol_per_sec). * CONSTANTS[39] is kappa_OC in component BG_parameters (fmol_per_sec). * CONSTANTS[40] is kappa_CCI in component BG_parameters (fmol_per_sec). * CONSTANTS[41] is kappa_CII in component BG_parameters (fmol_per_sec). * CONSTANTS[42] is kappa_IO in component BG_parameters (fmol_per_sec). * CONSTANTS[43] is kappa_SERCA_R1_2 in component BG_parameters (fmol_per_sec). * CONSTANTS[44] is kappa_SERCA_R2_4 in component BG_parameters (fmol_per_sec). * CONSTANTS[45] is kappa_SERCA_R2_2a in component BG_parameters (fmol_per_sec). * CONSTANTS[46] is kappa_SERCA_R4_5 in component BG_parameters (fmol_per_sec). * CONSTANTS[47] is kappa_SERCA_R5_6 in component BG_parameters (fmol_per_sec). * CONSTANTS[48] is kappa_SERCA_R6_8 in component BG_parameters (fmol_per_sec). * CONSTANTS[49] is kappa_SERCA_R8_9 in component BG_parameters (fmol_per_sec). * CONSTANTS[50] is kappa_SERCA_R9_10 in component BG_parameters (fmol_per_sec). * CONSTANTS[51] is kappa_SERCA_R10_1 in component BG_parameters (fmol_per_sec). * CONSTANTS[52] is K_TRPN in component BG_parameters (per_fmol). * CONSTANTS[53] is K_Cai in component BG_parameters (per_fmol). * CONSTANTS[54] is K_Ca_TRPN in component BG_parameters (per_fmol). * CONSTANTS[55] is K_B in component BG_parameters (per_fmol). * CONSTANTS[56] is K_U in component BG_parameters (per_fmol). * CONSTANTS[57] is K_W in component BG_parameters (per_fmol). * CONSTANTS[58] is K_S in component BG_parameters (per_fmol). * CONSTANTS[59] is K_Ca_i in component BG_parameters (per_fmol). * CONSTANTS[60] is K_Ca_o in component BG_parameters (per_fmol). * CONSTANTS[61] is K_K_i in component BG_parameters (per_fmol). * CONSTANTS[62] is K_K_o in component BG_parameters (per_fmol). * CONSTANTS[63] is K_000_LCC in component BG_parameters (per_fmol). * CONSTANTS[64] is K_010_LCC in component BG_parameters (per_fmol). * CONSTANTS[65] is K_020_LCC in component BG_parameters (per_fmol). * CONSTANTS[66] is K_100_LCC in component BG_parameters (per_fmol). * CONSTANTS[67] is K_110_LCC in component BG_parameters (per_fmol). * CONSTANTS[68] is K_120_LCC in component BG_parameters (per_fmol). * CONSTANTS[69] is K_001_LCC in component BG_parameters (per_fmol). * CONSTANTS[70] is K_011_LCC in component BG_parameters (per_fmol). * CONSTANTS[71] is K_021_LCC in component BG_parameters (per_fmol). * CONSTANTS[72] is K_101_LCC in component BG_parameters (per_fmol). * CONSTANTS[73] is K_111_LCC in component BG_parameters (per_fmol). * CONSTANTS[74] is K_121_LCC in component BG_parameters (per_fmol). * CONSTANTS[75] is K_PLB in component BG_parameters (per_fmol). * CONSTANTS[76] is K_PKACI in component BG_parameters (per_fmol). * CONSTANTS[77] is K_PLB_PKACI in component BG_parameters (per_fmol). * CONSTANTS[78] is K_PP1 in component BG_parameters (per_fmol). * CONSTANTS[79] is K_PLBp_PP1 in component BG_parameters (per_fmol). * CONSTANTS[80] is K_PLBp in component BG_parameters (per_fmol). * CONSTANTS[81] is K_Ip in component BG_parameters (per_fmol). * CONSTANTS[82] is K_Ip_PP1 in component BG_parameters (per_fmol). * CONSTANTS[83] is K_Ca_SR in component BG_parameters (per_fmol). * CONSTANTS[84] is K_C_RyR in component BG_parameters (per_fmol). * CONSTANTS[85] is K_CI_RyR in component BG_parameters (per_fmol). * CONSTANTS[86] is K_I_RyR in component BG_parameters (per_fmol). * CONSTANTS[87] is K_O_RyR in component BG_parameters (per_fmol). * CONSTANTS[88] is K_P1_SERCA in component BG_parameters (per_fmol). * CONSTANTS[89] is K_P2_SERCA in component BG_parameters (per_fmol). * CONSTANTS[90] is K_P2a_SERCA in component BG_parameters (per_fmol). * CONSTANTS[91] is K_P4_SERCA in component BG_parameters (per_fmol). * CONSTANTS[92] is K_P5_SERCA in component BG_parameters (per_fmol). * CONSTANTS[93] is K_P6_SERCA in component BG_parameters (per_fmol). * CONSTANTS[94] is K_P8_SERCA in component BG_parameters (per_fmol). * CONSTANTS[95] is K_P9_SERCA in component BG_parameters (per_fmol). * CONSTANTS[96] is K_P10_SERCA in component BG_parameters (per_fmol). * CONSTANTS[97] is K_H in component BG_parameters (per_fmol). * CONSTANTS[98] is K_MgATP in component BG_parameters (per_fmol). * CONSTANTS[99] is K_MgADP in component BG_parameters (per_fmol). * CONSTANTS[100] is K_Pi in component BG_parameters (per_fmol). * CONSTANTS[101] is z_rCa in component BG_parameters (dimensionless). * CONSTANTS[102] is z_fd in component BG_parameters (dimensionless). * CONSTANTS[103] is z_ff1 in component BG_parameters (dimensionless). * CONSTANTS[104] is z_ff2 in component BG_parameters (dimensionless). * CONSTANTS[105] is z_ff3 in component BG_parameters (dimensionless). * CONSTANTS[106] is z_rd in component BG_parameters (dimensionless). * CONSTANTS[107] is z_rf1 in component BG_parameters (dimensionless). * CONSTANTS[108] is z_rf2 in component BG_parameters (dimensionless). * CONSTANTS[109] is z_rf3 in component BG_parameters (dimensionless). * VOI is time in component environment (second). * CONSTANTS[110] is vol_myo in component environment (pL). * CONSTANTS[111] is freq in component environment (dimensionless). * CONSTANTS[112] is C_m in component environment (fF). * CONSTANTS[113] is F in component constants (C_per_mol). * STATES[0] is q_membrane in component environment (fC). * ALGEBRAIC[0] is V_m in component environment (volt). * STATES[1] is q_TRPN in component environment (fmol). * STATES[2] is q_Ca_TRPN in component environment (fmol). * STATES[3] is q_Cai in component environment (fmol). * STATES[4] is q_B in component environment (fmol). * STATES[5] is q_U in component environment (fmol). * STATES[6] is q_W in component environment (fmol). * STATES[7] is q_S in component environment (fmol). * STATES[8] is q_Cao in component environment (fmol). * STATES[9] is q_Ki in component environment (fmol). * STATES[10] is q_Ko in component environment (fmol). * STATES[11] is q_PLB in component environment (fmol). * STATES[12] is q_PKACI in component environment (fmol). * STATES[13] is q_PLB_PKACI in component environment (fmol). * STATES[14] is q_PP1 in component environment (fmol). * STATES[15] is q_PLBp_PP1 in component environment (fmol). * STATES[16] is q_PLBp in component environment (fmol). * STATES[17] is q_Ip in component environment (fmol). * STATES[18] is q_Ip_PP1 in component environment (fmol). * STATES[19] is q_Ca_SR in component environment (fmol). * STATES[20] is q_H in component environment (fmol). * STATES[21] is q_P in component environment (fmol). * STATES[22] is q_MgADP in component environment (fmol). * STATES[23] is q_MgATP in component environment (fmol). * STATES[24] is q_000_LCC in component environment (fmol). * STATES[25] is q_010_LCC in component environment (fmol). * STATES[26] is q_020_LCC in component environment (fmol). * STATES[27] is q_100_LCC in component environment (fmol). * STATES[28] is q_110_LCC in component environment (fmol). * STATES[29] is q_120_LCC in component environment (fmol). * STATES[30] is q_001_LCC in component environment (fmol). * STATES[31] is q_011_LCC in component environment (fmol). * STATES[32] is q_021_LCC in component environment (fmol). * STATES[33] is q_101_LCC in component environment (fmol). * STATES[34] is q_111_LCC in component environment (fmol). * STATES[35] is q_121_LCC in component environment (fmol). * STATES[36] is q_C_RyR in component environment (fmol). * STATES[37] is q_CI_RyR in component environment (fmol). * STATES[38] is q_I_RyR in component environment (fmol). * STATES[39] is q_O_RyR in component environment (fmol). * CONSTANTS[114] is q_P1_SERCA in component environment (fmol). * CONSTANTS[115] is q_P2_SERCA in component environment (fmol). * CONSTANTS[116] is q_P2a_SERCA in component environment (fmol). * CONSTANTS[117] is q_P4_SERCA in component environment (fmol). * CONSTANTS[118] is q_P5_SERCA in component environment (fmol). * CONSTANTS[119] is q_P6_SERCA in component environment (fmol). * CONSTANTS[120] is q_P8_SERCA in component environment (fmol). * CONSTANTS[121] is q_P9_SERCA in component environment (fmol). * CONSTANTS[122] is q_P10_SERCA in component environment (fmol). * ALGEBRAIC[232] is v_R_TRPNCa in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[43] is v_R_BU in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[46] is v_R_UW in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[49] is v_R_WS in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[52] is v_R_SU in component crossbridge_TRPN (fmol_per_sec). * ALGEBRAIC[229] is T_total in component crossbridge_TRPN (kPa). * STATES[40] is SL in component crossbridge_TRPN (metre). * ALGEBRAIC[163] is v_Ca_i_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[161] is v_Ca_o_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[167] is v_K_i_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[168] is v_K_o_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[169] is I_mem_LCC in component LCC (fA). * ALGEBRAIC[44] is v_PLBph1 in component PLB (fmol_per_sec). * ALGEBRAIC[47] is v_PLBph2 in component PLB (fmol_per_sec). * ALGEBRAIC[50] is v_PLBd1 in component PLB (fmol_per_sec). * ALGEBRAIC[53] is v_PLBd2 in component PLB (fmol_per_sec). * ALGEBRAIC[55] is v_Inh in component PLB (fmol_per_sec). * ALGEBRAIC[177] is v_RyR in component RyR (fmol_per_sec). * ALGEBRAIC[182] is v_RyRgate_Ca_i in component RyR (fmol_per_sec). * ALGEBRAIC[209] is v_Ca_i_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[219] is v_Ca_SR_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[226] is v_H_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[204] is v_MgADP_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[201] is v_MgATP_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[199] is v_Pi_SERCA in component SERCA (fmol_per_sec). * ALGEBRAIC[9] is I_pulse in component environment (fA). * CONSTANTS[123] is pulse_start in component environment (second). * CONSTANTS[124] is pulse_end in component environment (second). * CONSTANTS[125] is pulseMag in component environment (fA). * CONSTANTS[126] is pulseHolding in component environment (fA). * ALGEBRAIC[172] is sum_I in component environment (fA). * CONSTANTS[127] is zCa in component ion_valences (dimensionless). * CONSTANTS[128] is zK in component ion_valences (dimensionless). * ALGEBRAIC[2] is PLB_T in component environment (fmol). * ALGEBRAIC[3] is PKACI_T in component environment (fmol). * ALGEBRAIC[4] is Ip_T in component environment (fmol). * ALGEBRAIC[5] is K_i_T in component environment (fmol). * ALGEBRAIC[6] is K_o_T in component environment (fmol). * CONSTANTS[129] is R in component constants (J_per_K_per_mol). * CONSTANTS[130] is T in component constants (kelvin). * CONSTANTS[131] is zNa in component ion_valences (dimensionless). * CONSTANTS[132] is zCl in component ion_valences (dimensionless). * CONSTANTS[133] is n_Tm in component crossbridge_TRPN (dimensionless). * ALGEBRAIC[24] is mu_TRPN in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[28] is mu_Cai in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[30] is mu_Ca_TRPN in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[31] is mu_B in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[34] is mu_U in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[37] is mu_W in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[40] is mu_S in component crossbridge_TRPN (J_per_mol). * ALGEBRAIC[230] is tension in component crossbridge_TRPN (N_per_mm2). * ALGEBRAIC[231] is mu_tension in component crossbridge_TRPN (J_per_mol). * CONSTANTS[134] is kf_coeff in component crossbridge_TRPN (dimensionless). * CONSTANTS[160] is n in component crossbridge_TRPN (J_per_mol). * CONSTANTS[135] is hh in component crossbridge_TRPN (mm2_per_N). * CONSTANTS[136] is SL_0 in component crossbridge_TRPN (metre). * CONSTANTS[137] is q_MS in component crossbridge_TRPN (fmol). * CONSTANTS[138] is r_s in component crossbridge_TRPN (dimensionless). * CONSTANTS[139] is r_w in component crossbridge_TRPN (dimensionless). * CONSTANTS[140] is A_eff in component crossbridge_TRPN (dimensionless). * CONSTANTS[141] is phi in component crossbridge_TRPN (dimensionless). * CONSTANTS[142] is k_uw in component crossbridge_TRPN (per_sec). * CONSTANTS[143] is k_ws in component crossbridge_TRPN (per_sec). * STATES[41] is G_w in component crossbridge_TRPN (metre). * STATES[42] is G_s in component crossbridge_TRPN (metre). * CONSTANTS[157] is c_w in component crossbridge_TRPN (per_sec). * CONSTANTS[158] is c_s in component crossbridge_TRPN (per_sec). * CONSTANTS[159] is A_w in component crossbridge_TRPN (dimensionless). * CONSTANTS[163] is A_s in component crossbridge_TRPN (dimensionless). * ALGEBRAIC[22] is mu_1 in component crossbridge_TRPN (J_per_m). * CONSTANTS[164] is mu_2 in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[7] is mu_3 in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[23] is mu_4 in component crossbridge_TRPN (J_per_m). * CONSTANTS[167] is mu_5 in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[8] is mu_6 in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[26] is v_1 in component crossbridge_TRPN (m_per_s). * CONSTANTS[144] is v_2 in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[27] is v_4 in component crossbridge_TRPN (m_per_s). * CONSTANTS[145] is v_to_mu in component crossbridge_TRPN (Js_per_m2). * ALGEBRAIC[42] is T_active in component crossbridge_TRPN (kPa). * CONSTANTS[146] is T_ref in component crossbridge_TRPN (kPa). * ALGEBRAIC[39] is mu_T_a in component crossbridge_TRPN (kPa). * ALGEBRAIC[33] is mu_T_S in component crossbridge_TRPN (kPa). * ALGEBRAIC[36] is mu_T_W in component crossbridge_TRPN (kPa). * STATES[43] is Cdd in component crossbridge_TRPN (metre). * ALGEBRAIC[212] is eta in component crossbridge_TRPN (per_sec). * CONSTANTS[147] is eta_l in component crossbridge_TRPN (per_sec). * CONSTANTS[148] is eta_s in component crossbridge_TRPN (per_sec). * CONSTANTS[149] is k in component crossbridge_TRPN (dimensionless). * CONSTANTS[150] is alpha in component crossbridge_TRPN (kPa). * ALGEBRAIC[218] is v_Cdd in component crossbridge_TRPN (m_per_s). * ALGEBRAIC[213] is mu_d in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[214] is mu_k in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[215] is mu_f in component crossbridge_TRPN (J_per_m). * ALGEBRAIC[222] is mu_T_passive in component crossbridge_TRPN (kPa). * ALGEBRAIC[225] is T_passive in component crossbridge_TRPN (kPa). * ALGEBRAIC[48] is mu_Ca_o in component LCC (J_per_mol). * ALGEBRAIC[51] is mu_Ca_i in component LCC (J_per_mol). * ALGEBRAIC[54] is mu_K_o in component LCC (J_per_mol). * ALGEBRAIC[56] is mu_K_i in component LCC (J_per_mol). * ALGEBRAIC[45] is V_mem in component LCC (volt). * ALGEBRAIC[65] is Af_LCC_Ca1 in component LCC (J_per_mol). * ALGEBRAIC[97] is Ar_LCC_Ca1 in component LCC (J_per_mol). * ALGEBRAIC[93] is Am_LCC_Ca1 in component LCC (J_per_mol). * ALGEBRAIC[158] is v_LCC_Ca1 in component LCC (fmol_per_sec). * ALGEBRAIC[66] is Af_LCC_Ca2 in component LCC (J_per_mol). * ALGEBRAIC[98] is Ar_LCC_Ca2 in component LCC (J_per_mol). * ALGEBRAIC[94] is Am_LCC_Ca2 in component LCC (J_per_mol). * ALGEBRAIC[160] is v_LCC_Ca2 in component LCC (fmol_per_sec). * ALGEBRAIC[67] is Af_LCC_K1 in component LCC (J_per_mol). * ALGEBRAIC[99] is Ar_LCC_K1 in component LCC (J_per_mol). * ALGEBRAIC[95] is Am_LCC_K1 in component LCC (J_per_mol). * ALGEBRAIC[162] is v_LCC_K1 in component LCC (fmol_per_sec). * ALGEBRAIC[68] is Af_LCC_K2 in component LCC (J_per_mol). * ALGEBRAIC[100] is Ar_LCC_K2 in component LCC (J_per_mol). * ALGEBRAIC[96] is Am_LCC_K2 in component LCC (J_per_mol). * ALGEBRAIC[164] is v_LCC_K2 in component LCC (fmol_per_sec). * ALGEBRAIC[87] is Af_fCa00 in component LCC (J_per_mol). * ALGEBRAIC[137] is Ar_fCa00 in component LCC (J_per_mol). * ALGEBRAIC[138] is v_fCa00 in component LCC (fmol_per_sec). * ALGEBRAIC[88] is Af_fCa01 in component LCC (J_per_mol). * ALGEBRAIC[139] is Ar_fCa01 in component LCC (J_per_mol). * ALGEBRAIC[142] is v_fCa01 in component LCC (fmol_per_sec). * ALGEBRAIC[89] is Af_fCa02 in component LCC (J_per_mol). * ALGEBRAIC[143] is Ar_fCa02 in component LCC (J_per_mol). * ALGEBRAIC[146] is v_fCa02 in component LCC (fmol_per_sec). * ALGEBRAIC[90] is Af_fCa10 in component LCC (J_per_mol). * ALGEBRAIC[147] is Ar_fCa10 in component LCC (J_per_mol). * ALGEBRAIC[150] is v_fCa10 in component LCC (fmol_per_sec). * ALGEBRAIC[91] is Af_fCa11 in component LCC (J_per_mol). * ALGEBRAIC[151] is Ar_fCa11 in component LCC (J_per_mol). * ALGEBRAIC[154] is v_fCa11 in component LCC (fmol_per_sec). * ALGEBRAIC[92] is Af_fCa12 in component LCC (J_per_mol). * ALGEBRAIC[155] is Ar_fCa12 in component LCC (J_per_mol). * ALGEBRAIC[157] is v_fCa12 in component LCC (fmol_per_sec). * ALGEBRAIC[57] is mu_S000_LCC in component LCC (J_per_mol). * ALGEBRAIC[140] is v_S000_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[58] is mu_S010_LCC in component LCC (J_per_mol). * ALGEBRAIC[144] is v_S010_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[59] is mu_S020_LCC in component LCC (J_per_mol). * ALGEBRAIC[148] is v_S020_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[60] is mu_S100_LCC in component LCC (J_per_mol). * ALGEBRAIC[152] is v_S100_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[61] is mu_S110_LCC in component LCC (J_per_mol). * ALGEBRAIC[156] is v_S110_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[63] is mu_S120_LCC in component LCC (J_per_mol). * ALGEBRAIC[159] is v_S120_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[10] is mu_S001_LCC in component LCC (J_per_mol). * ALGEBRAIC[141] is v_S001_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[11] is mu_S011_LCC in component LCC (J_per_mol). * ALGEBRAIC[145] is v_S011_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[12] is mu_S021_LCC in component LCC (J_per_mol). * ALGEBRAIC[149] is v_S021_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[13] is mu_S101_LCC in component LCC (J_per_mol). * ALGEBRAIC[153] is v_S101_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[62] is mu_S111_LCC in component LCC (J_per_mol). * ALGEBRAIC[166] is v_S111_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[64] is mu_S121_LCC in component LCC (J_per_mol). * ALGEBRAIC[171] is v_S121_LCC in component LCC (fmol_per_sec). * ALGEBRAIC[69] is Af_d000 in component LCC (J_per_mol). * ALGEBRAIC[101] is Ar_d000 in component LCC (J_per_mol). * ALGEBRAIC[102] is v_d000 in component LCC (fmol_per_sec). * ALGEBRAIC[71] is Af_d010 in component LCC (J_per_mol). * ALGEBRAIC[105] is Ar_d010 in component LCC (J_per_mol). * ALGEBRAIC[106] is v_d010 in component LCC (fmol_per_sec). * ALGEBRAIC[73] is Af_d020 in component LCC (J_per_mol). * ALGEBRAIC[109] is Ar_d020 in component LCC (J_per_mol). * ALGEBRAIC[110] is v_d020 in component LCC (fmol_per_sec). * ALGEBRAIC[75] is Af_f1_000 in component LCC (J_per_mol). * ALGEBRAIC[113] is Ar_f1_000 in component LCC (J_per_mol). * ALGEBRAIC[114] is v_f1_000 in component LCC (fmol_per_sec). * ALGEBRAIC[79] is Af_f2_000 in component LCC (J_per_mol). * ALGEBRAIC[121] is Ar_f2_000 in component LCC (J_per_mol). * ALGEBRAIC[122] is v_f2_000 in component LCC (fmol_per_sec). * ALGEBRAIC[83] is Af_f3_010 in component LCC (J_per_mol). * ALGEBRAIC[129] is Ar_f3_010 in component LCC (J_per_mol). * ALGEBRAIC[130] is v_f3_010 in component LCC (fmol_per_sec). * ALGEBRAIC[77] is Af_f1_100 in component LCC (J_per_mol). * ALGEBRAIC[117] is Ar_f1_100 in component LCC (J_per_mol). * ALGEBRAIC[118] is v_f1_100 in component LCC (fmol_per_sec). * ALGEBRAIC[81] is Af_f2_100 in component LCC (J_per_mol). * ALGEBRAIC[125] is Ar_f2_100 in component LCC (J_per_mol). * ALGEBRAIC[126] is v_f2_100 in component LCC (fmol_per_sec). * ALGEBRAIC[85] is Af_f3_110 in component LCC (J_per_mol). * ALGEBRAIC[133] is Ar_f3_110 in component LCC (J_per_mol). * ALGEBRAIC[134] is v_f3_110 in component LCC (fmol_per_sec). * ALGEBRAIC[70] is Af_d001 in component LCC (J_per_mol). * ALGEBRAIC[103] is Ar_d001 in component LCC (J_per_mol). * ALGEBRAIC[104] is v_d001 in component LCC (fmol_per_sec). * ALGEBRAIC[72] is Af_d011 in component LCC (J_per_mol). * ALGEBRAIC[107] is Ar_d011 in component LCC (J_per_mol). * ALGEBRAIC[108] is v_d011 in component LCC (fmol_per_sec). * ALGEBRAIC[74] is Af_d021 in component LCC (J_per_mol). * ALGEBRAIC[111] is Ar_d021 in component LCC (J_per_mol). * ALGEBRAIC[112] is v_d021 in component LCC (fmol_per_sec). * ALGEBRAIC[76] is Af_f1_001 in component LCC (J_per_mol). * ALGEBRAIC[115] is Ar_f1_001 in component LCC (J_per_mol). * ALGEBRAIC[116] is v_f1_001 in component LCC (fmol_per_sec). * ALGEBRAIC[80] is Af_f2_001 in component LCC (J_per_mol). * ALGEBRAIC[123] is Ar_f2_001 in component LCC (J_per_mol). * ALGEBRAIC[124] is v_f2_001 in component LCC (fmol_per_sec). * ALGEBRAIC[84] is Af_f3_011 in component LCC (J_per_mol). * ALGEBRAIC[131] is Ar_f3_011 in component LCC (J_per_mol). * ALGEBRAIC[132] is v_f3_011 in component LCC (fmol_per_sec). * ALGEBRAIC[78] is Af_f1_101 in component LCC (J_per_mol). * ALGEBRAIC[119] is Ar_f1_101 in component LCC (J_per_mol). * ALGEBRAIC[120] is v_f1_101 in component LCC (fmol_per_sec). * ALGEBRAIC[82] is Af_f2_101 in component LCC (J_per_mol). * ALGEBRAIC[127] is Ar_f2_101 in component LCC (J_per_mol). * ALGEBRAIC[128] is v_f2_101 in component LCC (fmol_per_sec). * ALGEBRAIC[86] is Af_f3_111 in component LCC (J_per_mol). * ALGEBRAIC[135] is Ar_f3_111 in component LCC (J_per_mol). * ALGEBRAIC[136] is v_f3_111 in component LCC (fmol_per_sec). * ALGEBRAIC[1] is mu_PLB in component PLB (J_per_mol). * ALGEBRAIC[21] is mu_PKACI in component PLB (J_per_mol). * ALGEBRAIC[25] is mu_PLB_PKACI in component PLB (J_per_mol). * ALGEBRAIC[29] is mu_PP1 in component PLB (J_per_mol). * ALGEBRAIC[32] is mu_PLBp_PP1 in component PLB (J_per_mol). * ALGEBRAIC[35] is mu_PLBp in component PLB (J_per_mol). * ALGEBRAIC[38] is mu_Ip in component PLB (J_per_mol). * ALGEBRAIC[41] is mu_Ip_PP1 in component PLB (J_per_mol). * CONSTANTS[151] is nCa_1 in component RyR (dimensionless). * CONSTANTS[152] is nCa_2 in component RyR (dimensionless). * ALGEBRAIC[165] is mu_Ca_SR in component RyR (J_per_mol). * ALGEBRAIC[170] is mu_Ca_i in component RyR (J_per_mol). * ALGEBRAIC[14] is q_Ca_gate_complexes in component RyR (fmol). * ALGEBRAIC[173] is mu_C_RyR in component RyR (J_per_mol). * ALGEBRAIC[174] is mu_CI_RyR in component RyR (J_per_mol). * ALGEBRAIC[175] is mu_I_RyR in component RyR (J_per_mol). * ALGEBRAIC[176] is mu_O_RyR in component RyR (J_per_mol). * ALGEBRAIC[178] is v_OC in component RyR (fmol_per_sec). * ALGEBRAIC[179] is v_CCI in component RyR (fmol_per_sec). * ALGEBRAIC[180] is v_CII in component RyR (fmol_per_sec). * ALGEBRAIC[181] is v_IO in component RyR (fmol_per_sec). * CONSTANTS[153] is n_Ca_i in component SERCA (dimensionless). * CONSTANTS[154] is n_Ca_SR in component SERCA (dimensionless). * CONSTANTS[155] is n_H in component SERCA (dimensionless). * CONSTANTS[161] is q_SERCA_Ca_complexes in component SERCA (fmol). * ALGEBRAIC[15] is c_Ca_i in component SERCA (mM). * ALGEBRAIC[16] is c_Ca_SR in component SERCA (mM). * ALGEBRAIC[17] is c_H in component SERCA (mM). * ALGEBRAIC[18] is c_MgADP in component SERCA (mM). * ALGEBRAIC[19] is c_MgATP in component SERCA (mM). * ALGEBRAIC[20] is c_Pi in component SERCA (mM). * CONSTANTS[156] is vol_i in component SERCA (pL). * CONSTANTS[162] is vol_sr in component SERCA (pL). * CONSTANTS[166] is vol_isr in component SERCA (pL). * ALGEBRAIC[183] is mu_Ca_i in component SERCA (J_per_mol). * ALGEBRAIC[184] is mu_Ca_SR in component SERCA (J_per_mol). * ALGEBRAIC[185] is mu_H in component SERCA (J_per_mol). * ALGEBRAIC[186] is mu_MgADP in component SERCA (J_per_mol). * ALGEBRAIC[187] is mu_MgATP in component SERCA (J_per_mol). * CONSTANTS[165] is mu_P1 in component SERCA (J_per_mol). * ALGEBRAIC[203] is v_P1 in component SERCA (fmol_per_sec). * ALGEBRAIC[188] is mu_Pi in component SERCA (J_per_mol). * CONSTANTS[168] is mu_P2 in component SERCA (J_per_mol). * ALGEBRAIC[210] is v_P2 in component SERCA (fmol_per_sec). * CONSTANTS[169] is mu_P2a in component SERCA (J_per_mol). * ALGEBRAIC[207] is v_P2a in component SERCA (fmol_per_sec). * CONSTANTS[170] is mu_P4 in component SERCA (J_per_mol). * ALGEBRAIC[216] is v_P4 in component SERCA (fmol_per_sec). * CONSTANTS[171] is mu_P5 in component SERCA (J_per_mol). * ALGEBRAIC[217] is v_P5 in component SERCA (fmol_per_sec). * CONSTANTS[173] is mu_P6 in component SERCA (J_per_mol). * ALGEBRAIC[221] is v_P6 in component SERCA (fmol_per_sec). * CONSTANTS[174] is mu_P8 in component SERCA (J_per_mol). * ALGEBRAIC[224] is v_P8 in component SERCA (fmol_per_sec). * CONSTANTS[175] is mu_P9 in component SERCA (J_per_mol). * ALGEBRAIC[227] is v_P9 in component SERCA (fmol_per_sec). * CONSTANTS[172] is mu_P10 in component SERCA (J_per_mol). * ALGEBRAIC[228] is v_P10 in component SERCA (fmol_per_sec). * ALGEBRAIC[190] is Af_R1_2 in component SERCA (J_per_mol). * CONSTANTS[177] is Ar_R1_2 in component SERCA (J_per_mol). * ALGEBRAIC[200] is v_SERCA_R1_2 in component SERCA (fmol_per_sec). * CONSTANTS[178] is Af_R5_6 in component SERCA (J_per_mol). * ALGEBRAIC[191] is Ar_R5_6 in component SERCA (J_per_mol). * ALGEBRAIC[202] is v_SERCA_R5_6 in component SERCA (fmol_per_sec). * ALGEBRAIC[192] is Af_R2_2a in component SERCA (J_per_mol). * CONSTANTS[179] is Ar_R2_2a in component SERCA (J_per_mol). * ALGEBRAIC[205] is v_SERCA_R2_2a in component SERCA (fmol_per_sec). * ALGEBRAIC[193] is Af_R2_4 in component SERCA (J_per_mol). * CONSTANTS[180] is Ar_R2_4 in component SERCA (J_per_mol). * ALGEBRAIC[206] is v_SERCA_R2_4 in component SERCA (fmol_per_sec). * CONSTANTS[181] is Af_R4_5 in component SERCA (J_per_mol). * ALGEBRAIC[194] is Ar_R4_5 in component SERCA (J_per_mol). * ALGEBRAIC[208] is v_SERCA_R4_5 in component SERCA (fmol_per_sec). * CONSTANTS[182] is Af_R6_8 in component SERCA (J_per_mol). * ALGEBRAIC[195] is Ar_R6_8 in component SERCA (J_per_mol). * ALGEBRAIC[211] is v_SERCA_R6_8 in component SERCA (fmol_per_sec). * ALGEBRAIC[196] is Af_R8_9 in component SERCA (J_per_mol). * CONSTANTS[183] is Ar_R8_9 in component SERCA (J_per_mol). * ALGEBRAIC[220] is v_SERCA_R8_9 in component SERCA (fmol_per_sec). * CONSTANTS[184] is Af_R9_10 in component SERCA (J_per_mol). * ALGEBRAIC[197] is Ar_R9_10 in component SERCA (J_per_mol). * ALGEBRAIC[223] is v_SERCA_R9_10 in component SERCA (fmol_per_sec). * CONSTANTS[176] is Af_R10_1 in component SERCA (J_per_mol). * ALGEBRAIC[189] is Ar_R10_1 in component SERCA (J_per_mol). * ALGEBRAIC[198] is v_SERCA_R10_1 in component SERCA (fmol_per_sec). * RATES[3] is d/dt q_Cai in component environment (fmol). * RATES[1] is d/dt q_TRPN in component environment (fmol). * RATES[2] is d/dt q_Ca_TRPN in component environment (fmol). * RATES[4] is d/dt q_B in component environment (fmol). * RATES[5] is d/dt q_U in component environment (fmol). * RATES[6] is d/dt q_W in component environment (fmol). * RATES[7] is d/dt q_S in component environment (fmol). * RATES[8] is d/dt q_Cao in component environment (fmol). * RATES[9] is d/dt q_Ki in component environment (fmol). * RATES[10] is d/dt q_Ko in component environment (fmol). * RATES[11] is d/dt q_PLB in component environment (fmol). * RATES[12] is d/dt q_PKACI in component environment (fmol). * RATES[13] is d/dt q_PLB_PKACI in component environment (fmol). * RATES[14] is d/dt q_PP1 in component environment (fmol). * RATES[15] is d/dt q_PLBp_PP1 in component environment (fmol). * RATES[16] is d/dt q_PLBp in component environment (fmol). * RATES[17] is d/dt q_Ip in component environment (fmol). * RATES[18] is d/dt q_Ip_PP1 in component environment (fmol). * RATES[19] is d/dt q_Ca_SR in component environment (fmol). * RATES[20] is d/dt q_H in component environment (fmol). * RATES[22] is d/dt q_MgADP in component environment (fmol). * RATES[23] is d/dt q_MgATP in component environment (fmol). * RATES[21] is d/dt q_P in component environment (fmol). * RATES[0] is d/dt q_membrane in component environment (fC). * RATES[41] is d/dt G_w in component crossbridge_TRPN (metre). * RATES[40] is d/dt SL in component crossbridge_TRPN (metre). * RATES[42] is d/dt G_s in component crossbridge_TRPN (metre). * RATES[43] is d/dt Cdd in component crossbridge_TRPN (metre). * RATES[24] is d/dt q_000_LCC in component environment (fmol). * RATES[25] is d/dt q_010_LCC in component environment (fmol). * RATES[26] is d/dt q_020_LCC in component environment (fmol). * RATES[27] is d/dt q_100_LCC in component environment (fmol). * RATES[28] is d/dt q_110_LCC in component environment (fmol). * RATES[29] is d/dt q_120_LCC in component environment (fmol). * RATES[30] is d/dt q_001_LCC in component environment (fmol). * RATES[31] is d/dt q_011_LCC in component environment (fmol). * RATES[32] is d/dt q_021_LCC in component environment (fmol). * RATES[33] is d/dt q_101_LCC in component environment (fmol). * RATES[34] is d/dt q_111_LCC in component environment (fmol). * RATES[35] is d/dt q_121_LCC in component environment (fmol). * RATES[39] is d/dt q_O_RyR in component environment (fmol). * RATES[36] is d/dt q_C_RyR in component environment (fmol). * RATES[37] is d/dt q_CI_RyR in component environment (fmol). * RATES[38] is d/dt q_I_RyR in component environment (fmol). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 52.0453; CONSTANTS[1] = 2.90875; CONSTANTS[2] = 0.332501; CONSTANTS[3] = 0.00782356; CONSTANTS[4] = 140824; CONSTANTS[5] = 0.756518; CONSTANTS[6] = 0.996063; CONSTANTS[7] = 0.0247377; CONSTANTS[8] = 0.0325707; CONSTANTS[9] = 329.349; CONSTANTS[10] = 16.6321; CONSTANTS[11] = 21.8985; CONSTANTS[12] = 0.141736; CONSTANTS[13] = 0.00715766; CONSTANTS[14] = 0.00942408; CONSTANTS[15] = 31.7244; CONSTANTS[16] = 157.535; CONSTANTS[17] = 0.0136527; CONSTANTS[18] = 0.0677958; CONSTANTS[19] = 2.06197; CONSTANTS[20] = 10.2392; CONSTANTS[21] = 0.000887375; CONSTANTS[22] = 0.00440648; CONSTANTS[23] = 4498.81; CONSTANTS[24] = 22339.9; CONSTANTS[25] = 1.93608; CONSTANTS[26] = 9.61406; CONSTANTS[27] = 48530.3; CONSTANTS[28] = 240989; CONSTANTS[29] = 2450.77; CONSTANTS[30] = 12169.9; CONSTANTS[31] = 3226.79; CONSTANTS[32] = 16023.4; CONSTANTS[33] = 45.5263; CONSTANTS[34] = 6.55904; CONSTANTS[35] = 0.386674; CONSTANTS[36] = 1.21269; CONSTANTS[37] = 431.435; CONSTANTS[38] = 2.85549e+06; CONSTANTS[39] = 8.83262; CONSTANTS[40] = 0.12618; CONSTANTS[41] = 883.262; CONSTANTS[42] = 73.6052; CONSTANTS[43] = 1.57289e-05; CONSTANTS[44] = 175018; CONSTANTS[45] = 2.44679e+06; CONSTANTS[46] = 175018; CONSTANTS[47] = 3487.81; CONSTANTS[48] = 4.39557e+06; CONSTANTS[49] = 1.01676e+07; CONSTANTS[50] = 1.01676e+07; CONSTANTS[51] = 0.00441447; CONSTANTS[52] = 1.27424; CONSTANTS[53] = 1.27424; CONSTANTS[54] = 0.11171; CONSTANTS[55] = 6.06164; CONSTANTS[56] = 2.27312; CONSTANTS[57] = 14.8627; CONSTANTS[58] = 3.71567e-06; CONSTANTS[59] = 0.0343245; CONSTANTS[60] = 0.227859; CONSTANTS[61] = 0.000825154; CONSTANTS[62] = 0.00547767; CONSTANTS[63] = 1.47795; CONSTANTS[64] = 29.2665; CONSTANTS[65] = 22.2281; CONSTANTS[66] = 0.297629; CONSTANTS[67] = 5.89367; CONSTANTS[68] = 4.47629; CONSTANTS[69] = 3434.28; CONSTANTS[70] = 68005.7; CONSTANTS[71] = 51650.9; CONSTANTS[72] = 691.593; CONSTANTS[73] = 13695; CONSTANTS[74] = 10401.4; CONSTANTS[75] = 0.00235741; CONSTANTS[76] = 0.395191; CONSTANTS[77] = 0.638527; CONSTANTS[78] = 0.361989; CONSTANTS[79] = 0.203757; CONSTANTS[80] = 0.326014; CONSTANTS[81] = 5.41093; CONSTANTS[82] = 0.0673793; CONSTANTS[83] = 0.980701; CONSTANTS[84] = 3962.58; CONSTANTS[85] = 39.6258; CONSTANTS[86] = 0.06793; CONSTANTS[87] = 6.793; CONSTANTS[88] = 38162.5; CONSTANTS[89] = 128177; CONSTANTS[90] = 10585.3; CONSTANTS[91] = 147985; CONSTANTS[92] = 0.728251; CONSTANTS[93] = 29565.2; CONSTANTS[94] = 4226.28; CONSTANTS[95] = 2547.31; CONSTANTS[96] = 4643.82; CONSTANTS[97] = 212.798; CONSTANTS[98] = 1254.32; CONSTANTS[99] = 1.89443e-05; CONSTANTS[100] = 0.0257106; CONSTANTS[101] = 2; CONSTANTS[102] = 2.1404; CONSTANTS[103] = -1.1495; CONSTANTS[104] = 0.72162; CONSTANTS[105] = 4.2933; CONSTANTS[106] = -2.1404; CONSTANTS[107] = 1.8993; CONSTANTS[108] = -0.52288; CONSTANTS[109] = 0; CONSTANTS[110] = 34.4; CONSTANTS[111] = 500; CONSTANTS[112] = 1.381e5; CONSTANTS[113] = 96485; STATES[0] = -13039; STATES[1] = 2.57; STATES[2] = 1e-18; STATES[3] = 6.82e-3; STATES[4] = 1; STATES[5] = 0; STATES[6] = 0; STATES[7] = 0; STATES[8] = 6.84; STATES[9] = 5.51E+03; STATES[10] = 2.05E+02; STATES[11] = 4.028E+00; STATES[12] = 2.234E-03; STATES[13] = 1e-18; STATES[14] = 3.382E-02; STATES[15] = 1e-18; STATES[16] = 1e-18; STATES[17] = 1.999E-03; STATES[18] = 1e-18; STATES[19] = 0.641; STATES[20] = 0.004028; STATES[21] = 570; STATES[22] = 1.3794; STATES[23] = 3.8; STATES[24] = 1.3798e-05; STATES[25] = 1.3798e-05; STATES[26] = 1.3798e-05; STATES[27] = 1.3798e-05; STATES[28] = 1.3798e-05; STATES[29] = 1.3798e-05; STATES[30] = 1.3798e-05; STATES[31] = 1.3798e-05; STATES[32] = 1.3798e-05; STATES[33] = 1.3798e-05; STATES[34] = 1.3798e-05; STATES[35] = 1.3798e-05; STATES[36] = 4.13941e-05; STATES[37] = 4.13941e-05; STATES[38] = 4.13941e-05; STATES[39] = 4.13941e-05; CONSTANTS[114] = 2.83036e-14; CONSTANTS[115] = 2.83036e-14; CONSTANTS[116] = 2.83036e-14; CONSTANTS[117] = 2.83036e-14; CONSTANTS[118] = 2.83036e-14; CONSTANTS[119] = 2.83036e-14; CONSTANTS[120] = 2.83036e-14; CONSTANTS[121] = 2.83036e-14; CONSTANTS[122] = 2.83036e-14; STATES[40] = 1e-6; CONSTANTS[123] = 4e-2; CONSTANTS[124] = 4.1e-2; CONSTANTS[125] = 1e8; CONSTANTS[126] = 0; CONSTANTS[127] = 2; CONSTANTS[128] = 1; CONSTANTS[129] = 8.31; CONSTANTS[130] = 310; CONSTANTS[131] = 1; CONSTANTS[132] = -1; CONSTANTS[133] = 1.1; CONSTANTS[134] = -0.0118; CONSTANTS[135] = 1; CONSTANTS[136] = 2e-6; CONSTANTS[137] = 1e-6; CONSTANTS[138] = 0.25; CONSTANTS[139] = 0.5; CONSTANTS[140] = 25; CONSTANTS[141] = 2.23; CONSTANTS[142] = 26; CONSTANTS[143] = 4; STATES[41] = 1e-6; STATES[42] = 1e-6; CONSTANTS[144] = 0; CONSTANTS[145] = 1; CONSTANTS[146] = 40.5; STATES[43] = 0; CONSTANTS[147] = 200e3; CONSTANTS[148] = 200e3; CONSTANTS[149] = 7; CONSTANTS[150] = 2.1; CONSTANTS[151] = 1; CONSTANTS[152] = 2; CONSTANTS[153] = 2; CONSTANTS[154] = 2; CONSTANTS[155] = 2; CONSTANTS[156] = 34.0; CONSTANTS[157] = ( CONSTANTS[141]*CONSTANTS[142]*(1.00000 - CONSTANTS[139]))/CONSTANTS[139]; CONSTANTS[158] = ( CONSTANTS[141]*CONSTANTS[143]*(1.00000 - CONSTANTS[138])*CONSTANTS[139])/CONSTANTS[138]; CONSTANTS[159] = ( CONSTANTS[140]*CONSTANTS[138])/( (1.00000 - CONSTANTS[138])*CONSTANTS[139]+CONSTANTS[138]); CONSTANTS[160] = CONSTANTS[134]*CONSTANTS[129]*CONSTANTS[130]; CONSTANTS[161] = CONSTANTS[153]*CONSTANTS[117]+ CONSTANTS[153]*CONSTANTS[118]+ CONSTANTS[154]*CONSTANTS[119]; CONSTANTS[162] = CONSTANTS[156]*0.0350000; CONSTANTS[185] = CONSTANTS[144]; CONSTANTS[163] = CONSTANTS[159]; CONSTANTS[164] = - CONSTANTS[159]*CONSTANTS[144]*CONSTANTS[145]; CONSTANTS[165] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[88]*CONSTANTS[114]); CONSTANTS[166] = CONSTANTS[156]+CONSTANTS[162]; CONSTANTS[167] = - CONSTANTS[163]*CONSTANTS[144]*CONSTANTS[145]; CONSTANTS[168] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[89]*CONSTANTS[115]); CONSTANTS[169] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[90]*CONSTANTS[116]); CONSTANTS[170] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[91]*CONSTANTS[117]); CONSTANTS[171] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[92]*CONSTANTS[118]); CONSTANTS[172] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[96]*CONSTANTS[122]); CONSTANTS[173] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[93]*CONSTANTS[119]); CONSTANTS[174] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[94]*CONSTANTS[120]); CONSTANTS[175] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[95]*CONSTANTS[121]); CONSTANTS[176] = CONSTANTS[172]; CONSTANTS[177] = CONSTANTS[168]; CONSTANTS[178] = CONSTANTS[171]; CONSTANTS[179] = CONSTANTS[169]; CONSTANTS[180] = CONSTANTS[170]; CONSTANTS[181] = CONSTANTS[170]; CONSTANTS[182] = CONSTANTS[173]; CONSTANTS[183] = CONSTANTS[175]; CONSTANTS[184] = CONSTANTS[175]; } void computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { RATES[40] = CONSTANTS[185]; ALGEBRAIC[7] = CONSTANTS[157]*STATES[41]*CONSTANTS[145]; ALGEBRAIC[22] = - CONSTANTS[164] - ALGEBRAIC[7]; ALGEBRAIC[26] = ALGEBRAIC[22]/CONSTANTS[145]; RATES[41] = ALGEBRAIC[26]; ALGEBRAIC[8] = CONSTANTS[158]*STATES[42]*CONSTANTS[145]; ALGEBRAIC[23] = - CONSTANTS[167] - ALGEBRAIC[8]; ALGEBRAIC[27] = ALGEBRAIC[23]/CONSTANTS[145]; RATES[42] = ALGEBRAIC[27]; ALGEBRAIC[30] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[54]*STATES[2]); ALGEBRAIC[31] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[55]*STATES[4]); ALGEBRAIC[34] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[56]*STATES[5]); ALGEBRAIC[43] = CONSTANTS[1]*(exp((ALGEBRAIC[31]+ CONSTANTS[133]*ALGEBRAIC[30])/( CONSTANTS[129]*CONSTANTS[130])) - exp((ALGEBRAIC[34]+ CONSTANTS[133]*ALGEBRAIC[30])/( CONSTANTS[129]*CONSTANTS[130]))); RATES[4] = - ALGEBRAIC[43]; ALGEBRAIC[1] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[75]*STATES[11]); ALGEBRAIC[21] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[76]*STATES[12]); ALGEBRAIC[25] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[77]*STATES[13]); ALGEBRAIC[44] = CONSTANTS[33]*(exp((ALGEBRAIC[1]+ALGEBRAIC[21])/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[25]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[35] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[80]*STATES[16]); ALGEBRAIC[47] = CONSTANTS[34]*(exp(ALGEBRAIC[25]/( CONSTANTS[129]*CONSTANTS[130])) - exp((ALGEBRAIC[35]+ALGEBRAIC[21])/( CONSTANTS[129]*CONSTANTS[130]))); RATES[12] = ALGEBRAIC[47] - ALGEBRAIC[44]; RATES[13] = ALGEBRAIC[44] - ALGEBRAIC[47]; ALGEBRAIC[37] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[57]*STATES[6]); ALGEBRAIC[46] = CONSTANTS[2]*(exp(ALGEBRAIC[34]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[37]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[40] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[58]*STATES[7]); ALGEBRAIC[49] = CONSTANTS[3]*(exp(ALGEBRAIC[37]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[40]/( CONSTANTS[129]*CONSTANTS[130]))); RATES[6] = ALGEBRAIC[46] - ALGEBRAIC[49]; ALGEBRAIC[29] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[78]*STATES[14]); ALGEBRAIC[32] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[79]*STATES[15]); ALGEBRAIC[50] = CONSTANTS[35]*(exp((ALGEBRAIC[35]+ALGEBRAIC[29])/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[32]/( CONSTANTS[129]*CONSTANTS[130]))); RATES[16] = ALGEBRAIC[47] - ALGEBRAIC[50]; ALGEBRAIC[52] = CONSTANTS[4]*(exp(ALGEBRAIC[40]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[34]/( CONSTANTS[129]*CONSTANTS[130]))); RATES[5] = (ALGEBRAIC[43] - ALGEBRAIC[46])+ALGEBRAIC[52]; RATES[7] = ALGEBRAIC[49] - ALGEBRAIC[52]; ALGEBRAIC[53] = CONSTANTS[36]*(exp(ALGEBRAIC[32]/( CONSTANTS[129]*CONSTANTS[130])) - exp((ALGEBRAIC[1]+ALGEBRAIC[29])/( CONSTANTS[129]*CONSTANTS[130]))); RATES[11] = ALGEBRAIC[53] - ALGEBRAIC[44]; RATES[15] = ALGEBRAIC[50] - ALGEBRAIC[53]; ALGEBRAIC[38] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[81]*STATES[17]); ALGEBRAIC[41] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[82]*STATES[18]); ALGEBRAIC[55] = CONSTANTS[37]*(exp((ALGEBRAIC[29]+ALGEBRAIC[38])/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[41]/( CONSTANTS[129]*CONSTANTS[130]))); RATES[14] = (ALGEBRAIC[53] - ALGEBRAIC[50]) - ALGEBRAIC[55]; RATES[17] = - ALGEBRAIC[55]; RATES[18] = ALGEBRAIC[55]; ALGEBRAIC[57] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[63]*STATES[24]); ALGEBRAIC[87] = ALGEBRAIC[57]; ALGEBRAIC[51] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[59]*STATES[3]); ALGEBRAIC[137] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[69]*STATES[30])+ CONSTANTS[101]*ALGEBRAIC[51]; ALGEBRAIC[138] = CONSTANTS[27]*(exp(ALGEBRAIC[87]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[137]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[45] = STATES[0]/CONSTANTS[112]; ALGEBRAIC[69] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[63]*STATES[24])+ CONSTANTS[102]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[101] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[66]*STATES[27])+ CONSTANTS[106]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[102] = CONSTANTS[9]*(exp(ALGEBRAIC[69]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[101]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[75] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[63]*STATES[24])+ CONSTANTS[103]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[113] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[64]*STATES[25])+ CONSTANTS[107]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[114] = CONSTANTS[15]*(exp(ALGEBRAIC[75]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[113]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[79] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[63]*STATES[24])+ CONSTANTS[104]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[121] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[65]*STATES[26])+ CONSTANTS[108]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[122] = CONSTANTS[19]*(exp(ALGEBRAIC[79]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[121]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[140] = ((- ALGEBRAIC[102] - ALGEBRAIC[114]) - ALGEBRAIC[122]) - ALGEBRAIC[138]; RATES[24] = ALGEBRAIC[140]; ALGEBRAIC[70] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[69]*STATES[30])+ CONSTANTS[102]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[103] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[72]*STATES[33])+ CONSTANTS[106]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[104] = CONSTANTS[12]*(exp(ALGEBRAIC[70]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[103]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[76] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[69]*STATES[30])+ CONSTANTS[103]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[115] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[70]*STATES[31])+ CONSTANTS[107]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[116] = CONSTANTS[17]*(exp(ALGEBRAIC[76]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[115]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[80] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[69]*STATES[30])+ CONSTANTS[104]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[123] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[71]*STATES[32])+ CONSTANTS[108]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[124] = CONSTANTS[21]*(exp(ALGEBRAIC[80]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[123]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[141] = ((- ALGEBRAIC[104] - ALGEBRAIC[116]) - ALGEBRAIC[124]) - - ALGEBRAIC[138]; RATES[30] = ALGEBRAIC[141]; ALGEBRAIC[58] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[64]*STATES[25]); ALGEBRAIC[88] = ALGEBRAIC[58]; ALGEBRAIC[139] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[70]*STATES[31])+ CONSTANTS[101]*ALGEBRAIC[51]; ALGEBRAIC[142] = CONSTANTS[29]*(exp(ALGEBRAIC[88]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[139]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[71] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[64]*STATES[25])+ CONSTANTS[102]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[105] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[67]*STATES[28])+ CONSTANTS[106]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[106] = CONSTANTS[10]*(exp(ALGEBRAIC[71]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[105]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[83] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[64]*STATES[25])+ CONSTANTS[105]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[129] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[65]*STATES[26])+ CONSTANTS[109]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[130] = CONSTANTS[23]*(exp(ALGEBRAIC[83]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[129]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[144] = ((- ALGEBRAIC[106]+ALGEBRAIC[114]) - ALGEBRAIC[142]) - ALGEBRAIC[130]; RATES[25] = ALGEBRAIC[144]; ALGEBRAIC[72] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[70]*STATES[31])+ CONSTANTS[102]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[107] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[73]*STATES[34])+ CONSTANTS[106]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[108] = CONSTANTS[13]*(exp(ALGEBRAIC[72]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[107]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[84] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[70]*STATES[31])+ CONSTANTS[105]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[131] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[71]*STATES[32])+ CONSTANTS[109]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[132] = CONSTANTS[25]*(exp(ALGEBRAIC[84]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[131]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[145] = ((- ALGEBRAIC[108]+ALGEBRAIC[116]) - - ALGEBRAIC[142]) - ALGEBRAIC[132]; RATES[31] = ALGEBRAIC[145]; ALGEBRAIC[59] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[65]*STATES[26]); ALGEBRAIC[89] = ALGEBRAIC[59]; ALGEBRAIC[143] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[71]*STATES[32])+ CONSTANTS[101]*ALGEBRAIC[51]; ALGEBRAIC[146] = CONSTANTS[31]*(exp(ALGEBRAIC[89]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[143]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[73] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[65]*STATES[26])+ CONSTANTS[102]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[109] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[68]*STATES[29])+ CONSTANTS[106]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[110] = CONSTANTS[11]*(exp(ALGEBRAIC[73]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[109]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[148] = ((- ALGEBRAIC[110]+ALGEBRAIC[122]) - ALGEBRAIC[146])+ALGEBRAIC[130]; RATES[26] = ALGEBRAIC[148]; ALGEBRAIC[74] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[71]*STATES[32])+ CONSTANTS[102]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[111] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[74]*STATES[35])+ CONSTANTS[106]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[112] = CONSTANTS[14]*(exp(ALGEBRAIC[74]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[111]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[149] = ((- ALGEBRAIC[112]+ALGEBRAIC[124]) - - ALGEBRAIC[146])+ALGEBRAIC[132]; RATES[32] = ALGEBRAIC[149]; ALGEBRAIC[60] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[66]*STATES[27]); ALGEBRAIC[90] = ALGEBRAIC[60]; ALGEBRAIC[147] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[72]*STATES[33])+ CONSTANTS[101]*ALGEBRAIC[51]; ALGEBRAIC[150] = CONSTANTS[28]*(exp(ALGEBRAIC[90]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[147]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[77] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[66]*STATES[27])+ CONSTANTS[103]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[117] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[67]*STATES[28])+ CONSTANTS[107]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[118] = CONSTANTS[16]*(exp(ALGEBRAIC[77]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[117]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[81] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[66]*STATES[27])+ CONSTANTS[104]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[125] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[68]*STATES[29])+ CONSTANTS[108]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[126] = CONSTANTS[20]*(exp(ALGEBRAIC[81]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[125]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[152] = ((- ALGEBRAIC[126]+ALGEBRAIC[102]) - ALGEBRAIC[118]) - ALGEBRAIC[150]; RATES[27] = ALGEBRAIC[152]; ALGEBRAIC[78] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[72]*STATES[33])+ CONSTANTS[103]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[119] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[73]*STATES[34])+ CONSTANTS[107]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[120] = CONSTANTS[18]*(exp(ALGEBRAIC[78]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[119]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[82] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[72]*STATES[33])+ CONSTANTS[104]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[127] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[74]*STATES[35])+ CONSTANTS[108]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[128] = CONSTANTS[22]*(exp(ALGEBRAIC[82]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[127]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[153] = ((- ALGEBRAIC[128]+ALGEBRAIC[104]) - ALGEBRAIC[120]) - - ALGEBRAIC[150]; RATES[33] = ALGEBRAIC[153]; ALGEBRAIC[61] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[67]*STATES[28]); ALGEBRAIC[91] = ALGEBRAIC[61]; ALGEBRAIC[151] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[73]*STATES[34])+ CONSTANTS[101]*ALGEBRAIC[51]; ALGEBRAIC[154] = CONSTANTS[30]*(exp(ALGEBRAIC[91]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[151]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[85] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[67]*STATES[28])+ CONSTANTS[105]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[133] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[68]*STATES[29])+ CONSTANTS[109]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[134] = CONSTANTS[24]*(exp(ALGEBRAIC[85]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[133]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[156] = ((ALGEBRAIC[118]+ALGEBRAIC[106]) - ALGEBRAIC[154]) - ALGEBRAIC[134]; RATES[28] = ALGEBRAIC[156]; ALGEBRAIC[63] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[68]*STATES[29]); ALGEBRAIC[92] = ALGEBRAIC[63]; ALGEBRAIC[155] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[74]*STATES[35])+ CONSTANTS[101]*ALGEBRAIC[51]; ALGEBRAIC[157] = CONSTANTS[32]*(exp(ALGEBRAIC[92]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[155]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[159] = ((ALGEBRAIC[126]+ALGEBRAIC[110]) - ALGEBRAIC[157])+ALGEBRAIC[134]; RATES[29] = ALGEBRAIC[159]; ALGEBRAIC[62] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[73]*STATES[34]); ALGEBRAIC[65] = ALGEBRAIC[51]+ CONSTANTS[127]*CONSTANTS[113]*ALGEBRAIC[45]+ALGEBRAIC[62]; ALGEBRAIC[48] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[60]*STATES[8]); ALGEBRAIC[97] = ALGEBRAIC[48]+ALGEBRAIC[62]; ALGEBRAIC[93] = CONSTANTS[127]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[158] = (ALGEBRAIC[93]==0.00000 ? CONSTANTS[5]*(exp(ALGEBRAIC[65]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[97]/( CONSTANTS[129]*CONSTANTS[130]))) : ((( CONSTANTS[5]*ALGEBRAIC[93])/( CONSTANTS[129]*CONSTANTS[130]))/(exp(ALGEBRAIC[93]/( CONSTANTS[129]*CONSTANTS[130])) - 1.00000))*(exp(ALGEBRAIC[65]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[97]/( CONSTANTS[129]*CONSTANTS[130])))); ALGEBRAIC[64] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[74]*STATES[35]); ALGEBRAIC[66] = ALGEBRAIC[51]+ CONSTANTS[127]*CONSTANTS[113]*ALGEBRAIC[45]+ALGEBRAIC[64]; ALGEBRAIC[98] = ALGEBRAIC[48]+ALGEBRAIC[64]; ALGEBRAIC[94] = CONSTANTS[127]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[160] = (ALGEBRAIC[94]==0.00000 ? CONSTANTS[6]*(exp(ALGEBRAIC[66]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[98]/( CONSTANTS[129]*CONSTANTS[130]))) : ((( CONSTANTS[6]*ALGEBRAIC[94])/( CONSTANTS[129]*CONSTANTS[130]))/(exp(ALGEBRAIC[94]/( CONSTANTS[129]*CONSTANTS[130])) - 1.00000))*(exp(ALGEBRAIC[66]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[98]/( CONSTANTS[129]*CONSTANTS[130])))); ALGEBRAIC[161] = ALGEBRAIC[160]+ALGEBRAIC[158]; RATES[8] = ALGEBRAIC[161]; ALGEBRAIC[56] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[61]*STATES[9]); ALGEBRAIC[67] = ALGEBRAIC[56]+ CONSTANTS[128]*CONSTANTS[113]*ALGEBRAIC[45]+ALGEBRAIC[62]; ALGEBRAIC[54] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[62]*STATES[10]); ALGEBRAIC[99] = ALGEBRAIC[54]+ALGEBRAIC[62]; ALGEBRAIC[95] = CONSTANTS[128]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[162] = (ALGEBRAIC[95]==0.00000 ? CONSTANTS[7]*(exp(ALGEBRAIC[67]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[99]/( CONSTANTS[129]*CONSTANTS[130]))) : ((( CONSTANTS[7]*ALGEBRAIC[95])/( CONSTANTS[129]*CONSTANTS[130]))/(exp(ALGEBRAIC[95]/( CONSTANTS[129]*CONSTANTS[130])) - 1.00000))*(exp(ALGEBRAIC[67]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[99]/( CONSTANTS[129]*CONSTANTS[130])))); ALGEBRAIC[86] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[73]*STATES[34])+ CONSTANTS[105]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[135] = CONSTANTS[129]*CONSTANTS[130]*log( CONSTANTS[74]*STATES[35])+ CONSTANTS[109]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[136] = CONSTANTS[26]*(exp(ALGEBRAIC[86]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[135]/( CONSTANTS[129]*CONSTANTS[130]))); ALGEBRAIC[166] = ((ALGEBRAIC[120]+ALGEBRAIC[108]) - ((((- ALGEBRAIC[154]+ALGEBRAIC[158]) - ALGEBRAIC[158])+ALGEBRAIC[162]) - ALGEBRAIC[162])) - ALGEBRAIC[136]; RATES[34] = ALGEBRAIC[166]; ALGEBRAIC[68] = ALGEBRAIC[56]+ CONSTANTS[128]*CONSTANTS[113]*ALGEBRAIC[45]+ALGEBRAIC[64]; ALGEBRAIC[100] = ALGEBRAIC[54]+ALGEBRAIC[64]; ALGEBRAIC[96] = CONSTANTS[128]*CONSTANTS[113]*ALGEBRAIC[45]; ALGEBRAIC[164] = (ALGEBRAIC[96]==0.00000 ? CONSTANTS[8]*(exp(ALGEBRAIC[68]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[100]/( CONSTANTS[129]*CONSTANTS[130]))) : ((( CONSTANTS[8]*ALGEBRAIC[96])/( CONSTANTS[129]*CONSTANTS[130]))/(exp(ALGEBRAIC[96]/( CONSTANTS[129]*CONSTANTS[130])) - 1.00000))*(exp(ALGEBRAIC[68]/( CONSTANTS[129]*CONSTANTS[130])) - exp(ALGEBRAIC[100]/( CONSTANTS[129]*CONSTANTS[130])))); ALGEBRAIC[167] = - ALGEBRAIC[162] - ALGEBRAIC[164]; RATES[9] = ALGEBRAIC[167]; ALGEBRAIC[168] = ALGEBRAIC[164]+ALGEBRAIC[162]; RATES[10] = ALGEBRAIC[168]; ALGEBRAIC[171] = ((ALGEBRAIC[128]+ALGEBRAIC[112]) - ((((- ALGEBRAIC[157]+ALGEBRAIC[160]) - ALGEBRAIC[160])+ALGEBRAIC[164]) - ALGEBRAIC[164]))+ALGEBRAIC[136]; RATES[35] = ALGEBRAIC[171]; ALGEBRAIC[169] = CONSTANTS[113]*((((((((((((((((((((((((((((((((((((( CONSTANTS[108]*ALGEBRAIC[124] - CONSTANTS[104]*ALGEBRAIC[124])+ CONSTANTS[109]*ALGEBRAIC[132]) - CONSTANTS[105]*ALGEBRAIC[132])+ CONSTANTS[107]*ALGEBRAIC[120]) - CONSTANTS[103]*ALGEBRAIC[120])+ CONSTANTS[108]*ALGEBRAIC[128]) - CONSTANTS[104]*ALGEBRAIC[128])+ CONSTANTS[109]*ALGEBRAIC[136]) - CONSTANTS[105]*ALGEBRAIC[136])+ CONSTANTS[107]*ALGEBRAIC[116]) - CONSTANTS[103]*ALGEBRAIC[116])+ CONSTANTS[106]*ALGEBRAIC[112]) - CONSTANTS[102]*ALGEBRAIC[112])+ CONSTANTS[106]*ALGEBRAIC[108]) - CONSTANTS[102]*ALGEBRAIC[108])+ CONSTANTS[106]*ALGEBRAIC[104]) - CONSTANTS[102]*ALGEBRAIC[104])+ CONSTANTS[109]*ALGEBRAIC[134]) - CONSTANTS[105]*ALGEBRAIC[134])+ CONSTANTS[108]*ALGEBRAIC[126]) - CONSTANTS[104]*ALGEBRAIC[126])+ CONSTANTS[107]*ALGEBRAIC[118]) - CONSTANTS[103]*ALGEBRAIC[118])+ CONSTANTS[109]*ALGEBRAIC[130]) - CONSTANTS[105]*ALGEBRAIC[130])+ CONSTANTS[108]*ALGEBRAIC[122]) - CONSTANTS[104]*ALGEBRAIC[122])+ CONSTANTS[107]*ALGEBRAIC[114]) - CONSTANTS[103]*ALGEBRAIC[114])+ CONSTANTS[106]*ALGEBRAIC[110]) - CONSTANTS[102]*ALGEBRAIC[110])+ CONSTANTS[106]*ALGEBRAIC[106]) - CONSTANTS[102]*ALGEBRAIC[106])+ CONSTANTS[106]*ALGEBRAIC[102]) - CONSTANTS[102]*ALGEBRAIC[102]) - ( CONSTANTS[127]*ALGEBRAIC[158]+ CONSTANTS[127]*ALGEBRAIC[160])) - ( CONSTANTS[128]*ALGEBRAIC[162]+ CONSTANTS[128]*ALGEBRAIC[164])); ALGEBRAIC[9] = (VOI>CONSTANTS[123]&&VOICONSTANTS[123]&&VOIaVOI #define CONSTANTS rfi->aCONSTANTS #define RATES rfi->aRATES #define STATES rfi->aSTATES #define ALGEBRAIC rfi->aALGEBRAIC #define pret rfi->aPRET ALGEBRAIC[212] = p[0]; ALGEBRAIC[213] = p[1]; ALGEBRAIC[214] = p[2]; ALGEBRAIC[215] = p[3]; hx[0] = ALGEBRAIC[212] - (ALGEBRAIC[213]>0.00000 ? CONSTANTS[147] : CONSTANTS[148]); hx[1] = ALGEBRAIC[214] - (CONSTANTS[149]/ALGEBRAIC[212])*STATES[43]*CONSTANTS[145]; hx[2] = ALGEBRAIC[215] - (CONSTANTS[149]/ALGEBRAIC[212])*(STATES[40] - CONSTANTS[136])*CONSTANTS[145]; hx[3] = ALGEBRAIC[213] - (- ALGEBRAIC[214]+ALGEBRAIC[215]); #undef VOI #undef CONSTANTS #undef RATES #undef STATES #undef ALGEBRAIC #undef pret } void rootfind_0(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC, int* pret) { static double p[4] = {0.1,0.1,0.1,0.1}; double bp[4], work[LM_DIF_WORKSZ(4, 4)]; struct rootfind_info rfi; rfi.aVOI = VOI; rfi.aCONSTANTS = CONSTANTS; rfi.aRATES = RATES; rfi.aSTATES = STATES; rfi.aALGEBRAIC = ALGEBRAIC; rfi.aPRET = pret; do_levmar(objfunc_0, p, bp, work, pret, 4, &rfi); ALGEBRAIC[212] = p[0]; ALGEBRAIC[213] = p[1]; ALGEBRAIC[214] = p[2]; ALGEBRAIC[215] = p[3]; }