<p>A novel Legendre polynomial-based inverse fast Fourier transform (IFFT) technique, termed the Gauss-Legendre-IFFT, is presented for gravity forward modeling. This method combines the high-precision numerical integration of Gauss-Legendre quadrature with the computational efficiency of the IFFT, significantly reducing the truncation and aliasing errors inherent in conventional FFT-based approaches. Numerical experiments, including an infinite horizontal cylinder model and a right rectangular prism model, demonstrate that the proposed method achieves significantly higher accuracy compared to conventional frequency-domain methods while maintaining competitive computational efficiency. The Gauss-Legendre-IFFT scheme offers an efficient and reliable tool for high-resolution gravity forward modeling, facilitating advanced interpretation and inversion in exploration geophysics.</p>