function f = logisticmodel(t,p,kmax,pmax)
%
% Calculates rhs of the following ODE
% dp/dt = kmax*(1-p/pmax)*p
%
f = kmax*(1-p/pmax)*p;
