A numerical method for solving the inverse scattering problem with fixed-energy phase shifts