(* Input *)
(*kin = 20.;*)
U0 = -40.;
mass= 939.;
range = 3;

kIn[ee1_] := (Sqrt[2*mass*ee1])

kOut[ee1_] := (Sqrt[2 mass (ee1- U0)])

phiIn[A_,kk_,rr_] := A [((Sin[kk rr])/(kk rr))- Cos[kk rr]]

phiOut[del_,qq_,rr_] := ((Sin[qq*rr-del])/(qq rr))-Cos[qq*rr-del]

phiIn1[A_,kk_,rr_] := A  [((Cos[kk rr])/rr)-((Sin[kk rr])/(kk (rr^2)))+kk Sin[kk rr]]

phiOut1[del_,qq_,rr_] := ((Cos[qq rr-del])/rr)-((Sin[qq rr-del])/(qq (rr^2)))+qq Cos[qq rr-del]

Print[kIn[ee],"  ",phiIn[aa,kIn[ee],range]]

Print[kOut[ee],"  ",phiOut[aa,kOut[ee],range]]

Solve[phiIn[aa,kIn[ee],range]==phiOut[dd,kOut[ee],range]&&phiIn1[aa,kIn[ee],range]==phiOut1[dd,kOut[ee],range],{aa,dd}]

Print[((phiIn1[A,kIn[ee1],(range/197.327)])/(phiIn[A,kIn[ee1],(range/197.327)]))]

Print[((phiOut1[A,kOut[ee1],(range/197.327)])/(phiOut[A,kOut[ee1],(range/197.327)]))]

LeftS[ee1_]:=Simplify[((phiIn1[A,kIn[ee1],(range/197.327)])/(phiIn[A,kIn[ee1],(range/197.327)]))]

RightS[ee1_]:=Simplify[((phiOut1[A,kOut[ee1],(range/197.327)])/(phiOut[A,kOut[ee1],(range/197.327)]))]

roots = Solve[LeftS[ee1]==RightS[ee1],ee1,Assumptions->ee1>0]

Print[roots[[1]]]

CrossSection[ee_]:=((4 Pi)/((kIn[ee1])^2))((Sin[roots[[1]]])^2)

Plot[CrossSection[ee],{ee,-10,30},PlotRange->All]