L  =  0

mulead  =  {m -> 206.771,  Z -> 82,  
        Rad -> 0.0001343,  
        r -> 0.0001343, ϵ -> 
          1000000 ϵ / 27.21139565}

Rinterior  =  
      r^L E^( - 
                4 ϵ m r^2 β^2) Hypergeometric1F1[-γ,  L  +  3/2,  
          8 ϵ m β^2 r^2]
β  =  ((Z/(64 m Rad^3 ϵ^2)))^0.25
Λ  =  ((3 Z)/(8 Rad ϵ)) -  
        (1/4)
γ  =  (1/2) ((Λ/(2 β^2)) -  
            (1/2) -  (L  +  1))

Rexterior  =  W / r
W  =  (E^(- z / 2)) 
        z^(L  +  1) HypergeometricU[a, c, 
          z]
a  =  1  +  L  -  η
c  =  2 L  +  2
z  =  (Sqrt[8 m ϵ]) r
η  =  Z (Sqrt[(m/(2 ϵ))])

F  =  
      D[Rinterior, r]/Rinterior  -  D[Rexterior, r]/Rexterior

F  =  F /. mulead

Plot[F, {ϵ, 1.5, 12.5}, PlotRange -> {-45000, 45000}, 
      AxesLabel -> {"ϵ(MeV)", "F(ϵ)"}]