with(DEtools): de:=diff(x(t),t)=x(t)*(1-x(t)): vars:=x(t): wind:=t=0..5,x=-0.2..3: ic:=[[x(0)=1.2],[x(0)=1],[x(0)=0],[x(0)=0.5]]: opts:=stepsize=.2,title=`Verhulst model`,arrows=MEDIUM: DEplot(de,vars,wind,ic,opts);