Analytical DF replaced by a numerical one

This numerical DF is calculated by means of a numerical integration,
which will enable us to add another component to the density profile,
corresponding to the gas, in the future, while mantaining the same
rejection algorithm. The only thing that will have to be changed is the
integrand in the formula for the DF.