galpy.df.quasiisothermaldf.density

quasiisothermaldf.density(R, z, nsigma=None, mc=False, nmc=10000, gl=True, ngl=10, **kwargs)

NAME:

density

PURPOSE:

calculate the density at R,z by marginalizing over velocity

INPUT:

R - radius at which to calculate the density

z - height at which to calculate the density

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)