SOR-Like New Iterative Method for Solving the Epidemic Model and the Prey and Predator Problem

Our aim in this paper is to propose an SOR-like new iterative method by introducing a relaxation parameter ω to improve the new iterative method proposed by Daftardar-Gejji and Jafari (NIM) [J. Math. Anal. Appl. 316 (2006) 753–763] in order to solve two problems.)e 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.)e proposedmethod is not limited to these two problems but can be applicable to a wide range of systems of nonlinear functional problem. )e results, for different values of ω, show that we found some knownmethods and ourmethod compared tomethods using the calculation of special polynomials and derivatives like the Adomian decompositionmethod (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 ω.


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 S, infected I, and recovered R. e 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][3][4]. At time t, suppose the population consists of S(t) susceptible population: those so far uninfected and therefore liable to infection I(t) infective population: those who have the disease and are still at large R(t) 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. en, in time δt, δS of the susceptible population becomes infective, where where β is a positive constant. If c > 0 is the rate at which current infective population becomes isolated, then e number of new isolated population, δR, is given by Now let δt ⟶ 0. en, the following system determines the progress of the disease: dS dt � −βS(t)I(t), with the following initial conditions: S(0) � n 1 , I(0) � n 2 , and R(0) � n 3 .
ere 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 x(t) and y(t) are, respectively, the populations of the rabbits and the foxes at time t. 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 AL-Jawary [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, Daftardar-Gejji 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][18][19][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 SOR-like new iterative method and show that our proposed approach works well when the other methods fails and gives accurate results compared to the order-four Runge-Kutta (RK4) method. e motivation of this paper is to introduce a relaxation parameter ω ∈ (0, 2) and propose an SOR-like 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. e main attractive features of the current method are being derivative-free, 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. e paper is organized as follows. Section 2 is devoted to present an SOR-like 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 SOR-like method [22] to solve a system of nonlinear differential algebraic equations, and we give a numerical example for comparison with our proposed method. e numerical results show that our method is more accurate, converges rapidly, and is derivative-free. In Section 3, the mathematical model of the epidemic model and the problem of the prey and predator is illustrated. e two problems are solved by the SOR-like 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.

SOR-Like 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: . . , u n � 0, ⋮ du n dt + F n u 1 , . . . , u i , . . . , u n � 0, u i (0), i � 1, 2, . . . , n initial conditions, where F is the nonlinear functional operator. Without loss of generality, one can obtain u i from (7) as follows: 2 Discrete Dynamics in Nature and Society where f is a given function and N 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 ω ∈ (0, 2) and define the SOR-like new iterative method as follows.

Case 1: e 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 4 Discrete Dynamics in Nature and Society 2.2. Case 2: e 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 k th iteration (k � 2, 3, · · ·), is section is devoted to introduce the basic approach of the multistage principle: e SOR-like 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. is drawback prevents its application to general problems. To overcome this difficulty, the multistage strategy is applied to the SOR-like new iterative method as follows: Step Step 2: at each stage, a set of differential equations will be solved using the SOR-like 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 t 0 � 0. erefore, knowledge of a process at the initial time t 0 � 0 is used to calculate approximate solutions in the first subinterval [t 0 , t 1 ], and this solution serves as the initial condition for the next subinterval [t 1 , t 2 ]. 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 I � ∪ M−1 j�0 I j such that Using the initial conditions given in (7), let u ⋆ 0 � (u ⋆ 1,0 , . . . , u ⋆ n,0 ) be the solution of system (8) in the subinterval I 0 given by (11). e value at t � t 1 is used as the initial condition for the next time step. For j � 1, . . . , M − 1, the initial condition for the next time step is written as Discrete Dynamics in Nature and Society erefore, to obtain the approximated values of the u ⋆ at any grid point, we define the characteristic function χ I j on I j by en, for t ∈ [t 0 , T], the solution u ⋆ of (8) is given by

Comparison with Nonlinear SOR-Like Method.
In this section, we review the nonlinear SOR-like method defined in [22] to solve the nonlinear algebraic equations in R n or C n , which can simply denoted by with F � (F 1 , . . . , F n ) t and z � (z 1 , . . . , z n ) defined by solving the equation with respect to z i and setting z k+1 e SOR-like 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 (32) en, the corresponding nonlinear SOR-like method is defined by Example 1. Consider the following nonlinear algebraic equation: 15x + y 2 − 4z − 13 � 0 Equation (34) is converted to In Tables 1 and 2, we compare the number of iterations taken by both the SOR-like method (33) and the SOR-like 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 SOR-like method (33).

Remark 1.
e SOR-like new iterative method is derivativefree and converges rapidly to accurate solutions. It may be observed that our method is an extension of the SOR-like 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.

Applications
e existence of the solution is not easy to find for the nonlinear equations, so we use well-known approximation methods used as referenced approximation solutions.
us, for all numerical comparisons, the fourth-order Runge-Kutta (RK4) approximation solution will be considered as a benchmark for the considered problems. 6 Discrete Dynamics in Nature and Society In order to apply the SOR-like 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 d(.)/dt, which is the integration operator t 0 (.)dt to each equation in system (8) as follows.

(39)
Continuing in this manner, we can find the rest of the components which are calculated but not listed.

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 AL-Jawary [7] using NIM, the following values, for parameters, are considered: ree-, four-, five-, and six-term approximations for S(t), I(t), and R(t) 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 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 four-term approximations of the of both SOR-like and the NIM and VIM is the same as the results of six-term approximations of the ADM. In fact, the results obtained by the SOR-Like 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 order-four Runge-Kutta (RK4) and the SOR-like 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.

e Prey and Predator Problem.
e prey and predator problem is converted to By starting with initial approximations x 0 � x 0 and y 0 � y 0 , the k-term approximate solution of (6) by the SORlike new iterative method is represented by A few first term approximations can be calculated by the SOR-like method: (43) Continuing in this manner, we can find the rest of the components which are calculated but not listed.

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 AL-Jawary [7] using the NIM, the following values, for parameters, are considered: 10 Discrete Dynamics in Nature and Society A comparison between the results obtained by the SOR-like method and the NIM, ADM, and VIM in cases 1-2 is completely the same for ω � 0 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 Runge-Kutta (RK4) method. e following plots show the relations between the number of foxes and the rabbits versus time (Figures 3 and 4).

Comparison with the Multistage SOR-Like 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 SOR-like new iterative method and show that our proposed approach works well and gives accurate results compared to the order-four Runge-Kutta (RK4) method as illustrated in Figure 5.

Conclusion
In this paper, we propose an improvement of a new iterative method called the SOR-like new iterative method to solve the nonlinear system of differential equations. e method introduces a relaxation parameter which in some cases makes it possible to find the well-known method. Compared to the nonlinear SOR-like 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. e advantages of the proposed method are being derivative-free, 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. x(t) RK4 y(t) RK4 x(t) SOR y(t) SOR (b) Figure 5: Comparison of multistage three term SOR-like and RK4 for x(t), y(t), and ω � 1.6.

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
e authors declare that they have no conflicts of interest.