Wrote function which takes the 'old particle data', the cooling data, and a timestep as arguments, and returns a new internal energy. This function would have to be different depending on the cooling model used. For now we just have a constant cooling rate for all particles.
and try it :) compare to rgb's attempt
Another functions takes the 'new internal energy' and updates the entropy.