function v=xdot(x,t)
global r

v=x.*(r-x);