function pred_pray
% Модель "Хищник-жертва"
alpha= 1; beta= 1; gamma= 1;delta = 1; % параметры
% правая часть системы
rp= @(t,x)[alpha*x(1)-beta*x(1)*x(2);delta*x(1)*x(2)-gamma*x(2)];
close all
op= odeset('MaxStep',0.005,'Events',@eventf); % максимальный шаг
x0= [0.5 0.5]; % начальные условия
[t,x]= ode45(rp,[0 32],x0,op);
plot(t,x(:,1),t,x(:,2)); % графики
disp(['T= ' num2str(t(end))]);


figure
% фазовые траектории
comet(x(:,1),x(:,2)); % графики
grid on

function [value,isterminal,direction]= eventf(t,xy)
    value= xy(2)-x0(2);
    isterminal= t>0.01;
    direction= -1;
end

end
