t=0:1:6;
N=[1.99 2.68 3.63 4.89 6.63 8.93 12.10];

razlika=diff(N); %promena kolicine jedinki izmedju 2 vremenska trenutka

plot(t,N,'o'); %da li prati eksponencijalni model?
plot(t,log(N),'o'); % ako prati eksponencijalni model, ovo treba da je (skoro) prava
plot(N(1:end-1),razlika,'o'); %isto kao i red iznad

p=razlika./N(1:end-1); % p ce biti koeficijent prave
%p =0.3467    0.3545    0.3471    0.3558    0.3469    0.3550

% Posto prava nije idealna, treba nam najmanje odstupanje (srednja
% vrednost) od zadatih vrednosti
p=mean(p) %p = 0.3510

plot(t,N,'o',t,1.99*exp(0.351*t)); % nece se poklopiti --> nije dobro p

p=[0 razlika./N(1:end-1)]; % dopunimo sa 0 na pocetku (prava da sadrzi koordinatni pocetak)
p=mean(p)  % p = 0.3009
plot(t,N,'o',t,1.99*exp(0.3009*t)) % sada je ok