Generated Code

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