clear all
N0=4.5;  % pocetni broj jedinki/komaraca u populaciji
K0=5; %pocetna velicina stanista/bare
t0=0;   % pocetno t0
tmax=50; % kraj intervala

p=0.3; % stopa prirastaja komaraca
r=p;    % r=p(0) u trenutku t=0

% koeficijenti kojima je opisana periodicna promena velicine bare po
% jedinici vremena (a*cos(bt+c))
a=-0.5;
b=2*pi/12;  % da bi period bio 12 meseci
c=0;

% Periodicna promena okruzenja
f=@(t,x) [r*x(1)*(1-x(1)/x(2)); a*cos(b*t+c)];
[t,x]=ode45(f, [t0,tmax], [N0,K0]);

plot(t,x(:,1),t,x(:,2))
xlabel ('Vreme')
ylabel('Broj jedinki')
legend('N(t)','K(t)')
