On the Stochastic Modeling of COVID-19 under the Environmental White Noise

Faculty of Informatics and Computing, Universiti Sultan Zainal Abidin (UniSZA), Besut Campus, Terengganu, Malaysia Department of Mathematics, Shaheed Benazir Bhutto University, Sheringal, Dir Upper PO, 18000 Khyber Pakhtunkhwa, Pakistan Department of Biotechnology, Shaheed Benazir Bhutto University, Sheringal, Dir Upper PO, 18000 Khyber Pakhtunkhwa, Pakistan Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz, Iran Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan Gofa Camp, Near Gofa Industrial College and German Adebabay, Nifas Silk-Lafto, 26649 Addis Ababa, Ethiopia Institute of Mathematical Sciences, Faculty of Science, University of Malaya, Kuala Lumpur 50603, Malaysia


Introduction
A new virus known as "corona" was claimed to be wreaking havoc on the Chinese city of Wuhan in December 2019. This virus and its subsequent pandemic initially struck Wuhan and then spread to practically the entire world [1,2]. Hundreds of thousands of people died all across the world. It is difficult to come to a single conclusion on the virus's genesis. It could be the result of the migration patterns of people or a seafood market exchange. It is also possible that it is related to human-to-human relationships [3,4]. So far, the virus has wreaked havoc on practically everything on the planet. Every aspect of human life has been badly impacted, including social life, health, the economy, and education. In battling the fatal pandemic, health researchers, government politicians, and health care authorities are perplexed. They all have a different perspective on the problem and are working hard to reduce the number of people who die as a result of the outbreak. For detail, one can see the useful work in [5,6].
The coronavirus pandemic causes moderate respiratory symptoms in those who are affected. The symptoms of this illness include a fever, a dry cough, a throat infection, and exhaustion. Nasal infection, pains, and a sore throat are some of the other symptoms that people may experience.
In order to describe the pandemic of infectious illnesses, mathematical modeling is an important technique. The goal of mathematical modeling is to use mathematical language to depict many types of real-world situations. The COVID-19 was presented with the help of mathematical models in several recently published articles including a paper regarding coronavirus with two categories of asymptomatic and symptomatic individuals [7], a paper regarding the effect of quarantine on the spread of this disease [8], a paper in relation to this pandemic with deterministic structure [9], and a paper on the COVID-19 model with Haar wavelet collocation method [10]. Even, some researchers modeled this disease via different operators such as nonsingular derivatives [11] or Caputo singular derivatives [12].
A neural network is a type of nonparametric machine learning technique that is commonly used for prediction. New improved mathematical methods for predicting and analyzing pandemic casualties have been established by well-known researchers including Tutsoy et al. in [13,14]. The nature of nonlinear dynamics is determined by the countries' health care infrastructure, as well as the qualities of enforced pharmacological or nonpharmacological policies. In the context of nonlinear dynamics, the neural network technique clearly improves the batch type leastsquares technique. In the case of linear dynamics, batch type least squares typically outperforms. Only when the casualties have less diversity, such as in intensive care and death casualties, the neural network produces more accurate findings [15]. In future, we will analyze our biological model with the help of batch type least squares and neural network approach which will provide us more real results only when the causalities have less variation.
In probability theory and associated fields, a stochastic or random process is a mathematical object, usually characterized as a group of random variables. In the 1930s, the first mathematical definition of a stochastic process as a family of random variables indexed by means of the change of the real line was given by Aleksandr Khinchin. Khinchin, Andrey Kolmogorov, Joseph Doob, and William Feller furthered the probability concept and stochastic modeling.
A stochastic model means to recognize the random nature of the input components. By contrast, in a stochastic model, the output is random in nature. The output is only a snapshot or an estimate of the characteristic of the model for a given set of inputs. A stochastic model will be run many times with the same output and gives distributions of the relevant results for a distribution of scenarios. The results of stochastic models require Monte Carlo simulations, although some stochastic models can have analytical solutions. Correlation can be important in stochastic models as they indicate the behavior of one variable is associated with that of the other. Stochastic models are more complex and difficult to interpret, so stochastic models require more expertise, expense, and computer power. A stochastic model has the capacity to handle uncertainties in the inputs applied. Stochastic models possess some inherent randomness-the same set of parametric values and initial conditions will lead to an ensemble of different outputs. In mathematics, the theory of stochastic processes is a signifi-cant aspect of probability theory and a subject of research for both theory and applications [16,17]. The word stochastic is used to portray different terms and objects in mathematics. For example, this includes a stochastic matrix, which describes a stochastic process known as a Markov process [18], and stochastic calculus [19,20], which includes differential equations and integrals depending on stochastic processes such as the Wiener process, also called the Brownian motion (white noise) process [21][22][23].
In bioinformatics (in particular, in biochemical systems), there are two types of important white noises. These are intrinsic which could be seen as an inherent stochastic biochemical process and extrinsic white noise as cell-to-cell variability and external noise as environmental fluctuations, which in both cases can be viewed mathematically as white noise [24] Also, stochastic modeling can be studied in various fields of science and engineering. In the financial markets, stochastic models are used to express stocks, relative stock price, commodities, and interest rates. In biological systems, adding stochastic "noise" to internal links for balance and other auditory interactions has been shown to assist boost. Manufacturing processes are thought to be random. For both continuous and batch manufacturing processes, this assumption is essentially valid. A process control chart illustrates a process control parameter over time and is used to record the monitoring and testing of the existing process. In the same way that events are interactions of systems, stochastic social science theory emphasizes unconscious processes. Because of the large number of variables involved, every event creates its own set of possibilities, making it unpredictable. Stochastic modeling of the dynamical problems is very much popular in the community of scientists. In their work, they have different aspects of their new models including the existence, extinction, and numerical simulations. For a detailed study of the stochastic modeling and their theoretical plus numerical analysis, we recommend the readers to read [25][26][27] and references therein.
In this article, we consider that the fluctuation of the β and γ is random and β = β + σ 1 dB 1 and γ = γ + σ 2 dB 2 . With these considerations, we have the following new stochastic model for the COVID-19 spread.
Here, the symbol S * represents the class of susceptible compartment, I * represents the class of infected compartment, and R * represents the class of recovered compartment from COVID-19. Also, α 0 is the isolation rate, y is the natural death rate, γ is the recovery rate, β is the protection rate, b is the birth rate, k is the constant rate, and d 0 represents death due to corona. We have the following deterministic model [28]. 2 Journal of Function Spaces In this article, we will present an analytical and computational study of the stochastic epidemic model (1) and will compare our results with the deterministic model (2) numerically for different parametric values for the prediction of future. Thus, this study can help the community in the awareness about the spread of the disease.
Partition of the paper. In this article, a stochastic COVID-19 model based on the environmental white noise is developed. In the second section, the essential requirement for the existence criteria, extinction, and persistence with the help of Lyapunov function is analyzed. A numerical approach is developed for computational and graphical analysis of the mentioned model which is included for the numerical interpretation of the stochastic model with the help of available parametric and initial values for the model in the third section.

