galpy.df.quasiisothermaldf.density¶
-
quasiisothermaldf.
density
(R, z, nsigma=None, mc=False, nmc=10000, gl=True, ngl=10, **kwargs)[source]¶ NAME:
densityPURPOSE:
calculate the density at R,z by marginalizing over velocityINPUT:
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)