% Run the rabbit-fox model for various values of
% the encounter parameter alpha, plotting each
% solution.
global alpha
for i=2:-1:0,
alpha = 10^(-i)
[t,y] = ode45('rabfox',[0:.1:2], [20,10]);
plot(t,y(:,1),'r',t,y(:,2),'b');
legend('rabbits','foxes')
title(sprintf('alpha = %f',alpha));
pause
end