function ShowRabbits()
close all
% Model parameter
alfa = .01;
% Initial Population Densities
R(1) = 300;
F(1) = 150;
% Duration of simulation in years...
Y = 10;
% Time step length...
deltaT = 1/365;
% Length of result vectors..
N = 365*Y+1;
% The simulation...
for k=1:N-1
   % The populations at the current time step...
   Rc = R(k);
   Fc = F(k);
   % Model the rates of change...
   deltaR = 2*Rc - alfa*Rc*Fc;
   deltaF =  -Fc + alfa*Rc*Fc;
   % Estimate the populations at the next time step...
   R(k+1) = Rc + deltaT*deltaR;
   F(k+1) = Fc + deltaT*deltaF;
end
plot(1:N,R,'m',1:N,F,'k')
title(sprintf('Predator-Prey Simulation with alfa = %5.3f',alfa),'fontsize',16)
legend('Rabbits','Foxes')
shg