  Schwarzschild models of a Plummer sphere were constructed and compared with
analytical results to test the robustness of Schwarzschild's method and its
efficiency in probing the degeneracy of the solution space. The method was then
used to construct and study the nonuniqueness of triaxial equilibria with
central density cusps, for which no analytical solutions are known. The
stability of some models was also tested via N-body simulations.
