In this work we develop the Gaussian quadrature rule for weight functions involving fractional powers, exponentials and Bessel functions of the first kind. Besides the computation based on the use of the standard and the modied Chebyshev algorithm, here we present a very stable algorithm based on the preconditioning of the moment matrix. Numerical experiments are provided and a geophysical application is considered.
