function f = fdot1d(x,t); global a e f=x*(1-x)-a*(1+e*sin(2*pi*t));