clear all % figures for logisitic growth K=100; r=0.8; tim=[0:0.1:20]; N0=100; for kk=1:length(tim) N(kk)=K/(1+(K/N0-1)*exp(-r*tim(kk))); end plot(tim,N,'k')