Monte Carlo method is used to calculate an equilibrium distribution of interacting dipole hard-spheres in the presence of gravitational field