function xdot=f(x,t)
  xdot=x.*(1-x);