function [] = LogisticDiscrete(N0,Nmax,kappa,beta)

n=linspace(0,Nmax,Nmax+1);

N(1)=N0;

for i=1:Nmax,
    N(i+1)=N(i)+beta*(kappa-N(i))*N(i);
end

plot(n,N,'o-');
axis([0,Nmax,0,2*kappa])
end

