Papers

This paper proposes an efficient quadrature rule of precision eleven for the numerical evaluation of improper integrals with endpoint singularities. The proposed method combines Gaussian and Newton–Cotes quadrature rules of lower precision within a unified framework, incorporating cubic Newton divided differences and Romberg extrapolation to enhance computational efficiency and accuracy. An error bound is derived using a mixed quadrature approach, providing a theoretical basis for the performance of the proposed scheme. Numerical experiments and convergence analyses demonstrate its accuracy, stability, and efficiency. The proposed method is further validated through comparisons with existing quadrature techniques and available analytical solutions. Moreover, the Mellin transform is employed to accelerate the numerical computations and improve the accuracy of the resulting approximations. The results indicate that the proposed efficient quadrature scheme provides an effective and reliable approach for the numerical treatment of singular improper integrals