How to avoid an integrator in Matlab to step out the domain of valid solutions?

Question

I use ode45 to solve a system of ODE\'s. The system describes a planar flow, for which I intend to draw trajectories. The flow covers a convex shape on a plane and theoretically