demo_velocity_distribution_2.py

import thermosim

b = thermosim.Box.generic(N=500)
b.show('velocities')
b.dt = 1.
b.v[:, 1] = 0.
b.v[:, 0] = 1.
b._init()
b.set_fig_position(0, 0, 1286, 1150)
b.run(100000)