better implementation of estimate_hz
better implementation of estimate_hz
use Gauss-Legendre integration to calcuate the surface density