Research Article  Open Access
Noufe H. Aljahdaly, Manar A. Alqudah, "Analytical Solutions of a Modified PredatorPrey Model through a New Ecological Interaction", Computational and Mathematical Methods in Medicine, vol. 2019, Article ID 4849393, 7 pages, 2019. https://doi.org/10.1155/2019/4849393
Analytical Solutions of a Modified PredatorPrey Model through a New Ecological Interaction
Abstract
The predatorprey model is a common tool that researchers develop continuously to predict the dynamics of the animal population within a certain phenomenon. Due to the sexual interaction of the predator in the mating period, the males and females feed together on one or more preys. This scenario describes the ecological interaction between two predators and one prey. In this study, the nonlinear diffusive predatorprey model is presented where this type of interaction is accounted for. The influence of this interaction on the population of predators and preys is predicted through analytical solutions of the dynamical system. The solutions are obtained by using two reliable and simple methods and are presented in terms of hyperbolic functions. In addition, the biological relevance of the solutions is discussed.
1. Introduction
Many phenomena in natural science, biology, physics, or engineering are studied by developing a mathematical model that consists of partial differential equations. A preypredator system is a wellknown mathematical model for studying the species’ population and density [1, 2]. Many different interactions in this model are very significant phenomena in nature’s population [3, 4]. Alqudah [5] modified the diffusive predatorprey model in which two predators interact with one or more preys in the mating period as follows:
The biological meaning of each term is presented in Table 1. We assume that for all are positive constants and , which means that the predator is moving toward the prey.

