Abstract
This work proposes a numerical solution of a predator-prey mathematical model in presence of a new ecological interaction which between two predators (male and female) and one prey. This interaction arises in the predator-mating period when the male stays with female and feed together on same prey. In this work, we neglect the diffusion term and obtain a system of non-linear ordinary differential equations (ODES). The numerical solutions are computed using the differential transform method (DTM), adomian decomposition method (ADM), 4th Runge-Kutta method and the multistage differential transform method (MsDTM). As a results, we find that the MsDTM is a reliable and fast method to approximate the solution of the aforementioned model. In addition, the associated ecological interpretation of the results is provided.