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