This is a simple N-body code, written in python. The aim is to work out the merger rate of black holes distributed in small clusters or halos.