Action-angle formalism for axisymmetric potentials using Binney (2012)’s Staeckel approximation, grid-based interpolation
pot= potential or list of potentials
delta= focus of prolate confocal coordinate system
Rmax = Rmax for building grids
nE=, npsi=, nLz= grid size
numcores= number of cpus to use to parallellize
+scipy.integrate.quad keywords
OUTPUT: HISTORY:
2012-11-29 - Written - Bovy (IAS)