color = ['red', 'green', 'blue', 'yellow', 'magenta']
plot = plt.figure()
for i, v in enumerate(vect0):
    sol = odeint(system, v, t)
plt.quiver(sol[:-1, (0)], sol[:-1, (1)], sol[1:, (0)] - sol[:-1, (0)], sol[
    1:, (1)] - sol[:-1, (1)], scale_units='xy', angles='xy', scale=1, color
    =color[i])
plt.show(plot)
Drawing phase space trajectories with arrows in matplotlib