p=197.327;
m=939;
r0=3;

j1[r_]:=((Sin[r])/r)-Cos[r];

y1[r_]:=Sin[r]+((Cos[r])/r);

f1[r_,e_,U_]:=j1[(Sqrt[2 m (e-U)]) (r/p)];

F1[r_,e_]:=j1[(Sqrt[2 m e]) (r/p)];

G1[r_,e_]:=y1[(Sqrt[2 m e]) (r/p)];

dj1[r_]:=((Cos[r] r-Sin[r])/(r^2))+Sin[r];

dy1[r_]:=Cos[r]+((-Sin[r] r -Cos[r])/(r^2));

df1[r_,e_,U_]:=(( Cos[(Sqrt[2 m (e-U)]) (r/p)] r-(1/((Sqrt[2 m (e-U)])/p))Sin[(Sqrt[2 m (e-U)]) (r/p)])/(r^2))+((Sqrt[2 m (e-U)])/p)  Sin[(Sqrt[2 m (e-U)]) (r/p)];

dF1[r_,e_]:=(( Cos[(Sqrt[2 m e]) (r/p)] r-(1/((Sqrt[2 m e])/p))Sin[(Sqrt[2 m e]) (r/p)])/( (r^2)))+((Sqrt[2 m e])/p)  Sin[(Sqrt[2 m  e]) (r/p)];

dG1[r_,e_]:=((Sqrt[2 m  e])/p)Cos[((Sqrt[2 m  e])/p)r]+((-Sin[((Sqrt[2 m  e])/p) r] r -(1/((Sqrt[2 m e])/p))Cos[((Sqrt[2 m  e])/p) r])/( (r^2)));

δ[e_,U_]:=ArcTan[((f1[r0,e,U] dF1[r0,e]-df1[r0,e,U] F1[r0,e])/(df1[r0,e,U] G1[r0,e]-dG1[r0,e] f1[r0,e,U]))];

TotalSec[k_,E_,U_]:=((2 π)/(((m )/p)(E/p) 100)) (2 k+1) (((Sin[δ[E,U]]))^2);

k4=1;
Plot[{TotalSec[k4,ee,-17.],TotalSec[k4,ee,-18],TotalSec[k4,ee,-19],TotalSec[k4,ee,-20],TotalSec[k4,ee,-22.8],TotalSec[k4,ee,-23]},{ee,0,10},PlotRange->All,PlotStyle->sps]