We present a new method for constructing equilibrium phase models for stellar systems. Applications of the iterative method include both modelling of observational data and the construction of initial condition for N-body simulations.