Added a main in which we show how a system can equilibrates, and how many iterations are needed to equilibrate