Simulation of N spheres collidiing on an closed vessel

This codes models N polystyrene spheres colliding on an enclosed vessel, with an airflow, through Maxwell-Boltzmann distribution.

This codes models N polystyrene spheres colliding on an enclosed vessel, with an airflow. the first code is where all the motion happened and the positions of the spheres, in every step of time, is saved, the second one animate the movement of the spheres and the third one calculate the velocity in each point of the box using Helmholtz theorem with the approximation of irrotational and incompressible fluid.