function [Ydot] = predprey(t,Y)
%r=2; K=1.4; s=1; c=0.5; a=2; d=0.1;
global time IC r K s c a d;
Ydot(1) = r*Y(1)*(1-Y(1)/K)-s*Y(1)*Y(2)/(a*Y(1)+1);
Ydot(2) = Y(2)*(c*s*Y(1)/(a*Y(1)+1)-d);
Ydot=Ydot';