As the ordinary differential system of equations is considerably used to explain mathematical properties, they are used to represent the real-world problems arising in the field of science such as biology, physics, engineering, finance, or sociology. It is well-known that many models in these areas can be described by using a system of ordinary differential equations. Analytically, it is generally very hard or impossible to solve these models. Moreover, as the exact solutions of these models cannot be found, we need to use the numerical methods to find approximate solutions. Many researchers used methods for solving such equations. Hojjati et al. [1] investigated the multistep method. Mastorakis [2] presented the collocation method. Shawagfed et al. [3] proposed the Adomian decomposition method. While Yüzbaşı ve Karaçayır [4] handled the Galerkin method, Yüzbaşı [5] proposed the exponential collocation method. Al-Omari et al. [6] used the Galerkin finite element method to observe numerical properties of models.
In modern times, many experts studied the epidemic models. Iwa et al. studied malaria and cholera models by using the Caputo-Fabrizio derivative of order
This paper considers the SIR model being an epidemic model. This model is used to study the spread of infectious diseases in a population. Unlike the SIR model, there is an additional factor for vaccination. In a vaccination-inclusive model, controlling the epidemic can be achieved by forecasting its future direction. Because biologists can use it to deal with natural initial conditions to design and conduct experiments in observing the spreading of diseases, with the help of these experiments and modeling, the way of controlling the spreading of the epidemic can be learned. Getting exact solutions for the problems representing such natural events is extremely difficult. Searching for convenient methods is an excellent mission for the science community.
The SIR model is composed of three members: those who are susceptible to disease, those who might be infected with the disease, and those who recovered from or become immune to it. Kermack and McKendrick earlier [13] discovered the SIR model.
Rachah and Torres improved the SIR model [14]. They expressed the susceptible-infectious-recovery (SIR) model with vaccination and the initial conditions as follows,
A person can only be in one of the three groups in any time.
Recently, some researchers have analyzed the SIR model using various methods to obtain an approximate solution. Argub and El-Ajou [16] applied the Homotopy analysis method for different parameter values for the SIR model. Biazar [17] presented the Adomian decomposition method. While Rafai et al. [18] studied the homotopy perturbation method, İbrahim and Ismail [19] used a differential transformation approach. This model was also solved using the Laplace-Adomian decomposition method [20]. Harman ve Johnston [21] solved by using the stochastic Galerkin method. This epidemic model was solved by Kousar et al. [22] by using the 4th-order Runge Kutta method. It was also solved by the Hussain et al. [23] Runge Kutta-2 ve Runge Kutta-4 methods. Kousar solved this problem via the Euler method [22].
In this paper, the solutions of the SIR model have been found using the Morgan-Voyce collocation method (M-VCM). Many researchers have recently utilized the collocation method to solve this differential equation. Taylor [24,25,26,27,28], Pell-Lucas [ 29], Morgan-Voyce [ 30,31,32,33,34], Dickson [35], ortho-exponential [36], and Bernoulli [37] matrix methods have been used to solve the linear differential, integro-differential equations, fractional differential equations [38], nonlinear differential equations, multi-pantograph [39], delay differential equations and difference and Fredholm integro-difference equation [40]. This method can significantly aid in reducing the complexity of solutions for large-scale epidemic models and can contribute to obtaining results that closely approximate the exact solutions of ordinary differential equation systems.
This study is organized as follows. Section 2 describes the method. We give the error estimation of the method in Section 3. Several applications are provided in Section 4. Finally, we introduce the results of this paper in Conclusion 5.
This section presents the M-VCM in clear. This method is used to obtain the approximate solution of the equation (1). Equation (1) is expressed in the truncated Morgan-Voyce series form as follows,
Here,
To write equation (1) in matrix form, we can first write the approximate solutions as follows,
Therefore, we can write the following equations,
The relationship between matrix
From equations (4) and (5), we can get the following equations,
Thus, we can obtain the matrices
Equation (1) can be written in matrix form as follows,
Now, to determine the unknown coefficients of
Substituting the collocation points in equation (6), we get the following system,
So, equation (1) can be written in the following matrix form utilizing the upper matrices,
Similarly, the following matrix form is obtained by substituting the collocation points in
The fundamental matrix equation from equations (7)–(10) is as follows,
Briefly, equation (11) can be written as,
Equation (12) refers to a nonlinear system of 3(
Thus, these matrix forms can be written as,
Finally, we replace the row matrices in equation (13) with any three rows of the matrix in equation (12), and we get a solution for equation (1) under initial conditions. Finally, we get the augmented matrix as,
First of all, this system is solved and the unknown coefficients
In this section, we investigate the accuracy of the presented method.
This section supports the accuracy and efficiency of the presented method. We offer two numerical examples, one comparing the errors and solutions with the other methods. In the second sample, the results for different vaccination rates are compared. The tables and figures show the comparisons at specific points within the given interval. Numerical computations are performed using MATLABR2021a software.
Comparisons of the values of
0.2 | 19.39842557 | 19.39842556 | 19.39842557 | 19.3984255735 |
0.4 | 18.79461232 | 18.79461227 | 18.79461232 | 18.7946122770 |
0.6 | 18.18993727 | 18.18993677 | 18.18993727 | 18.1899367871 |
0.8 | 17.58578366 | 17.58578115 | 17.58578366 | 17.5857811555 |
1.0 | 16.98352883 | 16.98352001 | 16.98352883 | 16.9835200246 |
Comparisons of the residual errors
0.2 | 2.2416×10−8 | 1.0080×10−9 | 2.2417×10−8 | 3.8045×10−9 |
0.4 | 6.6425×10−7 | 1.0201×10−9 | 6.6427×10−7 | 4.4120×10−9 |
0.6 | 4.6388×10−6 | 1.0363×10−9 | 4.6389×10−6 | 5.2320×10−9 |
0.8 | 1.7832×10−5 | 1.0566×10−9 | 1.7833×10−5 | 6.3120×10−9 |
1.0 | 4.9171×10−5 | 1.9918×10−7 | 4.9172×10−5 | 1.9256×10−7 |
Comparisons of the values of
0.2 | 15.54049369 | 15.54049370 | 15.54049369 | 15.5404936937 |
0.4 | 16.08106363 | 16.08106367 | 16.08106363 | 16.081063672 |
0.6 | 16.62033527 | 16.62033570 | 16.62033527 | 16.6203357015 |
0.8 | 17.15693346 | 17.15693567 | 17.15693345 | 17.156935671 |
1.0 | 17.68949465 | 17.68950236 | 17.68949464 | 17.6895023739 |
Comparisons of the residual errors
0.2 | 2.0373×10−8 | 1.0080×10−9 | 2.0375×10−8 | 8.9114×10−9 |
0.4 | 5.9888 ×10−7 | 1.0201×10−9 | 5.9892×10−7 | 9.0230×10−9 |
0.6 | 4.1424×10−6 | 1.0363×10−9 | 4.1426×10−6 | 9.0493×10−9 |
0.8 | 1.5741×10−5 | 1.0566×10−9 | 1.5741×10−5 | 8.9952×10−9 |
1.0 | 4.2787×10−5 | 1.9918×10−7 | 4.2789×10−5 | 4.4706×10−7 |
Comparisons of the values of
0.2 | 10.06108073 | 10.06108073 | 10.06108073 | 10.0610807329 |
0.4 | 10.12432405 | 10.12432405 | 10.12432405 | 10.1243240513 |
0.6 | 10.18972750 | 10.18972751 | 10.18972750 | 10.1897275117 |
0.8 | 10.25728304 | 10.25728317 | 10.25728304 | 10.2572831742 |
1.0 | 10.32697698 | 10.32697760 | 10.32697698 | 10.3269776025 |
Comparisons of the residual errors
0.2 | 1.5573×10−9 | 1.7802×10−12 | 1.5573×10−9 | 1.9663×10−9 |
0.4 | 7.7657×10−9 | 3.5556×10−12 | 7.7656×10−9 | 2.7950×10−9 |
0.6 | 2.0477×10−7 | 5.3304×10−12 | 2.0477×10−7 | 3.7913×10−9 |
0.8 | 1.1701×10−6 | 7.1090×10−12 | 1.1701×10−6 | 4.9600×10−9 |
1.0 | 4.1334×10−6 | 4.5486 ×10−7 | 4.1334×10−6 | 2.4930×10−7 |
The residual errors for equation (1),
0.1 | 2.0017×10−7 | 2.0805×10−8 | 1.0544×10−8 |
0.2 | 2.4419×10−7 | 2.5143×10−8 | 1.1137×10−8 |
0.3 | 2.9299×10−7 | 2.9958×10−8 | 1.1797×10−8 |
0.4 | 3.4678×10−7 | 3.5273×10−8 | 1.2529×10−8 |
0.5 | 4.0577×10−7 | 4.1114×10−8 | 1.3337×10−8 |
0.6 | 4.7020×10−7 | 4.7505×10−8 | 1.4226×10−8 |
0.7 | 5.4027×10−7 | 5.4471×10−8 | 1.5200×10−8 |
0.8 | 6.1622×10−7 | 6.2036×10−8 | 1.6262×10−8 |
0.9 | 6.9824×10−7 | 7.0221×10−8 | 1.7416×10−8 |
1 | 7.8657×10−7 | 7.9046×10−8 | 1.8661×10−8 |
The residual errors for equation (1),
0.1 | 5.1591×10−8 | 9.3905×10−8 | 1.7736×10−9 |
0.2 | 5.5453×10−8 | 1.1294×10−7 | 2.5089×10−9 |
0.3 | 5.9708×10−8 | 1.3399×10−7 | 3.4104×10−9 |
0.4 | 6.4395×10−8 | 1.5712×10−7 | 4.4980×10−9 |
0.5 | 6.9554×10−8 | 1.8243×10−7 | 5.7929×10−9 |
0.6 | 7.5227×10−8 | 2.1000×10−7 | 7.3172×10−9 |
0.7 | 8.1458×10−8 | 2.3993×10−7 | 9.0941×10−9 |
0.8 | 8.8295×10−8 | 2.7229×10−7 | 1.1148×10−8 |
0.9 | 9.5794×10−8 | 3.0717×10−7 | 1.3503×10−8 |
1 | 1.0402×10−7 | 3.4464×10−7 | 1.6184×10−8 |
The residual errors for equation (1),
0.1 | 3.1084×10−7 | 1.0500×10−9 | 2.4970×10−7 |
0.2 | 4.0778×10−7 | 2.0300×10−9 | 2.9760×10−7 |
0.3 | 5.2230×10−7 | 3.2950×10−9 | 3.5110×10−7 |
0.4 | 6.5640×10−7 | 4.0000×10−9 | 4.1066×10−7 |
0.5 | 8.1100×10−7 | 6.8300×10−9 | 4.7650×10−7 |
0.6 | 9.9065×10−7 | 9.2000×10−9 | 5.4925×10−7 |
0.7 | 1.1948×10−7 | 1.2043×10−8 | 6.2904×10−7 |
0.8 | 1.4266×10−7 | 1.5405×10−8 | 7.1630×10−7 |
0.9 | 1.6880×10−6 | 1.9000×10−8 | 8.1140×10−7 |
1 | 1.9819×10−6 | 2.4010×10−8 | 9.1494×10−7 |
The residual errors for equation (1),
0.1 | 2.2100×10−8 | 3.7200×10−8 | 6.9600×10−10 |
0.2 | 2.7380×10−8 | 4.9700×10−8 | 1.7400×10−9 |
0.3 | 3.3490×10−8 | 6.4701×10−8 | 3.0619×10−9 |
0.4 | 4.0500×10−8 | 8.2200×10−8 | 4.6670×10−9 |
0.5 | 4.8700×10−8 | 1.0270×10−7 | 6.6010×10−9 |
0.6 | 5.7000×10−8 | 1.2640×10−7 | 8.9000×10−9 |
0.7 | 6.8390×10−8 | 1.5347×10−7 | 1.1605×10−8 |
0.8 | 8.0117×10−8 | 1.8420×10−7 | 1.4750×10−8 |
0.9 | 9.3220×10−8 | 2.1914×10−7 | 1.8380×10−8 |
1 | 1.0780×10−7 | 2.5836×10−7 | 2.2531×10−8 |
The residual errors for equation (1),
0.1 | 8.5647×10−7 | 5.4000×10−8 | 2.0762×10−7 |
0.2 | 9.5890×10−7 | 5.3700×10−8 | 2.4940×10−7 |
0.3 | 1.0660×10−6 | 5.2000×10−8 | 2.9568×10−7 |
0.4 | 1.1792×10−6 | 4.9930×10−8 | 3.4662×10−7 |
0.5 | 1.2972×10−6 | 4.6100×10−8 | 4.0240×10−7 |
0.6 | 1.4202×10−6 | 4.0964×10−8 | 4.6333×10−7 |
0.7 | 1.5484×10−6 | 3.4200×10−8 | 5.2955×10−7 |
0.8 | 1.6819×10−6 | 2.5800×10−8 | 6.0120×10−7 |
0.9 | 1.8205×10−6 | 1.5672×10−8 | 6.7870×10−7 |
1 | 1.9641×10−6 | 3.7710×10−9 | 7.6213×10−7 |
In this paper, we used the M-VCM to find the approximate solution of the governing model. Several applications were presented to show the efficiency and accuracy of the proposed method. To display the accuracy of the approximate solutions, these solutions are once again substituted into equation (1) using MATLAB software. Thus, this work offers a different measure for the reliability of the approximate solutions. The obtained approximate solutions and errors have been compared by the HPM [18] and the LADM [20]. These comparisons have shown that the presented method to find approximate solutions to the SIR model from epidemic models is more efficient and valuable. Tables 2, 4, and 6 illustrated that the numerical solutions of HPM [18] and LADM [20] are nearly the same. As can be seen from the figures and tables, minimal errors were obtained in this study, indicating that the results are very close to the exact solutions. In particular, in the second example, we can accurately determine the vaccination rate needed to control the epidemic. In conclusion, all these findings show that our method is efficient and feasible for solving the nonlinear ODE systems. The most significant advantage of the proposed method is that all calculations can be easily and quickly applied with the help of MATLAB code. In addition, the fact that the SIR model is a real-life problem increases the importance of our study. One of the most significant advantages of the method is changing the system equation (1) into an easily solved nonlinear algebraic equation system. As a further study, this method can be applied to other epidemic models for longer time periods.
The authors declare no conflict of interest.
Not applicable.
Ö.İ.-Methodology, Supervision, Resources, and Methodology, Validation, Conceptualization and Formal analysis; G.Ş.-Data Curation, Writing-Original Draft, Investigation and Visualization. All authors read and approved the final submitted version of this manuscript.
Not applicable.
All data that support the findings of this study are included within the article.
The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.