Scheme 2 : dX= dB + beta/2*forces*dt
Zoom: Final time=0.1
initial condition = triangular lattice