function[dX]=logistic_rhs(t,X,b)
dX = b(2)*X*(1-X/b(1));