function dydt = predprey(t,y)
% The Lotka-Volterra predator-prey model
%   y1d = (1 - alpha*y2)*y1;
%   y2d = (-1 + beta*y1)*y2;
% To run the model try: 
% ode45(@predprey,[0 25],[20;20]);
alpha = 0.02; beta = 0.01;
y1d = (1 - alpha*y(2))*y(1);
y2d = (-1 + beta*y(1))*y(2);
dydt = [y1d;y2d];
% ----- end of function M-file file predprey.m -----


