function[X]=logistic(b,tvec)

[t,X]=ode45(@logistic_rhs,tvec,b(1),[],b(2:3));