function yprime = rabfox(t,y);
% rabbit and fox preditor-prey model - FMM P. 6.4
% return a column vector
yprime = [ 2*y(1)-y(2)*y(1)/100;  -y(2)+y(2)*y(1)/100];