For solving various equations of the form $v=f+N (v)$, the new iterative method and the new algorithm proposed by V. Daftardar–Gejji et al. [Daftardar–Gejji V., Jafari H. J. Math. Anal. Appl. 316 (2), 753–763 (2006); Kumar M., Jhinga A., Daftardar–Gejji V. Int. J. Appl. Comp. Math.
This study proposes and explores a prey–predator model that presents a functional response to group behavior of prey–predator harvesting. We study a non-linear model of prey–predator growths in two species. The proposed model is supported by theoretical and numerical results. Some numerical descriptions are provided to help our analytical and theoretical conclusions. For all possible parameter values occurring in a prey–predator system, we solved it by using both VIM (variational iteration method) and HPM (homotopy perturbation method). We also used MATLAB coding to compare our approxi