clear all

a= 2;   % prirodni prirastaj plena
b=1.1;  % stopa umanjenja plena usled grabljivaca
c=1;    % natalitet grabljivaca
m=0.9;  % mortalitet grabljivaca
tp = 16; % u trenutku t_p lovci ulove 50% od obe vrste
P0=2;
G0=1;

%P0=m/c; % pocetna populacija plena (stac.)
%G0=a/b; % pocetna populacija grabljivaca (stac.)


f=@(t,x) [a*x(1)-b*x(1)*x(2); c*x(1)*x(2) - m*x(2)];
[t,x]=ode45(f, [0,50], [P0,G0]);
%[t,x]=ode45(f, [0,tp], [P0,G0]);
%[t1,x1]=ode45(f, [tp,50], [x(end,1)*0.5,x(end,2)*0.5]);

%t=[t; t1];
%x=[x;x1]

figure(1)
plot(t,x(:,1),t,x(:,2));
xlabel('vreme')
ylabel('populacija')
legend('plen','grabljivice')

figure(2)
plot(x(:,1),x(:,2))
%plot(x(:,1),x(:,2),'o')
xlabel('plen')
ylabel('grabljivice')
