Schroedinger Equation for Yukawa Potential

Clear["Global`*"]

v[x_] = 
      - g ((E^(-x/r))/x) + ((l (l + 1))/(2 m (x^2)))

v[x] ψ[x] -(1/2) m (((ψ'))^(′ ))[ x]==en ψ[x]

eq  =  %

m = 1;

r = 1 ;

g = 2;

l = 0;

en = -0.18;

en1 = -0.16;

en2 =  -0.20;

eq

ulim  =  6;

x0 =  0.001;

plotpot =  Plot[ v[x],  {x, x0, ulim}];

NDSolve[ {eq, ψ[x0]==1, ψ'[x0]==1}, ψ, {x,x0,ulim} ] ;

plot1 = Plot[ Evaluate[ ψ[x] /. % ], {x,x0,ulim} ] ;

en  =  en1;

NDSolve[ {eq, ψ[x0]==1, ψ'[x0]==1}, ψ, {x,x0,ulim} ] ;

\(\(plot2\  = \ 
      Plot[\ Evaluate[\ ψ[x]\  /. \ %\ ], \ {x, x0, ulim}, 

        \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
        PlotStyle\  -> \ Dashing[{0.01, 0.01}]\ ]\ ; \)\)

en  =  en2;

NDSolve[ {eq,  ψ[x0] == 1,  ψ'[x0] == 1},  ψ, 
       {x, x0, ulim} ] ;

\(\(plot3\  = \ 
      Plot[\ Evaluate[\ ψ[x]\  /. \ %\ ], \ {x, x0, ulim}, 

        \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
        PlotStyle\  -> \ Dashing[{0.01, 0.01}]\ ]\ ; \)\)

Show[plotpot, plot1, plot2, plot3];