(* Input *)
(*kin = 20.;*)
mass= 936.;
range = 3;

sps={{Dashing[None],RGBColor[0,0,0]}, {Dashing[None],RGBColor[1,0,0]},{Dashing[None],RGBColor[0,0,1]},{Dashing[None],RGBColor[0,1,0]},{Dashing[None],RGBColor[0,1,1]}};

kIn[e_] := ((Sqrt[2*mass*e])/197.327)

kOut[e_,U0_] := ((Sqrt[2 mass (e- U0)])/197.327)

ang = 9;
Plot[{e*SphericalBesselJ[ang,e],Sin[e-ang*(Pi/2)]},{e,150,200},PlotStyle->sps]

δ[e_,ang_,r_,u0_]:= ((kIn[e] SphericalBesselJ[ang,kOut[e,u0]r] D[SphericalBesselJ[ang,kIn[e]r],r] - kOut[e,u0]D[SphericalBesselJ[ang,kOut[e,u0]r],r] SphericalBesselJ[ang,kIn[e]r])/(-kOut[e,u0]D[SphericalBesselJ[ang,kOut[e,u0]r],r]SphericalBesselY[ang,kIn[e] r]+kIn[e]SphericalBesselJ[ang,kOut[e,u0]r]D[SphericalBesselY[ang,kIn[e]r],r]))

myDelta[e_,ang_,u0_]:= Simplify[δ[e,ang,r,u0]/.r->3]

myU = -20;
Plot[{myDelta[e,0,myU],myDelta[e,1,myU],myDelta[e,2,myU],myDelta[e,3,myU]},{e,0,20},PlotRange->All,PlotStyle->sps]

CS[e_,ang_,u0_]:=((4 π )/(2m e 100)) (2ang +1) (((Sin[myDelta[e,ang,u0]]))^2)

myAng = 1;
Plot[{CS[ee,myAng,-11],CS[ee,myAng,-10.5],CS[ee,myAng,-11.5],CS[ee,myAng,-12]},{ee,0,1},PlotRange->All,PlotStyle->sps]