“…To perform the multiple derivations, Jagiełło et al have used a virial expansion in order to fit the experimental isotherms. Some examples of distribution functions calculated according to this method are reported in the literature. ,, The high flexibility of this fitting equation, using a limited number of parameters, brings in the well-known default of polynomial fittings, i.e., adding residual oscillations to the experimental data that are greatly amplified by the multiple derivation procedure leading to artifacts …”