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.