function[X]=num_logistic(b,tvec,yvec)

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