The STARS code was originally developed in the 1970's by Peter Eggleton. Over the years, the code has been gradually updated to it's current form.
Sverre Aarseth's research into numerical simulations of many-body (N-body) gravitational interactions has developed into a set of FORTRAN codes which describe the dynamics very closely, and freely downloadable.