New investigation of bats-hosts-reservoir-people coronavirus model and application to 2019-nCoV system

According to the report presented by the World Health Organization, a new member of viruses, namely, coronavirus, shortly 2019-nCoV, which arised in Wuhan, China, on January 7, 2020, has been introduced to the literature. The main aim of this paper is investigating and finding the optimal values for better understanding the mathematical model of the transfer of 2019-nCoV from the reservoir to people. This model, named Bats-Hosts-Reservoir-People coronavirus (BHRPC) model, is based on bats as essential animal beings. By using a powerful numerical method we obtain simulations of its spreading under suitably chosen parameters. Whereas the obtained results show the effectiveness of the theoretical method considered for the governing system, the results also present much light on the dynamic behavior of the Bats-Hosts-Reservoir-People transmission network coronavirus model.

On December 2019, the symptom of 2019-nCoV infected patients was identified as fever, cough, breathing difficulties, and some other. Due to the long incubation period and mild symptoms, the suspicious infected people need to be observed for around 14 days. To reduce population flow and restrict the spread of the virus, corresponding virus dissemination control policies and relevant actions are being carried out at different levels. On January 23, Wuhan was locked down by strict restrictions of transport, and soon some other provinces announced to lockdown as well [3]. Chinese people were suggested to stay at home and avoid gathering, assembling, celebrations, visiting, and so on to reduce the virus dissemination. People were required to wear respirators in public areas. Some researchers pointed out some effective ways to control the spread of infectious virus, including school closure, case isolation, household quarantine, internal travel restrictions, and border control, which were proved to be helpful for the delay or reduction of virus infections [4][5][6]. Almost all the economic activities in some countries have been paused, which caused countless damages to humans' lives, development, and also a large amount of financial pressure on these countries. Globally, Australia and New Zealand firstly released regulations to ban travellers who had been to China in the past 14 days. From December 2019 up to now, 2019-nCoV virus suddenly out-broke in other countries in flood. Italy, Japan, South Korea, Spain, and other countries were in the clouds of the virus as well, and some of them started to lockdown cities, restrict transports, close schools, and so on. The whole world is in damage by the 2019-nCoV and in the combat against coronavirus. Some other researchers devoted themselves to the drug and vaccine developments, however, there is no way to effectively eliminate the virus in the human's body.
A large number of researchers have studied the 2019-nCoV from a wide range of perspectives, including diverse infectious diseases, microbiology, virology, respiratory system, biochemistry molecular biology, immunology, public environmental occupational health, genetics heredity, veterinary sciences, environmental sciences ecology, and pathology. Most of them are conducted by the USA and China, followed by Saudi Arabia, South Korea, and Germany. According to Tian [7] and others, the differences of RBD between SARS-CoV and 2019-nCoV is of significance for the cross-reactivity of neutralizing antibodies, and a SARS-CoV-specific human monoclonal antibody CR3022 could bind potently with 2019-nCoV RBD (KD of 6.3 nM). As for the origins of 2019-nCoV virus, 2019-nCoV, Benvenuto et al. [8] held that 2019nCoV could be considered a coronavirus distinct from SARS virus, probably transmitted from bats or another host, where mutations conferred upon it the ability to infect humans, and they also proposed a preliminary evolutionary and molecular epidemiological analysis of this new virus, considering high genetic similarity between 2019-nCoV and the Severe Acute Respiratory Syndrome coronavirus (SARS-CoV), and leveraging existing immunological studies of SARS-CoV, Ahmed, Quadeer, and McKay [9] devoted themselves to seeking for gaining insights for vaccine design against 2019-nCoV.
With the increasing of virus spread and the ongoing related research both at domestic and international, one question still hinders human's knowledge of the 2019-nCoV: What is the original source of such a virus and how can it transmit to human. In this paper, we intend to study a mathematical model called the Bats-Hosts-Reservoir-People coronavirus (BHRPC) model for the transfer of 2019-nCoV from the reservoir to people. By using a powerful numerical method we gain its spreading simulations under suitably chosen parameters. The obtained results show the effectiveness of the theoretical method considering the governing system and also present much light on the dynamic behavior of the Bats-Hosts-Reservoir-People transmission network coronavirus model. In this regards, more recently, some experts have investigated some important nonlinear models arising in real-world problems [10][11][12][13][14][15][16].
One of such problems has been mathematically developed by Chen et al. for simulating the phase-based transmissibility of a novel coronavirus as 2019-nCoV [17] defined by where n p , m p , b p , κ, b w , δ p , ε, c are real nonzero constants. The initial conditions for this system are shortly given by Equation (1) is used to describe the phase-based transmissibility of a novel coronavirus from source to people. In Eq. (1), u is the susceptible people, v is used to symbolize exposed people, y is the symptomatic infected people, f is asymptomatic infected people, r is removed people (recovered and died people), n p is the birth rate, m p is the death rate of people, w is the reservoir (the seafood area), 1/ω B is the incubation period of bat infection, and 1/γ B is the infectious period of bat infection [17]. Khan et al. [10] have investigated the endemic equilibria, stability, and global sensitivity of system (1). If r 0 < 1, then system (1) is locally asymptotically stable [10,18], and the outbreak will fade away [18]. When r 0 > 1, the outbreak will occur [18], and it is not stable. The used data for system (1) are for Wuhan, China, [17,19]. Moreover, mathematical analysis and applications of dengue fever outbreak and epidemiology in the sense of fractional have been investigated [20,21]. Recently, a numerical scheme based on the Newton polynomial has been applied successfully to observe important properties of the spread of COVID-19 with new fractal-fractional operators [22]. Some close relationships of COVID-19 with HIV have been presented in [23]. Many applications of fractional-or integer-order mathematical models explaining more detailed informations about the real-world problems have been presented in a detailed manner . In this paper, we investigate the numerical distributions of 2019-nCoV according to time with the help of several approaching terms of VIM.

