n0 10 r 0.5 k 100 r1 1.5 nt n0 nt*exp(r*(1-nt/k)) qt n0 qt*r1