  We present a significantly improved version of our numerical code JASMINE,
that can now solve the Jeans equations for axisymmetric models of stellar
systems, composed of an arbitrary number of stellar populations, a Dark Matter
halo, and a central Black Hole. The stellar components can have different
structural (density profile, flattening, mass, scale length), dynamical
(rotational support, velocity dispersion anisotropy), and population (age,
metallicity, Initial Mass Function, mass-to-light ratio) properties. These
models, when combined with observations, will allow to investigate important
issues, such as quantifying the systematic effects of IMF variations, of
mass-to-light ratio gradients, and of different stellar kinematic components
(e.g. counter rotating disks, kinematically decoupled cores) on
luminosity-weighted properties. The developed analytical and numerical
framework aims at modeling Early-Type Galaxies, but it can also be applied to
dwarf Spheroidal galaxies and Globular Clusters.
