/* 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[0] 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[1] is ADP in component model (dimensionless). * STATES[0] is ATPi in component model (dimensionless). * ALGEBRAIC[2] is inv_ATPi in component model (dimensionless). * ALGEBRAIC[3] is VF_over_RT in component model (dimensionless). * STATES[1] is V in component model (dimensionless). * ALGEBRAIC[4] is exp_VF_over_RT in component model (dimensionless). * ALGEBRAIC[5] is VFsq_over_RT in component model (dimensionless). * ALGEBRAIC[6] is exp2VFRT in component model (dimensionless). * ALGEBRAIC[7] is O1_RyR in component model (dimensionless). * STATES[2] is C1_RyR in component model (dimensionless). * STATES[3] is C2_RyR in component model (dimensionless). * STATES[4] is O2_RyR in component model (dimensionless). * ALGEBRAIC[8] is V_30 in component model (dimensionless). * ALGEBRAIC[124] is V_E_K in component model (dimensionless). * ALGEBRAIC[125] is INa in component model (dimensionless). * STATES[5] is mNa in component model (dimensionless). * STATES[6] is hNa in component model (dimensionless). * STATES[7] is jNa in component model (dimensionless). * ALGEBRAIC[126] is IKs in component model (dimensionless). * STATES[8] is xKs in component model (dimensionless). * ALGEBRAIC[127] is K1Alpha in component model (dimensionless). * ALGEBRAIC[128] is K1Beta in component model (dimensionless). * ALGEBRAIC[129] is K1_inf in component model (dimensionless). * ALGEBRAIC[130] is IK1 in component model (dimensionless). * ALGEBRAIC[131] is INab in component model (dimensionless). * ALGEBRAIC[132] is IKp in component model (dimensionless). * ALGEBRAIC[9] is ICamax in component model (dimensionless). * ALGEBRAIC[10] is ICaK in component model (dimensionless). * STATES[9] is Open in component model (dimensionless). * STATES[10] is OCa in component model (dimensionless). * STATES[11] is yCa in component model (dimensionless). * STATES[12] is Ki in component model (dimensionless). * ALGEBRAIC[11] is ICa in component model (dimensionless). * ALGEBRAIC[12] is NaiP1p5 in component model (dimensionless). * STATES[13] is Nai in component model (dimensionless). * ALGEBRAIC[13] is INaK in component model (dimensionless). * ALGEBRAIC[14] is exp_eta_VF_over_RT in component model (dimensionless). * ALGEBRAIC[15] is exp_eta1_VF_over_RT in component model (dimensionless). * ALGEBRAIC[16] is INaCa in component model (dimensionless). * STATES[14] is Cai in component model (dimensionless). * ALGEBRAIC[133] is ICab in component model (dimensionless). * ALGEBRAIC[17] is IpCa in component model (dimensionless). * ALGEBRAIC[18] is CaiP3 in component model (dimensionless). * ALGEBRAIC[19] is common in component model (dimensionless). * ALGEBRAIC[20] is InsNa in component model (dimensionless). * ALGEBRAIC[21] is InsK in component model (dimensionless). * ALGEBRAIC[22] is InsCa in component model (dimensionless). * ALGEBRAIC[23] is V_AM in component model (dimensionless). * ALGEBRAIC[24] is ATPm in component model (dimensionless). * STATES[15] is ADPm in component model (dimensionless). * ALGEBRAIC[25] is DmuH in component model (dimensionless). * STATES[16] is Dpsi in component model (dimensionless). * ALGEBRAIC[26] is NAD in component model (dimensionless). * STATES[17] is NADH in component model (dimensionless). * ALGEBRAIC[27] is KmIDNAD_NAD in component model (dimensionless). * ALGEBRAIC[28] is exp_FRT_6_g_DmuH in component model (dimensionless). * ALGEBRAIC[29] is FRT2_Dpsi in component model (dimensionless). * ALGEBRAIC[30] is VCS in component model (dimensionless). * STATES[18] is Oaa in component model (dimensionless). * ALGEBRAIC[31] is VACO in component model (dimensionless). * STATES[19] is AKG in component model (dimensionless). * STATES[20] is SCoA in component model (dimensionless). * STATES[21] is Succ in component model (dimensionless). * STATES[22] is FUM in component model (dimensionless). * STATES[23] is MAL in component model (dimensionless). * STATES[24] is ISOC in component model (dimensionless). * ALGEBRAIC[32] is Fa in component model (dimensionless). * STATES[25] is Cam in component model (dimensionless). * ALGEBRAIC[33] is Fi in component model (dimensionless). * ALGEBRAIC[34] is VIDH in component model (dimensionless). * ALGEBRAIC[35] is a__1 in component model (dimensionless). * ALGEBRAIC[36] is VKGDH in component model (dimensionless). * ALGEBRAIC[37] is VSL in component model (dimensionless). * ALGEBRAIC[38] is VSDH in component model (dimensionless). * ALGEBRAIC[39] is VFH in component model (dimensionless). * ALGEBRAIC[40] is VMDH in component model (dimensionless). * ALGEBRAIC[41] is VAAT in component model (dimensionless). * ALGEBRAIC[42] is AREN in component model (dimensionless). * ALGEBRAIC[43] is denominator1 in component model (dimensionless). * ALGEBRAIC[46] is denominator2 in component model (dimensionless). * ALGEBRAIC[50] is denominator3 in component model (dimensionless). * ALGEBRAIC[44] is VNO in component model (dimensionless). * ALGEBRAIC[45] is VHNe in component model (dimensionless). * ALGEBRAIC[47] is VHFe in component model (dimensionless). * ALGEBRAIC[48] is exp_3FRT_DmuH in component model (dimensionless). * ALGEBRAIC[49] is AF1 in component model (dimensionless). * ALGEBRAIC[51] is VATPase in component model (dimensionless). * ALGEBRAIC[52] is Vhu in component model (dimensionless). * ALGEBRAIC[53] is ATPi_ADP in component model (dimensionless). * ALGEBRAIC[54] is ADPm_ATPm in component model (dimensionless). * ALGEBRAIC[55] is VANT in component model (dimensionless). * ALGEBRAIC[56] is Vhleak in component model (dimensionless). * ALGEBRAIC[57] is MAlpha in component model (dimensionless). * ALGEBRAIC[58] is MBeta in component model (dimensionless). * ALGEBRAIC[59] is inv_MBeta_MAlpha in component model (dimensionless). * ALGEBRAIC[60] is tmNa in component model (dimensionless). * ALGEBRAIC[61] is HAlpha in component model (dimensionless). * ALGEBRAIC[62] is HBeta in component model (dimensionless). * ALGEBRAIC[63] is JAlpha in component model (dimensionless). * ALGEBRAIC[64] is JBeta in component model (dimensionless). * ALGEBRAIC[65] is fb in component model (dimensionless). * ALGEBRAIC[66] is rb in component model (per_ms). * STATES[26] is CaNSR in component model (dimensionless). * ALGEBRAIC[67] is Jup in component model (dimensionless). * ALGEBRAIC[68] is Jrel in component model (dimensionless). * STATES[27] is CaJSR in component model (dimensionless). * STATES[28] is CaSS in component model (dimensionless). * ALGEBRAIC[69] is Jtr in component model (dimensionless). * ALGEBRAIC[70] is Jxfer in component model (dimensionless). * ALGEBRAIC[71] is P1_N1_P2_P3 in component model (dimensionless). * STATES[29] is N1 in component model (dimensionless). * STATES[30] is P0 in component model (dimensionless). * STATES[31] is P1 in component model (dimensionless). * STATES[32] is P2 in component model (dimensionless). * STATES[33] is P3 in component model (dimensionless). * ALGEBRAIC[80] is N1_differential in component model (dimensionless). * ALGEBRAIC[76] is P0_differential in component model (dimensionless). * ALGEBRAIC[77] is P1_differential in component model (dimensionless). * ALGEBRAIC[78] is P2_differential in component model (dimensionless). * ALGEBRAIC[79] is P3_differential in component model (dimensionless). * ALGEBRAIC[72] is FN_Ca in component model (dimensionless). * ALGEBRAIC[73] is force_norm in component model (dimensionless). * ALGEBRAIC[74] is force in component model (dimensionless). * ALGEBRAIC[75] is kTrop_np in component model (dimensionless). * STATES[34] is LTRPNCa in component model (dimensionless). * STATES[35] is HTRPNCa in component model (dimensionless). * ALGEBRAIC[81] is LTRPNCa_differential in component model (dimensionless). * ALGEBRAIC[82] is HTRPNCa_differential in component model (dimensionless). * ALGEBRAIC[83] is Jtrpn in component model (dimensionless). * ALGEBRAIC[84] is beta_SS in component model (dimensionless). * ALGEBRAIC[85] is beta_JSR in component model (dimensionless). * ALGEBRAIC[86] is beta_i in component model (dimensionless). * ALGEBRAIC[87] is Cai_ktrans_plus1 in component model (dimensionless). * ALGEBRAIC[88] is Cai_ktrans_plus1_p3 in component model (dimensionless). * ALGEBRAIC[89] is Vuni in component model (dimensionless). * ALGEBRAIC[90] is VnaCa in component model (dimensionless). * ALGEBRAIC[91] is alpha in component model (dimensionless). * ALGEBRAIC[92] is beta in component model (dimensionless). * ALGEBRAIC[93] is alpha_prime in component model (dimensionless). * ALGEBRAIC[94] is beta_prime in component model (dimensionless). * ALGEBRAIC[95] is C0_to_C1 in component model (dimensionless). * ALGEBRAIC[96] is C1_to_C2 in component model (dimensionless). * ALGEBRAIC[97] is C2_to_C3 in component model (dimensionless). * ALGEBRAIC[98] is C3_to_C4 in component model (dimensionless). * ALGEBRAIC[99] is CCa0_to_CCa1 in component model (dimensionless). * ALGEBRAIC[100] is CCa1_to_CCa2 in component model (dimensionless). * ALGEBRAIC[101] is CCa2_to_CCa3 in component model (dimensionless). * ALGEBRAIC[102] is CCa3_to_CCa4 in component model (dimensionless). * ALGEBRAIC[103] is C1_to_C0 in component model (dimensionless). * ALGEBRAIC[104] is C2_to_C1 in component model (dimensionless). * ALGEBRAIC[105] is C3_to_C2 in component model (dimensionless). * ALGEBRAIC[106] is C4_to_C3 in component model (dimensionless). * ALGEBRAIC[107] is CCa1_to_CCa0 in component model (dimensionless). * ALGEBRAIC[108] is CCa2_to_CCa1 in component model (dimensionless). * ALGEBRAIC[109] is CCa3_to_CCa2 in component model (dimensionless). * ALGEBRAIC[110] is CCa4_to_CCa3 in component model (dimensionless). * ALGEBRAIC[111] is gamma in component model (dimensionless). * ALGEBRAIC[112] is C0_to_CCa0 in component model (dimensionless). * ALGEBRAIC[113] is C1_to_CCa1 in component model (dimensionless). * ALGEBRAIC[114] is C2_to_CCa2 in component model (dimensionless). * ALGEBRAIC[115] is C3_to_CCa3 in component model (dimensionless). * ALGEBRAIC[116] is C4_to_CCa4 in component model (dimensionless). * CONSTANTS[311] is CCa0_to_C0 in component model (dimensionless). * CONSTANTS[315] is CCa1_to_C1 in component model (dimensionless). * CONSTANTS[320] is CCa2_to_C2 in component model (dimensionless). * CONSTANTS[325] is CCa3_to_C3 in component model (dimensionless). * CONSTANTS[328] is CCa4_to_C4 in component model (dimensionless). * STATES[36] is C1 in component model (dimensionless). * STATES[37] is CCa0 in component model (dimensionless). * STATES[38] is C0 in component model (dimensionless). * STATES[39] is C2 in component model (dimensionless). * STATES[40] is CCa1 in component model (dimensionless). * STATES[41] is C3 in component model (dimensionless). * STATES[42] is CCa2 in component model (dimensionless). * STATES[43] is C4 in component model (dimensionless). * STATES[44] is CCa3 in component model (dimensionless). * STATES[45] is CCa4 in component model (dimensionless). * ALGEBRAIC[117] is Vt_CRP2 in component model (dimensionless). * STATES[46] is CrPi_mito in component model (dimensionless). * STATES[47] is CrPi_cyto in component model (dimensionless). * ALGEBRAIC[118] is VCK_cyto in component model (dimensionless). * STATES[48] is ATPi_cyto in component model (dimensionless). * ALGEBRAIC[119] is VCK_mito in component model (dimensionless). * STATES[49] is ASP in component model (dimensionless). * STATES[50] is N0 in component model (dimensionless). * ALGEBRAIC[121] is reversal_potential in component one_ion (dimensionless). * ALGEBRAIC[120] is reversal_potential in component one_ion (dimensionless). * ALGEBRAIC[123] is reversal_potential in component two_ions (dimensionless). * ALGEBRAIC[122] is reversal_potential in component one_ion (dimensionless). * RATES[5] is d/dt mNa in component model (dimensionless). * RATES[6] is d/dt hNa in component model (dimensionless). * RATES[7] is d/dt jNa in component model (dimensionless). * RATES[8] is d/dt xKs in component model (dimensionless). * RATES[30] is d/dt P0 in component model (dimensionless). * RATES[31] is d/dt P1 in component model (dimensionless). * RATES[32] is d/dt P2 in component model (dimensionless). * RATES[33] is d/dt P3 in component model (dimensionless). * RATES[29] is d/dt N1 in component model (dimensionless). * RATES[50] is d/dt N0 in component model (dimensionless). * RATES[34] is d/dt LTRPNCa in component model (dimensionless). * RATES[35] is d/dt HTRPNCa in component model (dimensionless). * RATES[13] is d/dt Nai in component model (dimensionless). * RATES[12] is d/dt Ki in component model (dimensionless). * RATES[14] is d/dt Cai in component model (dimensionless). * RATES[28] is d/dt CaSS in component model (dimensionless). * RATES[27] is d/dt CaJSR in component model (dimensionless). * RATES[26] is d/dt CaNSR in component model (dimensionless). * RATES[1] is d/dt V in component model (dimensionless). * RATES[2] is d/dt C1_RyR in component model (dimensionless). * RATES[4] is d/dt O2_RyR in component model (dimensionless). * RATES[3] is d/dt C2_RyR in component model (dimensionless). * RATES[38] is d/dt C0 in component model (dimensionless). * RATES[36] is d/dt C1 in component model (dimensionless). * RATES[39] is d/dt C2 in component model (dimensionless). * RATES[41] is d/dt C3 in component model (dimensionless). * RATES[43] is d/dt C4 in component model (dimensionless). * RATES[9] is d/dt Open in component model (dimensionless). * RATES[37] is d/dt CCa0 in component model (dimensionless). * RATES[40] is d/dt CCa1 in component model (dimensionless). * RATES[42] is d/dt CCa2 in component model (dimensionless). * RATES[44] is d/dt CCa3 in component model (dimensionless). * RATES[45] is d/dt CCa4 in component model (dimensionless). * RATES[11] is d/dt yCa in component model (dimensionless). * RATES[10] is d/dt OCa in component model (dimensionless). * RATES[46] is d/dt CrPi_mito in component model (dimensionless). * RATES[47] is d/dt CrPi_cyto in component model (dimensionless). * RATES[0] is d/dt ATPi in component model (dimensionless). * RATES[48] is d/dt ATPi_cyto in component model (dimensionless). * RATES[25] is d/dt Cam in component model (dimensionless). * RATES[15] is d/dt ADPm in component model (dimensionless). * RATES[16] is d/dt Dpsi in component model (dimensionless). * RATES[17] is d/dt NADH in component model (dimensionless). * RATES[24] is d/dt ISOC in component model (dimensionless). * RATES[19] is d/dt AKG in component model (dimensionless). * RATES[20] is d/dt SCoA in component model (dimensionless). * RATES[21] is d/dt Succ in component model (dimensionless). * RATES[22] is d/dt FUM in component model (dimensionless). * RATES[23] is d/dt MAL in component model (dimensionless). * RATES[18] is d/dt Oaa in component model (dimensionless). * RATES[49] is d/dt ASP in component model (dimensionless). * There are a total of 9 condition variables. */ 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; 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]; STATES[0] = CONSTANTS[196]; STATES[1] = CONSTANTS[197]; STATES[2] = CONSTANTS[198]; STATES[3] = CONSTANTS[199]; STATES[4] = CONSTANTS[200]; STATES[5] = CONSTANTS[201]; STATES[6] = CONSTANTS[202]; STATES[7] = CONSTANTS[203]; STATES[8] = CONSTANTS[204]; STATES[9] = CONSTANTS[205]; STATES[10] = CONSTANTS[206]; STATES[11] = CONSTANTS[207]; STATES[12] = CONSTANTS[208]; STATES[13] = CONSTANTS[209]; STATES[14] = CONSTANTS[210]; STATES[15] = CONSTANTS[211]; STATES[16] = CONSTANTS[212]; STATES[17] = CONSTANTS[213]; STATES[18] = CONSTANTS[214]; STATES[19] = CONSTANTS[215]; STATES[20] = CONSTANTS[216]; STATES[21] = CONSTANTS[217]; STATES[22] = CONSTANTS[218]; STATES[23] = CONSTANTS[219]; STATES[24] = CONSTANTS[220]; STATES[25] = CONSTANTS[221]; STATES[26] = CONSTANTS[222]; STATES[27] = CONSTANTS[223]; STATES[28] = CONSTANTS[224]; STATES[29] = CONSTANTS[225]; STATES[30] = CONSTANTS[226]; STATES[31] = CONSTANTS[227]; STATES[32] = CONSTANTS[228]; STATES[33] = CONSTANTS[229]; STATES[34] = CONSTANTS[230]; STATES[35] = CONSTANTS[231]; STATES[36] = CONSTANTS[232]; STATES[37] = CONSTANTS[233]; STATES[38] = CONSTANTS[234]; STATES[39] = CONSTANTS[235]; STATES[40] = CONSTANTS[236]; STATES[41] = CONSTANTS[237]; STATES[42] = CONSTANTS[238]; STATES[43] = CONSTANTS[239]; STATES[44] = CONSTANTS[240]; STATES[45] = CONSTANTS[241]; STATES[46] = CONSTANTS[242]; STATES[47] = CONSTANTS[243]; STATES[48] = CONSTANTS[244]; STATES[49] = CONSTANTS[245]; STATES[50] = CONSTANTS[246]; RATES[5] = 0.1001; RATES[6] = 0.1001; RATES[7] = 0.1001; RATES[8] = 0.1001; RATES[30] = 0.1001; RATES[31] = 0.1001; RATES[32] = 0.1001; RATES[33] = 0.1001; RATES[29] = 0.1001; RATES[50] = 0.1001; RATES[34] = 0.1001; RATES[35] = 0.1001; RATES[13] = 0.1001; RATES[12] = 0.1001; RATES[14] = 0.1001; RATES[28] = 0.1001; RATES[27] = 0.1001; RATES[26] = 0.1001; RATES[1] = 0.1001; RATES[2] = 0.1001; RATES[4] = 0.1001; RATES[3] = 0.1001; RATES[38] = 0.1001; RATES[36] = 0.1001; RATES[39] = 0.1001; RATES[41] = 0.1001; RATES[43] = 0.1001; RATES[9] = 0.1001; RATES[37] = 0.1001; RATES[40] = 0.1001; RATES[42] = 0.1001; RATES[44] = 0.1001; RATES[45] = 0.1001; RATES[11] = 0.1001; RATES[10] = 0.1001; RATES[46] = 0.1001; RATES[47] = 0.1001; RATES[0] = 0.1001; RATES[48] = 0.1001; RATES[25] = 0.1001; RATES[15] = 0.1001; RATES[16] = 0.1001; RATES[17] = 0.1001; RATES[24] = 0.1001; RATES[19] = 0.1001; RATES[20] = 0.1001; RATES[21] = 0.1001; RATES[22] = 0.1001; RATES[23] = 0.1001; RATES[18] = 0.1001; RATES[49] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[5] - ALGEBRAIC[57]*(1.00000 - ALGEBRAIC[60]) - ALGEBRAIC[58]*ALGEBRAIC[60]; resid[1] = RATES[6] - ALGEBRAIC[61]*(1.00000 - STATES[6]) - ALGEBRAIC[62]*STATES[6]; resid[2] = RATES[7] - ALGEBRAIC[63]*(1.00000 - STATES[7]) - ALGEBRAIC[64]*STATES[7]; resid[3] = RATES[8] - 7.19000e-05*(ALGEBRAIC[8]/(1.00000 - exp( -0.148000*ALGEBRAIC[8])))*(1.00000 - STATES[8]) - 0.000131000*(ALGEBRAIC[8]/(exp( 0.0687000*ALGEBRAIC[8]) - 1.00000))*STATES[8]; resid[4] = RATES[30] - ALGEBRAIC[76]; resid[5] = RATES[31] - ALGEBRAIC[77]; resid[6] = RATES[32] - ALGEBRAIC[78]; resid[7] = RATES[33] - ALGEBRAIC[79]; resid[8] = RATES[29] - ALGEBRAIC[80]; resid[9] = RATES[50] - - ALGEBRAIC[76] - (ALGEBRAIC[77]+ALGEBRAIC[78]+ALGEBRAIC[79]+ALGEBRAIC[80]); resid[10] = RATES[34] - ALGEBRAIC[81]; resid[11] = RATES[35] - ALGEBRAIC[82]; resid[12] = RATES[13] - - ( (ALGEBRAIC[125]+ALGEBRAIC[131]+ALGEBRAIC[20]+ 3.00000*(ALGEBRAIC[16]+ALGEBRAIC[13]))*CONSTANTS[379]) - ALGEBRAIC[90]*0.615000; resid[13] = RATES[12] - - ( (ALGEBRAIC[21]+ALGEBRAIC[126]+ALGEBRAIC[130]+ALGEBRAIC[132]+ALGEBRAIC[10]+(ALGEBRAIC[0] - 2.00000*ALGEBRAIC[13]))*CONSTANTS[379]); resid[14] = RATES[14] - ALGEBRAIC[86]*((ALGEBRAIC[70] - (ALGEBRAIC[67]+ALGEBRAIC[83]+ 0.250000*CONSTANTS[379]*((ALGEBRAIC[133] - 2.00000*ALGEBRAIC[16])+ALGEBRAIC[17])))+ (ALGEBRAIC[90] - ALGEBRAIC[89])*0.615000); resid[15] = RATES[28] - ALGEBRAIC[84]*( ALGEBRAIC[68]*CONSTANTS[381] - ( ALGEBRAIC[70]*CONSTANTS[382]+ ALGEBRAIC[11]*CONSTANTS[380])); resid[16] = RATES[27] - ALGEBRAIC[85]*(ALGEBRAIC[69] - ALGEBRAIC[68]); resid[17] = RATES[26] - ALGEBRAIC[67]*CONSTANTS[383] - ALGEBRAIC[69]*CONSTANTS[384]; resid[18] = RATES[1] - - ( CONSTANTS[385]*(ALGEBRAIC[125]+ALGEBRAIC[11]+ALGEBRAIC[10]+ALGEBRAIC[126]+ALGEBRAIC[130]+ALGEBRAIC[132]+ALGEBRAIC[16]+ALGEBRAIC[13]+ALGEBRAIC[22]+ALGEBRAIC[17]+ALGEBRAIC[133]+ALGEBRAIC[131]+ALGEBRAIC[0])); resid[19] = RATES[2] - - ( CONSTANTS[64]*pow(STATES[28], CONSTANTS[145])*STATES[2])+ CONSTANTS[63]*ALGEBRAIC[7]; resid[20] = RATES[4] - CONSTANTS[66]*pow(STATES[28], CONSTANTS[140])*ALGEBRAIC[7] - CONSTANTS[65]*STATES[4]; resid[21] = RATES[3] - CONSTANTS[70]*ALGEBRAIC[7] - CONSTANTS[68]*STATES[3]; resid[22] = RATES[38] - ALGEBRAIC[103]*STATES[36]+( CONSTANTS[311]*STATES[37] - (ALGEBRAIC[95]+ALGEBRAIC[112])*STATES[38]); resid[23] = RATES[36] - ALGEBRAIC[95]*STATES[38]+ ALGEBRAIC[104]*STATES[39]+( CONSTANTS[315]*STATES[40] - (ALGEBRAIC[103]+ALGEBRAIC[96]+ALGEBRAIC[113])*STATES[36]); resid[24] = RATES[39] - ALGEBRAIC[96]*STATES[36]+ ALGEBRAIC[105]*STATES[41]+( CONSTANTS[320]*STATES[42] - (ALGEBRAIC[104]+ALGEBRAIC[97]+ALGEBRAIC[114])*STATES[39]); resid[25] = RATES[41] - ALGEBRAIC[97]*STATES[39]+ ALGEBRAIC[106]*STATES[43]+( CONSTANTS[325]*STATES[44] - (ALGEBRAIC[105]+ALGEBRAIC[98]+ALGEBRAIC[115])*STATES[41]); resid[26] = RATES[43] - ALGEBRAIC[98]*STATES[41]+ CONSTANTS[49]*STATES[9]+( CONSTANTS[328]*STATES[45] - (ALGEBRAIC[106]+CONSTANTS[40]+ALGEBRAIC[116])*STATES[43]); resid[27] = RATES[9] - CONSTANTS[40]*STATES[43] - CONSTANTS[49]*STATES[9]; resid[28] = RATES[37] - ALGEBRAIC[107]*STATES[40]+( ALGEBRAIC[112]*STATES[38] - (ALGEBRAIC[99]+CONSTANTS[311])*STATES[37]); resid[29] = RATES[40] - ALGEBRAIC[99]*STATES[37]+ ALGEBRAIC[108]*STATES[42]+( ALGEBRAIC[113]*STATES[36] - (ALGEBRAIC[107]+ALGEBRAIC[100]+CONSTANTS[315])*STATES[40]); resid[30] = RATES[42] - ALGEBRAIC[100]*STATES[40]+ ALGEBRAIC[109]*STATES[44]+( ALGEBRAIC[114]*STATES[39] - (ALGEBRAIC[108]+ALGEBRAIC[101]+CONSTANTS[320])*STATES[42]); resid[31] = RATES[44] - ALGEBRAIC[101]*STATES[42]+ ALGEBRAIC[110]*STATES[45]+( ALGEBRAIC[115]*STATES[41] - (ALGEBRAIC[109]+ALGEBRAIC[102]+CONSTANTS[325])*STATES[44]); resid[32] = RATES[45] - ALGEBRAIC[102]*STATES[44]+ CONSTANTS[51]*STATES[10]+( ALGEBRAIC[116]*STATES[43] - (ALGEBRAIC[110]+CONSTANTS[42]+CONSTANTS[328])*STATES[45]); resid[33] = 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)))); resid[34] = RATES[10] - CONSTANTS[42]*STATES[45] - CONSTANTS[51]*STATES[10]; resid[35] = RATES[46] - ALGEBRAIC[119] - ALGEBRAIC[117]; resid[36] = RATES[47] - ALGEBRAIC[117]+ALGEBRAIC[118]; resid[37] = RATES[0] - (1.00000 - CONSTANTS[3])*( 0.615000*ALGEBRAIC[55] - (ALGEBRAIC[23]+ 0.500000*ALGEBRAIC[67]+ 6.37100e-05*(ALGEBRAIC[13]+ALGEBRAIC[17])+ALGEBRAIC[119])); resid[38] = RATES[48] - (1.00000 - CONSTANTS[4])*(- ALGEBRAIC[118] - CONSTANTS[15]); resid[39] = RATES[25] - CONSTANTS[41]*(ALGEBRAIC[89] - ALGEBRAIC[90]); resid[40] = RATES[15] - ALGEBRAIC[55] - (ALGEBRAIC[51]+ALGEBRAIC[37]); resid[41] = RATES[16] - - ( ((- ALGEBRAIC[45] - ALGEBRAIC[47])+ALGEBRAIC[52]+ALGEBRAIC[55]+ALGEBRAIC[56]+ CONSTANTS[387]*ALGEBRAIC[90]+ 2.00000*ALGEBRAIC[89])*CONSTANTS[386]); resid[42] = RATES[17] - - ALGEBRAIC[44]+ALGEBRAIC[34]+ALGEBRAIC[36]+ALGEBRAIC[40]; resid[43] = RATES[24] - ALGEBRAIC[31] - ALGEBRAIC[34]; resid[44] = RATES[19] - ALGEBRAIC[34]+(ALGEBRAIC[41] - ALGEBRAIC[36]); resid[45] = RATES[20] - ALGEBRAIC[36] - ALGEBRAIC[37]; resid[46] = RATES[21] - ALGEBRAIC[37] - ALGEBRAIC[38]; resid[47] = RATES[22] - ALGEBRAIC[38] - ALGEBRAIC[39]; resid[48] = RATES[23] - ALGEBRAIC[39] - ALGEBRAIC[40]; resid[49] = RATES[18] - ALGEBRAIC[40] - (ALGEBRAIC[30]+ALGEBRAIC[41]); resid[50] = RATES[49] - ALGEBRAIC[41] - CONSTANTS[69]*STATES[49]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[73] = CONSTANTS[366]*(ALGEBRAIC[71]+STATES[32]+STATES[33]+STATES[33]); ALGEBRAIC[74] = CONSTANTS[368]*(ALGEBRAIC[71]+STATES[32]+STATES[33]+STATES[33]); } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[0] = (CONDVAR[0]<0.00000 ? CONSTANTS[2] : 0.00000); ALGEBRAIC[7] = 1.00000 - (STATES[2]+STATES[3]+STATES[4]); ALGEBRAIC[8] = STATES[1]+30.0000; ALGEBRAIC[3] = STATES[1]*CONSTANTS[319]; ALGEBRAIC[4] = exp(ALGEBRAIC[3]); ALGEBRAIC[6] = ALGEBRAIC[4]*ALGEBRAIC[4]; ALGEBRAIC[5] = CONSTANTS[314]*ALGEBRAIC[3]; ALGEBRAIC[9] = (CONDVAR[1]<0.00000 ? CONSTANTS[256]*(ALGEBRAIC[6] - CONSTANTS[254])*(0.500000 - 0.0200000*STATES[1]) : CONSTANTS[256]*ALGEBRAIC[5]*((ALGEBRAIC[6] - CONSTANTS[254])/(ALGEBRAIC[6] - 1.00000))); ALGEBRAIC[10] = (CONDVAR[2]<0.00000 ? CONSTANTS[317]*(STATES[9]+STATES[10])*STATES[11]*( STATES[12]*ALGEBRAIC[6] - CONSTANTS[127])*((0.500000 - 0.0200000*STATES[1])/(1.00000+ ALGEBRAIC[9]*CONSTANTS[257])) : CONSTANTS[317]*(STATES[9]+STATES[10])*STATES[11]*( STATES[12]*ALGEBRAIC[6] - CONSTANTS[127])*(ALGEBRAIC[3]/( (ALGEBRAIC[6] - 1.00000)*(1.00000+ ALGEBRAIC[9]*CONSTANTS[257])))); ALGEBRAIC[11] = 6.00000*ALGEBRAIC[9]*STATES[11]*STATES[9]; ALGEBRAIC[1] = 8.00000 - STATES[0]; ALGEBRAIC[2] = 1.00000/STATES[0]; ALGEBRAIC[12] = pow(( STATES[13]*STATES[13]*STATES[13]), 1.0 / 2); ALGEBRAIC[13] = CONSTANTS[327]*(ALGEBRAIC[12]/( (ALGEBRAIC[12]+ pow(( CONSTANTS[119]*CONSTANTS[119]*CONSTANTS[119]), 1.0 / 2))*(1.00000+ 0.124500*exp( -0.100000*ALGEBRAIC[3])+CONSTANTS[324]/ALGEBRAIC[4])*(1.00000+ ( CONSTANTS[100]*ALGEBRAIC[2])*(1.00000+ ALGEBRAIC[1]*CONSTANTS[331])))); ALGEBRAIC[14] = exp( CONSTANTS[32]*ALGEBRAIC[3]); ALGEBRAIC[15] = ALGEBRAIC[14]/ALGEBRAIC[4]; ALGEBRAIC[16] = ( ALGEBRAIC[14]*STATES[13]*STATES[13]*STATES[13] - ALGEBRAIC[15]*CONSTANTS[260]*STATES[14])/(CONSTANTS[261]+ CONSTANTS[306]*ALGEBRAIC[15]); ALGEBRAIC[17] = CONSTANTS[57]*(STATES[14]/(CONSTANTS[122]+STATES[14]))*(1.00000/(1.00000+ ( CONSTANTS[101]*ALGEBRAIC[2])*(1.00000+ ALGEBRAIC[1]*CONSTANTS[262]))+1.00000/(1.00000+ CONSTANTS[102]*ALGEBRAIC[2])); ALGEBRAIC[18] = STATES[14]*STATES[14]*STATES[14]; ALGEBRAIC[19] = (CONDVAR[3]<0.00000 ? 0.750000*ALGEBRAIC[18]*((1.00000 - 0.0200000*STATES[1])/(ALGEBRAIC[18]+CONSTANTS[263])) : 0.750000*ALGEBRAIC[18]*(ALGEBRAIC[5]/( (ALGEBRAIC[4] - 1.00000)*(ALGEBRAIC[18]+CONSTANTS[263])))); ALGEBRAIC[20] = CONSTANTS[162]*ALGEBRAIC[19]*( STATES[13]*ALGEBRAIC[4] - CONSTANTS[144]); ALGEBRAIC[21] = CONSTANTS[161]*ALGEBRAIC[19]*( STATES[12]*ALGEBRAIC[4] - CONSTANTS[127]); ALGEBRAIC[22] = ALGEBRAIC[20]+ALGEBRAIC[21]; ALGEBRAIC[23] = CONSTANTS[318]*(( CONSTANTS[248]*STATES[30]+ CONSTANTS[304]*STATES[31]+ CONSTANTS[312]*STATES[32])/(1.00000+ ALGEBRAIC[2]*(CONSTANTS[105]+ CONSTANTS[264]*ALGEBRAIC[1]))); ALGEBRAIC[30] = CONSTANTS[265]*(STATES[18]/(STATES[18]+CONSTANTS[121])); ALGEBRAIC[31] = CONSTANTS[74]*(CONSTANTS[23] - (STATES[19]+STATES[20]+STATES[21]+STATES[22]+STATES[23]+STATES[18]+ STATES[24]*CONSTANTS[266])); ALGEBRAIC[26] = CONSTANTS[27] - STATES[17]; ALGEBRAIC[27] = CONSTANTS[112]/ALGEBRAIC[26]; ALGEBRAIC[32] = 1.00000/( (1.00000+ STATES[15]*CONSTANTS[269])*(1.00000+ STATES[25]*CONSTANTS[270])); ALGEBRAIC[33] = 1.00000+ STATES[17]*CONSTANTS[271]; ALGEBRAIC[34] = CONSTANTS[268]/(CONSTANTS[267]+ ALGEBRAIC[27]*ALGEBRAIC[33]+ pow(CONSTANTS[113]/STATES[24], CONSTANTS[147])*ALGEBRAIC[32]*(1.00000+ ALGEBRAIC[27]*ALGEBRAIC[33])); ALGEBRAIC[35] = CONSTANTS[273]+ CONSTANTS[308]*STATES[25]; ALGEBRAIC[36] = CONSTANTS[274]*(ALGEBRAIC[35]/(ALGEBRAIC[35]+pow(CONSTANTS[114]/STATES[19], CONSTANTS[148])+ CONSTANTS[272]*ALGEBRAIC[27])); ALGEBRAIC[24] = CONSTANTS[24] - STATES[15]; ALGEBRAIC[37] = CONSTANTS[78]*( STATES[20]*STATES[15] - CONSTANTS[275]*STATES[21]*ALGEBRAIC[24]); ALGEBRAIC[38] = CONSTANTS[276]*(STATES[21]/(STATES[21]+ (CONSTANTS[123]+ CONSTANTS[277]*STATES[22])*(1.00000+ CONSTANTS[278]*STATES[18]))); ALGEBRAIC[39] = CONSTANTS[76]*STATES[22] - CONSTANTS[279]*STATES[23]; ALGEBRAIC[40] = CONSTANTS[280]*STATES[23]*(ALGEBRAIC[26]/( (STATES[23]+CONSTANTS[104]+ STATES[18]*CONSTANTS[281])*(CONSTANTS[117]+ALGEBRAIC[26]))); ALGEBRAIC[41] = CONSTANTS[282]*(STATES[18]/(CONSTANTS[283]+STATES[19])); ALGEBRAIC[25] = CONSTANTS[307]+STATES[16]; ALGEBRAIC[28] = exp( CONSTANTS[342]*ALGEBRAIC[25]); ALGEBRAIC[42] = pow(( STATES[17]*CONSTANTS[286]*ALGEBRAIC[27]), 1.0 / 2); ALGEBRAIC[43] = 1.00000/((CONSTANTS[337]+ CONSTANTS[346]*ALGEBRAIC[42])+ (CONSTANTS[164]+ CONSTANTS[165]*ALGEBRAIC[42])*ALGEBRAIC[28]); ALGEBRAIC[44] = ( (CONSTANTS[348]+ CONSTANTS[350]*ALGEBRAIC[28])*ALGEBRAIC[42] - CONSTANTS[352]*ALGEBRAIC[28])*ALGEBRAIC[43]; ALGEBRAIC[45] = ( CONSTANTS[354]*ALGEBRAIC[42] - CONSTANTS[356]*ALGEBRAIC[28])*ALGEBRAIC[43]; ALGEBRAIC[46] = CONSTANTS[171]/(CONSTANTS[369]+ CONSTANTS[370]*ALGEBRAIC[28]); ALGEBRAIC[47] = (CONSTANTS[364] - CONSTANTS[367]*ALGEBRAIC[28])*ALGEBRAIC[46]; ALGEBRAIC[48] = exp( CONSTANTS[372]*ALGEBRAIC[25]); ALGEBRAIC[49] = CONSTANTS[287]*(ALGEBRAIC[24]/STATES[15]); ALGEBRAIC[50] = - (CONSTANTS[170]/(CONSTANTS[371]+ CONSTANTS[376]*ALGEBRAIC[49]+ (CONSTANTS[152]+ CONSTANTS[153]*ALGEBRAIC[49])*ALGEBRAIC[48])); ALGEBRAIC[51] = ( (CONSTANTS[373]+ CONSTANTS[157]*ALGEBRAIC[48])*ALGEBRAIC[49] - CONSTANTS[154]*ALGEBRAIC[48])*ALGEBRAIC[50]; ALGEBRAIC[52] = (CONSTANTS[375]+( CONSTANTS[375]*ALGEBRAIC[49] - CONSTANTS[374]*ALGEBRAIC[48]))*ALGEBRAIC[50]; ALGEBRAIC[53] = STATES[0]/ALGEBRAIC[1]; ALGEBRAIC[54] = STATES[15]/ALGEBRAIC[24]; ALGEBRAIC[55] = (CONSTANTS[288] - CONSTANTS[289]*ALGEBRAIC[53]*ALGEBRAIC[54]*exp(- ( CONSTANTS[319]*STATES[16])))/( (1.00000+ (10.0000/9.00000)*ALGEBRAIC[53]*exp(- ( CONSTANTS[322]*STATES[16])))*(1.00000+ 18.0000*ALGEBRAIC[54])); ALGEBRAIC[56] = CONSTANTS[48]*ALGEBRAIC[25]; ALGEBRAIC[57] = (STATES[1]==-47.1300 ? 3.20000 : 0.320000*((STATES[1]+47.1300)/(1.00000 - exp( -0.100000*(STATES[1]+47.1300))))); ALGEBRAIC[58] = 0.0800000*exp(- ( STATES[1]*(1.00000/11.0000))); ALGEBRAIC[59] = 1.00000/(ALGEBRAIC[58]+ALGEBRAIC[57]); ALGEBRAIC[60] = (CONDVAR[4]<0.00000 ? ALGEBRAIC[57]*ALGEBRAIC[59] : STATES[5]); ALGEBRAIC[61] = (CONDVAR[5]<0.00000 ? ( 0.135000*exp(-80.0000/6.80000))*exp( (-1.00000/6.80000)*STATES[1]) : 0.00000); ALGEBRAIC[62] = (CONDVAR[6]<0.00000 ? 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)))); ALGEBRAIC[63] = (CONDVAR[7]<0.00000 ? ( -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[64] = (CONDVAR[8]<0.00000 ? 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)))); ALGEBRAIC[65] = pow( STATES[14]*CONSTANTS[290], CONSTANTS[146]); ALGEBRAIC[66] = pow( STATES[26]*CONSTANTS[291], CONSTANTS[149]); ALGEBRAIC[67] = CONSTANTS[135]*(( CONSTANTS[179]*ALGEBRAIC[65] - CONSTANTS[180]*ALGEBRAIC[66])/( (1.00000+ALGEBRAIC[65]+ALGEBRAIC[66])*( ALGEBRAIC[2]*(CONSTANTS[106]+ ALGEBRAIC[1]*CONSTANTS[377])+(1.00000+ ALGEBRAIC[1]*CONSTANTS[378])))); ALGEBRAIC[68] = CONSTANTS[177]*(ALGEBRAIC[7]+STATES[4])*(STATES[27] - STATES[28]); ALGEBRAIC[69] = (STATES[26] - STATES[27])*CONSTANTS[292]; ALGEBRAIC[70] = (STATES[28] - STATES[14])*CONSTANTS[293]; ALGEBRAIC[75] = CONSTANTS[137]*pow( STATES[34]*CONSTANTS[313], CONSTANTS[251]); ALGEBRAIC[76] = CONSTANTS[309]*STATES[30]+ ALGEBRAIC[75]*STATES[50]+ CONSTANTS[333]*STATES[31]; ALGEBRAIC[77] = CONSTANTS[340]*STATES[31]+ ALGEBRAIC[75]*STATES[29]+ CONSTANTS[248]*STATES[30]+ CONSTANTS[334]*STATES[32]; ALGEBRAIC[78] = CONSTANTS[341]*STATES[32]+ CONSTANTS[304]*STATES[31]+ CONSTANTS[335]*STATES[33]; ALGEBRAIC[79] = - ( CONSTANTS[335]*STATES[33])+ CONSTANTS[312]*STATES[32]; ALGEBRAIC[80] = CONSTANTS[137]*STATES[31] - (ALGEBRAIC[75]+CONSTANTS[339])*STATES[29]; ALGEBRAIC[71] = STATES[31]+STATES[29]+STATES[32]+STATES[33]; ALGEBRAIC[72] = CONSTANTS[365]*ALGEBRAIC[71]; ALGEBRAIC[81] = CONSTANTS[98]*STATES[14]*(CONSTANTS[139] - STATES[34]) - ( CONSTANTS[99]*STATES[34])*(1.00000 - (2.00000/3.00000)*ALGEBRAIC[72]); ALGEBRAIC[82] = CONSTANTS[86]*STATES[14]*(CONSTANTS[54] - STATES[35]) - CONSTANTS[85]*STATES[35]; ALGEBRAIC[83] = ALGEBRAIC[81]+ALGEBRAIC[82]; ALGEBRAIC[84] = 1.00000/(1.00000+CONSTANTS[294]/( (STATES[28]+CONSTANTS[107])*(STATES[28]+CONSTANTS[107]))); ALGEBRAIC[85] = 1.00000/(1.00000+CONSTANTS[295]/( (STATES[27]+CONSTANTS[109])*(STATES[27]+CONSTANTS[109]))); ALGEBRAIC[86] = 1.00000/(1.00000+CONSTANTS[294]/( (STATES[14]+CONSTANTS[107])*(STATES[14]+CONSTANTS[107]))); ALGEBRAIC[29] = CONSTANTS[323]*(STATES[16] - 91.0000); ALGEBRAIC[87] = 1.00000+ STATES[14]*CONSTANTS[296]; ALGEBRAIC[88] = ALGEBRAIC[87]*ALGEBRAIC[87]*ALGEBRAIC[87]; ALGEBRAIC[89] = CONSTANTS[298]*STATES[14]*ALGEBRAIC[29]*(ALGEBRAIC[88]/( ( ALGEBRAIC[88]*ALGEBRAIC[87]+CONSTANTS[138]/pow(1.00000+ STATES[14]*CONSTANTS[297], CONSTANTS[143]))*(1.00000 - exp(- ALGEBRAIC[29])))); ALGEBRAIC[90] = CONSTANTS[182]*exp( CONSTANTS[299]*ALGEBRAIC[29])*(STATES[25]/( STATES[14]*pow(1.00000+CONSTANTS[124]/STATES[13], CONSTANTS[142])*(1.00000+CONSTANTS[126]/STATES[25]))); ALGEBRAIC[91] = 0.400000*exp( (STATES[1]+2.00000)*0.100000); ALGEBRAIC[95] = 4.00000*ALGEBRAIC[91]; ALGEBRAIC[96] = 3.00000*ALGEBRAIC[91]; ALGEBRAIC[97] = 2.00000*ALGEBRAIC[91]; ALGEBRAIC[98] = ALGEBRAIC[91]; ALGEBRAIC[93] = CONSTANTS[18]*ALGEBRAIC[91]; ALGEBRAIC[99] = 4.00000*ALGEBRAIC[93]; ALGEBRAIC[100] = 3.00000*ALGEBRAIC[93]; ALGEBRAIC[101] = 2.00000*ALGEBRAIC[93]; ALGEBRAIC[102] = ALGEBRAIC[93]; ALGEBRAIC[92] = 0.0500000*exp( (STATES[1]+2.00000)*(-1.00000/13.0000)); ALGEBRAIC[103] = ALGEBRAIC[92]; ALGEBRAIC[104] = 2.00000*ALGEBRAIC[92]; ALGEBRAIC[105] = 3.00000*ALGEBRAIC[92]; ALGEBRAIC[106] = 4.00000*ALGEBRAIC[92]; ALGEBRAIC[94] = ALGEBRAIC[92]*CONSTANTS[302]; ALGEBRAIC[107] = ALGEBRAIC[94]; ALGEBRAIC[108] = 2.00000*ALGEBRAIC[94]; ALGEBRAIC[109] = 3.00000*ALGEBRAIC[94]; ALGEBRAIC[110] = 4.00000*ALGEBRAIC[94]; ALGEBRAIC[111] = 0.187500*STATES[28]; ALGEBRAIC[112] = ALGEBRAIC[111]; ALGEBRAIC[113] = CONSTANTS[18]*ALGEBRAIC[112]; ALGEBRAIC[114] = CONSTANTS[18]*ALGEBRAIC[113]; ALGEBRAIC[115] = CONSTANTS[18]*ALGEBRAIC[114]; ALGEBRAIC[116] = CONSTANTS[18]*ALGEBRAIC[115]; ALGEBRAIC[117] = CONSTANTS[9]*(STATES[46] - STATES[47]); ALGEBRAIC[118] = CONSTANTS[10]*( (CONSTANTS[13] - STATES[47])*STATES[48] - STATES[47]*(8.00000 - STATES[48])*CONSTANTS[300]); ALGEBRAIC[119] = CONSTANTS[11]*( (CONSTANTS[14] - STATES[46])*STATES[0] - STATES[46]*ALGEBRAIC[1]*CONSTANTS[300]); ALGEBRAIC[120] = (( CONSTANTS[6]*CONSTANTS[7])/( CONSTANTS[191]*CONSTANTS[5]))*log(CONSTANTS[144]/STATES[13]); ALGEBRAIC[125] = CONSTANTS[46]*STATES[5]*STATES[5]*STATES[5]*STATES[6]*STATES[7]*(STATES[1] - ALGEBRAIC[120]); ALGEBRAIC[123] = (( CONSTANTS[6]*CONSTANTS[7])/( CONSTANTS[192]*CONSTANTS[5]))*log(( CONSTANTS[194]*CONSTANTS[127]+ CONSTANTS[195]*CONSTANTS[144])/( CONSTANTS[194]*STATES[12]+ CONSTANTS[195]*STATES[13])); ALGEBRAIC[126] = CONSTANTS[252]*STATES[8]*STATES[8]*((STATES[1] - ALGEBRAIC[123])/(1.00000+exp((STATES[1] - 40.0000)/40.0000))); ALGEBRAIC[121] = (( CONSTANTS[6]*CONSTANTS[7])/( CONSTANTS[192]*CONSTANTS[5]))*log(CONSTANTS[127]/STATES[12]); ALGEBRAIC[124] = STATES[1] - ALGEBRAIC[121]; ALGEBRAIC[127] = 1.02000/(1.00000+exp( 0.238500*(ALGEBRAIC[124] - 59.2150))); ALGEBRAIC[128] = ( 0.491200*exp( 0.0803200*(ALGEBRAIC[124]+5.47600))+exp( 0.0617500*(ALGEBRAIC[124] - 594.310)))/(1.00000+exp( -0.514300*((STATES[1] - ALGEBRAIC[121])+4.75300))); ALGEBRAIC[129] = ALGEBRAIC[127]/(ALGEBRAIC[127]+ALGEBRAIC[128]); ALGEBRAIC[130] = CONSTANTS[310]*ALGEBRAIC[129]*ALGEBRAIC[124]; ALGEBRAIC[131] = CONSTANTS[47]*(STATES[1] - ALGEBRAIC[120]); ALGEBRAIC[132] = CONSTANTS[45]*(ALGEBRAIC[124]/(1.00000+exp( (7.48800 - STATES[1])*CONSTANTS[253]))); ALGEBRAIC[122] = (( CONSTANTS[6]*CONSTANTS[7])/( CONSTANTS[193]*CONSTANTS[5]))*log(CONSTANTS[22]/STATES[14]); ALGEBRAIC[133] = CONSTANTS[44]*(STATES[1] - ALGEBRAIC[122]); } void getStateInformation(double* SI) { SI[0] = 1.0; SI[1] = 1.0; SI[2] = 1.0; SI[3] = 1.0; SI[4] = 1.0; SI[5] = 1.0; SI[6] = 1.0; SI[7] = 1.0; SI[8] = 1.0; SI[9] = 1.0; SI[10] = 1.0; SI[11] = 1.0; SI[12] = 1.0; SI[13] = 1.0; SI[14] = 1.0; SI[15] = 1.0; SI[16] = 1.0; SI[17] = 1.0; SI[18] = 1.0; SI[19] = 1.0; SI[20] = 1.0; SI[21] = 1.0; SI[22] = 1.0; SI[23] = 1.0; SI[24] = 1.0; SI[25] = 1.0; SI[26] = 1.0; SI[27] = 1.0; SI[28] = 1.0; SI[29] = 1.0; SI[30] = 1.0; SI[31] = 1.0; SI[32] = 1.0; SI[33] = 1.0; SI[34] = 1.0; SI[35] = 1.0; SI[36] = 1.0; SI[37] = 1.0; SI[38] = 1.0; SI[39] = 1.0; SI[40] = 1.0; SI[41] = 1.0; SI[42] = 1.0; SI[43] = 1.0; SI[44] = 1.0; SI[45] = 1.0; SI[46] = 1.0; SI[47] = 1.0; SI[48] = 1.0; SI[49] = 1.0; SI[50] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = (int)(VOI) % (int)(CONSTANTS[0]) - CONSTANTS[1]; CONDVAR[1] = fabs(STATES[1]) - CONSTANTS[8]; CONDVAR[2] = fabs(STATES[1]) - CONSTANTS[8]; CONDVAR[3] = fabs(STATES[1]) - CONSTANTS[8]; CONDVAR[4] = ALGEBRAIC[59] - 0.0300000; CONDVAR[5] = STATES[1] - -40.0000; CONDVAR[6] = STATES[1] - -40.0000; CONDVAR[7] = STATES[1] - -40.0000; CONDVAR[8] = STATES[1] - -40.0000; }