This paper treats the interactions among four population species. The system includes one mutuality prey, one harvested prey and two predators. The four species interaction can be described as a food chain, where the first prey helps the second harvested prey. The first and the second predator attack the first and the second prey, respectively, according to Lotka-Volterra type functional responses. The model is formulated using differential equations. One equilibrium point of the model is found and analysed to reveal a threshold that will allow the coexistence of all species. All other equilibrium points of the system are located, with their local and global stability being assessed. To back up the conclusions of the mathematical analysis, a numerical simulation examination of the model is carried out. The system's coexistence can be achieved as long as the harvesting rate of the prey population is lower than its intrinsic growth rate.
