function xdot=f(x,t) %Logistic vector field xdot=x.*(1-x);