Some important properties of VIM
In 1999, VIM, one of the most powerful numerical methods, was firstly developed by He [51][52][53][54] for numerical investigation and exceeds the difficulties of the perturbation or Adomian functions. Later, Wazwaz has applied VIM for investigating linear and nonlinear wave equations along with wave-like equations [55] and Laplace equation [56]. Moreover, many applications of VIM have been observed for various models [14][15][16][17][51][52][53][54][55][56][57][58][59][60]. We consider differential equations of the form where L and N are linear and nonlinear operators, respectively, [55], and f (x) is a source inhomogeneous term. According to basic concepts of VIM presented by He, we construct the following iteration formula for Eq. (3) [51][52][53][54]: where the parameter λ is a general Lagrange multiplier, which can be optimally identified via the variational theory, the subscript n denotes the nth-order approximation, and u n is considered as a restricted variation, which means δ u n = 0. Clearly, the main steps of VIM first require the determination of the Lagrange multiplier λ, which needs to be optimally identified. Once λ is determined, the successive approximations u n+1 , n ≥ 0, of the solution u are obtained upon using a suitably selected function u 0 , which satisfies the boundary conditions. Then the solution is given by

Application of VIM to COVID-19 system
In this subsection, by using VIM we numerically investigate the Bats-Hosts-Reservoir-People coronavirus model. According to VIM iteration structure, we can write Eq. (1) in the following form: where k = 0, 1, 2, 3, . . . . It produces the stationary condition λ = -1.
Substituting Eq. (7) into Eq. (6), we find the following iteration equation: With the help of some computational software algorithm, by considering Eq. (2) with initial values we get the first approximate components of the system for k = 0: Thus we get the first approach of u n , v n , y n , f n , r n , w n as follows: f 1 (t) = β 4 + τ 4 t, r 1 (t) = β 5 + τ 5 t, where, for simplicity, we have taken Now we obtain the second components of the variables u n , v n , y n , f n , r n , w n for k = 1: y 2 (t) = β 3 + (τ 3 + τ 13 )t + τ 14 2 t 2 , f 2 (t) = β 4 + (τ 4 + τ 15 )t + τ 16 2 t 2 , r 2 (t) = β 5 + (τ 5 + τ 17 )t + τ 18 2 t 2 , w 2 (t) = β 10 + (τ 6 + τ 19 )t + τ 20 2 Figure 1 Two-dimensional surfaces of system (10) where τ 14 = +τ 2 (-1 + δ p )ω p + τ 3 (γ p + m p ), 10 , The remaining components of the iteration formula (5) can be found in the same manner using a similar algorithm via various computational schemes. In this work, we observe the spreading rate of the Bats-Hosts-Reservoir-People coronavirus model from reservoir to people under suitably chosen values of the parameters, reported by experts in Wuhan area of China. We obtain two-dimensional simulations of the second terms u 2 , v 2 , y 2 , f 2 , r 2 , w 2 of u n , v n , y n , f n , r n , w n as in Fig. 1.

Conclusions
In this paper, we have successfully applied VIM to numerical investigation of the 2019-nCoV model. This method is based on a series solution terms of iteration. Only in the second terms of iteration, we have obtained numerical results for system (10). Under suit- Figure 2 Two-dimensional surfaces of susceptible people and exposed people of system (10) Figure 3 Two-dimensional surfaces of symptomatic infected people and asymptomatic infected people of system (10) Figure 4 Two-dimensional surfaces of removed people (recovered or died people) and reservoir (the seafood area) of system (10) ably chosen values of the parameters, reported by WHO, we have plotted the numerical results. According to Figs. 1, 2, 3, 4, we observed that VIM produces similar distributions for susceptible people, which increase exponentially, when we compare the simulated data with Fig. 2 in [17]. Moreover, we can say that each class of system (1) has also simulated estimated behaviors from them. Furthermore, the spread of the 2019-nCoV with susceptible people u 2 is faster than the others, such as v 2 , y 2 , f 2 , r 2 , w 2 . Finally, from the second terms of the proposed algorithm we can observe that susceptible people will affect more and more people from all over the world. As a future direction of this concept, application of powerful projected tools may be studied. These produce more comprehensive results on the mathematical system of 2019-nCoV.