Research Article  Open Access
Atika Radid, Karim Rhofir, "SORLike New Iterative Method for Solving the Epidemic Model and the Prey and Predator Problem", Discrete Dynamics in Nature and Society, vol. 2020, Article ID 9053754, 13 pages, 2020. https://doi.org/10.1155/2020/9053754
SORLike New Iterative Method for Solving the Epidemic Model and the Prey and Predator Problem
Abstract
Our aim in this paper is to propose an SORlike new iterative method by introducing a relaxation parameter to improve the new iterative method proposed by DaftardarGejji 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 .
1. Introduction
Mathematical models are an essential tool for understanding, analysing, and describing results from natural phenomena modeling. One of the important topics studied in the field of mathematical modeling is epidemiology, because it can provide very useful results and help public health agencies better manage epidemics.
Kermack [1], a pioneer in epidemiological modeling, proposed a first epidemic model, called the SIR model, describing the spread of an infectious disease following a division of population into three distinct categories: susceptible , infected , and recovered .
The problem of spreading of a nonfatal disease in a population which is assumed to have constant size over the period of the epidemic is considered in [2–4]. At time t, suppose the population consists of susceptible population: those so far uninfected and therefore liable to infection infective population: those who have the disease and are still at large isolated population, or who have recovered and are therefore immune
Assume that there is a steady constant rate between susceptible population and infective population and that a constant proportion of these constants results in transmission. Then, in time , of the susceptible population becomes infective, wherewhere is a positive constant. If is the rate at which current infective population becomes isolated, then
The number of new isolated population, , is given by
Now let . Then, the following system determines the progress of the disease:with the following initial conditions:
To solve this problem, Biazar [2, 5] used the Adomian decomposition method (ADM), Rafei et al. [3–6] used the variational iterative method (VIM) and the homotopy perturbation method (HPM), and ALJawary [7] used the new iterative method (NIM).
There are some rabbits and foxes living together. Foxes eat the rabbits and rabbits eat clover, and there is an increase and decrease in the number of foxes and rabbits. For more details on the mathematical modeling leading to the following system of nonlinear equations governing the problem, we refer the readers to [6, 8]:where and are, respectively, the populations of the rabbits and the foxes at time . To solve the prey and predator problem, Biazar et al. [2, 5] and Chowdhury et al. [9] used the Adomian decomposition and the power series methods, Rafei et al. [3] used the variational iterative method (VIM), and ALJawary [7] used the new iterative method (NIM).
Recently, many attempts have been made to develop analytic and approximate methods to solve the epidemic model. Although such methods have been successfully applied, some difficulties have appeared, for examples, constructing a homotopy as in the HPM [10, 11] and solving the corresponding the algebraic equations in calculating Adomian polynomials to handle the nonlinear terms in the ADM [12, 13] and calculate the Lagrange multiplier as in the VIM [14, 15], respectively. To overcome these difficulties, DaftardarGejji and Jafari [16] have proposed a new technique for solving linear/nonlinear functional equations, namely, the new iterative method (NIM) or the DJ method [17–20].
We noticed that for solving the prey and predator problem, the methods mentioned above give the accurate results but in a small time interval. Once the time interval becomes wide, these methods no longer give good results or fail. To overcome this difficulty, we introduce the multistage approach to the SORlike new iterative method and show that our proposed approach works well when the other methods fails and gives accurate results compared to the orderfour RungeKutta (RK4) method.
The motivation of this paper is to introduce a relaxation parameter and propose an SORlike new iterative method which improves the DJ method. As special cases for certain values of , we find some wellknown methods like the new iterative method [16] and the revised new iterative method [21]. It has been used to solve effectively, easily, and accurately a large class of nonlinear problems with approximations. The main attractive features of the current method are being derivativefree, overcoming the difficulty in some existing techniques, simple to understand, and easy to implement, and the numerical results show that the approximations converge rapidly to accurate solutions.
The paper is organized as follows. Section 2 is devoted to present an SORlike new iterative method, and for some cases for particular values of , we find some known methods. In the case of a large interval, we present a multistage approach applied to our method. We recall a nonlinear SORlike method [22] to solve a system of nonlinear differential algebraic equations, and we give a numerical example for comparison with our proposed method. The numerical results show that our method is more accurate, converges rapidly, and is derivativefree. In Section 3, the mathematical model of the epidemic model and the problem of the prey and predator is illustrated. The two problems are solved by the SORlike iterative method, and the numerical results for the approximate solutions are discussed in comparison with some existing techniques, and finally in Section 4, the conclusion is presented with some possible perspectives.
2. SORLike New Iterative Method
A variety of problems in physics, chemistry, biology, and engineering can be formulated in terms of the system of nonlinear functional equations as follows:where is the nonlinear functional operator. Without loss of generality, one can obtain from (7) as follows:where is a given function and is the nonlinear operator. Equation (8) represents integral equations, ordinary differential equations (ODEs), partial differential equations (PDEs), and differential equations involving fractional order systems of ODE/PDE.
In order to solve system (8), let introduce a real parameter and define the SORlike new iterative method as follows.
For ,
For ,
For ,
Suppose that system (8) has a solution , then we are looking for having the series form:
By (10) and (11), we have for and
Then,
As , we get
For ,
Then,
As , we get
Thus, we have for ,
Hence, is a solution of (8).
2.1. Case 1: The New Iterative Method
We replace in (11) by 0, and we find the new iterative method, also called the DG method, which is defined in [16, 20] by
2.2. Case 2: The Revised New Iterative Method
We replace in (11) by 1, and we find the revised new iterative method, which is defined in [21] as follows.
Initial step:
First iteration:
For iteration (),
2.3. Multistage SORLike New Iterative Method
This section is devoted to introduce the basic approach of the multistage principle: The SORlike new iterative method as other semianalytical methods can provide valid approximate solutions only in the neighbourhood of initial time. Sometime, those methods have a slow convergence or fail when using the large domain. This drawback prevents its application to general problems. To overcome this difficulty, the multistage strategy is applied to the SORlike new iterative method as follows: Step 1: divide the time interval into subintervals, i.e., , , , with constant step size . Step 2: at each stage, a set of differential equations will be solved using the SORlike new iterative method with initial conditions. In order to compute approximate solutions in each stage, we need to know initial conditions at the left end point of every subinterval. Usually, we do not have these information except at the initial time . Therefore, knowledge of a process at the initial time is used to calculate approximate solutions in the first subinterval , and this solution serves as the initial condition for the next subinterval . Likewise, we can obtain approximate solutions in every subinterval, and this procedure will continue until a desired time span is reached. Step 3: formulation. Let such that , , where and . Using the initial conditions given in (7), let be the solution of system (8) in the subinterval given by (11).
The value at is used as the initial condition for the next time step. For , the initial condition for the next time step is written as
Therefore, to obtain the approximated values of the at any grid point, we define the characteristic function on by
Then, for , the solution of (8) is given by
2.4. Comparison with Nonlinear SORLike Method
In this section, we review the nonlinear SORlike method defined in [22] to solve the nonlinear algebraic equationsin or , which can simply denoted bywith and defined by solving the equationwith respect to and setting , , .
The SORlike iteration is written in the following form:where
In order to compare this method with our proposed method, we transform the nonlinear algebraic equation (28) into the form defined in (8), by
Then, the corresponding nonlinear SORlike method is defined by
Example 1. Consider the following nonlinear algebraic equation:Equation (34) is converted toIn Tables 1 and 2, we compare the number of iterations taken by both the SORlike method (33) and the SORlike new iterative method ((10) and (13)) for equations (34) and (35) as follows.
In Tables 1 and 2, we observe that a few terms of the proposed method (10) and (11) are enough to give an accurate solution approximation compared to the SORlike method (33).


