function   y=xdot(x,t)

global e r

y=x.*(1-x)-r*(1+e*sin(2*pi*t));
