If your aim is to visualize the solution trajectory with phase plane plots, you can plot multiple phase plane plots each at a different time `t` value. You can use quiver function for phase plane plots.
You can also animate these plots in order of increasing time using getframe and movies function.
Hope it helps.
