code: def density(self,R,z,nsigma=None,mc=False,nmc=10000,
                  gl=True,ngl=_DEFAULTNGL,**kwargs):
          
          return self._vmomentdensity(R,z,0.,0.,0.,
                                     nsigma=nsigma,mc=mc,nmc=nmc,
                                     gl=gl,ngl=ngl,
                                     **kwargs)
docs: NAME:
  
             density
  
          PURPOSE:
  
             calculate the density at R,z by marginalizing over velocity
  
          INPUT:
  
             R - radius at which to calculate the density (can be Quantity)
  
             z - height at which to calculate the density (can be Quantity)
  
          OPTIONAL INPUT:
  
             nsigma - number of sigma to integrate the velocities over
  
             scipy.integrate.tplquad kwargs epsabs and epsrel
  
             mc= if True, calculate using Monte Carlo integration
  
             nmc= if mc, use nmc samples
  
             gl= if True, calculate using Gauss-Legendre integration
  
             ngl= if gl, use ngl-th order Gauss-Legendre integration for each dimension
  
          OUTPUT:
  
             density at (R,z)
  
          HISTORY:
  
             2012-07-26 - Written - Bovy (IAS@MPIA)