function xdot = PP(t,x)
global eps beta delta mu
xdot = [x(1).*(1-x(1)- x(2)./(beta+x(1)))/eps; 
        x(2).*(x(1)./(beta+x(1))-delta-mu*x(2))];
        
    if x(1)<0
        x(1)=0;
        xdot(1)=0;
    end
    if x(2)<0
        x(2)=0;
        xdot(2)=0;
    end
   