galpy.df.quasiisothermaldf.sigmaz2

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

NAME:

sigmaz2

PURPOSE:

calculate sigma_z^2 by marginalizing over velocity

INPUT:

R - radius at which to calculate this

z - height at which to calculate this

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_z^2

HISTORY:

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