def density(self,R,z,nsigma=None,mc=False,nmc=10000,
                gl=True,ngl=_DEFAULTNGL,**kwargs):
        """
        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)

        """
        return self._vmomentdensity(R,z,0.,0.,0.,
                                   nsigma=nsigma,mc=mc,nmc=nmc,
                                   gl=gl,ngl=ngl,
                                   **kwargs)
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 )
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)