function xdot=fishdot(x,t)

global h e

xdot=x.*(1-x)-h*(1+e*sin(t));

