clear all
N0=10;  % pocetni broj jedinki u populaciji
n=3;    % stopa nataliteta
m=2.5;    % stopa mortaliteta
K=200; % maksimalni broj jedinki
t0=0;   % pocetno t0
tmax=9; % kraj intervala
h=0.1;  % korak
alpha=0.4; 

p=n-m; % stopa prirastaja
r=p;    % r=p(0)

i=1;
for t=t0:h:tmax
    Ne(i)=N0*exp(p*t); %eksponencijalni (Malthusov)
    N(i)=K/(1+(K/N0-1)*exp(-r*t)); % logisticki (Verhulstov)
    Ng(i)=N0*exp(r*(1-exp(-alpha*t)/alpha)); % Gompertz
    i=i+1;    
end
t=t0:h:tmax;
hold on
plot(t,Ne)
plot(t,N)
plot(t,Ng)
xlabel ('Vreme')
ylabel('Broj jedinki')
legend('eksponencijalna','Verhulst','Gombertz')