We use numerical simulations to investigate the cusp at the centre of elliptical galaxies, due to the slow growth of a super-massive black hole. We study this problem for axisymmetric models of galaxies, with or without rotation. The numerical simulations are based on the `Perturbation Particles' method, and use GRAPEs to compute the force due to the cusp. We study how the density cusp is affected by the initial flattening of the model, as well as the role played by initial rotation. The logarithmic slope of the density cusp is found to be very much insensitive to flattening; as a consequence, we deduce that tangential velocity anisotropy -which supports the flattening- is also of little influence on the final cusp. We investigate via two different kinds of rotating models the efficiency with which a rotation velocity component builds within the cusp. A cusp in rotation develops only for models where a net rotation component is initially present at high energy levels. The eventual observation of a central rotational velocity peak in E galaxies has therefore some implications for the galaxy dynamical history.

We use a set of Stackel potentials to obtain a local approximation for an effective third integral in axisymmetric systems. We present a study on the feasibility and effectiveness of this approach. We have applied it to three trial potentials of various flattenings, corresponding to nearly ellipsoidal, disky and boxy density isophotes. In all three cases, a good fit to the potential requires only a small set of Stackel potentials, and the associated Stackel third integral provides a very satisfactory, yet analytically simple, approximation to the trial potentials effective third integral.