function dpop=beach(t,pop,b,T)
% I assume pop contains: females, searching males, busy males
dpop=[-b*pop(1)*pop(2); -b*pop(1)*pop(2)+pop(3)/T; b*pop(1)*pop(2)-pop(3)/T];