Vlasov-Poisson方程组（VP）描述了在集体场势影响下无碰撞粒子分布的演化。VP是冷暗物质（DM）宇宙学密度扰动引力不稳定性研究的基础，其应用范围还扩展到其他领域，如等离子体物理。对于冷暗物质，每个流体单元（或粒子）只对应一个速度，初始条件存在刚性间断。这会导致当采用基于网格的方法求解VP时出现扩散或负分布函数等问题。本文旨在强调这一问题，并聚焦于该现象的技术细节。通过比较不同的有限体积方法和谱方法，我们观察到虽然所有积分方案都能保持系统的不变量（例如能量），但感兴趣的物理可观测量——密度却未被正确再现。因此，我们将不同欧拉积分方案得到的密度与参考N体方法的结果进行对比。我们指出，求解自引力系统VP系统的最适方法是谱方法。

The Vlasov-Poisson systems of equations (VP) describes the evolution of a distribution of collisionless particles under the effect of a collective-field potential. VP is at the basis of the study of the gravitational instability of cosmological density perturbations in Dark-Matter (DM), but its range of application extends to other fields, such as plasma physics. In the case of Cold Dark Matter, a single velocity is associated with each fluid-element (or particle) , the initial condition presents a stiff discontinuity. This creates problems such as diffusion or negative distribution function when a grid based method is used to solve VP. In this work we want to highlight this problem, focusing on the technical aspects of this phenomenon. By comparing different finite volume methods and a spectral method we observe that, while all integration schemes preserve the invariants of the system (e.g, energy), the physical observable of interest, i.e., the density, is not correctly reproduced. We thus compare the density obtained with the different Eulerian integration schemes with the result obtained from a reference N-body method. We point out that the most suitable method to solve the VP system for a self-gravitating system is a spectral method.