A simple mathematical model for Coronavirus (COVID-19) ====================================================== * Said Melliani * Abdelati El Allaoui * Lalla Saadia Chadli ## Abstract A novel coronavirus (COVID-19) was identified in Wuhan, China in the end of 2019, it causing an outbreak of viral pneumonia. It caused to the death rate of 4.63% among 571, 678 confirmed cases around the world to the March 28th, 2020. In this brief currentstudy, we will present a simple mathematical model where we show how the probability of successfully getting infected when coming into contact with an infected individual and the per-capita contact rate affect the healthy and infected population with time. The proposed model is used to offer predictions about the behavior of COVID-19 for a shorter period of time. Keywords * COVID-19 * Epidemic * Mathematical modelling ## 1 Introduction Early January 2020, a novel coronavirus (COVID-19) was identified as the infectious agent causing an outbreak of viral pneumonia in Wuhan, China, despite, the first cases had their symptom onset in December 2019. It has caused at least 571, 678 confirmed cases with 26, 494 deaths and 158, 700 recovered in the world, from the start of this pandemic until March 28th, 2020, according to World Health Organization [2]. The ongoing outbreak of coronavirus disease is affecting 199 countries. The most infected countries distributed in descending order of infected cases are: USA, Italy, Spain, China, Germany, Iran and France and other less infected countries. The seriousness of this epidemic is evident in the speed of its spread with a proportion of 4.63% of deaths among the infected worldwide, which is a low proportion compared to those of other epedemics. Since the beginning of COVID-19, Several researchers attacked the study of this pandemic, each of his field [5, 13–15]. Also, Mathematicians contributed their part to these studies. In [10], Qianying Lin et al studied the transmission of COVID-19 in Wuhan, China with more details. In [11], Chen et al. treated a mathematical model for simulating the phase-based transmissibility of a novel coronavirus, the goal of this study was to calculate the transmissibility of the virus by using a developed mathematical model. Other authors have studied this new disease mathematically (see [5, 11, 15]). The purpose of this paper is to consider a mathematical model to study predictions of the outbreak of coronavirus in the most infected countries in the world such as: USA, Italy, Spain, Germany, Iran and France (China is excluded from these studies because several works have been devoted only to China [4–6,10,12–15]). And, how to reduce the seriousness of Coronavirus which is the speed of its outbreak. For that, based on [1], we consider the following system of non-linear differential equations: ![Formula][1] Where the meaning of the parameters is given as follows: *S*(*t*) the susceptible at time *t, I*(*t*) the infectious at time *t, R*(*t*) the removed at time *t, C*(*t*) the exposed at time *t*, *N* the population studied (total population), *α* is an average of the number of infected people by contact with one symptomatic individual who has been infected with Corona virus, *γ* is the per-capita infectious rate, *σ* is the per-capita death rate. We shall assume that the population *N* is constant, and *N* = *S*(*t*) + *E*(*t*) + *I*(*t*) + *R*(*t*). we are interested in predictions of the outbreak of the novel coronavirus disease (COVID-19) for the next 15 weeks from March 29th, 2020. At the same time, we are interested in how the spread of this infectious disease. We have followed the Corona disease outbreak news from the start of this pandemic, January 10th to March 28th, 2020. The table 1 shows the numbers of confirmed cases of Corona virus for the most infected countries in the world: View this table: [Table 1:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/T1) Table 1: Numbers of confirmed cases of COVID-19 outbreak as they appear in official data obtained from WHO from January 10th to March 28th, 2020 In this paper, we will consider the following parameters: Coronavirus is twice as contagious as the flu. Research signals a person with the flu infects an average of 1.28 other people, CNN Chief Medical Correspondent Dr. Sanjay Gupta said. But with coronavirus, “it’s likely between two and three” other people”. For these reasons and based on the paper [10], we suggest that transmission rate *α* varies between 0.5944 and 1.68. Coronavirus can be transmitted by people who are infected by this pandemic. The chance of getting infected for an infectious individual is very higher. which means that the factor *β* satisfies 0 *≤ β ≤* 1. In the following figures we will show the prediction of our model given by the proposed system of non linear differential equations in a shorter interval of time [0 15], where 0 is the day the predictions start, it’s March 29th, 2020, and 15 designates the 15th week, i.e. from March 29th to July 11th, 2020. Figures 1.1, 1.2, 1.3, 1.4, 1.5 and 1.6 show cumulative predictions of infected in USA, Italy, Spain, Germany, France and Iran respectively from March 29th to July 11th, 2020. Where we assumed that: *α* = 0.6 and *β* = 0.037, with the other parameters given in the table 2. To analyze these figures, for example from figure 1.1, we can notice that the growth of the infected is faster than the other countries. We can see that the predictions of the number of infected with COVID-19 towards the end of the 15th week may reach: 1.944 *×* 106 cases in USA, 8.1168 *×* 105 in Italy, 8.1144 *×* 105 in Spain, 8.825 *×* 10 in Germany, 4.4216 *×* 10 in France and 3.979 *×* 105 in Iran. View this table: [Table 2:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/T2) Table 2: Summary table of the parameters used in model (1.1) ![Figure 1.1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F1.medium.gif) [Figure 1.1:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F1) Figure 1.1: Model prediction for USA ![Figure 1.2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F2.medium.gif) [Figure 1.2:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F2) Figure 1.2: Model prediction for Italy ![Figure 1.3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F3.medium.gif) [Figure 1.3:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F3) Figure 1.3: Model prediction for Spain ![Figure 1.4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F4.medium.gif) [Figure 1.4:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F4) Figure 1.4: Model prediction for Germany ![Figure 1.5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F5.medium.gif) [Figure 1.5:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F5) Figure 1.5: Model prediction for France ![Figure 1.6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F6.medium.gif) [Figure 1.6:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F6) Figure 1.6: Model prediction for Iran In the following figures, we will present the number of infected for USA and Italy most infected in the world for different values of *α*: From figure 1.7, we can see that, the higher the coefficient, the higher the number of infected, likewise for the case of Italy in figure 1.8. So to have a minimum possible of infected, we must have a minimum value of *α*. Since, the parameter is defined by the product of the probability of successfully getting infected when coming into contact with an infected individual and the per-capita contact rate, for that, it is necessary to decrease these two parameters which depend both on “the contact”. That is why reducing contact remains among the reasons for reducing infections. ![Figure 1.7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F7.medium.gif) [Figure 1.7:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F7) Figure 1.7: Model prediction for USA with different values of *α* ![Figure 1.8:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/04/28/2020.04.23.20076919/F8.medium.gif) [Figure 1.8:](http://medrxiv.org/content/early/2020/04/28/2020.04.23.20076919/F8) Figure 1.8: Model prediction for Italy with different values of *α* To evaluate the equilibrium points, we solve the following equations: *S*′(*t*) = *I*′(*t*) = *R*′(*t*) = *C*′(*t*) = 0. That is, ![Formula][2] Which is equivalent to *I*(*t*) = *C*(*t*) = 0. The Jacobian is given by ![Formula][3] The Jacobian Calculated at the equilibrium point, i.e. we take into account *I*(*t*) = *C*(*t*) = 0 will be given by ![Formula][4] Whose eigenvalues are given by ![Formula][5] and ![Formula][6] since, the parameters of the studied model verify ![Formula][7] For the case of USA, we get ![Formula][8] and ![Formula][9] With small calculations, we can easily verify that ![Formula][10] We have *β ∈* [0 1], so, a sufficient condition for *λ*3 to be negative is *S*(0) *≥* 5, 453, 335, while in this case *λ*4 *≥* 0. ## Data Availability Te data used to support the findings of this study are available from the corresponding author upon request [https://www.researchgate.net/profile/Abdelati\_El\_Allaoui2](https://www.researchgate.net/profile/Abdelati_El_Allaoui2) ## Data Availability Te data used to support the findings of this study are available from the corresponding author upon request. ## Competing interests The authors declare that they have no competing interests. ## Author’s contributions All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript. ## Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. ## Acknowledgements The authors express their sincere thanks to the anonymous referees for numerous helpful and constructive suggestions which have improved the manuscript. * Received April 23, 2020. * Revision received April 23, 2020. * Accepted April 28, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].Anderson, R, May, R: Infectious Diseases of Humans: Dynamics and Control. Oxford Science Publications, Oxford (1992). 2. [2].World Health Organization, Coronavirus disease 2019 (COVID-19) Situation Report-68. 3. [3].Bogoch II, Watts A, Thomas-Bachli A, Huber C, Kraemer MUG, Khan K. Pneumonia of unknown etiology in Wuhan, China: potential for international spread via commercial air travel. J Travel Med 2020;, doi:[http://dx.doi.org/10.1056/NEJMoa2002032](http://dx.doi.org/10.1056/NEJMoa2002032). 4. [4].Bogoch II, A. Watts, A. Thomas-Bachli, C. Huber, M.U.G. Kraemer, K. Khan, Pneumonia of unknown etiology in Wuhan, China: potential for international spread via commercial air travel, J Travel Med (2020), doi:10.1056/NEJMoa2002032. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2002032&link_type=DOI) 5. [5].Chen T, Rui J, Wang Q, Zhao Z, Cui J-A, Yin L. A mathematical model for simulating the transmission of Wuhan novel Coronavirus. bioRxiv. 2020:2020.2001.2019.911669. accessed 13 Feb 2020. 6. [6]. W. Guan, Z. Ni, Y. Hu, W. Liang, C. Ou, J. He, et al. Clinical characteristics of 2019 novel coronavirus infection in China NEJM (2020), doi:10.1056/NEJMoa2002032. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2002032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32109013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F28%2F2020.04.23.20076919.atom) 7. [7]. N. Imai, A. Cori, I. Dorigatti, M. Baguelin, C.A. Donnelly, S. Riley, et al. Report 3: transmissibility of 2019-nCoV British Medical Journal Publishing Group (2020). 8. [8].Li Q., Guan X., Wu P., Wang X., Zhou L., Tong Y., et al. Early transmission dynamics in Wuhan, China, of novel Coronavirus–infected pneumonia. N Engl J Med 2020. 9. [9].Riou J, Althaus CL. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020. Eurosurveillance 2020;25(4). 10. [10]. Qianying Lin, Shi Zhao, Daozhou Gao, Yijun Lou, Shu Yang, Salihu S. Musa, Maggie H. Wang, Yongli Cai, Weiming Wang, Lin Yang, Daihai He, A conceptual model for the coronavirus disease 2019 (COVID - 19) outbreak in Wuhan, China with individual reaction and governmental action, International Journal of Infectious Diseases, Volume 93, 2020, Pages 211–216, [https://doi.org/10.1016/j.ijid.2020.02.058](https://doi.org/10.1016/j.ijid.2020.02.058). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.058&link_type=DOI) 11. [11]. Tian-Mu Chen, Jia Ruil, Qiu-Peng Wang, Ze- Yu Zhao, Jing-An Cui and Ling Yin, A mathematical model for simulating the phase-based transmissibility of a novel coronavirus, Infectious Diseases of Poverty (2020) 9:24. 12. [12].Wu Z., McGoogan J. M. Characteristics of and important lessons from the Coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72,314 cases from the Chinese center for disease control and prevention. JAMA 2020;, doi: [http://dx.doi.org/10.1001/jama.2020.2648](http://dx.doi.org/10.1001/jama.2020.2648). 13. [13].Wu J. T., Leung K, Leung G. M. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV out- break originating in Wuhan, China: a modelling study. Lancet 2020. 14. [14].Wu P., Hao X., Lau EHY, Wong J. Y., Leung KSM, Wu J. T., et al. Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22 January 2020. Euro- surveillance 2020;25. 15. [15].Zhao S, Lin Q, Ran J, Musa SS, Yang G, Wang W, et al. Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak. Int J Infect Dis. 2020. [https://doi.org/10.1016/j.ijid.2020.01.050](https://doi.org/10.1016/j.ijid.2020.01.050). [1]: /embed/graphic-1.gif [2]: /embed/graphic-12.gif [3]: /embed/graphic-13.gif [4]: /embed/graphic-14.gif [5]: /embed/graphic-15.gif [6]: /embed/graphic-16.gif [7]: /embed/graphic-17.gif [8]: /embed/graphic-18.gif [9]: /embed/graphic-19.gif [10]: /embed/graphic-20.gif