function f=dxdt_1d(x,t);

global a e

f=x*(1-x)-a*(1-e*sin(2*pi*t));