In this paper, we use the Nyström method based on a quadrature formula associated with non-uniform quasi-interpolation for solving Hammerstein integral equations. This method results in a system of nonlinear equations for determining approximate values of the exact solution of the integral equation. Numerical examples are given for illustrating the performance of the method.