galpy.df.quasiisothermaldf.jmomentdensity

quasiisothermaldf.jmomentdensity(R, z, n, m, o, nsigma=None, mc=True, nmc=10000, _returnmc=False, _vrs=None, _vts=None, _vzs=None, **kwargs)

NAME:

jmomentdensity

PURPOSE:

calculate the an arbitrary moment of an action of the velocity distribution at R times the surfacmass

INPUT:

R - radius at which to calculate the moment(/ro)

n - jr^n

m - lz^m

o - jz^o

OPTIONAL INPUT:

nsigma - number of sigma to integrate the velocities over (when doing explicit numerical integral)

mc= if True, calculate using Monte Carlo integration

nmc= if mc, use nmc samples

OUTPUT:

<jr^n lz^m jz^o x density> at R

HISTORY:

2012-08-09 - Written - Bovy (IAS@MPIA)