% Master
% X - a vector with two components telling me the population of
%     rabbits and coyotes
% X(1) - rabbits
% X(2) - coyotes

%
% The initial population of rabbits = 100
% The initial population of coyotes = 5
%
X0 = [100; 5];

%
% The time points. We chose evenly spaces time points from t=0 to t=10
%
T = linspace(0, 10, 1000);

%
% A call to lsode. Produces the solution in the form of
% a sequence of column vectors of height 2, one per time point.
%
y=lsode ("rabbits2", X0, T);

% Phase plot
plot(y(:,1),y(:,2));


pause(5);

% Time plot
hold
% Rabbits
plot(T, y(:,1));		

% Coyotes
plot(T, y(:,2));		


pause(5);
