  A global approximation method of Nystr\"om type is explored for the numerical
solution of a class of nonlinear integral equations of the second kind. The
cases of smooth and weakly singular kernels are both considered. In the first
occurrence, the method uses a Gauss-Legendre rule whereas in the second one
resorts to a product rule based on Legendre nodes. Stability and convergence
are proved in functional spaces equipped with the uniform norm and several
numerical tests are given to show the good performance of the proposed method.
An application to the interior Neumann problem for the Laplace equation with
nonlinear boundary conditions is also considered.
