galpy.df.quasiisothermaldf.sigmaRz

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

NAME:

sigmaRz

PURPOSE:

calculate sigma_RZ^2 by marginalizing over velocity

INPUT:

R - radius at which to calculate this (can be Quantity)

z - height at which to calculate this (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:

sigma_Rz^2

HISTORY:

2012-07-30 - Written - Bovy (IAS@MPIA)