function xdot=fish(x,t)

global a e

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

%xdot=-a*x+e*sin(2*pi*t);