This model focuses on the term which refers to a positive interaction with respect to the predator species, like sexual interaction in the mating period, to proliferate the predators. Thus, the males and females are together and feed on the same prey or more, but this term is negative with respect to the prey. In [5], the existence and the uniqueness of the solution of model (1) are studied, and the scaling is introduced to get the dimensionless parameters. Thus, model (1) can be expressed aswhere and .
One way to understand the applications and dynamic system is finding the analytical or approximate solutions. In the field of mathematical computation, several methods have been developed to find the analytical solutions such as the method [6], the Riccati equation expansion method [7], the ()expansion method [8, 9], the algebraic method [10], and so on. In this paper, the exact and stable solutions of the system are obtained using two sufficient methods: (i) the ()expansion method [11] and (ii) the generalized auxiliary equation method [12].
This paper is organized as follows. Section 2 is the description of two employed methods. Section 3 presents the analytical solutions of the considered dynamical system. Section 4 is the discussion of the results, including the biological interpretation. The last section is the conclusion of the work and the results.
2. Description of Algorithms
This section shows briefly the steps of applying the functionexpansion method for finding analytical solutions. We chose two methods that are simple and reliable using the asymptotic software (MATHEMATICA), which are the ()expansion method and the generalized auxiliary equation method. Let us consider the nonlinear evaluation equation of the formwhere and are polynomials in and their derivatives and and are two unknown functions. The methods that are considered in this article have the common following steps:(i)Transforming the function and to and , respectively, by applying the transformation , where c is the wave speed. Thus, the system of equations (3) is reduced to the ordinary differential equation (ODE)(ii)Assuming the function and are expressed by the following polynomials: where and are the degrees of the polynomial of and , respectively. We obtain the value of n and m by the homogeneous balance theory [13–15].(iii)Substituting the polynomials (5) into equation (4) to obtain polynomials of the function .(iv)Equating the confections of (or m) to zero to obtain a system in terms of variables and c.(v)Solving the system to find the value of and c.(vi)Substituting the obtained variables and the value of the function into equation (5) to construct the solutions of U and V.
The main difference between the methods utilized in this paper is the value of the function .(1)The ()expansion method where satisfies and G satisfies the second order ODE: The function is the general solution of the equation (6) as follows:(2)The generalized auxiliary equation method where F satisfies the auxiliary equation The general solutions of the auxiliary equation (9) have several types of solutions depending on the value of , . The following are some types of the solution and the reader is referred to [10, 12] to find more cases.(I)If , the solution of equation (9) is expressed as follows:
where .(II)If ,(III)If ,
3. Analytical Solutions of the Diffusive PredatorPrey System
In this section, the analytical solutions of model (1) will be found by the methods that are described in the previous section. First, introducing the traveling wave solution reduces model (1) aswhere , , and c is the wave speed. Then, the equations (13a) and (13b) are combined as follows:
Following the steps of the algorithm in Section 2, the solutions of and can be expressed as follows:
Substituting equation (15) into equation (14) to obtain
3.1. The ()Expansion Method
In order to utilize the ()expansion method, the definition of equation (6) is applied to equation (16), and then the following algebraic system is obtained:
Hence, two cases of the value of and are realized:(i)Case I(ii)Case II
The solutions within both cases exist if . As we see in Figure 1, for . Therefore, the solution of the system is in terms of the hyperbolic function and parameters and . Assuming yields which is not a reasonable case in actual situations. Therefore, the solution is considered only for , (). The following is the solution by the ()expansion method which is a kink soliton solution for U and V as we see in Figure 2:
(a)
(b)
3.2. The Generalized Auxiliary Equation Method
This subsection presents the solution of the considered system by applying the auxiliary equation method. Thus, we use the definition of (9) into equation (16) which yields to the following algebraic systems:
The solutions of the algebraic system aforementioned give , and and are arbitrary constants. Hence, the appropriate solution of the auxiliary equation is in terms of the hyperbolic function. Assuming yields the solution of the problem in the following expression:
The computed solution by the generalized auxiliary equation method exists if and is a soliton solution in terms of and (see Figure 3).
4. Discussion
4.1. Biological Implication
The solution is obtained by the ()expansion method for and with the condition . Thus, this solution is obtained when the decay rate of the prey population is greater due to the interaction between two predators and one prey () than the decay rate of the prey population due to the competition on the food supply over the growth rate of the prey (). The solution in Figure 2 shows the predator population increases because of plentiful prey and sexual interaction while the prey population decline because of its high consumption during the mating period of the predator in a close environmental area. Ultimately, the predator population will become dominant in the area. This situation is expected when the predators overgraze in the mating period where the prey is plentiful.
Figure 3 presents the solution by applying the generalized auxiliary equation method for . This solution depicts that in certain environmental area, the prey population u will grow due to the absence of the predator . In the mating period, the predators graze in where the prey density (u) is large and as a result u will decay. However, the obtained solutions by both considered methods remain constant away from places of grazing predator during the predatormating period and over time.
4.2. Connection with the Previous Studies
The diffusive predatorprey model has been solved numerically and analytically in some papers. The analytical solutions of the model when both species having the same diffusivities and the prey density being of the Allee type are computed by the ()expansion method [16], and the results are either single structure or periodic structure. Also, it was computed by the improved Riccati equation mapping method [17], and the solutions are presented in three structures: single, periodic, and kink. The same model was solved analytically by utilizing the expexpansion method [18], and the obtained solutions were kink solutions, singular kink solutions, dark soliton solutions, bright soliton solutions, soliton solutions, singular soliton solutions, multiple solitonlike solutions, and triangular periodic solutions, but all these solutions were not discussed in terms of the actual situation.
The numerical solution of a system of nonlinear Volterra differential equations governing on the problem of prey and predator was solved numerically by the Adomian decomposition method [19], Runge–Kutta–Fehlberg method and Laplace Adomian decomposition method [20], and the results state that the number of predator increases, as the number of prey decreases.
The present paper studied the predatorprey model with adding the term and in absence of diffusion among preys . We obtained that predator’s population increases, as the prey’s population decreases in a certain interval, and the solutions are constant out of this interval. Comparing the results obtained in the present study with different modifications of the diffusive predatorprey model, it can be concluded that obtained results are new and different.
5. Conclusion
In this paper, both the ()expansion method and the generalized auxiliary equation method have been used to find new stable analytical solutions of the modified nonlinear diffusive predatorprey model by considering the term . To the best of the researchers’ knowledge, this modified model has not been solved analytically or numerically in previous research. The result showed the effect of the term in the diffusive predatorprey dynamic system which leads to a decline in the prey population and an increase in the predator population in the mating period at predator’s density points. However, the idea of this modified system can be used in epidemiological researches to construct several models. For example, if we consider u as the infected human body by any disease, in many cases we use two types of medicine which is considered as . Thus, the interaction term refers to combatting the epidemic or the disease through two types of medicine. In addition, the preypredator system has been studied widely in the literature regarding the time delay [21–23]. In our future work, we will study time delay for our new modified system with the effect of fluctuation of the biological parameters or fuzzy parameters.
Data Availability
The data supporting this research are from previously reported studies, which have been cited.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Authors’ Contributions
Noufe Aljahdaly implemented the analytical solutions with software and analyzed the data. Manar Alqudah proposed the mathematical model of this paper. All authors wrote this paper, read, and approved the final manuscript.
Acknowledgments
This research was funded by the Deanship of Scientific Research at Princess Nourah Bint Abdulrahman University through the FastTrack Research Funding Program.
References
 A. J. Lotka, “Elements of physical biology,” Science Progress in the Twentieth Century, vol. 21, no. 82, pp. 341–343, 1926. View at: Google Scholar
 V. Volterra, “Variations and Fluctuations of the Number of Individuals in Animal Species Living Together,” Animal Ecology, pp. 409–448, 1926. View at: Google Scholar
 O. P. Misra, P. Sinha, and C. S. Sikarwar, “Dynamical study of a preypredator system with a commensal species competing with prey species: effect of time lag and alternative food source,” Computational and Applied Mathematics, vol. 34, no. 1, pp. 343–361, 2015. View at: Publisher Site  Google Scholar
 Q. Wang, Z. Liu, X. Zhang, and R. A. Cheke, “Incorporating prey refuge into a predatorprey system with imprecise parameter estimates,” Computational and Applied Mathematics, vol. 36, no. 2, pp. 1067–1084, 2017. View at: Publisher Site  Google Scholar
 M. A. Al Qudah, “Existence, uniqueness solution of a modified predatorprey model,” Nonlinear Analysis and Differential Equations, vol. 4, pp. 669–677, 2016. View at: Publisher Site  Google Scholar
 S. A. Elwakil, S. K., M. A. Ellabany, R. Sabry, and R. Sabry, “Modified extended tanhfunction method for solving nonlinear partial differential equations,” Physics Letters A, vol. 299, no. 23, pp. 179–188, 2002. View at: Publisher Site  Google Scholar
 C. Yong, L. Biao, and Z. HongQing, “Symbolic computation and construction of solitonlike solutions to the (2+1)dimensional breaking soliton equation,” Communications in Theoretical Physics, vol. 40, no. 2, pp. 137–142, 2003. View at: Publisher Site  Google Scholar
 K. Zhouzheng, “()expansion solutions to MBBM and OBBM equations,” Journal of Partial Differential Equations, vol. 28, no. 2, pp. 158–166, 2015. View at: Publisher Site  Google Scholar
 N. H. Aljahdaly, “Some applications of the modified ()expansion method in mathematical physics,” Results in Physics, vol. 13, Article ID 102272, 2019. View at: Publisher Site  Google Scholar
 E. Fan, “Multiple travelling wave solutions of nonlinear evolution equations using a unified algebraic method,” Journal of Physics A: Mathematical and General, vol. 35, no. 32, pp. 6853–6872, 2002. View at: Publisher Site  Google Scholar
 M. Wang, X. Li, and J. Zhang, “The ()expansion method and travelling wave solutions of nonlinear evolution equations in mathematical physics,” Physics Letters A, vol. 372, no. 4, pp. 417–423, 2008. View at: Publisher Site  Google Scholar
 E. Yomba, “A generalized auxiliary equation method and its application to nonlinear KleinGordon and generalized nonlinear CamassaHolm equations,” Physics Letters A, vol. 372, no. 7, pp. 1048–1060, 2007. View at: Publisher Site  Google Scholar
 M. Wang, “Exact solutions for a compound KdVBurgers equation,” Physics Letters A, vol. 213, no. 56, pp. 279–287, 1996. View at: Publisher Site  Google Scholar
 M. Wang, “Solitary wave solutions for variant Boussinesq equations,” Physics Letters A, vol. 199, no. 34, pp. 169–172, 1995. View at: Publisher Site  Google Scholar
 M. Wang, Y. Zhou, and Z. Li, “Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics,” Physics Letters A, vol. 216, no. 15, pp. 67–75, 1996. View at: Publisher Site  Google Scholar
 R. A. Kraenkel, K. Manikandan, and M. Senthilvelan, “On certain new exact solutions of a diffusive predatorprey system,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 5, pp. 1269–1274, 2013. View at: Publisher Site  Google Scholar
 H. Kim and J. H. Choi, “Exact solutions of a diffusive predatorprey system by the generalized Riccati equation,” Bulletin of the Malaysian Mathematical Sciences Society, vol. 39, no. 3, pp. 1125–1143, 2016. View at: Publisher Site  Google Scholar
 M. N. Alam and C. Tunc, “An analytical method for solving exact solutions of the nonlinear Bogoyavlenskii equation and the nonlinear diffusive predatorprey system,” Alexandria Engineering Journal, vol. 55, no. 2, pp. 1855–1865, 2016. View at: Publisher Site  Google Scholar
 J. Biazar and R. Montazeri, “A computational method for solution of the prey and predator problem,” Applied Mathematics and Computation, vol. 163, no. 2, pp. 841–847, 2005. View at: Publisher Site  Google Scholar
 S. Paul, S. P. Mondal, and P. Bhattacharya, “Numerical solution of Lotka Volterra prey predator model by using RungeKuttaFehlberg method and Laplace Adomian decomposition method,” Alexandria Engineering Journal, vol. 55, no. 1, pp. 613–617, 2016. View at: Publisher Site  Google Scholar
 D. Pal and G. S. Mahapatra, “Effect of toxic substance on delayed competitive allelopathic phytoplankton system with varying parameters through stability and bifurcation analysis,” Chaos, Solitons & Fractals, vol. 87, pp. 109–124, 2016. View at: Publisher Site  Google Scholar
 D. Pal, G. S. Mahapatra, and G. P. Samanta, “New approach for stability and bifurcation analysis on predatorprey harvesting model for interval biological parameters with time delays,” Computational and Applied Mathematics, vol. 37, no. 3, pp. 3145–3171, 2018. View at: Publisher Site  Google Scholar
 D. Pal, G. S. Mahapatra, and G. P. Samanta, “A study of bifurcation of preypredator model with time delay and harvesting using fuzzy parameters,” Journal of Biological Systems, vol. 26, no. 2, pp. 339–372, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Noufe H. Aljahdaly and Manar A. Alqudah. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.