Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
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;
}