Accuracy analysis of a Bessel spectral function method for the solution of scattering equations