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