Self-consistent axisymmetric Sridhar-Touma models

We construct phase-space distribution functions for the oblate, cuspy mass models of Sridhar & Touma, which may contain a central point mass (black hole) and have potentials of St"ackel form in parabolic coordinates. The density in the ST models is proportional to a power $r^{-\gamma }$ of the radius, with $0<\gamma <1$. We derive distribution functions $f(E,L_z)$ for the scale-free ST models (no black hole) using a power series of the energy $E$ and the component $L_z$ of the angular momentum parallel to the symmetry axis. We use the contour integral method of Hunter & Qian to construct $f(E,L_z)$ for ST models with central black holes, and employ the scheme introduced by Dejonghe & de Zeeuw to derive more general distribution functions which depend on $E$, $L_z$ and the exact third integral $I_3$. We find that self-consistent two- and three-integral distribution functions exist for all values $0<\gamma <1$.