I implemented a very simple molecular dynamics engine of hard-sphere particles bouncing around inside a box. With this engine, I demonstrate that particle speeds approach a Maxwell-Boltzmann distribution.
Example of simulation: https://www.youtube.com/watch?v=KZ4s24x_FTQ&feature=youtu.be
