Stars are formed in groups or clusters. We perform N-body simulations following the evolution of such systems. Close encounters between stars in clusters transfer energy, leading to the ejection of some of the stars. Binaries and triples (three stars bound together) may also form through dynamical encounters. Here we wish to have software which looks at the output of N-body runs of a small number of stars (10-30) and automatically identifies the bound, sub-systems of stars and the multiple stars (binaries, triples) within each sub-system.