Click to view movie
The initial parameters are x,y and z average
position and velocity,gaussian scale lenght and
velocity dispersion value,time step and interval.
There is also and obtion for the type of calculation,
mass model and nature of gravitational field
There is also an obtion to include drag

The runs I did used the Dehnen and Binney gravitational
field model and no drag.