function rhs=shoot2(xspan,x,dummy,beta,n0)

rhs=[x(2)
     (beta-n0)*x(1)];