  We present a new method to solve nonlinear Hammerstein equations with weakly
singular kernels. The process to approximate the solution, followed usually,
consists in adapting the discretization scheme from the linear case in order to
obtain a nonlinear system in a finite dimensional space and solve it by any
linearization method. In this paper, we propose to first linearize, via Newton
method, the nonlinear operator equation and only then to discretize the
obtained linear equations by the product integration method. We prove that the
iterates, issued from our method, tends to the exact solution of the nonlinear
Hammerstein equation when the number of Newton iterations tends to infinity,
whatever the discretization parameter can be. This is not the case when the
discretization is done first: in this case, the accuracy of the approximation
is limited by the mesh size discretization. A Numerical example is given to
confirm the theorical result.
