function Y=dyrepop(X)
a=1.01; b=0.00001; c=0.99; d=0.0000075; 
x=X(1); y=X(2);
Y(1)=a*x-b*x*y;
Y(2)=c*y+d*x*y;
