The current research approximates the unknown function based on the normalized Müntz−Legendre polynomials (NMLPs) in conjunction with a spectral method for the solution of nonlinear Fredholm and Volterra integral equations. In this method, by using operational matrices, a system of algebraic equations is derived that can be readily handled through the use of the Newton scheme. The stability, error bound, and convergence analysis of the method are discussed in detail by preparing some theorems. Several illustrative examples are provided formally to show the efficiency of the proposed method. Manuscript profile