Our aim in this paper is to propose an SOR-like new iterative method by introducing a relaxation parameter ω to improve the new iterative method proposed by Daftardar-Gejji and Jafari (NIM) [J. Math. Anal. Appl. 316 (2006) 753–763] in order to solve two problems. The first one is the problem of the spread of a nonfatal disease in a population which is assumed to have constant size over the period of the epidemic, and the other one is the problem of prey and predator. The proposed method is not limited to these two problems but can be applicable to a wide range of systems of nonlinear functional problem. The results, for different values of ω, show that we found some known methods and our method compared to methods using the calculation of special polynomials and derivatives like the Adomian decomposition method (ADM), the calculation of the Lagrange multiplier as in the variational iterative method (VIM), or the construction of a homotopy as in the homotopy perturbation method (HPM) has several advantages, such as very effective and very simple to implement. Unfortunately, these methods do not guarantee a valid approximation in large time interval. To overcome this, we applied our method for approximating the solution of the problems in a sequence of time intervals as a multistage approach. Some numerical results are presented with plots according to the parameter ω.