“…Now we consider the real HCl reaction (40). For the simulation, we used the information from [43][44][45]: the equilibrium point was c eq = (0.198, 0.004, 0.1995, 0.001, 0.6) and the reaction rate constants of the direct reactions were k + = (10 16 , 10 16 , 1.7 × 10 11 , 1.59 × 10 8 ). This reaction system is very stiff and the equilibrium point is almost on the edge between the vertices (0, 0, 0, 0, 1) and (0.5, 0, 0.5, 0, 0).…”