A new and robust numerical method for Hammerstein functional integral equations is proposed. The convergence and the numerical stability of the method are mathematically proved and tested on some examples.