function dy = PredatorPrey_eqn(t,y)
global a b c d
dy = zeros(2,1); % need to be a column vector
dy(1) = a*y(1)-b*y(1)*y(2);
dy(2) = -c*y(2)+d*y(1)*y(2);