AbstractIn this paper, we consider the computation of the modified moments for the system of Laguerre polynomials on the real semiaxis with the Hermite weight. These moments can be used for the computation of integrals with the Hermite weight on the real semiaxis via product rules. We propose a new computational method based on the construction of the null-space of a rectangular matrix derived from the three-term recurrence relation of the system of orthonormal Laguerre polynomials. It is shown that the proposed algorithm computes the modified moments with high relative accuracy and linear complexity. Numerical examples illustrate the effectiveness of the proposed method.