Existence of Solution
For the existence of solution of our problem (1), we consider a general 3-dimensional stochastic ODE.
with Uð0Þ = U 0 ∈ ℝ 3 + and define an operator L by By operating L on V ∈ C 2,1 ðℝ d × ½t 0 ,+∞ i ℝÞ, then in In the light of the Ito's formula, if UðtÞ ∈ ℝ 3 , we get To prove the global positivity of solution for (1), we use Lyapunov analysis method. For this, consider Ω given by The following lemma gives the proof for the existence of solution of model (1).

Journal of Function Spaces
In the light of the Ito's formula, we get where Therefore, we have Integrating (16) from 0 to τ∧℘ k and taking expectation, we can get Set Ω k = f℘ k ≤ tg for k ≥ k 1 where P ðΩ k Þ ≥ ζ; that is, ∀ω ∈ Ω k , and at least one of ðS * ð℘ k , ωÞ, I * ð℘ k , ωÞ, R * ð℘ k , ωÞÞ exists so that equals to k or 1/k. Consequently, in which a∧b stands for the minimum of a and b. In view of (17) and (18), we have in which I * Ω k stands for the indicator function of Ω k . By k ⟶ ∞, we arrive at the following contradiction.
Therefore, ℘ ∞ = ∞ is valid. In view of system (1), we have Solving the above equation, we obtain that which implies that On the other hand, This implies that is, Solving the above inequality, For the larger time, we have and the proof is ended.

Theorem 2.
Consider ðS * ðtÞ, I * ðtÞ, R * ðtÞÞ as the solution of the stochastic epidemic model (1). Then, the disease will die out exponentially with probability 1. In other words, as t ⟶ ∞.

Journal of Function Spaces
Proof. In the light of (1), we have Integrating (1) implies that This gives as t ⟶ ∞, and we have Q ⟶ 0. By putting I * = 0 and R * = 0 in (1), we have Solving this equation, we have which yields and so This ends the argument.

Theorem 3.
For each arbitrary initial value ðS * ð0Þ, I * ð0Þ, R * ð0ÞÞ ∈ ℝ 3 + , the disease is persistent and Proof. From system (1), we have Integrating the above equation, we have Also, integrating the third equation of (1), we obtain which implies From (27) and (28), and also, We construct a real-valued C 2 -function V on ℝ by the following formulation. Thus, By integrating on both sides of above relation, we get This scheme is utilized for the numerical results of the stochastic model (1). The numerical results are presented via three graphs which show a comparative analysis of the deterministic and stochastic models.

Conclusions
We provided and analyzed a novel stochastic mathematical model for the COVID-19 spread with white noise in this work. The model was mathematically examined for the exis-tence of a solution and other factors such as disease persistence and disease extinction. We presented the numerical scheme for the computational analysis of the model after the mathematical analysis, and the scheme was then tested for the available parametric and initial data in the literature.  The computed findings clarified us with more realistic facts. Several diagrams were used to illustrate the numerical analysis. We mostly used two parameters in the graphs. We gave a graphical explanation of model (1) for different values of the white noise σ 1 , σ 2 . Also, We plotted all classes for different orders. The numerical simulations are presented with the help of the graphs. In the graphs, we have numerical data for the different values of the white noise σ 1 , σ 2 . In Figure 1, we have the susceptible class of model (1) for σ 1 = σ 2 = 0:0,0:02, 0:025, 0:03. Figure 2 represents numerical data for the infection rate in the model as per given initial and parametric values while varying the white noise for σ 1 = σ 2 = 0:0,0:02, 0:025, 0:03. Figure 3 shows simulations for the recovered class of model (1), changing the white noise for σ 1 = σ 2 = 0:0,0:02, 0:025, 0:03. For future work as a continuation of the study, we aim to analyze model (1) in the fractional order settings.

Data Availability
No data were used to support this study.

Conflicts of Interest
The authors declare that they have no competing interests.