Abstract: The analytical continuation method is used for evaluating efficiently and accuratyintegrals, such as∫₀^xr²χ₁(kr)χ₁,(k'r)χ_L(pr)dr, where χ₁(x) can be a spherical Bessel function j_l(x), or a spherical Neumanm function n_j(x). In particular, a practical program based on the method mintioned above has been made and its results show that the present method is very efficient and highly accurate.