Remark 1. The SORlike new iterative method is derivativefree and converges rapidly to accurate solutions. It may be observed that our method is an extension of the SORlike method to solve a general nonlinear equation of the form (8).
One remark is that the proposed method can be considered as a relaxation of the Picard method for this type of the problem but also in the resolution of a wide class of nonlinear equations.
3. Applications
The existence of the solution is not easy to find for the nonlinear equations, so we use wellknown approximation methods used as referenced approximation solutions.
Thus, for all numerical comparisons, the fourthorder RungeKutta (RK4) approximation solution will be considered as a benchmark for the considered problems.
In order to apply the SORlike new iterative method to solve the two systems of differential equations (4) and (6), we need to convert them into the corresponding system of the integral equation, and this can be done by applying the inverse of the operator , which is the integration operator to each equation in system (8) as follows.
3.1. The Epidemic Model
The epidemic model (4) is converted towhere , , and are the initial conditions.
The term approximate solution of (4) by the SORlike new iterative method is represented bywhere the components , , and are recursively computed by (11) with being the integral operator.
A few first approximation terms are calculated by the SORlike new iterative method:
Continuing in this manner, we can find the rest of the components which are calculated but not listed.
3.1.1. Numerical Results and Comparison with the ADM, HPM, NIM, and RK4
For comparison with the results obtained by Biazar [2] using the ADM, the results obtained by Rafei et al. [3], and the results obtained by ALJawary [7] using NIM, the following values, for parameters, are considered:
Three, four, five, and sixterm approximations for , , and are calculated and are plotted in Figure 1. As the plots show while the number of susceptible population increases, the population who are infective decreases in the period of the epidemic, and meanwhile the number of immune population increases. But the size of the population over the period of the epidemic is constant.
(a)
(b)
(c)
(d)
A comparison between the results obtained by the SORlike new iterative method and the NIM, ADM, and VIM in [2, 3, 7] shows that the results of fourterm approximations of the of both SORlike and the NIM and VIM is the same as the results of sixterm approximations of the ADM. In fact, the results obtained by the SORLike new iterative method are completely the same with the VIM and NIM, but it appears the results obtained by ADM in [2] are not accurate enough. Other results from comparison with orderfour RungeKutta (RK4) and the SORlike new iterative method for different values of the parameter are given in Figure 2. We notice that these plots show that our method gives better results and approaches the solution better.
(a)
(b)
(c)
3.2. The Prey and Predator Problem
The prey and predator problem is converted to
By starting with initial approximations and , the term approximate solution of (6) by the SORlike new iterative method is represented by
A few first term approximations can be calculated by the SORlike method:
Continuing in this manner, we can find the rest of the components which are calculated but not listed.
3.2.1. Numerical Results and Comparison with the ADM, HPM, NIM, and RK4
For comparison with the results obtained by Biazar [2] using the ADM, the results obtained by Rafei et al. [3], and the results obtained by ALJawary [7] using the NIM, the following values, for parameters, are considered:
A comparison between the results obtained by the SORlike method and the NIM, ADM, and VIM in cases 1–2 is completely the same for and will not be plotted (see [2, 3, 7]).
But when we vary the values of , we clearly notice that our proposed method gives better results compared to the NIM method if we compare them to the RungeKutta (RK4) method. The following plots show the relations between the number of foxes and the rabbits versus time (Figures 3 and 4).
(a)
(b)
(a)
(b)
3.2.2. Comparison with the Multistage SORLike New Iterative Method and RK4
When solving the prey and predator problem, all the methods ADM, VIM, and NIM give the accurate results but in a small time interval [0, 0.16] (Figure 4). Once the time interval becomes wide (as example [0, 1]), these methods fail as shown in Figure 3. To overcome this difficulty, we use the multistage SORlike new iterative method and show that our proposed approach works well and gives accurate results compared to the orderfour RungeKutta (RK4) method as illustrated in Figure 5.
(a)
(b)
4. Conclusion
In this paper, we propose an improvement of a new iterative method called the SORlike new iterative method to solve the nonlinear system of differential equations. The method introduces a relaxation parameter which in some cases makes it possible to find the wellknown method. Compared to the nonlinear SORlike method to solve the nonlinear algebraic system of equations, our method gives an accurate and fast approximation solution. We applied the proposed method to solve the epidemic model and the prey and predator problem successfully. When the time interval is wide, we introduce the multistage. The advantages of the proposed method are being derivativefree, simple to understand, and easy to implement and can be easily comprehended with only a basic knowledge of calculus. Compared to the Adomian decomposition method, the homotopy perturbation method, the variational iteration method, and the new iterative method, there is less computation needed. As perspectives, we will apply our method to solve the system of nonlinear stiff differential equations. For computations and plots, the MATLAB R2012 is used.
Data Availability
All algorithms and data used to support the findings of this study are included within the article and available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 W. O. Kermack, A. G. McKendrick, and G. T. Walker, “A contribution to the mathematical theory of epidemics,” Proceedings of the Royal Society, vol. 115, pp. 700–721, 1927. View at: Publisher Site  Google Scholar
 J. Biazar, “Solution of the epidemic model by Adomian decomposition method,” Applied Mathematics and Computation, vol. 173, no. 2, pp. 1101–1106, 2006. View at: Publisher Site  Google Scholar
 M. Rafei, H. Daniali, and D. D. Ganji, “Variational iteration method for solving the epidemic model and the prey and predator problem,” Applied Mathematics and Computation, vol. 186, no. 2, pp. 1701–1709, 2007. View at: Publisher Site  Google Scholar
 Y. M. Rangkuti, S. Side, S. Side, and M. S. M. Noorani, “Numerical analytic solution of SIR model of dengue fever disease in South Sulawesi using homotopy perturbation method and variational iteration method,” Journal of Mathematical and Fundamental Sciences, vol. 46, no. 1, pp. 91–105, 2014. 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
 M. Rafei, H. Daniali, D. D. Ganji, and H. Pashaei, “Solution of the prey and predator problem by homotopy perturbation method,” Applied Mathematics and Computation, vol. 188, no. 2, pp. 1419–1425, 2007. View at: Publisher Site  Google Scholar
 M. ALJawary, “A reliable iterative method for solving the epidemic model and the prey and predator problems,” International Journal of Basic and Applied Sciences, vol. 3, no. 4, pp. 441–450, 2014. View at: Publisher Site  Google Scholar
 G. F. Simmons, Differential Equations with Applications and Historical Notes, McGrawHill, New York, NY, USA, 1972.
 M. S. H. Chowdhury, I. Hashim, and S. Mawa, “Solution of prey and predator problem by numericanalytic technique,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, pp. 1008–1012, 2009. View at: Google Scholar
 J.H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004. View at: Publisher Site  Google Scholar
 J.H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004. View at: Publisher Site  Google Scholar
 G. Adomian and G. E. Adomian, “A global method for solution of complex systems,” Mathematical Modelling, vol. 5, no. 4, pp. 251–263, 1984. View at: Publisher Site  Google Scholar
 G. Adomian and R. Rach, “On the solution of Algebraic equations by the decomposition method,” Journal of Mathematical Analysis and Applications, vol. 105, no. 1, pp. 141–166, 1985. View at: Publisher Site  Google Scholar
 J.H. He, “Variational iteration method  a kind of nonlinear analytical technique: some examples,” International Journal of Nonlinear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Publisher Site  Google Scholar
 J. H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 23, pp. 115–123, 2000. View at: Publisher Site  Google Scholar
 V. DaftardarGejji and H. Jafari, “An iterative method for solving nonlinear functional equations,” Journal of Mathematical Analysis and Applications, vol. 316, no. 2, pp. 753–763, 2006. View at: Publisher Site  Google Scholar
 S. Bhalekar and V. DaftardarGejji, “New iterative method: application to partial differential equations,” Applied Mathematics and Computation, vol. 203, no. 2, pp. 778–783, 2008. View at: Publisher Site  Google Scholar
 S. Bhalekar and V. DaftardarGejji, “Solving evolution equations using a new iterative method,” Numerical Methods for Partial Differential Equations, vol. 26, pp. 906–916, 2010. View at: Google Scholar
 V. DaftardarGejji and S. Bhalekar, “Solving fractional boundary value problems with Dirichlet boundary conditions using a new iterative method,” Computers & Mathematics with Applications, vol. 59, no. 5, pp. 1801–1809, 2010. View at: Publisher Site  Google Scholar
 V. DaftardarGejji and H. Jafari, “Convergence of the new iterative method,” International Journal of Differential Equations, vol. 2011, Article ID 989065, 10 pages, 2011. View at: Google Scholar
 S. Bhalekar and V. DaftardarGejji, “Solving a system of nonlinear functional equations using revised new iterative method,” International Journal of Mathematical and Computational Sciences, vol. 6, no. 8, pp. 968–972, 2012. View at: Google Scholar
 K. Ishihara, Y. Muroya, and T. Yamamoto, “On nonlinear SORlike methods, II  convergence of the SORNewton method for mildly nonlinear equations,” Japan Journal of Industrial and Applied Mathematics, vol. 14, no. 1, pp. 99–110, 1997. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Atika Radid and Karim Rhofir. 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.