  We consider a macroscopic model describing a system of self-gravitating
particles. We study the existence and uniqueness of non-negative stationary
solutions and allude the differences to results obtained from classical
gravitational models. The problem is considered on a bounded domain up to three
space dimension, subject to Neumann boundary condition for the particle
density, and Dirichlet boundary condition for the self-interacting potential.
Finally, we show numerical simulations that affirm our findings.
