In the present work, we develop a calculational method of solving the scattering equations for spherically symmetric potentials by expanding the solutions on Coulomb functions. We utilize a multistep integration scheme together with the standard partial wave analysis in a region where the potential term dominates. The method applies to any physical problem expressed as [∇
2 + V(r) + k
2]ψ(r) = 0, while the extension of the method to more general scattering problems is briefly discussed. At present, we demonstrate a two-step Coulomb-fitted integration scheme by calculating the short-range scattering phase shifts for various potentials V (r).