function [gstop,isterminal,direction]=g(t,y,a)
gstop=y(1)-a(2);     %Measures when y(1)=Initial population
isterminal=0;   %Should the solver stop? 0 means no
direction=-1;   %What direction should the function be going for event finder to look for y(1)=0?
