function rhs=shoot2(xspan,phi,dummy,n0,beta)
% Right hand side for the schrodinger equation
rhs = [phi(2) ; (beta-n0)*phi(1)];