Stream Line Plots
var('x, y, t') F=(0.3*x-0.2*x*y, -0.2*y+0.1*x*y) ans=desolve_odeint(F,[1,1],srange(0,40,0.1),[x,y],ivar=t) pp=streamline_plot(F, (x, 0, 4), (y, 0, 4),density=2, color='gray') soln=line(ans,thickness=3) pp+soln
streamline_plot(0.5*y*(1-y/5)*(y-10),(x,0,10),(y,0,15))