/* There are a total of 134 entries in the algebraic variable array. There are a total of 51 entries in each of the rate and state variable arrays. There are a total of 388 entries in the constant variable array. */ /* * VOI is time in component time (dimensionless). * ALGEBRAIC[1] is Istim in component stimulus_current (dimensionless). * CONSTANTS[0] is stimPeriod in component stimulus_current (dimensionless). * CONSTANTS[1] is stimDuration in component stimulus_current (dimensionless). * CONSTANTS[2] is stimCurrent in component stimulus_current (dimensionless). * CONSTANTS[3] is clamp_ATPi in component parameters (dimensionless). * CONSTANTS[4] is clamp_ATPi_cyto in component parameters (dimensionless). * CONSTANTS[5] is Faraday in component parameters (C_per_mmol). * CONSTANTS[6] is R in component parameters (J_per_mol_per_K). * CONSTANTS[7] is T in component parameters (K). * CONSTANTS[8] is LHospitalThreshold in component parameters (mV). * CONSTANTS[9] is kt_2 in component parameters (dimensionless). * CONSTANTS[10] is kf_2 in component parameters (dimensionless). * CONSTANTS[11] is kf_3 in component parameters (dimensionless). * CONSTANTS[12] is keq in component parameters (dimensionless). * CONSTANTS[13] is CRT_cyto in component parameters (mM). * CONSTANTS[14] is CRT_mito in component parameters (mM). * CONSTANTS[15] is VATPase_cyto in component parameters (dimensionless). * CONSTANTS[16] is Acap in component parameters (cm2). * CONSTANTS[17] is AcCoA in component parameters (mM). * CONSTANTS[18] is aL in component parameters (dimensionless). * CONSTANTS[19] is b in component parameters (dimensionless). * CONSTANTS[20] is bL in component parameters (dimensionless). * CONSTANTS[21] is C_m in component parameters (uF_per_cm2). * CONSTANTS[22] is Cao in component parameters (mM). * CONSTANTS[23] is CIK in component parameters (mM). * CONSTANTS[24] is Cm in component parameters (mM). * CONSTANTS[25] is CMDNtot in component parameters (mM). * CONSTANTS[26] is Cmito in component parameters (mM_per_mV). * CONSTANTS[27] is CPN in component parameters (mM). * CONSTANTS[28] is CoA in component parameters (mM). * CONSTANTS[29] is CSQNtot in component parameters (mM). * CONSTANTS[30] is DpH in component parameters (pH_units). * CONSTANTS[31] is Dpsio in component parameters (mV). * CONSTANTS[32] is eta in component parameters (dimensionless). * CONSTANTS[33] is EtCS in component parameters (mM). * CONSTANTS[34] is EtID in component parameters (mM). * CONSTANTS[35] is EtKG in component parameters (mM). * CONSTANTS[36] is EtMD in component parameters (mM). * CONSTANTS[37] is EtSDH in component parameters (mM). * CONSTANTS[38] is FAD in component parameters (mM). * CONSTANTS[39] is FADH2 in component parameters (mM). * CONSTANTS[40] is fL in component parameters (per_ms). * CONSTANTS[41] is fm in component parameters (dimensionless). * CONSTANTS[42] is fprime in component parameters (per_ms). * CONSTANTS[43] is g in component parameters (dimensionless). * CONSTANTS[44] is G_Cab in component parameters (mS_per_cm2). * CONSTANTS[45] is G_Kp in component parameters (mS_per_cm2). * CONSTANTS[46] is G_Na in component parameters (mS_per_cm2). * CONSTANTS[47] is G_Nab in component parameters (mS_per_cm2). * CONSTANTS[48] is gh in component parameters (mM_per_ms_per_mV). * CONSTANTS[49] is gL in component parameters (per_ms). * CONSTANTS[50] is GLU in component parameters (mM). * CONSTANTS[51] is gprime in component parameters (per_ms). * CONSTANTS[52] is H in component parameters (mM). * CONSTANTS[53] is hm in component parameters (dimensionless). * CONSTANTS[54] is HTRPNtot in component parameters (mM). * CONSTANTS[55] is ICahalf in component parameters (uA_per_cm2). * CONSTANTS[56] is INaKmax in component parameters (uA_per_cm2). * CONSTANTS[57] is IpCamax in component parameters (uA_per_cm2). * CONSTANTS[58] is KAATeq in component parameters (dimensionless). * CONSTANTS[59] is KaCa in component parameters (mM). * CONSTANTS[60] is KACOeq in component parameters (dimensionless). * CONSTANTS[61] is kact in component parameters (mM). * CONSTANTS[62] is KADP in component parameters (mM). * CONSTANTS[63] is kaminus in component parameters (per_ms). * CONSTANTS[64] is kaplus in component parameters (per_mM4_per_ms). * CONSTANTS[65] is kbminus in component parameters (per_ms). * CONSTANTS[66] is kbplus in component parameters (per_mM3_per_ms). * CONSTANTS[67] is Kca in component parameters (mM). * CONSTANTS[68] is kcminus in component parameters (per_ms). * CONSTANTS[69] is kcnsASP in component parameters (per_ms). * CONSTANTS[70] is kcplus in component parameters (per_ms). * CONSTANTS[71] is KCS in component parameters (per_ms). * CONSTANTS[72] is kf1 in component parameters (dimensionless). * CONSTANTS[73] is kfAAT in component parameters (per_ms). * CONSTANTS[74] is kfACO in component parameters (per_ms). * CONSTANTS[75] is Kfb in component parameters (mM). * CONSTANTS[76] is kfFH in component parameters (per_ms). * CONSTANTS[77] is KFHeq in component parameters (dimensionless). * CONSTANTS[78] is kfSL in component parameters (per_mM_per_ms). * CONSTANTS[79] is kh_1 in component parameters (mM). * CONSTANTS[80] is kh_2 in component parameters (mM). * CONSTANTS[81] is Kh1 in component parameters (mM). * CONSTANTS[82] is Kh2 in component parameters (mM). * CONSTANTS[83] is Kh3 in component parameters (mM). * CONSTANTS[84] is Kh4 in component parameters (mM). * CONSTANTS[85] is khtrpn_minus in component parameters (per_ms). * CONSTANTS[86] is khtrpn_plus in component parameters (per_mM_per_ms). * CONSTANTS[87] is Ki_AM in component parameters (mM). * CONSTANTS[88] is Ki_prime_SR in component parameters (mM). * CONSTANTS[89] is Ki_SR in component parameters (mM). * CONSTANTS[90] is Ki1AD_NaK in component parameters (mM). * CONSTANTS[91] is KiADP_CaP in component parameters (mM). * CONSTANTS[92] is kIDH in component parameters (per_ms). * CONSTANTS[93] is KidhNADH in component parameters (mM). * CONSTANTS[94] is KiFUM in component parameters (mM). * CONSTANTS[95] is Kioaa in component parameters (mM). * CONSTANTS[96] is KiOxaa in component parameters (mM). * CONSTANTS[97] is kKGDH in component parameters (per_ms). * CONSTANTS[98] is kltrpn_plus in component parameters (per_mM_per_ms). * CONSTANTS[99] is kltrpn_minus in component parameters (per_ms). * CONSTANTS[100] is Km1AT_NaK in component parameters (mM). * CONSTANTS[101] is Km1ATP_CaP in component parameters (mM). * CONSTANTS[102] is Km2ATP_CaP in component parameters (mM). * CONSTANTS[103] is KmAcCoA in component parameters (mM). * CONSTANTS[104] is Kmal in component parameters (mM). * CONSTANTS[105] is KmATP_AM in component parameters (mM). * CONSTANTS[106] is KmATP_SR in component parameters (mM). * CONSTANTS[107] is KmCMDN in component parameters (mM). * CONSTANTS[108] is KmCa in component parameters (mM). * CONSTANTS[109] is KmCSQN in component parameters (mM). * CONSTANTS[110] is kMDH in component parameters (per_ms). * CONSTANTS[111] is Kmg in component parameters (mM). * CONSTANTS[112] is KmIDNAD in component parameters (mM). * CONSTANTS[113] is Kmiso in component parameters (mM). * CONSTANTS[114] is KmKG in component parameters (mM). * CONSTANTS[115] is KmKGNAD in component parameters (mM). * CONSTANTS[116] is KmKo in component parameters (mM). * CONSTANTS[117] is KmmNAD in component parameters (mM). * CONSTANTS[118] is KmNa in component parameters (mM). * CONSTANTS[119] is KmNai in component parameters (mM). * CONSTANTS[120] is KmnsCa in component parameters (mM). * CONSTANTS[121] is KmOaa in component parameters (mM). * CONSTANTS[122] is KmpCa in component parameters (mM). * CONSTANTS[123] is KmSucc in component parameters (mM). * CONSTANTS[124] is Kna in component parameters (mM). * CONSTANTS[125] is kNaCa in component parameters (uA_per_cm2). * CONSTANTS[126] is Knca in component parameters (mM). * CONSTANTS[127] is Ko in component parameters (mM). * CONSTANTS[128] is Koff in component parameters (dimensionless). * CONSTANTS[129] is Krb in component parameters (mM). * CONSTANTS[130] is kres in component parameters (dimensionless). * CONSTANTS[131] is kresf in component parameters (dimensionless). * CONSTANTS[132] is ksat in component parameters (dimensionless). * CONSTANTS[133] is kSDH in component parameters (per_ms). * CONSTANTS[134] is KSLeq in component parameters (dimensionless). * CONSTANTS[135] is KSR in component parameters (dimensionless). * CONSTANTS[136] is ktrans in component parameters (mM). * CONSTANTS[137] is kTrop_pn in component parameters (per_ms). * CONSTANTS[138] is L in component parameters (dimensionless). * CONSTANTS[139] is LTRPNtot in component parameters (mM). * CONSTANTS[140] is mcoop in component parameters (dimensionless). * CONSTANTS[141] is Mg in component parameters (mM). * CONSTANTS[142] is n in component parameters (dimensionless). * CONSTANTS[143] is na in component parameters (dimensionless). * CONSTANTS[144] is Nao in component parameters (mM). * CONSTANTS[145] is ncoop in component parameters (dimensionless). * CONSTANTS[146] is Nfb in component parameters (dimensionless). * CONSTANTS[147] is nID in component parameters (dimensionless). * CONSTANTS[148] is nKG in component parameters (dimensionless). * CONSTANTS[149] is Nrb in component parameters (dimensionless). * CONSTANTS[150] is omega in component parameters (per_ms). * CONSTANTS[151] is p1 in component parameters (dimensionless). * CONSTANTS[152] is p2 in component parameters (dimensionless). * CONSTANTS[153] is p3 in component parameters (dimensionless). * CONSTANTS[154] is pa in component parameters (per_ms). * CONSTANTS[155] is pb in component parameters (per_ms). * CONSTANTS[156] is pc1 in component parameters (per_ms). * CONSTANTS[157] is pc2 in component parameters (per_ms). * CONSTANTS[158] is PCa in component parameters (cm_per_s). * CONSTANTS[159] is Pi in component parameters (mM). * CONSTANTS[160] is PK in component parameters (cm_per_s). * CONSTANTS[161] is PnsK in component parameters (cm_per_s). * CONSTANTS[162] is PnsNa in component parameters (cm_per_s). * CONSTANTS[163] is r1 in component parameters (dimensionless). * CONSTANTS[164] is r2 in component parameters (dimensionless). * CONSTANTS[165] is r3 in component parameters (dimensionless). * CONSTANTS[166] is ra in component parameters (per_ms). * CONSTANTS[167] is rb in component parameters (per_ms). * CONSTANTS[168] is rc1 in component parameters (per_ms). * CONSTANTS[169] is rc2 in component parameters (per_ms). * CONSTANTS[170] is rhoF1 in component parameters (mM). * CONSTANTS[171] is rhoREF in component parameters (mM). * CONSTANTS[172] is rhoREN in component parameters (mM). * CONSTANTS[173] is tautr in component parameters (ms). * CONSTANTS[174] is tauxfer in component parameters (ms). * CONSTANTS[175] is V_AM_scaler in component parameters (dimensionless). * CONSTANTS[176] is V_AM_max in component parameters (mM_per_ms). * CONSTANTS[177] is v1 in component parameters (per_ms). * CONSTANTS[178] is VJSR in component parameters (uL). * CONSTANTS[179] is vmaxf in component parameters (per_ms). * CONSTANTS[180] is vmaxr in component parameters (per_ms). * CONSTANTS[181] is VmDT in component parameters (mM_per_ms). * CONSTANTS[182] is VmNC in component parameters (mM_per_ms). * CONSTANTS[183] is Vmuni in component parameters (mM_per_ms). * CONSTANTS[184] is Vmyo in component parameters (uL). * CONSTANTS[185] is VNSR in component parameters (uL). * CONSTANTS[186] is VSS in component parameters (uL). * CONSTANTS[187] is zeta in component parameters (N_per_mm2). * CONSTANTS[188] is f_xb in component parameters (per_ms). * CONSTANTS[189] is SL in component parameters (um). * CONSTANTS[190] is gmin_xb in component parameters (per_ms). * CONSTANTS[191] is zNa in component parameters (dimensionless). * CONSTANTS[192] is zK in component parameters (dimensionless). * CONSTANTS[193] is zCa in component parameters (dimensionless). * CONSTANTS[194] is one in component parameters (dimensionless). * CONSTANTS[195] is p_KNa in component parameters (dimensionless). * CONSTANTS[196] is ATPi_initial in component initial_values (dimensionless). * CONSTANTS[197] is V_initial in component initial_values (dimensionless). * CONSTANTS[198] is C1_RyR_initial in component initial_values (dimensionless). * CONSTANTS[199] is C2_RyR_initial in component initial_values (dimensionless). * CONSTANTS[200] is O2_RyR_initial in component initial_values (dimensionless). * CONSTANTS[201] is mNa_initial in component initial_values (dimensionless). * CONSTANTS[202] is hNa_initial in component initial_values (dimensionless). * CONSTANTS[203] is jNa_initial in component initial_values (dimensionless). * CONSTANTS[204] is xKs_initial in component initial_values (dimensionless). * CONSTANTS[205] is Open_initial in component initial_values (dimensionless). * CONSTANTS[206] is OCa_initial in component initial_values (dimensionless). * CONSTANTS[207] is yCa_initial in component initial_values (dimensionless). * CONSTANTS[208] is Ki_initial in component initial_values (dimensionless). * CONSTANTS[209] is Nai_initial in component initial_values (dimensionless). * CONSTANTS[210] is Cai_initial in component initial_values (dimensionless). * CONSTANTS[211] is ADPm_initial in component initial_values (dimensionless). * CONSTANTS[212] is Dpsi_initial in component initial_values (dimensionless). * CONSTANTS[213] is NADH_initial in component initial_values (dimensionless). * CONSTANTS[214] is Oaa_initial in component initial_values (dimensionless). * CONSTANTS[215] is AKG_initial in component initial_values (dimensionless). * CONSTANTS[216] is SCoA_initial in component initial_values (dimensionless). * CONSTANTS[217] is Succ_initial in component initial_values (dimensionless). * CONSTANTS[218] is FUM_initial in component initial_values (dimensionless). * CONSTANTS[219] is MAL_initial in component initial_values (dimensionless). * CONSTANTS[220] is ISOC_initial in component initial_values (dimensionless). * CONSTANTS[221] is Cam_initial in component initial_values (dimensionless). * CONSTANTS[222] is CaNSR_initial in component initial_values (dimensionless). * CONSTANTS[223] is CaJSR_initial in component initial_values (dimensionless). * CONSTANTS[224] is CaSS_initial in component initial_values (dimensionless). * CONSTANTS[225] is N1_initial in component initial_values (dimensionless). * CONSTANTS[226] is P0_initial in component initial_values (dimensionless). * CONSTANTS[227] is P1_initial in component initial_values (dimensionless). * CONSTANTS[228] is P2_initial in component initial_values (dimensionless). * CONSTANTS[229] is P3_initial in component initial_values (dimensionless). * CONSTANTS[230] is LTRPNCa_initial in component initial_values (dimensionless). * CONSTANTS[231] is HTRPNCa_initial in component initial_values (dimensionless). * CONSTANTS[232] is C1_initial in component initial_values (dimensionless). * CONSTANTS[233] is CCa0_initial in component initial_values (dimensionless). * CONSTANTS[234] is C0_initial in component initial_values (dimensionless). * CONSTANTS[235] is C2_initial in component initial_values (dimensionless). * CONSTANTS[236] is CCa1_initial in component initial_values (dimensionless). * CONSTANTS[237] is C3_initial in component initial_values (dimensionless). * CONSTANTS[238] is CCa2_initial in component initial_values (dimensionless). * CONSTANTS[239] is C4_initial in component initial_values (dimensionless). * CONSTANTS[240] is CCa3_initial in component initial_values (dimensionless). * CONSTANTS[241] is CCa4_initial in component initial_values (dimensionless). * CONSTANTS[242] is CrPi_mito_initial in component initial_values (dimensionless). * CONSTANTS[243] is CrPi_cyto_initial in component initial_values (dimensionless). * CONSTANTS[244] is ATPi_cyto_initial in component initial_values (dimensionless). * CONSTANTS[245] is ASP_initial in component initial_values (dimensionless). * CONSTANTS[246] is N0_initial in component initial_values (dimensionless). * CONSTANTS[247] is Vtotal in component computed_constants (uL). * CONSTANTS[303] is Vmito in component computed_constants (uL). * CONSTANTS[248] is f_01 in component computed_constants (per_ms). * CONSTANTS[304] is f_12 in component computed_constants (per_ms). * CONSTANTS[312] is f_23 in component computed_constants (per_ms). * CONSTANTS[316] is g0_01 in component computed_constants (per_ms). * CONSTANTS[321] is g0_12 in component computed_constants (per_ms). * CONSTANTS[326] is g0_23 in component computed_constants (per_ms). * CONSTANTS[329] is paths in component computed_constants (per_ms3). * CONSTANTS[332] is P1max in component computed_constants (dimensionless). * CONSTANTS[338] is P2max in component computed_constants (dimensionless). * CONSTANTS[343] is P3max in component computed_constants (dimensionless). * CONSTANTS[345] is Fmax in component computed_constants (dimensionless). * CONSTANTS[347] is fnormmax in component computed_constants (dimensionless). * CONSTANTS[249] is SLnorm in component computed_constants (dimensionless). * CONSTANTS[250] is Ktrop_Ca in component computed_constants (mM). * CONSTANTS[305] is Ktrop_half in component computed_constants (dimensionless). * CONSTANTS[251] is Ntrop in component computed_constants (dimensionless). * CONSTANTS[349] is fnormmax2 in component computed_constants (dimensionless). * CONSTANTS[351] is La in component computed_constants (um). * CONSTANTS[353] is Lm_prime in component computed_constants (um). * CONSTANTS[355] is Lz in component computed_constants (um). * CONSTANTS[357] is Lb in component computed_constants (um). * CONSTANTS[359] is Lm in component computed_constants (um). * CONSTANTS[330] is mod_factor in component computed_constants (dimensionless). * CONSTANTS[333] is g_01_mod in component computed_constants (per_ms). * CONSTANTS[334] is g_12_mod in component computed_constants (per_ms). * CONSTANTS[335] is g_23_mod in component computed_constants (per_ms). * CONSTANTS[336] is g_01_off in component computed_constants (per_ms). * CONSTANTS[339] is g_01_off_mod in component computed_constants (per_ms). * CONSTANTS[301] is RT_over_F in component computed_constants (mV). * CONSTANTS[252] is G_Ks in component computed_constants (mS_per_cm2). * CONSTANTS[310] is G_K1 in component computed_constants (mS_per_cm2). * CONSTANTS[253] is inv_5p98 in component computed_constants (dimensionless). * CONSTANTS[314] is FaradayE3 in component computed_constants (C_per_mol). * CONSTANTS[254] is Cao_341 in component computed_constants (mM). * CONSTANTS[255] is ICamax_LHospital in component computed_constants (uA_per_cm2). * CONSTANTS[256] is Pca_4En3 in component computed_constants (cm_per_ms). * CONSTANTS[319] is F_over_RT in component computed_constants (per_mV). * CONSTANTS[257] is inv_ICahalf in component computed_constants (cm2_per_uA). * CONSTANTS[317] is PKFe3 in component computed_constants (C_cm_per_s_per_mol). * CONSTANTS[324] is sigma in component computed_constants (dimensionless). * CONSTANTS[258] is inv_KmNai in component computed_constants (per_mM). * CONSTANTS[327] is INaKmax_Ko_Ko_KmKo in component computed_constants (uA_per_cm2). * CONSTANTS[331] is inv_Ki1AD_NaK in component computed_constants (per_mM). * CONSTANTS[259] is eta_1 in component computed_constants (dimensionless). * CONSTANTS[260] is Nao_p3 in component computed_constants (mM2). * CONSTANTS[261] is KmCa_Cao in component computed_constants (mM3_cm2_per_uA). * CONSTANTS[306] is KmCa_Cao_ksat in component computed_constants (mM3_cm2_per_uA). * CONSTANTS[262] is inv_KiADP_CaP in component computed_constants (mM). * CONSTANTS[263] is KmnsCa_p3 in component computed_constants (mM3). * CONSTANTS[318] is V_AM_scaler_max_1_f_01_12_23 in component computed_constants (mM_per_ms4). * CONSTANTS[264] is KmATP_AM_Ki_AM in component computed_constants (dimensionless). * CONSTANTS[307] is DmuH_Constant in component computed_constants (mV). * CONSTANTS[265] is VCS_C1 in component computed_constants (mM_per_ms). * CONSTANTS[266] is one_inv_KACOeq in component computed_constants (dimensionless). * CONSTANTS[267] is VIDH_Constant in component computed_constants (dimensionless). * CONSTANTS[268] is kIDH_EtID in component computed_constants (mM_per_ms). * CONSTANTS[269] is inv_KADP in component computed_constants (per_mM). * CONSTANTS[270] is inv_KaCa in component computed_constants (per_mM). * CONSTANTS[271] is inv_KidhNADH in component computed_constants (per_mM). * CONSTANTS[272] is KmKGNAD_KmIDNAD in component computed_constants (dimensionless). * CONSTANTS[273] is Mg_Kmg_1 in component computed_constants (dimensionless). * CONSTANTS[308] is Mg_Kmg_1_Kca in component computed_constants (mM). * CONSTANTS[274] is kKGDH_EtKG in component computed_constants (mM_per_ms). * CONSTANTS[275] is CoA_KSLeq in component computed_constants (mM). * CONSTANTS[276] is kSDH_EtSDH in component computed_constants (mM_per_ms). * CONSTANTS[277] is KmSucc_KiFUM in component computed_constants (dimensionless). * CONSTANTS[278] is inv_KiOxaa in component computed_constants (per_mM). * CONSTANTS[279] is kfFH_KFHeq in component computed_constants (per_ms). * CONSTANTS[280] is kMDH_Fh_EtMD in component computed_constants (mM_per_ms). * CONSTANTS[281] is Kmal_Kioaa in component computed_constants (dimensionless). * CONSTANTS[282] is VAAT_Constant in component computed_constants (mM_per_ms). * CONSTANTS[283] is kcnsASP_KAATeq_kfAAT in component computed_constants (dimensionless). * CONSTANTS[284] is KfAAT_GLU in component computed_constants (per_ms_per_mM). * CONSTANTS[285] is KfAAT_KAATeq in component computed_constants (per_ms). * CONSTANTS[286] is kres_sq_KmIDNAD in component computed_constants (per_mM). * CONSTANTS[337] is exp_6_FRT_Dpsio in component computed_constants (dimensionless). * CONSTANTS[342] is FRT_6_g in component computed_constants (per_mV). * CONSTANTS[344] is ra_rc1_exp_6_FRT_Dpsio in component computed_constants (per_ms). * CONSTANTS[346] is r1_exp_6_FRT_Dpsio in component computed_constants (dimensionless). * CONSTANTS[348] is rhoREN_ra_rc1_exp_6_FRT_Dpsio in component computed_constants (mM_per_ms). * CONSTANTS[350] is rhoREN_rc2 in component computed_constants (mM_per_ms). * CONSTANTS[352] is rhoREN_ra in component computed_constants (mM_per_ms). * CONSTANTS[354] is rhoRen_6_ra in component computed_constants (mM_per_ms). * CONSTANTS[356] is rhoRen_6_ra_rb in component computed_constants (mM_per_ms). * CONSTANTS[358] is AREF in component computed_constants (mV). * CONSTANTS[360] is exp_AREF_FRT in component computed_constants (dimensionless). * CONSTANTS[361] is ra_rc2_exp_AREF_FRT in component computed_constants (per_ms). * CONSTANTS[362] is VFO_C1 in component computed_constants (per_ms). * CONSTANTS[364] is ra_exp_AREF_FRT in component computed_constants (per_ms). * CONSTANTS[367] is ra_rb in component computed_constants (per_ms). * CONSTANTS[369] is VFO_VHFe_C1 in component computed_constants (per_ms). * CONSTANTS[370] is r2_r3_exp_AREF_FRT in component computed_constants (dimensionless). * CONSTANTS[371] is exp_3_FRT_Dpsio in component computed_constants (dimensionless). * CONSTANTS[372] is FRT_3 in component computed_constants (per_mV). * CONSTANTS[287] is kf1_Pi in component computed_constants (mM). * CONSTANTS[373] is VATPase_C1 in component computed_constants (per_ms). * CONSTANTS[374] is pa_pb_3 in component computed_constants (per_ms). * CONSTANTS[375] is pa_300 in component computed_constants (per_ms). * CONSTANTS[376] is p1_exp_3_FRT_Dpsio in component computed_constants (dimensionless). * CONSTANTS[322] is hm_F_over_RT in component computed_constants (per_mV). * CONSTANTS[288] is VmDT_75 in component computed_constants (mM_per_ms). * CONSTANTS[289] is VmDT_20 in component computed_constants (mM_per_ms). * CONSTANTS[290] is inv_Kfb in component computed_constants (per_mM). * CONSTANTS[291] is inv_Krb in component computed_constants (per_mM). * CONSTANTS[292] is inv_tautr in component computed_constants (per_ms). * CONSTANTS[293] is inv_tauxfer in component computed_constants (per_ms). * CONSTANTS[377] is KmATP_SR_Ki_SR in component computed_constants (dimensionless). * CONSTANTS[378] is inv_Ki_prime_SR in component computed_constants (per_mM). * CONSTANTS[365] is alpha_SL_fnormmax2 in component computed_constants (dimensionless). * CONSTANTS[366] is alpha_SL_fnormmax in component computed_constants (dimensionless). * CONSTANTS[313] is inv_LTRPNtot_Ktrop_half in component computed_constants (per_mM). * CONSTANTS[309] is kTrop_pn_f_01 in component computed_constants (per_ms). * CONSTANTS[340] is kTrop_pn_f_12_g_01_mod in component computed_constants (per_ms). * CONSTANTS[341] is f_23_g_12_mod in component computed_constants (per_ms). * CONSTANTS[294] is CMDNtot_KmCMDN in component computed_constants (mM2). * CONSTANTS[295] is CSQNtot_KmCSQN in component computed_constants (mM2). * CONSTANTS[296] is inv_ktrans in component computed_constants (per_mM). * CONSTANTS[297] is inv_kact in component computed_constants (per_mM). * CONSTANTS[298] is Vmuni_ktrans in component computed_constants (per_ms). * CONSTANTS[323] is FRT2 in component computed_constants (per_mV). * CONSTANTS[299] is b_05 in component computed_constants (dimensionless). * CONSTANTS[379] is Acap_Vmyo_F in component computed_constants (mmol_per_C_per_cm). * CONSTANTS[380] is Acap_VSS_F in component computed_constants (mmol_per_C_per_cm). * CONSTANTS[381] is VJSR_VSS in component computed_constants (dimensionless). * CONSTANTS[382] is Vmyo_VSS in component computed_constants (dimensionless). * CONSTANTS[383] is Vmyo_VNSR in component computed_constants (dimensionless). * CONSTANTS[384] is VJSR_VNSR in component computed_constants (dimensionless). * CONSTANTS[385] is inv_C_m in component computed_constants (cm2_per_uF). * CONSTANTS[302] is inv_bL in component computed_constants (dimensionless). * CONSTANTS[386] is inv_Cmito in component computed_constants (mV_per_mM). * CONSTANTS[387] is two_b in component computed_constants (dimensionless). * CONSTANTS[300] is inv_keq in component computed_constants (dimensionless). * CONSTANTS[368] is zeta_alpha_SL_fnormmax in component computed_constants (N_per_mm2). * CONSTANTS[363] is alpha_SL in component computed_constants (dimensionless). * ALGEBRAIC[7] is ADP in component model_2 (dimensionless). * STATES[0] is ATPi in component model_2 (dimensionless). * ALGEBRAIC[15] is inv_ATPi in component model_2 (dimensionless). * ALGEBRAIC[19] is VF_over_RT in component model_2 (dimensionless). * STATES[1] is V in component model_2 (dimensionless). * ALGEBRAIC[22] is exp_VF_over_RT in component model_2 (dimensionless). * ALGEBRAIC[25] is VFsq_over_RT in component model_2 (dimensionless). * ALGEBRAIC[27] is exp2VFRT in component model_2 (dimensionless). * ALGEBRAIC[29] is O1_RyR in component model_2 (dimensionless). * STATES[2] is C1_RyR in component model_2 (dimensionless). * STATES[3] is C2_RyR in component model_2 (dimensionless). * STATES[4] is O2_RyR in component model_2 (dimensionless). * ALGEBRAIC[0] is V_30 in component model_2 (dimensionless). * ALGEBRAIC[124] is V_E_K in component model_2 (dimensionless). * ALGEBRAIC[121] is INa in component model_2 (dimensionless). * STATES[5] is mNa in component model_2 (dimensionless). * STATES[6] is hNa in component model_2 (dimensionless). * STATES[7] is jNa in component model_2 (dimensionless). * ALGEBRAIC[133] is IKs in component model_2 (dimensionless). * STATES[8] is xKs in component model_2 (dimensionless). * ALGEBRAIC[125] is K1Alpha in component model_2 (dimensionless). * ALGEBRAIC[126] is K1Beta in component model_2 (dimensionless). * ALGEBRAIC[127] is K1_inf in component model_2 (dimensionless). * ALGEBRAIC[128] is IK1 in component model_2 (dimensionless). * ALGEBRAIC[122] is INab in component model_2 (dimensionless). * ALGEBRAIC[129] is IKp in component model_2 (dimensionless). * ALGEBRAIC[31] is ICamax in component model_2 (dimensionless). * ALGEBRAIC[33] is ICaK in component model_2 (dimensionless). * STATES[9] is Open in component model_2 (dimensionless). * STATES[10] is OCa in component model_2 (dimensionless). * STATES[11] is yCa in component model_2 (dimensionless). * STATES[12] is Ki in component model_2 (dimensionless). * ALGEBRAIC[35] is ICa in component model_2 (dimensionless). * ALGEBRAIC[37] is NaiP1p5 in component model_2 (dimensionless). * STATES[13] is Nai in component model_2 (dimensionless). * ALGEBRAIC[39] is INaK in component model_2 (dimensionless). * ALGEBRAIC[41] is exp_eta_VF_over_RT in component model_2 (dimensionless). * ALGEBRAIC[43] is exp_eta1_VF_over_RT in component model_2 (dimensionless). * ALGEBRAIC[45] is INaCa in component model_2 (dimensionless). * STATES[14] is Cai in component model_2 (dimensionless). * ALGEBRAIC[131] is ICab in component model_2 (dimensionless). * ALGEBRAIC[47] is IpCa in component model_2 (dimensionless). * ALGEBRAIC[49] is CaiP3 in component model_2 (dimensionless). * ALGEBRAIC[51] is common in component model_2 (dimensionless). * ALGEBRAIC[53] is InsNa in component model_2 (dimensionless). * ALGEBRAIC[55] is InsK in component model_2 (dimensionless). * ALGEBRAIC[57] is InsCa in component model_2 (dimensionless). * ALGEBRAIC[59] is V_AM in component model_2 (dimensionless). * ALGEBRAIC[61] is ATPm in component model_2 (dimensionless). * STATES[15] is ADPm in component model_2 (dimensionless). * ALGEBRAIC[63] is DmuH in component model_2 (dimensionless). * STATES[16] is Dpsi in component model_2 (dimensionless). * ALGEBRAIC[65] is NAD in component model_2 (dimensionless). * STATES[17] is NADH in component model_2 (dimensionless). * ALGEBRAIC[67] is KmIDNAD_NAD in component model_2 (dimensionless). * ALGEBRAIC[68] is exp_FRT_6_g_DmuH in component model_2 (dimensionless). * ALGEBRAIC[69] is FRT2_Dpsi in component model_2 (dimensionless). * ALGEBRAIC[70] is VCS in component model_2 (dimensionless). * STATES[18] is Oaa in component model_2 (dimensionless). * ALGEBRAIC[71] is VACO in component model_2 (dimensionless). * STATES[19] is AKG in component model_2 (dimensionless). * STATES[20] is SCoA in component model_2 (dimensionless). * STATES[21] is Succ in component model_2 (dimensionless). * STATES[22] is FUM in component model_2 (dimensionless). * STATES[23] is MAL in component model_2 (dimensionless). * STATES[24] is ISOC in component model_2 (dimensionless). * ALGEBRAIC[72] is Fa in component model_2 (dimensionless). * STATES[25] is Cam in component model_2 (dimensionless). * ALGEBRAIC[73] is Fi in component model_2 (dimensionless). * ALGEBRAIC[74] is VIDH in component model_2 (dimensionless). * ALGEBRAIC[75] is a__1 in component model_2 (dimensionless). * ALGEBRAIC[76] is VKGDH in component model_2 (dimensionless). * ALGEBRAIC[77] is VSL in component model_2 (dimensionless). * ALGEBRAIC[78] is VSDH in component model_2 (dimensionless). * ALGEBRAIC[79] is VFH in component model_2 (dimensionless). * ALGEBRAIC[80] is VMDH in component model_2 (dimensionless). * ALGEBRAIC[81] is VAAT in component model_2 (dimensionless). * ALGEBRAIC[82] is AREN in component model_2 (dimensionless). * ALGEBRAIC[83] is denominator1 in component model_2 (dimensionless). * ALGEBRAIC[86] is denominator2 in component model_2 (dimensionless). * ALGEBRAIC[90] is denominator3 in component model_2 (dimensionless). * ALGEBRAIC[84] is VNO in component model_2 (dimensionless). * ALGEBRAIC[85] is VHNe in component model_2 (dimensionless). * ALGEBRAIC[87] is VHFe in component model_2 (dimensionless). * ALGEBRAIC[88] is exp_3FRT_DmuH in component model_2 (dimensionless). * ALGEBRAIC[89] is AF1 in component model_2 (dimensionless). * ALGEBRAIC[91] is VATPase in component model_2 (dimensionless). * ALGEBRAIC[92] is Vhu in component model_2 (dimensionless). * ALGEBRAIC[93] is ATPi_ADP in component model_2 (dimensionless). * ALGEBRAIC[94] is ADPm_ATPm in component model_2 (dimensionless). * ALGEBRAIC[95] is VANT in component model_2 (dimensionless). * ALGEBRAIC[96] is Vhleak in component model_2 (dimensionless). * ALGEBRAIC[2] is MAlpha in component model_2 (dimensionless). * ALGEBRAIC[8] is MBeta in component model_2 (dimensionless). * ALGEBRAIC[13] is inv_MBeta_MAlpha in component model_2 (dimensionless). * ALGEBRAIC[17] is tmNa in component model_2 (dimensionless). * ALGEBRAIC[3] is HAlpha in component model_2 (dimensionless). * ALGEBRAIC[9] is HBeta in component model_2 (dimensionless). * ALGEBRAIC[4] is JAlpha in component model_2 (dimensionless). * ALGEBRAIC[10] is JBeta in component model_2 (dimensionless). * ALGEBRAIC[97] is fb in component model_2 (dimensionless). * ALGEBRAIC[98] is rb in component model_2 (per_ms). * STATES[26] is CaNSR in component model_2 (dimensionless). * ALGEBRAIC[99] is Jup in component model_2 (dimensionless). * ALGEBRAIC[101] is Jrel in component model_2 (dimensionless). * STATES[27] is CaJSR in component model_2 (dimensionless). * STATES[28] is CaSS in component model_2 (dimensionless). * ALGEBRAIC[102] is Jtr in component model_2 (dimensionless). * ALGEBRAIC[105] is Jxfer in component model_2 (dimensionless). * ALGEBRAIC[108] is P1_N1_P2_P3 in component model_2 (dimensionless). * STATES[29] is N1 in component model_2 (dimensionless). * STATES[30] is P0 in component model_2 (dimensionless). * STATES[31] is P1 in component model_2 (dimensionless). * STATES[32] is P2 in component model_2 (dimensionless). * STATES[33] is P3 in component model_2 (dimensionless). * ALGEBRAIC[24] is N1_differential in component model_2 (dimensionless). * ALGEBRAIC[11] is P0_differential in component model_2 (dimensionless). * ALGEBRAIC[14] is P1_differential in component model_2 (dimensionless). * ALGEBRAIC[18] is P2_differential in component model_2 (dimensionless). * ALGEBRAIC[21] is P3_differential in component model_2 (dimensionless). * ALGEBRAIC[111] is FN_Ca in component model_2 (dimensionless). * ALGEBRAIC[109] is force_norm in component model_2 (dimensionless). * ALGEBRAIC[110] is force in component model_2 (dimensionless). * ALGEBRAIC[5] is kTrop_np in component model_2 (dimensionless). * STATES[34] is LTRPNCa in component model_2 (dimensionless). * STATES[35] is HTRPNCa in component model_2 (dimensionless). * ALGEBRAIC[112] is LTRPNCa_differential in component model_2 (dimensionless). * ALGEBRAIC[113] is HTRPNCa_differential in component model_2 (dimensionless). * ALGEBRAIC[114] is Jtrpn in component model_2 (dimensionless). * ALGEBRAIC[106] is beta_SS in component model_2 (dimensionless). * ALGEBRAIC[107] is beta_JSR in component model_2 (dimensionless). * ALGEBRAIC[115] is beta_i in component model_2 (dimensionless). * ALGEBRAIC[116] is Cai_ktrans_plus1 in component model_2 (dimensionless). * ALGEBRAIC[117] is Cai_ktrans_plus1_p3 in component model_2 (dimensionless). * ALGEBRAIC[118] is Vuni in component model_2 (dimensionless). * ALGEBRAIC[119] is VnaCa in component model_2 (dimensionless). * ALGEBRAIC[6] is alpha in component model_2 (dimensionless). * ALGEBRAIC[12] is beta in component model_2 (dimensionless). * ALGEBRAIC[16] is alpha_prime in component model_2 (dimensionless). * ALGEBRAIC[20] is beta_prime in component model_2 (dimensionless). * ALGEBRAIC[23] is C0_to_C1 in component model_2 (dimensionless). * ALGEBRAIC[26] is C1_to_C2 in component model_2 (dimensionless). * ALGEBRAIC[28] is C2_to_C3 in component model_2 (dimensionless). * ALGEBRAIC[30] is C3_to_C4 in component model_2 (dimensionless). * ALGEBRAIC[32] is CCa0_to_CCa1 in component model_2 (dimensionless). * ALGEBRAIC[34] is CCa1_to_CCa2 in component model_2 (dimensionless). * ALGEBRAIC[36] is CCa2_to_CCa3 in component model_2 (dimensionless). * ALGEBRAIC[38] is CCa3_to_CCa4 in component model_2 (dimensionless). * ALGEBRAIC[40] is C1_to_C0 in component model_2 (dimensionless). * ALGEBRAIC[42] is C2_to_C1 in component model_2 (dimensionless). * ALGEBRAIC[44] is C3_to_C2 in component model_2 (dimensionless). * ALGEBRAIC[46] is C4_to_C3 in component model_2 (dimensionless). * ALGEBRAIC[48] is CCa1_to_CCa0 in component model_2 (dimensionless). * ALGEBRAIC[50] is CCa2_to_CCa1 in component model_2 (dimensionless). * ALGEBRAIC[52] is CCa3_to_CCa2 in component model_2 (dimensionless). * ALGEBRAIC[54] is CCa4_to_CCa3 in component model_2 (dimensionless). * ALGEBRAIC[56] is gamma in component model_2 (dimensionless). * ALGEBRAIC[58] is C0_to_CCa0 in component model_2 (dimensionless). * ALGEBRAIC[60] is C1_to_CCa1 in component model_2 (dimensionless). * ALGEBRAIC[62] is C2_to_CCa2 in component model_2 (dimensionless). * ALGEBRAIC[64] is C3_to_CCa3 in component model_2 (dimensionless). * ALGEBRAIC[66] is C4_to_CCa4 in component model_2 (dimensionless). * CONSTANTS[311] is CCa0_to_C0 in component model_2 (dimensionless). * CONSTANTS[315] is CCa1_to_C1 in component model_2 (dimensionless). * CONSTANTS[320] is CCa2_to_C2 in component model_2 (dimensionless). * CONSTANTS[325] is CCa3_to_C3 in component model_2 (dimensionless). * CONSTANTS[328] is CCa4_to_C4 in component model_2 (dimensionless). * STATES[36] is C1 in component model_2 (dimensionless). * STATES[37] is CCa0 in component model_2 (dimensionless). * STATES[38] is C0 in component model_2 (dimensionless). * STATES[39] is C2 in component model_2 (dimensionless). * STATES[40] is CCa1 in component model_2 (dimensionless). * STATES[41] is C3 in component model_2 (dimensionless). * STATES[42] is CCa2 in component model_2 (dimensionless). * STATES[43] is C4 in component model_2 (dimensionless). * STATES[44] is CCa3 in component model_2 (dimensionless). * STATES[45] is CCa4 in component model_2 (dimensionless). * ALGEBRAIC[100] is Vt_CRP2 in component model_2 (dimensionless). * STATES[46] is CrPi_mito in component model_2 (dimensionless). * STATES[47] is CrPi_cyto in component model_2 (dimensionless). * ALGEBRAIC[103] is VCK_cyto in component model_2 (dimensionless). * STATES[48] is ATPi_cyto in component model_2 (dimensionless). * ALGEBRAIC[104] is VCK_mito in component model_2 (dimensionless). * STATES[49] is ASP in component model_2 (dimensionless). * STATES[50] is N0 in component model_2 (dimensionless). * ALGEBRAIC[123] is reversal_potential in component EK (dimensionless). * ALGEBRAIC[120] is reversal_potential in component ENa (dimensionless). * ALGEBRAIC[132] is reversal_potential in component EKs (dimensionless). * ALGEBRAIC[130] is reversal_potential in component ECa (dimensionless). * RATES[5] is d/dt mNa in component model_2 (dimensionless). * RATES[6] is d/dt hNa in component model_2 (dimensionless). * RATES[7] is d/dt jNa in component model_2 (dimensionless). * RATES[8] is d/dt xKs in component model_2 (dimensionless). * RATES[30] is d/dt P0 in component model_2 (dimensionless). * RATES[31] is d/dt P1 in component model_2 (dimensionless). * RATES[32] is d/dt P2 in component model_2 (dimensionless). * RATES[33] is d/dt P3 in component model_2 (dimensionless). * RATES[29] is d/dt N1 in component model_2 (dimensionless). * RATES[50] is d/dt N0 in component model_2 (dimensionless). * RATES[34] is d/dt LTRPNCa in component model_2 (dimensionless). * RATES[35] is d/dt HTRPNCa in component model_2 (dimensionless). * RATES[13] is d/dt Nai in component model_2 (dimensionless). * RATES[12] is d/dt Ki in component model_2 (dimensionless). * RATES[14] is d/dt Cai in component model_2 (dimensionless). * RATES[28] is d/dt CaSS in component model_2 (dimensionless). * RATES[27] is d/dt CaJSR in component model_2 (dimensionless). * RATES[26] is d/dt CaNSR in component model_2 (dimensionless). * RATES[1] is d/dt V in component model_2 (dimensionless). * RATES[2] is d/dt C1_RyR in component model_2 (dimensionless). * RATES[4] is d/dt O2_RyR in component model_2 (dimensionless). * RATES[3] is d/dt C2_RyR in component model_2 (dimensionless). * RATES[38] is d/dt C0 in component model_2 (dimensionless). * RATES[36] is d/dt C1 in component model_2 (dimensionless). * RATES[39] is d/dt C2 in component model_2 (dimensionless). * RATES[41] is d/dt C3 in component model_2 (dimensionless). * RATES[43] is d/dt C4 in component model_2 (dimensionless). * RATES[9] is d/dt Open in component model_2 (dimensionless). * RATES[37] is d/dt CCa0 in component model_2 (dimensionless). * RATES[40] is d/dt CCa1 in component model_2 (dimensionless). * RATES[42] is d/dt CCa2 in component model_2 (dimensionless). * RATES[44] is d/dt CCa3 in component model_2 (dimensionless). * RATES[45] is d/dt CCa4 in component model_2 (dimensionless). * RATES[11] is d/dt yCa in component model_2 (dimensionless). * RATES[10] is d/dt OCa in component model_2 (dimensionless). * RATES[46] is d/dt CrPi_mito in component model_2 (dimensionless). * RATES[47] is d/dt CrPi_cyto in component model_2 (dimensionless). * RATES[0] is d/dt ATPi in component model_2 (dimensionless). * RATES[48] is d/dt ATPi_cyto in component model_2 (dimensionless). * RATES[25] is d/dt Cam in component model_2 (dimensionless). * RATES[15] is d/dt ADPm in component model_2 (dimensionless). * RATES[16] is d/dt Dpsi in component model_2 (dimensionless). * RATES[17] is d/dt NADH in component model_2 (dimensionless). * RATES[24] is d/dt ISOC in component model_2 (dimensionless). * RATES[19] is d/dt AKG in component model_2 (dimensionless). * RATES[20] is d/dt SCoA in component model_2 (dimensionless). * RATES[21] is d/dt Succ in component model_2 (dimensionless). * RATES[22] is d/dt FUM in component model_2 (dimensionless). * RATES[23] is d/dt MAL in component model_2 (dimensionless). * RATES[18] is d/dt Oaa in component model_2 (dimensionless). * RATES[49] is d/dt ASP in component model_2 (dimensionless). */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 4e3; CONSTANTS[1] = 0.5; CONSTANTS[2] = -100.0; CONSTANTS[3] = 0.0; CONSTANTS[4] = 0.0; CONSTANTS[5] = 96.5; CONSTANTS[6] = 8.314; CONSTANTS[7] = 310.0; CONSTANTS[8] = 1e-7; CONSTANTS[9] = 2.0e-3; CONSTANTS[10] = 1.4e-4; CONSTANTS[11] = 1.33e-6; CONSTANTS[12] = 0.0095; CONSTANTS[13] = 25.0; CONSTANTS[14] = 25.0; CONSTANTS[15] = 1.0E-5; CONSTANTS[16] = 1.534E-4; CONSTANTS[17] = 1.0; CONSTANTS[18] = 2.0; CONSTANTS[19] = 0.5; CONSTANTS[20] = 2.0; CONSTANTS[21] = 1.0; CONSTANTS[22] = 2.0; CONSTANTS[23] = 1.0; CONSTANTS[24] = 1.5; CONSTANTS[25] = 50E-3; CONSTANTS[26] = 1.812E-3; CONSTANTS[27] = 10.0; CONSTANTS[28] = 0.02; CONSTANTS[29] = 35.0; CONSTANTS[30] = -0.6; CONSTANTS[31] = 50.0; CONSTANTS[32] = 0.35; CONSTANTS[33] = 0.4; CONSTANTS[34] = 0.109; CONSTANTS[35] = 0.5; CONSTANTS[36] = 0.154; CONSTANTS[37] = 0.5; CONSTANTS[38] = 0.01; CONSTANTS[39] = 1.24; CONSTANTS[40] = 0.3; CONSTANTS[41] = 0.0003; CONSTANTS[42] = 0.0; CONSTANTS[43] = 0.85; CONSTANTS[44] = 0.003217; CONSTANTS[45] = 0.00828; CONSTANTS[46] = 12.8; CONSTANTS[47] = 0.000545; CONSTANTS[48] = 2.0E-7; CONSTANTS[49] = 2.0; CONSTANTS[50] = 10.0; CONSTANTS[51] = 0.0; CONSTANTS[52] = 2.5E-5; CONSTANTS[53] = 0.5; CONSTANTS[54] = 140E-3; CONSTANTS[55] = -0.4583; CONSTANTS[56] = 3.147; CONSTANTS[57] = 0.575; CONSTANTS[58] = 6.6; CONSTANTS[59] = 0.0005; CONSTANTS[60] = 2.22; CONSTANTS[61] = 0.00038; CONSTANTS[62] = 0.62; CONSTANTS[63] = 0.576; CONSTANTS[64] = 1.215E10; CONSTANTS[65] = 1.930; CONSTANTS[66] = 4.05E6; CONSTANTS[67] = 0.00127; CONSTANTS[68] = 0.0008; CONSTANTS[69] = 1.5E-6; CONSTANTS[70] = 0.10; CONSTANTS[71] = 0.5; CONSTANTS[72] = 1.71E6; CONSTANTS[73] = 6.44E-4; CONSTANTS[74] = 0.0125; CONSTANTS[75] = 0.00024; CONSTANTS[76] = 0.00332; CONSTANTS[77] = 1.0; CONSTANTS[78] = 0.005; CONSTANTS[79] = 8.1E-5; CONSTANTS[80] = 5.98E-5; CONSTANTS[81] = 1.131E-5; CONSTANTS[82] = 26.7; CONSTANTS[83] = 6.68E-9; CONSTANTS[84] = 5.62E-6; CONSTANTS[85] = 3.3E-4; CONSTANTS[86] = 100.0; CONSTANTS[87] = 0.26; CONSTANTS[88] = 5.1; CONSTANTS[89] = 0.14; CONSTANTS[90] = 0.1; CONSTANTS[91] = 1.0; CONSTANTS[92] = 0.05; CONSTANTS[93] = 0.19; CONSTANTS[94] = 1.3; CONSTANTS[95] = 0.0031; CONSTANTS[96] = 0.15; CONSTANTS[97] = 7.5E-2; CONSTANTS[98] = 100.0; CONSTANTS[99] = 40E-3; CONSTANTS[100] = 0.008; CONSTANTS[101] = 0.012; CONSTANTS[102] = 0.23; CONSTANTS[103] = 0.0126; CONSTANTS[104] = 1.493; CONSTANTS[105] = 0.03; CONSTANTS[106] = 0.01; CONSTANTS[107] = 2.38E-3; CONSTANTS[108] = 1.38; CONSTANTS[109] = 0.8; CONSTANTS[110] = 0.111; CONSTANTS[111] = 0.0308; CONSTANTS[112] = 0.923; CONSTANTS[113] = 1.52; CONSTANTS[114] = 1.94; CONSTANTS[115] = 38.7; CONSTANTS[116] = 1.5; CONSTANTS[117] = 0.2244; CONSTANTS[118] = 87.5; CONSTANTS[119] = 10.0; CONSTANTS[120] = 0.0012; CONSTANTS[121] = 0.00064; CONSTANTS[122] = 0.0005; CONSTANTS[123] = 0.03; CONSTANTS[124] = 9.4; CONSTANTS[125] = 9000.0; CONSTANTS[126] = 3.75E-4; CONSTANTS[127] = 5.4; CONSTANTS[128] = 0.0399; CONSTANTS[129] = 1.64269; CONSTANTS[130] = 1.35E18; CONSTANTS[131] = 5.765E13; CONSTANTS[132] = 0.1; CONSTANTS[133] = 0.005; CONSTANTS[134] = 3.115; CONSTANTS[135] = 2.5; CONSTANTS[136] = 0.019; CONSTANTS[137] = 0.04; CONSTANTS[138] = 110.0; CONSTANTS[139] = 70E-3; CONSTANTS[140] = 3.0; CONSTANTS[141] = 0.4; CONSTANTS[142] = 3.0; CONSTANTS[143] = 2.8; CONSTANTS[144] = 140.0; CONSTANTS[145] = 4.0; CONSTANTS[146] = 1.4; CONSTANTS[147] = 2.0; CONSTANTS[148] = 1.2; CONSTANTS[149] = 1.0; CONSTANTS[150] = 0.01; CONSTANTS[151] = 1.346E-8; CONSTANTS[152] = 7.739E-7; CONSTANTS[153] = 6.65E-15; CONSTANTS[154] = 1.656E-8; CONSTANTS[155] = 3.373E-10; CONSTANTS[156] = 9.651E-17; CONSTANTS[157] = 4.585E-17; CONSTANTS[158] = 0.00124; CONSTANTS[159] = 2.0; CONSTANTS[160] = 1.11E-11; CONSTANTS[161] = 0.0; CONSTANTS[162] = 1.75E-7; CONSTANTS[163] = 2.077E-18; CONSTANTS[164] = 1.728E-9; CONSTANTS[165] = 1.059E-26; CONSTANTS[166] = 6.394E-13; CONSTANTS[167] = 1.762E-16; CONSTANTS[168] = 2.656E-22; CONSTANTS[169] = 8.632E-30; CONSTANTS[170] = 0.05; CONSTANTS[171] = 3.75E-4; CONSTANTS[172] = 1.0E-1; CONSTANTS[173] = 0.574713; CONSTANTS[174] = 9.09091; CONSTANTS[175] = 15.0; CONSTANTS[176] = 0.00048; CONSTANTS[177] = 3.6; CONSTANTS[178] = 0.16E-6; CONSTANTS[179] = 0.29892E-3; CONSTANTS[180] = 0.31793E-3; CONSTANTS[181] = 0.015; CONSTANTS[182] = 0.625E-4; CONSTANTS[183] = 0.0275; CONSTANTS[184] = 25.84E-6; CONSTANTS[185] = 1.4E-6; CONSTANTS[186] = 0.495E-9; CONSTANTS[187] = 0.1; CONSTANTS[188] = 0.05; CONSTANTS[189] = 2.15; CONSTANTS[190] = 0.1; CONSTANTS[191] = 1; CONSTANTS[192] = 1; CONSTANTS[193] = 2; CONSTANTS[194] = 1; CONSTANTS[195] = 0.01833; CONSTANTS[196] = 7.98691; CONSTANTS[197] = -86.4154; CONSTANTS[198] = 0.748798; CONSTANTS[199] = 0.250904; CONSTANTS[200] = 3.5175e-10; CONSTANTS[201] = 0.0327969; CONSTANTS[202] = 0.988991; CONSTANTS[203] = 0.992918; CONSTANTS[204] = 0.000127552; CONSTANTS[205] = 6.96577e-24; CONSTANTS[206] = 0; CONSTANTS[207] = 0.941505; CONSTANTS[208] = 148.065; CONSTANTS[209] = 7.51215; CONSTANTS[210] = 6.71512e-05; CONSTANTS[211] = 0.0258454; CONSTANTS[212] = 164.073; CONSTANTS[213] = 7.76311; CONSTANTS[214] = 1.08962e-07; CONSTANTS[215] = 0.000610554; CONSTANTS[216] = 0.247364; CONSTANTS[217] = 0.000394005; CONSTANTS[218] = 0.0191924; CONSTANTS[219] = 0.00925885; CONSTANTS[220] = 0.496682; CONSTANTS[221] = 0.000137836; CONSTANTS[222] = 0.255288; CONSTANTS[223] = 0.255138; CONSTANTS[224] = 8.25461e-05; CONSTANTS[225] = 5.71121e-05; CONSTANTS[226] = 6.63724e-05; CONSTANTS[227] = 5.73114e-05; CONSTANTS[228] = 0.000106992; CONSTANTS[229] = 9.31884e-05; CONSTANTS[230] = 0.0100663; CONSTANTS[231] = 0.133453; CONSTANTS[232] = 1.04295e-05; CONSTANTS[233] = 0.00156783; CONSTANTS[234] = 0.998424; CONSTANTS[235] = 4.08551e-11; CONSTANTS[236] = 6.55103e-08; CONSTANTS[237] = 7.11287e-17; CONSTANTS[238] = 1.02648e-12; CONSTANTS[239] = 4.64382e-23; CONSTANTS[240] = 7.14846e-18; CONSTANTS[241] = 1.86682e-23; CONSTANTS[242] = 19.4385; CONSTANTS[243] = 19.4302; CONSTANTS[244] = 7.97824; CONSTANTS[245] = 0.0502614; CONSTANTS[246] = 0.999618; STATES[0] = 7.98691; STATES[1] = -86.4154; STATES[2] = 0.748798; STATES[3] = 0.250904; STATES[4] = 3.5175e-10; STATES[5] = 0.0327969; STATES[6] = 0.988991; STATES[7] = 0.992918; STATES[8] = 0.000127552; STATES[9] = 6.96577e-24; STATES[10] = 0; STATES[11] = 0.941505; STATES[12] = 148.065; STATES[13] = 7.51215; STATES[14] = 6.71512e-05; STATES[15] = 0.0258454; STATES[16] = 164.073; STATES[17] = 7.76311; STATES[18] = 1.08962e-07; STATES[19] = 0.000610554; STATES[20] = 0.247364; STATES[21] = 0.000394005; STATES[22] = 0.0191924; STATES[23] = 0.00925885; STATES[24] = 0.496682; STATES[25] = 0.000137836; STATES[26] = 0.255288; STATES[27] = 0.255138; STATES[28] = 8.25461e-05; STATES[29] = 5.71121e-05; STATES[30] = 6.63724e-05; STATES[31] = 5.73114e-05; STATES[32] = 0.000106992; STATES[33] = 9.31884e-05; STATES[34] = 0.0100663; STATES[35] = 0.133453; STATES[36] = 1.04295e-05; STATES[37] = 0.00156783; STATES[38] = 0.998424; STATES[39] = 4.08551e-11; STATES[40] = 6.55103e-08; STATES[41] = 7.11287e-17; STATES[42] = 1.02648e-12; STATES[43] = 4.64382e-23; STATES[44] = 7.14846e-18; STATES[45] = 1.86682e-23; STATES[46] = 19.4385; STATES[47] = 19.4302; STATES[48] = 7.97824; STATES[49] = 0.0502614; STATES[50] = 0.999618; CONSTANTS[247] = (CONSTANTS[184]+CONSTANTS[178]+CONSTANTS[185]+CONSTANTS[186])/0.640000; CONSTANTS[248] = 3.00000*CONSTANTS[188]; CONSTANTS[249] = (CONSTANTS[189] - 1.70000)/0.600000; CONSTANTS[250] = CONSTANTS[99]/CONSTANTS[98]; CONSTANTS[251] = 3.50000*CONSTANTS[189] - 2.00000; CONSTANTS[252] = 0.282000* pow((CONSTANTS[127]/5.40000), 1.0 / 2); CONSTANTS[253] = 1.00000/5.98000; CONSTANTS[254] = CONSTANTS[22]*341.000; CONSTANTS[255] = 2.00000*CONSTANTS[158]*1000.00*CONSTANTS[5]*(1.00000 - 341.000*CONSTANTS[22]); CONSTANTS[256] = 4.00000*0.00100000*CONSTANTS[158]; CONSTANTS[257] = 1.00000/CONSTANTS[55]; CONSTANTS[258] = 1.00000/CONSTANTS[119]; CONSTANTS[259] = CONSTANTS[32] - 1.00000; CONSTANTS[260] = pow(CONSTANTS[144], 3.00000)/CONSTANTS[22]; CONSTANTS[261] = (CONSTANTS[108]+CONSTANTS[22])*((pow(CONSTANTS[118], 3.00000)+pow(CONSTANTS[144], 3.00000))/( CONSTANTS[125]*CONSTANTS[22])); CONSTANTS[262] = 1.00000/CONSTANTS[91]; CONSTANTS[263] = pow(CONSTANTS[120], 3.00000); CONSTANTS[264] = CONSTANTS[105]/CONSTANTS[87]; CONSTANTS[265] = ( CONSTANTS[71]*CONSTANTS[33]*CONSTANTS[17])/(CONSTANTS[103]+CONSTANTS[17]); CONSTANTS[266] = 1.00000+1.00000/CONSTANTS[60]; CONSTANTS[267] = 1.00000+CONSTANTS[52]/CONSTANTS[79]+CONSTANTS[80]/CONSTANTS[52]; CONSTANTS[268] = CONSTANTS[92]*CONSTANTS[34]; CONSTANTS[269] = 1.00000/CONSTANTS[62]; CONSTANTS[270] = 1.00000/CONSTANTS[59]; CONSTANTS[271] = 1.00000/CONSTANTS[93]; CONSTANTS[272] = CONSTANTS[115]/CONSTANTS[112]; CONSTANTS[273] = CONSTANTS[141]/CONSTANTS[111]+1.00000; CONSTANTS[274] = CONSTANTS[97]*CONSTANTS[35]; CONSTANTS[275] = CONSTANTS[28]/CONSTANTS[134]; CONSTANTS[276] = CONSTANTS[133]*CONSTANTS[37]; CONSTANTS[277] = CONSTANTS[123]/CONSTANTS[94]; CONSTANTS[278] = 1.00000/CONSTANTS[96]; CONSTANTS[279] = CONSTANTS[76]/CONSTANTS[77]; CONSTANTS[280] = pow(1.00000/(1.00000+CONSTANTS[83]/CONSTANTS[52]+ CONSTANTS[83]*(CONSTANTS[84]/pow(CONSTANTS[52], 2.00000))), 2.00000)*(1.00000/(1.00000+CONSTANTS[52]/CONSTANTS[81]+pow(CONSTANTS[52], 2.00000)/( CONSTANTS[81]*CONSTANTS[82]))+CONSTANTS[128])*CONSTANTS[110]*CONSTANTS[36]; CONSTANTS[281] = CONSTANTS[104]/CONSTANTS[95]; CONSTANTS[282] = CONSTANTS[73]*CONSTANTS[50]*CONSTANTS[69]*(CONSTANTS[58]/CONSTANTS[73]); CONSTANTS[283] = CONSTANTS[69]*(CONSTANTS[58]/CONSTANTS[73]); CONSTANTS[284] = CONSTANTS[73]*CONSTANTS[50]; CONSTANTS[285] = CONSTANTS[73]/CONSTANTS[58]; CONSTANTS[286] = CONSTANTS[130]*(CONSTANTS[130]/CONSTANTS[112]); CONSTANTS[287] = CONSTANTS[72]/CONSTANTS[159]; CONSTANTS[288] = 0.750000*CONSTANTS[181]; CONSTANTS[289] = 20.0000*CONSTANTS[181]; CONSTANTS[290] = 1.00000/CONSTANTS[75]; CONSTANTS[291] = 1.00000/CONSTANTS[129]; CONSTANTS[292] = 1.00000/CONSTANTS[173]; CONSTANTS[293] = 1.00000/CONSTANTS[174]; CONSTANTS[294] = CONSTANTS[25]*CONSTANTS[107]; CONSTANTS[295] = CONSTANTS[29]*CONSTANTS[109]; CONSTANTS[296] = 1.00000/CONSTANTS[136]; CONSTANTS[297] = 1.00000/CONSTANTS[61]; CONSTANTS[298] = CONSTANTS[183]/CONSTANTS[136]; CONSTANTS[299] = CONSTANTS[19]*0.500000; CONSTANTS[300] = 1.00000/CONSTANTS[12]; CONSTANTS[301] = ( CONSTANTS[6]*CONSTANTS[7])/CONSTANTS[5]; CONSTANTS[302] = 1.00000/CONSTANTS[20]; CONSTANTS[303] = CONSTANTS[247]*0.360000; CONSTANTS[304] = 10.0000*CONSTANTS[188]; CONSTANTS[305] = 1.00000/(1.00000+CONSTANTS[250]/(1.70000/1000.00+ ((0.900000/1000.00 - 1.70000/1000.00)/(2.30000 - 1.70000))*(CONSTANTS[189] - 1.70000))); CONSTANTS[306] = CONSTANTS[261]*CONSTANTS[132]; CONSTANTS[307] = -2.30300*CONSTANTS[301]*CONSTANTS[30]; CONSTANTS[308] = CONSTANTS[273]/CONSTANTS[67]; CONSTANTS[309] = - CONSTANTS[137] - CONSTANTS[248]; CONSTANTS[310] = 0.750000* pow((CONSTANTS[127]/5.40000), 1.0 / 2); CONSTANTS[311] = CONSTANTS[150]; CONSTANTS[312] = 7.00000*CONSTANTS[188]; CONSTANTS[313] = 1.00000/( CONSTANTS[139]*CONSTANTS[305]); CONSTANTS[314] = 1000.00*CONSTANTS[5]; CONSTANTS[315] = CONSTANTS[311]*CONSTANTS[302]; CONSTANTS[316] = 1.00000*CONSTANTS[190]; CONSTANTS[317] = CONSTANTS[314]*CONSTANTS[160]; CONSTANTS[318] = CONSTANTS[175]*(CONSTANTS[176]/(CONSTANTS[248]+CONSTANTS[304]+CONSTANTS[312])); CONSTANTS[319] = 1.00000/CONSTANTS[301]; CONSTANTS[320] = CONSTANTS[315]*CONSTANTS[302]; CONSTANTS[321] = 2.00000*CONSTANTS[190]; CONSTANTS[322] = CONSTANTS[53]*CONSTANTS[319]; CONSTANTS[323] = 2.00000*CONSTANTS[319]; CONSTANTS[324] = 0.0365000*((exp(CONSTANTS[144]/67.3000) - 1.00000)/7.00000); CONSTANTS[325] = CONSTANTS[320]*CONSTANTS[302]; CONSTANTS[326] = 3.00000*CONSTANTS[190]; CONSTANTS[327] = CONSTANTS[56]*(CONSTANTS[127]/(CONSTANTS[127]+CONSTANTS[116])); CONSTANTS[328] = CONSTANTS[325]*CONSTANTS[302]; CONSTANTS[329] = CONSTANTS[316]*CONSTANTS[321]*CONSTANTS[326]+ CONSTANTS[248]*CONSTANTS[321]*CONSTANTS[326]+ CONSTANTS[248]*CONSTANTS[304]*CONSTANTS[326]+ CONSTANTS[248]*CONSTANTS[304]*CONSTANTS[312]; CONSTANTS[330] = 1.00000+(2.30000 - CONSTANTS[189])/pow((2.30000 - 1.70000)/1.00000, 1.60000); CONSTANTS[331] = 1.00000/CONSTANTS[90]; CONSTANTS[332] = ( CONSTANTS[248]*( 2.00000*CONSTANTS[190])*( 3.00000*CONSTANTS[190]))/CONSTANTS[329]; CONSTANTS[333] = CONSTANTS[316]*CONSTANTS[330]; CONSTANTS[334] = CONSTANTS[321]*CONSTANTS[330]; CONSTANTS[335] = CONSTANTS[326]*CONSTANTS[330]; CONSTANTS[336] = 30.0000/1000.00; CONSTANTS[337] = exp( 6.00000*CONSTANTS[31]*CONSTANTS[319]); CONSTANTS[338] = ( CONSTANTS[248]*CONSTANTS[304]*( 3.00000*CONSTANTS[190]))/CONSTANTS[329]; CONSTANTS[339] = CONSTANTS[336]*CONSTANTS[330]; CONSTANTS[340] = - (CONSTANTS[137]+CONSTANTS[304]+CONSTANTS[333]); CONSTANTS[341] = - (CONSTANTS[312]+CONSTANTS[334]); CONSTANTS[342] = 6.00000*CONSTANTS[43]*CONSTANTS[319]; CONSTANTS[343] = ( CONSTANTS[248]*CONSTANTS[304]*CONSTANTS[312])/CONSTANTS[329]; CONSTANTS[344] = CONSTANTS[166]+ CONSTANTS[168]*CONSTANTS[337]; CONSTANTS[345] = CONSTANTS[332]+ 2.00000*CONSTANTS[338]+ 3.00000*CONSTANTS[343]; CONSTANTS[346] = CONSTANTS[163]*CONSTANTS[337]; CONSTANTS[347] = CONSTANTS[345]/3.00000; CONSTANTS[348] = 0.500000*CONSTANTS[172]*CONSTANTS[344]; CONSTANTS[349] = CONSTANTS[332]+CONSTANTS[338]+CONSTANTS[343]; CONSTANTS[350] = 0.500000*CONSTANTS[172]*CONSTANTS[169]; CONSTANTS[351] = 1.00000; CONSTANTS[352] = 0.500000*CONSTANTS[172]*CONSTANTS[166]; CONSTANTS[353] = 1.50000; CONSTANTS[354] = 6.00000*CONSTANTS[172]*CONSTANTS[166]; CONSTANTS[355] = 0.100000; CONSTANTS[356] = 6.00000*CONSTANTS[172]*(CONSTANTS[166]+CONSTANTS[167]); CONSTANTS[357] = 0.100000; CONSTANTS[358] = CONSTANTS[301]*arbitrary_log( CONSTANTS[131]* pow((CONSTANTS[39]/CONSTANTS[38]), 1.0 / 2), 10); CONSTANTS[359] = CONSTANTS[353] - CONSTANTS[357]; CONSTANTS[360] = exp( CONSTANTS[358]*CONSTANTS[319]); CONSTANTS[361] = 0.500000*(CONSTANTS[166]+ CONSTANTS[169]*CONSTANTS[360]); CONSTANTS[362] = (CONSTANTS[166]+ CONSTANTS[168]*CONSTANTS[337])*CONSTANTS[360]*0.500000; CONSTANTS[363] = (CONSTANTS[189]<2.20000 ? multi_min(2, 1.00000, ((CONSTANTS[189] - 2.00000*CONSTANTS[351])+(CONSTANTS[353] - CONSTANTS[355]))/CONSTANTS[359]) : 1.00000 - (CONSTANTS[189] - 2.20000)/CONSTANTS[359]); CONSTANTS[364] = 4.00000*CONSTANTS[166]*CONSTANTS[360]; CONSTANTS[365] = CONSTANTS[363]/CONSTANTS[349]; CONSTANTS[366] = CONSTANTS[363]/( CONSTANTS[347]*3.00000); CONSTANTS[367] = 4.00000*(CONSTANTS[166]+CONSTANTS[167]); CONSTANTS[368] = CONSTANTS[187]*CONSTANTS[366]; CONSTANTS[369] = (1.00000+ CONSTANTS[163]*CONSTANTS[360])*CONSTANTS[337]; CONSTANTS[370] = CONSTANTS[164]+ CONSTANTS[165]*CONSTANTS[360]; CONSTANTS[371] = exp( 3.00000*CONSTANTS[31]*CONSTANTS[319]); CONSTANTS[372] = 3.00000*CONSTANTS[319]; CONSTANTS[373] = 100.000*CONSTANTS[154]+ CONSTANTS[156]*CONSTANTS[371]; CONSTANTS[374] = 3.00000*(CONSTANTS[154]+CONSTANTS[155]); CONSTANTS[375] = 300.000*CONSTANTS[154]; CONSTANTS[376] = CONSTANTS[151]*CONSTANTS[371]; CONSTANTS[377] = CONSTANTS[106]/CONSTANTS[89]; CONSTANTS[378] = 1.00000/CONSTANTS[88]; CONSTANTS[379] = CONSTANTS[16]/( CONSTANTS[184]*CONSTANTS[5]*1000.00); CONSTANTS[380] = CONSTANTS[16]/( 2.00000*CONSTANTS[186]*CONSTANTS[5]*1000.00); CONSTANTS[381] = CONSTANTS[178]/CONSTANTS[186]; CONSTANTS[382] = CONSTANTS[184]/CONSTANTS[186]; CONSTANTS[383] = CONSTANTS[184]/CONSTANTS[185]; CONSTANTS[384] = CONSTANTS[178]/CONSTANTS[185]; CONSTANTS[385] = 1.00000/CONSTANTS[21]; CONSTANTS[386] = 1.00000/CONSTANTS[26]; CONSTANTS[387] = 2.00000*CONSTANTS[19]; } void computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { RATES[9] = CONSTANTS[40]*STATES[43] - CONSTANTS[49]*STATES[9]; RATES[11] = (1.00000/(1.00000+exp( (STATES[1]+55.0000)*(1.00000/7.50000)))+(0.500000/(1.00000+exp( (21.0000 - STATES[1])*(1.00000/6.00000))) - STATES[11]))/(20.0000+600.000/(1.00000+exp( (STATES[1]+30.0000)*(1.00000/9.50000)))); RATES[10] = CONSTANTS[42]*STATES[45] - CONSTANTS[51]*STATES[10]; ALGEBRAIC[0] = STATES[1]+30.0000; RATES[8] = 7.19000e-05*(ALGEBRAIC[0]/(1.00000 - exp( -0.148000*ALGEBRAIC[0])))*(1.00000 - STATES[8]) - 0.000131000*(ALGEBRAIC[0]/(exp( 0.0687000*ALGEBRAIC[0]) - 1.00000))*STATES[8]; ALGEBRAIC[3] = (STATES[1]<-40.0000 ? ( 0.135000*exp(-80.0000/6.80000))*exp( (-1.00000/6.80000)*STATES[1]) : 0.00000); ALGEBRAIC[9] = (STATES[1]<-40.0000 ? 3.56000*exp( 0.0790000*STATES[1])+ 310000.*exp( 0.350000*STATES[1]) : 1.00000/(0.130000+ ( 0.130000*exp(- (10.6600/11.1000)))*exp( STATES[1]*(-1.00000/11.1000)))); RATES[6] = ALGEBRAIC[3]*(1.00000 - STATES[6]) - ALGEBRAIC[9]*STATES[6]; ALGEBRAIC[4] = (STATES[1]<-40.0000 ? ( -127140.*exp( 0.244400*STATES[1]) - 3.47400e-05*exp( -0.0439100*STATES[1]))*((STATES[1]+37.7800)/(1.00000+exp( 0.311000*(STATES[1]+79.2300)))) : 0.00000); ALGEBRAIC[10] = (STATES[1]<-40.0000 ? 0.121200*(exp( -0.0105200*STATES[1])/(1.00000+exp( -0.137800*(STATES[1]+40.1400)))) : 0.300000*(exp( -2.53500e-07*STATES[1])/(1.00000+exp( -0.100000*STATES[1] - 3.20000)))); RATES[7] = ALGEBRAIC[4]*(1.00000 - STATES[7]) - ALGEBRAIC[10]*STATES[7]; ALGEBRAIC[5] = CONSTANTS[137]*pow( STATES[34]*CONSTANTS[313], CONSTANTS[251]); ALGEBRAIC[11] = CONSTANTS[309]*STATES[30]+ ALGEBRAIC[5]*STATES[50]+ CONSTANTS[333]*STATES[31]; RATES[30] = ALGEBRAIC[11]; ALGEBRAIC[14] = CONSTANTS[340]*STATES[31]+ ALGEBRAIC[5]*STATES[29]+ CONSTANTS[248]*STATES[30]+ CONSTANTS[334]*STATES[32]; RATES[31] = ALGEBRAIC[14]; ALGEBRAIC[2] = (STATES[1]==-47.1300 ? 3.20000 : 0.320000*((STATES[1]+47.1300)/(1.00000 - exp( -0.100000*(STATES[1]+47.1300))))); ALGEBRAIC[8] = 0.0800000*exp(- ( STATES[1]*(1.00000/11.0000))); ALGEBRAIC[13] = 1.00000/(ALGEBRAIC[8]+ALGEBRAIC[2]); ALGEBRAIC[17] = (ALGEBRAIC[13]<0.0300000 ? ALGEBRAIC[2]*ALGEBRAIC[13] : STATES[5]); RATES[5] = ALGEBRAIC[2]*(1.00000 - ALGEBRAIC[17]) - ALGEBRAIC[8]*ALGEBRAIC[17]; ALGEBRAIC[18] = CONSTANTS[341]*STATES[32]+ CONSTANTS[304]*STATES[31]+ CONSTANTS[335]*STATES[33]; RATES[32] = ALGEBRAIC[18]; ALGEBRAIC[21] = - ( CONSTANTS[335]*STATES[33])+ CONSTANTS[312]*STATES[32]; RATES[33] = ALGEBRAIC[21]; ALGEBRAIC[24] = CONSTANTS[137]*STATES[31] - (ALGEBRAIC[5]+CONSTANTS[339])*STATES[29]; RATES[29] = ALGEBRAIC[24]; RATES[50] = - ALGEBRAIC[11] - (ALGEBRAIC[14]+ALGEBRAIC[18]+ALGEBRAIC[21]+ALGEBRAIC[24]); ALGEBRAIC[29] = 1.00000 - (STATES[2]+STATES[3]+STATES[4]); RATES[2] = - ( CONSTANTS[64]*pow(STATES[28], CONSTANTS[145])*STATES[2])+ CONSTANTS[63]*ALGEBRAIC[29]; RATES[4] = CONSTANTS[66]*pow(STATES[28], CONSTANTS[140])*ALGEBRAIC[29] - CONSTANTS[65]*STATES[4]; RATES[3] = CONSTANTS[70]*ALGEBRAIC[29] - CONSTANTS[68]*STATES[3]; ALGEBRAIC[6] = 0.400000*exp( (STATES[1]+2.00000)*0.100000); ALGEBRAIC[23] = 4.00000*ALGEBRAIC[6]; ALGEBRAIC[12] = 0.0500000*exp( (STATES[1]+2.00000)*(-1.00000/13.0000)); ALGEBRAIC[40] = ALGEBRAIC[12]; ALGEBRAIC[56] = 0.187500*STATES[28]; ALGEBRAIC[58] = ALGEBRAIC[56]; RATES[38] = ALGEBRAIC[40]*STATES[36]+( CONSTANTS[311]*STATES[37] - (ALGEBRAIC[23]+ALGEBRAIC[58])*STATES[38]); ALGEBRAIC[16] = CONSTANTS[18]*ALGEBRAIC[6]; ALGEBRAIC[32] = 4.00000*ALGEBRAIC[16]; ALGEBRAIC[20] = ALGEBRAIC[12]*CONSTANTS[302]; ALGEBRAIC[48] = ALGEBRAIC[20]; RATES[37] = ALGEBRAIC[48]*STATES[40]+( ALGEBRAIC[58]*STATES[38] - (ALGEBRAIC[32]+CONSTANTS[311])*STATES[37]); ALGEBRAIC[26] = 3.00000*ALGEBRAIC[6]; ALGEBRAIC[42] = 2.00000*ALGEBRAIC[12]; ALGEBRAIC[60] = CONSTANTS[18]*ALGEBRAIC[58]; RATES[36] = ALGEBRAIC[23]*STATES[38]+ ALGEBRAIC[42]*STATES[39]+( CONSTANTS[315]*STATES[40] - (ALGEBRAIC[40]+ALGEBRAIC[26]+ALGEBRAIC[60])*STATES[36]); ALGEBRAIC[34] = 3.00000*ALGEBRAIC[16]; ALGEBRAIC[50] = 2.00000*ALGEBRAIC[20]; RATES[40] = ALGEBRAIC[32]*STATES[37]+ ALGEBRAIC[50]*STATES[42]+( ALGEBRAIC[60]*STATES[36] - (ALGEBRAIC[48]+ALGEBRAIC[34]+CONSTANTS[315])*STATES[40]); ALGEBRAIC[28] = 2.00000*ALGEBRAIC[6]; ALGEBRAIC[44] = 3.00000*ALGEBRAIC[12]; ALGEBRAIC[62] = CONSTANTS[18]*ALGEBRAIC[60]; RATES[39] = ALGEBRAIC[26]*STATES[36]+ ALGEBRAIC[44]*STATES[41]+( CONSTANTS[320]*STATES[42] - (ALGEBRAIC[42]+ALGEBRAIC[28]+ALGEBRAIC[62])*STATES[39]); ALGEBRAIC[36] = 2.00000*ALGEBRAIC[16]; ALGEBRAIC[52] = 3.00000*ALGEBRAIC[20]; RATES[42] = ALGEBRAIC[34]*STATES[40]+ ALGEBRAIC[52]*STATES[44]+( ALGEBRAIC[62]*STATES[39] - (ALGEBRAIC[50]+ALGEBRAIC[36]+CONSTANTS[320])*STATES[42]); ALGEBRAIC[30] = ALGEBRAIC[6]; ALGEBRAIC[46] = 4.00000*ALGEBRAIC[12]; ALGEBRAIC[64] = CONSTANTS[18]*ALGEBRAIC[62]; RATES[41] = ALGEBRAIC[28]*STATES[39]+ ALGEBRAIC[46]*STATES[43]+( CONSTANTS[325]*STATES[44] - (ALGEBRAIC[44]+ALGEBRAIC[30]+ALGEBRAIC[64])*STATES[41]); ALGEBRAIC[38] = ALGEBRAIC[16]; ALGEBRAIC[54] = 4.00000*ALGEBRAIC[20]; RATES[44] = ALGEBRAIC[36]*STATES[42]+ ALGEBRAIC[54]*STATES[45]+( ALGEBRAIC[64]*STATES[41] - (ALGEBRAIC[52]+ALGEBRAIC[38]+CONSTANTS[325])*STATES[44]); ALGEBRAIC[66] = CONSTANTS[18]*ALGEBRAIC[64]; RATES[43] = ALGEBRAIC[30]*STATES[41]+ CONSTANTS[49]*STATES[9]+( CONSTANTS[328]*STATES[45] - (ALGEBRAIC[46]+CONSTANTS[40]+ALGEBRAIC[66])*STATES[43]); RATES[45] = ALGEBRAIC[38]*STATES[44]+ CONSTANTS[51]*STATES[10]+( ALGEBRAIC[66]*STATES[43] - (ALGEBRAIC[54]+CONSTANTS[42]+CONSTANTS[328])*STATES[45]); ALGEBRAIC[71] = CONSTANTS[74]*(CONSTANTS[23] - (STATES[19]+STATES[20]+STATES[21]+STATES[22]+STATES[23]+STATES[18]+ STATES[24]*CONSTANTS[266])); ALGEBRAIC[65] = CONSTANTS[27] - STATES[17]; ALGEBRAIC[67] = CONSTANTS[112]/ALGEBRAIC[65]; ALGEBRAIC[72] = 1.00000/( (1.00000+ STATES[15]*CONSTANTS[269])*(1.00000+ STATES[25]*CONSTANTS[270])); ALGEBRAIC[73] = 1.00000+ STATES[17]*CONSTANTS[271]; ALGEBRAIC[74] = CONSTANTS[268]/(CONSTANTS[267]+ ALGEBRAIC[67]*ALGEBRAIC[73]+ pow(CONSTANTS[113]/STATES[24], CONSTANTS[147])*ALGEBRAIC[72]*(1.00000+ ALGEBRAIC[67]*ALGEBRAIC[73])); RATES[24] = ALGEBRAIC[71] - ALGEBRAIC[74]; ALGEBRAIC[75] = CONSTANTS[273]+ CONSTANTS[308]*STATES[25]; ALGEBRAIC[76] = CONSTANTS[274]*(ALGEBRAIC[75]/(ALGEBRAIC[75]+pow(CONSTANTS[114]/STATES[19], CONSTANTS[148])+ CONSTANTS[272]*ALGEBRAIC[67])); ALGEBRAIC[61] = CONSTANTS[24] - STATES[15]; ALGEBRAIC[77] = CONSTANTS[78]*( STATES[20]*STATES[15] - CONSTANTS[275]*STATES[21]*ALGEBRAIC[61]); RATES[20] = ALGEBRAIC[76] - ALGEBRAIC[77]; ALGEBRAIC[78] = CONSTANTS[276]*(STATES[21]/(STATES[21]+ (CONSTANTS[123]+ CONSTANTS[277]*STATES[22])*(1.00000+ CONSTANTS[278]*STATES[18]))); RATES[21] = ALGEBRAIC[77] - ALGEBRAIC[78]; ALGEBRAIC[79] = CONSTANTS[76]*STATES[22] - CONSTANTS[279]*STATES[23]; RATES[22] = ALGEBRAIC[78] - ALGEBRAIC[79]; ALGEBRAIC[80] = CONSTANTS[280]*STATES[23]*(ALGEBRAIC[65]/( (STATES[23]+CONSTANTS[104]+ STATES[18]*CONSTANTS[281])*(CONSTANTS[117]+ALGEBRAIC[65]))); RATES[23] = ALGEBRAIC[79] - ALGEBRAIC[80]; ALGEBRAIC[81] = CONSTANTS[282]*(STATES[18]/(CONSTANTS[283]+STATES[19])); RATES[19] = ALGEBRAIC[74]+(ALGEBRAIC[81] - ALGEBRAIC[76]); ALGEBRAIC[70] = CONSTANTS[265]*(STATES[18]/(STATES[18]+CONSTANTS[121])); RATES[18] = ALGEBRAIC[80] - (ALGEBRAIC[70]+ALGEBRAIC[81]); RATES[49] = ALGEBRAIC[81] - CONSTANTS[69]*STATES[49]; ALGEBRAIC[63] = CONSTANTS[307]+STATES[16]; ALGEBRAIC[68] = exp( CONSTANTS[342]*ALGEBRAIC[63]); ALGEBRAIC[82] = pow(( STATES[17]*CONSTANTS[286]*ALGEBRAIC[67]), 1.0 / 2); ALGEBRAIC[83] = 1.00000/((CONSTANTS[337]+ CONSTANTS[346]*ALGEBRAIC[82])+ (CONSTANTS[164]+ CONSTANTS[165]*ALGEBRAIC[82])*ALGEBRAIC[68]); ALGEBRAIC[84] = ( (CONSTANTS[348]+ CONSTANTS[350]*ALGEBRAIC[68])*ALGEBRAIC[82] - CONSTANTS[352]*ALGEBRAIC[68])*ALGEBRAIC[83]; RATES[17] = - ALGEBRAIC[84]+ALGEBRAIC[74]+ALGEBRAIC[76]+ALGEBRAIC[80]; ALGEBRAIC[88] = exp( CONSTANTS[372]*ALGEBRAIC[63]); ALGEBRAIC[89] = CONSTANTS[287]*(ALGEBRAIC[61]/STATES[15]); ALGEBRAIC[90] = - (CONSTANTS[170]/(CONSTANTS[371]+ CONSTANTS[376]*ALGEBRAIC[89]+ (CONSTANTS[152]+ CONSTANTS[153]*ALGEBRAIC[89])*ALGEBRAIC[88])); ALGEBRAIC[91] = ( (CONSTANTS[373]+ CONSTANTS[157]*ALGEBRAIC[88])*ALGEBRAIC[89] - CONSTANTS[154]*ALGEBRAIC[88])*ALGEBRAIC[90]; ALGEBRAIC[7] = 8.00000 - STATES[0]; ALGEBRAIC[93] = STATES[0]/ALGEBRAIC[7]; ALGEBRAIC[94] = STATES[15]/ALGEBRAIC[61]; ALGEBRAIC[95] = (CONSTANTS[288] - CONSTANTS[289]*ALGEBRAIC[93]*ALGEBRAIC[94]*exp(- ( CONSTANTS[319]*STATES[16])))/( (1.00000+ (10.0000/9.00000)*ALGEBRAIC[93]*exp(- ( CONSTANTS[322]*STATES[16])))*(1.00000+ 18.0000*ALGEBRAIC[94])); RATES[15] = ALGEBRAIC[95] - (ALGEBRAIC[91]+ALGEBRAIC[77]); ALGEBRAIC[15] = 1.00000/STATES[0]; ALGEBRAIC[97] = pow( STATES[14]*CONSTANTS[290], CONSTANTS[146]); ALGEBRAIC[98] = pow( STATES[26]*CONSTANTS[291], CONSTANTS[149]); ALGEBRAIC[99] = CONSTANTS[135]*(( CONSTANTS[179]*ALGEBRAIC[97] - CONSTANTS[180]*ALGEBRAIC[98])/( (1.00000+ALGEBRAIC[97]+ALGEBRAIC[98])*( ALGEBRAIC[15]*(CONSTANTS[106]+ ALGEBRAIC[7]*CONSTANTS[377])+(1.00000+ ALGEBRAIC[7]*CONSTANTS[378])))); ALGEBRAIC[102] = (STATES[26] - STATES[27])*CONSTANTS[292]; RATES[26] = ALGEBRAIC[99]*CONSTANTS[383] - ALGEBRAIC[102]*CONSTANTS[384]; ALGEBRAIC[100] = CONSTANTS[9]*(STATES[46] - STATES[47]); ALGEBRAIC[104] = CONSTANTS[11]*( (CONSTANTS[14] - STATES[46])*STATES[0] - STATES[46]*ALGEBRAIC[7]*CONSTANTS[300]); RATES[46] = ALGEBRAIC[104] - ALGEBRAIC[100]; ALGEBRAIC[103] = CONSTANTS[10]*( (CONSTANTS[13] - STATES[47])*STATES[48] - STATES[47]*(8.00000 - STATES[48])*CONSTANTS[300]); RATES[47] = ALGEBRAIC[100]+ALGEBRAIC[103]; ALGEBRAIC[19] = STATES[1]*CONSTANTS[319]; ALGEBRAIC[22] = exp(ALGEBRAIC[19]); ALGEBRAIC[37] = pow(( STATES[13]*STATES[13]*STATES[13]), 1.0 / 2); ALGEBRAIC[39] = CONSTANTS[327]*(ALGEBRAIC[37]/( (ALGEBRAIC[37]+ pow(( CONSTANTS[119]*CONSTANTS[119]*CONSTANTS[119]), 1.0 / 2))*(1.00000+ 0.124500*exp( -0.100000*ALGEBRAIC[19])+CONSTANTS[324]/ALGEBRAIC[22])*(1.00000+ ( CONSTANTS[100]*ALGEBRAIC[15])*(1.00000+ ALGEBRAIC[7]*CONSTANTS[331])))); ALGEBRAIC[47] = CONSTANTS[57]*(STATES[14]/(CONSTANTS[122]+STATES[14]))*(1.00000/(1.00000+ ( CONSTANTS[101]*ALGEBRAIC[15])*(1.00000+ ALGEBRAIC[7]*CONSTANTS[262]))+1.00000/(1.00000+ CONSTANTS[102]*ALGEBRAIC[15])); ALGEBRAIC[59] = CONSTANTS[318]*(( CONSTANTS[248]*STATES[30]+ CONSTANTS[304]*STATES[31]+ CONSTANTS[312]*STATES[32])/(1.00000+ ALGEBRAIC[15]*(CONSTANTS[105]+ CONSTANTS[264]*ALGEBRAIC[7]))); RATES[0] = (1.00000 - CONSTANTS[3])*( 0.615000*ALGEBRAIC[95] - (ALGEBRAIC[59]+ 0.500000*ALGEBRAIC[99]+ 6.37100e-05*(ALGEBRAIC[39]+ALGEBRAIC[47])+ALGEBRAIC[104])); RATES[48] = (1.00000 - CONSTANTS[4])*(- ALGEBRAIC[103] - CONSTANTS[15]); ALGEBRAIC[25] = CONSTANTS[314]*ALGEBRAIC[19]; ALGEBRAIC[27] = ALGEBRAIC[22]*ALGEBRAIC[22]; ALGEBRAIC[31] = (fabs(STATES[1])