Toomre model to simulate galaxy interactions and collisions. Applies a second-order finite difference approximation.