I am trying to model some elastic collisions which i can then transform into an SIR model, however for some reason my vertical collisons with the boundary of the domain (in
