Theoretical and mathematical codynamics of nonlinear tuberculosis and COVID-19 model pertaining to fractional calculus and probabilistic approach
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel virus known as coronavirus 2 (SARS-CoV-2) that affects the pulmonary structure and results in the coronavirus illness 2019 (COVID-19). Tuberculosis (TB) and COVID-19 codynamics have been documented in numerous nations. Understanding the complexities of codynamics is now critically necessary as a consequence. The aim of this research is to construct a co-infection model of TB and COVID-19 in the context of fractional calculus operators, white noise and probability density functions, employing a rigorous biological investigation. By exhibiting that the system possesses non-negative and bounded global outcomes, it is shown that the approach is both mathematically and biologically practicable. The required conditions are derived, guaranteeing the eradication of the infection. Sensitivity analysis and bifurcation of the submodel are also investigated with system parameters. Furthermore, existence and uniqueness results are established, and the configuration is tested for the existence of an ergodic stationary distribution. For discovering the system's long-term behavior, a deterministic-probabilistic technique for modeling is designed and operated in MATLAB. By employing an extensive review, we hope that the previously mentioned approach improves and leads to mitigating the two diseases and their co-infections by examining a variety of behavioral trends, such as transitions to unpredictable procedures. In addition, the piecewise differential strategies are being outlined as having promising potential for scholars in a range of contexts because they empower them to include particular characteristics across multiple time frame phases. Such formulas can be strengthened via classical technique, power-law, exponential decay, generalized Mittag–Leffler kernels, probability density functions and random procedures. Furthermore, we get an accurate description of the probability density function encircling a quasi-equilibrium point if the effect of TB and COVID-19 minimizes the propagation of the codynamics. Consequently, scholars can obtain better outcomes when analyzing facts using random perturbations by implementing these strategies for challenging issues. Random perturbations in TB and COVID-19 co-infection are crucial in controlling the spread of an epidemic whenever the suggested circulation is steady and the amount of infection eliminated is closely correlated with the random perturbation level.
Keywords: Co-infection TB and COVID-19 model; Deterministic and probabilistic model; Fractional calculus; Global positive solution; Ergodic stationary distribution; Probability density function; Quasi-equilibrium
Introduction
The COVID-19 outbreak has posed novel obstacles to worldwide medical systems, resulting in enormous impacts on nations around the globe. Undoubtedly, the battle against COVID-19 has taken up much of the attention, but it is important to remember that TB has existed as a problem for quite a while. Mankind has been plagued by this extremely contagious sickness for ages. Ultimately, 2020 will likely go down in history as the year that the coronavirus ailments, or COVID-19, took center stage. The outbreak's causative agent, the SARS-CoV-2, first appeared in China in the second half of 2019[1],[2]. Even though COVID-19 continues to be a topic widely discussed in academic journals and news reports, it's crucial to remember about other infectious illnesses, such as TB[3],[4].
The COVID-19 outbreak has had a major effect on the TB treatment mechanism, resulting in a reduction in both detection and transmission. This is explained by the repercussions of TB care and limitations on accessibility for patients, which have led to an increase in TB-related mortality[5],[6]. In order to successfully combat both of these transmissible illnesses, this viewpoint assessment seeks to point out the overlap between COVID-19 and TB, emphasizing their combined menace and suggesting common approaches.
Furthermore, there are some notable clinical commonalities between the COVID-19 epidemic and TB. Since pulmonary secretions are an important way that these ailments are communicated, proximity and congested surroundings are favorable. Furthermore, COVID-19 and TB are especially dangerous for disadvantaged and underprivileged groups, such as the elderly, people with preexisting medical disorders, and people with compromised immunological capabilities. The COVID-19 epidemic has had a complex effect on TB. The increased challenge has caused a diversion to medical supplies, which has disrupted attempts to diagnose, address and regulate TB. Security measures, prohibitions on traveling and restricted availability of healthcare resources have made it more difficult to identify cases of TB and diagnose patients on time. The combination of these two contagious illnesses has produced a complicated scenario that needs prompt monitoring and all-encompassing approaches. Both COVID-19 and TB have a number of similarities, most notably the way in which their respective causal agents-mycobacterium TB and SARS-CoV-2-are transmitted[7]. Pulmonary system emissions are the route of transmission for both infections[8],[9]. Both COVID-19 and TB can spread via aerosols and droppings, with the respiratory tract being their usual site of infection. It is crucial to remember, though, that such illnesses may impact a variety of organs[10]. In addition, finding and evaluating interactions as well as safeguarding medical personnel and individuals at risk are essential elements of healthcare safety for these illnesses. To create comprehensive and inexpensive prevention and treatment strategies, it is essential to comprehend the channels and components impacting propagation. Numerous decades of therapeutic and laboratory research on TB have yielded a plethora of data that can be used to identify, prioritize and evaluate exposures[8]. It should come as no surprise that more research is needed to better understand how SARS-CoV-2 spreads, and there is ongoing debate regarding the distinct functions played by airborne particles, microbes and big pulmonary secretions[11]. In particular, excessive growth occurrences have been linked to the propagation of these two diseases[12],[13]. Figure 1 listed below shows a graphic that illustrates several of the prevalent therapeutic manifestations and multi-organ dysfunction.
Graph: Figure 1 Identical indications and multi-organ connection to TB and COVID-19.
Whereas the implantation time for TB can range from 2 weeks to many decades until the TB infection advances, that of COVID-19 is less lengthy, ranging from 1 to 14 days. The manifestations of COVID-19 include anemia, wheezing, throat irritation, diminished or absent perception, flavor loss, vomiting, muscular soreness, and exhaustion. Usually, these indications start off suddenly. On the other hand, TB causes a high temperature, perspiration at night, a persistently persistent cough, bleeding in the cough, decreased hunger, heartburn, and exhaustion. On the other hand, TB manifestations appear gradually and with a subtle beginning. When it comes to COVID-19, those with coexisting illnesses, including HIV, insulin resistance, being overweight, persistent lung disease, persistent cardiac problems, and impaired immune systems, are more likely to have extreme symptoms. These inherent medical issues may exacerbate the advancement of the sickness. Conversely, concomitant conditions, including type II diabetes, sickle cell syndrome, severe obstructive pulmonary ailments, HIV, and a weakened immune system, are recognized to escalate the likelihood and intensity of TB transmission. Combating and curing serious forms of both COVID-19 and TB require an understanding of and commitment to controlling these coexisting conditions[14],[15]. Figure 2 lists some of the prevalent danger indicators for both TB and COVID-19.
Graph: Figure 2 Basic danger signs for TB and COVID-19.
Meanwhile, the testing process facilities have been disrupted by COVID-19, resulting in decreased personnel objectives, longer evaluation processing times, and the unavailability of forensic equipment. The timely delivery of TB screening tests as well as their accessibility have been greatly impacted by these delays. Screening findings may take longer to reach people, which could postpone therapeutic beginnings and raise the danger of tuberculosis spreading throughout populations[16]. Screening TB infections and locating regions with widespread dissemination require efficient acquisition and inspection methods. The distribution of resources and focused treatments can be guided by observational reports. Effective use of statistical analysis and health monitoring networks can help with preventive choices and offer real-time information[17].
During the years, a great deal of mathematical concepts have been developed to help us understand the world in which we live. In order to regulate presentations involving considerable obstacles, powerful artificial intelligence algorithms have been constructed, and the concept of space and time modeling has been put into practice. Some of the algorithmic techniques that are particularly commonly applied in modeling and prediction involve the idea of differentiation. Differential equations (DEs) are scientific techniques that have been created using this concept. In the beginning endeavor, researchers suggested a number of algorithms via multifaceted associations. The variation in the compositions could include local (exchange rate, conformable derivative, and fractal derivative)[18]–[20]; nonlocal/singular kernel (Riemann–Liouville, Liouville–Caputo, and multiple expressions)[21]; local/non-singular kernel (Caputo–Fabrizio operators)[22]; and finally non-local/non-singular (Atangana–Baleanu–Caputo operators)[23]. For a variety of interpretations of differential derivatives or the individuals who structured the foundations, a number of academics suggest numerous novel approaches. Fractional-order (FO) calculus has a connection to realistic endeavors and finds extensive application in multiple domains such as atomic physics, optics, image encryption, nanotechnology, and infectious disease[18]–[20],[24].
Recently, a subfield of mathematical physics and comprehension known as fractional calculus uses FO derivatives to study how inventions and documentation operate. FO modeling, as opposed to integer-order settings, can employ reminiscence memory of the power, exponential decay, or generalized Mittag–Leffler (GML) formation kernel to capture non-local spatial–temporal interactions. The conceivable benefits of using the fractional approach by Atangana–Baleanu involve all non-localities that are inherent within the explanation, just like in all previous variations. However, the most significant characteristic is that it has a nonsingular and non-local kernel, represented by the GML functionality, which, from an empirical viewpoint, includes the clarification and advancement of competencies delineated by a series of privileges. Kumar et al.[25] contemplated a new investigation on fractional HBV models through Caputo and Atangana–Baleanu–Caputo derivatives. Mekkaoui et al.[26] presented the predictor–corrector for non-linear DEs and integral equations with fractional operators. Atangana and Araz[27] described a successive midpoint method for nonlinear DEs with integer and Caputo–Fabrizio derivatives. On the other hand, it has been shown that the previously mentioned approach precisely conveys the complex compositions of many practical representations[23],[27]. The piecewise derivative, which has recently gained prominence[28], was presented by Atangana and Araz[29] and distinguishes from every derivative by the fact that it may reprise the interconnected paths that comprise these fractional algorithms in a differentiation technique. Every aspect that happens demonstrates that while the prevalence of the codynamics of COVID-19 and TB is probabilistic instead of deterministic in nature, knowledgeable research is based on an empirical methodology. A number of academics investigated the real-world growth of viruses and bacteria using the fundamental concept of probabilistic modeling, as reported in Refs.[30],[31].
Certain probabilistic COVID-19 individuals via TB concurrent infection outbreak frameworks using theropetic representatives hydroxychloroquine, azithromycin, lopinavir/ritonavir, and darunavir/cobicistat conjunction systems have been successfully defined to examine the influence of probabilistic white noise and offer several efficient initiatives for governing infection interactions. These frameworks are founded on the randomly generated linear disruptive methodology, which assumes the biological nature of ambient white noise correlates to the dimension of every compartment. Moreover, it has been demonstrated experimentally that a probabilistic COVID-19 and TB system that includes immunological dysfunction affected by inherent and adapted resistance can prevent an epidemic of co-infection. Motivated by these findings[32],[33], we also presume that the random perturbation is closely connected to specific populations of evolution of TB and COVID-19 diagnostics. In order to illustrate the significant influence of a probabilistic framework condition mentioned in Ref.[34], we performed to create this paper. Additionally, we create a probabilistic mathematical structure utilizing piecewise fractional derivative expressions to analyze the co-infection process incorporating the positive immunomodulation against COVID-19, likely because of trained innate immunity and crossed heterologous immunity within predetermined time intervals. In order to achieve this, we separate the population into two groups: the incidence and occurrence of exacerbated immune dysregulation and decreased lymphocyte function, along with erroneous variations. The probability that the most recent COVID-19-infected TB will be engaged is represented by the proportion , whereas the unexplained component will not be implicated. In addition, we established the global positive solutions of the co-infection model with a unique ergodic and stationary distribution (ESD) technique to illustrate the biological properties and statistical viability of this structure. We also provide the precise definition of the probabilistic density function (P.D.F) at a quasi-equilibrium point that represents the probabilistic COVID-19 approach, which reflects significant spontaneous features in probabilistic relevance. The ESD and P.D.F surrounding the quasi-equilibrium point of the randomized multidimensional codynamics framework will be better understood as a result of this investigation. The intention of the investigation is to acquire an improved comprehension of how the infection persists over time in the probabilistic codynamics system. In general, fractional operators examine simulations conducted numerically of the proposed system that include crossover structures and white noise.
Codynamics model and preliminaries
The general population is divided into eight indistinguishable groups in this category, which are designated as susceptible people, , latent TB patients who do not exhibit TB-associated indications and are not pathogenic , influential TB-infected people , COVID-19-infested humans who do not exhibit indications but are transmissible , COVID-19-diagnosed people who exhibit scientific backing indications and are pathogenic , both inactive TB and COVID-19-contaminated people , current TB and COVID-19-contaminated humans , and retrieved people consisting of both TB and COVID-19. The underlying computational framework for the codynamics of TB and COVID-19 is developed in this portion. Considering such, all people at moment , represented by , are provided by
1
Graph
We hypothesized that acquisition increases the vulnerable community at an intensity of . Every person in every compartment experiences an inevitable mortality rate of . Equivalent to formula (1), vulnerable individuals contract TB via interaction with current TB individuals via agent transmission . The acceptable interaction rate for TB transmission is indicated by within this manifestation. It is believed that people with persistent TB are undiagnosed and cannot pass on the illness[35]. Likewise, those at risk contract COVID-19 at an intensity of transmission , which is determined as in formula (1), after effectively coming into proximity to COVID-19-infected people. The efficient interaction probability for COVID-19 infection is represented by in this case. Furthermore, we hypothesized that people in the hidden TB segment depart at an incidence of to segment , at an incidence of transmission of to the persistent TB as well as COVID-19 contaminated group, whilst certain of them recuperate at an intervention incidence of . Additionally, those in the TB-infected category recuperate due to the illness at an incidence of , with the surviving percentage either transferring to the transmission category at a pace of or dying at a speed of via TB-induced mortality.
The overall community in cohort potentially dies at COVID-19-induced mortality pace or advances at an intensity of to become contaminated category . As seen in Fig. 3, it is believed that the other people will be moved to the other cohort at a consistent multiplicity of . In other words, the general population classified as migrates at an intensity of to category , then at a pace of to compartment group, and finally recovers at a pace of . Additionally, we hypothesized that, although the codynamics-induced mortality prevalence is represented by , people in compartment depart for compartments or , correspondingly, at an intensity of or Furthermore, at an intensity of or the COVID-19 exposure people can choose to depart to compartment or , respectively. Comparably, the number of individuals in compartment is either moved to the codynamics cohort at an intensity of or restored at a steady pace of . represents the disease-induced fatality rate within this category. Figure 3 displays the suggested system's process layout.
Graph: Figure 3 Flow diagram for depicting the codynamics process of TB-COVID-19 model (2).
It leads to frameworks for the subsequent nonlinear DEs determined by the procedure illustration:
2
Graph
where and containing positive initial conditions (ICs)
Table 1 Description of model's parameters.
Symbols | Description |
---|
| People' spontaneous mortality rate |
| Rate of transmission of COVID-19 andTB exposure within the contaminated group |
| Transmission rate of infection among those inoculated to COVID-19 |
| Probability of recuperation for a COVID-19 influenced person |
| Percentage of TB exposure in people subjected to COVID-19 |
| Recuperation percentage of inactive TB infections |
| Percentage of COVID-19 afflicted people who have TB disease |
| COVID-19 contamination incidence among TB patients |
| Probability of recuperation for those subjected to COVID-19 |
| Probability at which people exit the group |
| TB healing rate for of long-term care residents |
| Percentage of patients recuperating with COVID-19 |
| Proportion at which people quit the affected group |
| patients' rate of TB recurrence |
| COVID-19 recuperation percentage among participants |
| Percentage of people suffering from TB who also get COVID-19 |
| Recoverability percentage of TB patients |
| Mortality rate from TB |
| Mortality rate from COVID-19 infection |
| Mortality incidence as a result of both infections co-occurring |
| Recruiting rate for those who are vulnerable |
| Prevalence of TB infection propagation |
| Rate of COVID-19 propagation |
| Percentage of people confronted with TB who get the disease |
| Intensity transmission for TB (the likelihood of contracting a virus from a TB sick person) |
| COVID-19 intensity of illness: the likelihood of contracting the virus from a person who has COVID-19 illness |
Table 1 provides a description of the system's characteristics.
To help readers that are acquainted with fractional calculus, we provide the related summary herein (see[21]–[23] comprehensive discussion on fractional calculus).
Graph
The index kernel is involved in the Caputo fractional derivative (CFD). Whenever experimenting with a particular integral transform, such as the Laplace transform[36],[37], the CFD accommodates regular ICs.
Graph
where indicates the normalization function
The non-singular kernel of the Caputo-Fabrizio fractional derivative (CFFD) operator has drawn the attention of numerous researchers. Furthermore, representing an assortment of prevalent issues that obey the exponential decay memory is best suited to utilize the CFFD operator[38]. With the passage of time, developing a mathematical model using the CFFD became a remarkable field of research. In recent times, several mathematicians have been busy with the development and simulation of CFFD DEs[39].
The ABC fractional derivative operator is described as follows:
Graph
where represents the normalization function.
The memory utilized in Atangana–Baleanu–Caputo fractional derivative (ABCFD) can be found intuitively within the index-law analogous for an extended period as well as exponential decay in a number of scientific concerns[40],[41]. The broad scope of the connection and the non-power-law nature of the underlying tendency are the driving forces behind the selection of this version. The impact of the kernel, considered crucial in the dynamic Baggs–Freedman framework, was fully produced by the GML function[42].
To far better perceive the propagation of TB and COVID-19, we indicate a dynamic mechanism (2) that includes the co-infection within the context of CFD, CFFD and ABCFD, respectively. This is because FO algorithms possess inherited properties that characterize the local/non-local and singular/non-singular dynamics of natural phenomena, presented as follows:
3
Graph
4
Graph
5
Graph
The arrangement of this article is as follows: In "Codynamics model and preliminaries" section, explanations for fractional calculus, along with several key notions and model (2) details, are provided. Moreover, a detailed analysis of the FO co-infection system's (3) equilibrium stability is presented in "Codynamics model and preliminaries" section. In "Stochastic configuration of codynamics of TB-COVID-19 model" section, a probabilistic form of the TB and COVID-19 models' (28) codynamics is proposed and a detailed description of the unique global positive solution for each positive initial requirement is presented. The dynamical characteristics of the mechanism's appropriate conditions for the presence of the distinctive stationary distribution are provided. The P.D.F enclosing a quasi-stable equilibrium of the probabilistic COVID-19 framework is presented in "Stochastic COVID-19 model without TB infection" section. Numerous numerical simulations in view of piecewise fractional derivative operators are presented in "Numerical solutions of co-dynamics model using random perturbations" section to validate the diagnostic findings we obtained in "Stochastic configuration of codynamics of TB-COVID-19 model" and "Stochastic COVID-19 model without TB infection" sections. In conclusion, we conceal our findings to conclude this study.
Positivity and boundedness
Since we interact with living communities, each approach ought to be constructive and centred on a workable area. We utilized the subsequent hypothesis that guarantees these.
Theorem 1
Assume that the set is a positive invariant set for the suggested FO model (3).
Proof
In order to demonstrate whether the solution to a set of equations (3) is positive, then (3) yields
6
Graph
Therefore, the outcomes related to the FO model (3) are positive. Finally, the variation in the entire community is described by
Graph
Addressing the variant previously mentioned, we get
Graph
Consequently, we derive the GML function's asymptotic operation[43] as
Graph
Adopting the same procedure for other systems of equations in the model (3), which indicates that the closed set is a positive invariant domain for the FO system (3).
- Assuming that every requirement is non-negative throughout time , we exhibit that the outcomes remain non-negative and bounded in the proposed region, . We'll look at the co-infection model (3) spreads in the domain, which is described as
- According to the afflicted categories in co-infection model (3), disease-free equilibrium (DFE) and endemic equilibrium (EE) are the biologically significant steady states of FO model (3). We establish the fractional derivative to get the immune-to-infection steady state as to zero of the FO model (3) have no infection, and get
-
Graph
- The dominating eigenvalue of the matrix correlates with the basic reproductive quantity of structure (3), in accordance with the next generation matrix approach[44]. Thus, we find
-
Graph
The next generation matrix at DEF can then be obtained by using the Jacobian of and examined at as
Graph
where
Graph
The fundamental reproducing quantity of the pairing system is shown by the highest spectral radius of the subsequent generation's matrix. It is evident that the matrix has four eigenvalues that are equivalent to zero. The truncated matrix yields the additional eigenvalues as
Graph
Consequently, by calculating the eigenvalues of , it is possible to simply determine that
Graph
where
Graph
Therefore, the co-dynamics structure's (3) fundamental reproductive quantity is provided by
Here, we shall then demonstrate how transmission persists in the FO mechanism. It explains how widespread the virus is within the framework. From the viewpoint of biology, the virus continues in the bloodstream if the infectious proportion is elevated for a sufficiently long time .
However, the linearization technique is used to examine the local stabilization of the codynamics algorithm's DFE state. At the DFE state the Jacobean matrix of system (3) is displayed as
7
Graph
The analysis of 's localized temporal equilibrium relies upon the eigenvalues' interpretation. Here, and are obtained by broadening the following polynomial . Moreover, we get the additional 's based on the simplified matrix's described as
Graph
where After simple computations, the characteristic polynomial of the above matrix is presented as
8
Graph
In other words, the outcomes to the are the eigenvalues:
9
Graph
where
10
Graph
Therefore, if the subsequent requirements apply, the roots of expression (10) exhibit negative real portions according to the Routh–Hurwitz stability specifications as
11
Graph
Graph: Figure 4 Evolution of the basic reproduction number R0CT with the aid of R0C and R0T.
Figure 4 is illustrated by depicting in 3D evolution of the threshold parameter of model (3) as a function of and
The forthcoming result is established thanks to Theorem 2 in[44].
Theorem 2
The DFE point of the FO codynamics model (3) is locally asymptotically stable if the prerequisite specified in formula (12) is satisfied.
Existence and uniqueness of solutions
This section shows that there is only one solution for the system (3). Now, we demonstrate that the framework's solution is distinctive. Initially, we construct framework (3) in the form of:
12
Graph
where
13
Graph
Integral transform applied to both sides of equations (14) yields
14
Graph
The kernels satisfies the Lipschitz condition and contraction, as demonstrated.
Theorem 3
satisfies the Lipschitz condition and contraction if the following condition holds:
Proof
For and we have
Graph
Suppose , where are a bounded functions. So, we have
15
Graph
After obtaining the Lipschitz criterion for , hence, is a contraction if .
In the same manner, satisfy the Lipschitz condition as follows:
Graph
where
For we find then are contractions. Assume the following recursive pattern, as suggested by system (15):
Graph
with
Throughout the above system, we compute the norm of its first equation and then
Graph
Therefore, (16) possesses Lipschitz's condition, then we have
Graph
Analogously, we find
16
Graph
As a consequence, we can write
Graph
Theorem 4
A system of solutions described by the codynamics model (3) exists if there exists such that
Proof
By means of (16) and (17), we have
Graph
Thus, the system is continuous and has a solution. Now we shall explain how the functions listed above may be used to construct a model solution (15). We make the assumption that
Graph
Therefore, we have
Graph
After recursive procedure, we have the following:
Graph
Thus,
Similarly, we may establish that
To examine the uniqueness of the solution, we assume that there is another solution of the system, such as Then
Graph
After taking norm, we get
Graph
Utilizing the Lipschitz condition, we have
Graph
Consequently, we have
17
Graph
Theorem 5
The codynamics model (3) has a unique solution, provided that
Proof
Assuming that condition (18) is vaild,
Graph
Then Hence, we have Similarly, we can prove that
Influence of TB on COVID-19
We started by describing the basic reproductive quantity, , by means of (and vice versa), in order to examine the effect of TB illness on COVID-19 (and vice versa)[45]. By interpreting the value of as a component of using the formula (2), we get
18
Graph
Now, we have
19
Graph
where Furthermore, the it also indicates that the COVID-19 outbreak is made worse by the spread of TB viruses.
Remark 1
The population's TB proliferation possesses no noticeable influence with the propagation of COVID-19 provided . On the other hand, the transmission of COVID-19 will be significantly adversely affected by the outbreak of TB if .
Furthermore, by quantifying in the context of and determining the meaning of the partial derivative of with regard to , the effect of COVID-19 of TB infections is able to be rectified.
Analysis of COVID-19
When the infection of TB is disregarded, the deterministic model (2) becomes the subsequent system:
20
Graph
The basic reproduction number of the model (21) is presented as
21
Graph
Sensitivity analysis
The sensitivity analysis of the model parameters for the COVID-19 submodel, as stated in (21), is carried out in this subsection. The sensitivity of a parameter, contemplate, is expressed as[46] and indicates how the framework behaves in response to a slight variation in a parameter value as
Graph
In our case, the sensitivity analysis of each parameters for (21) becomes:
22
Graph
Here, we notice that the dissemination of COVID-19 is boosted by the contact rate . Additionally, the transmission rate from the unprotected group to the afflicted group has a positive effect on the dissemination of the virus if In other words, the prevalence will rise as the values of these factors rise. The additional parameters and have adverse effects; therefore, raising their values will result in a drop in the frequency of COVID-19 infections. Nonetheless, the sensitivity analysis investigation does not take into account the immoral increase in the individual fatality rate as a means of controlling the spread of illness.
Bifurcation analysis
In what follows, we investigate the solution behavior of (21) by taking as the bifurcation parameter. Calculating the value of from , i.e, we have By replacing , we can determine the eigenvalues of the Jacobin matrix at the DFE point, as per the outcome provided in Ref.[44]. Thus, substituting in (7), it gives zero eigenvalue. This means that the Jacobean matrix in (7) at has a left eigenvector (associated with the zero eigenvalue) which is calculated from . Here, , for which and
Likewise, can be used to determine the right eigenvector linked to the zero eigenvalue when for which and
Now, suppose that represents the right-hand side of the equation in the COVID-19 submodel (21) and let denote the corresponding state variable for .
Introduce
23
Graph
The local dynamics of (21) near the bifurcation point are then calculated by the signs of two associated constants and with . Note that, in , the first zero corresponds to the DFE, , for (21). In other words, for if and only if the right-hand sides of (21) are equal to zero at
Moreover, from , we have when , which is the second zero component in . For the (21), the associated nonzero partial derivatives at the are
24
Graph
Next, with the aforementioned expressions for , it is evident that
25
Graph
It can be demonstrated that the corresponding non-vanishing partial derivatives for the corresponding sign of are
26
Graph
It is evident from the aforementioned statements as well that
From the bifurcation coefficient 's sign varies on the minimal value of recurrence that generates bistability , and we find that is always positive from the estimates of and . Therefore, a subsequent proposition is established by applying the result of Ref.[44].
Proposition 1
The system (21) has a forward bifurcation if the minimal value of the virus infection that causes bistability is smaller than unity.
We then conclude with the following theorem
Theorem 6
If then
- (i)The model (21) undergoes a backward bifurcation whenever
-
• .
- (ii)The model (21) undergoes a forward bifurcation whenever
-
• .
The bifurcation phenomenon is illustrated in Fig. 5a,b, where we have carried out a numerical simulation of the infection model (21). The system parameter values are presented in Table 2, the calculation gives and , the backward bifurcation condition is then satisfied and we obtain Fig. 5a. Also, the forward bifurcation Fig. 5b is obtained for Here, the parameters we have used may not all be epidemiologically realistic (see[47]).
Graph: Figure 5 Simulation of the codynamics model (3) to illustrate the occurrence of (a) forward (b) backward bifurcation.
Stochastic configuration of codynamics of TB-COVID-19 model
We examine how random interference affects the distinctiveness and presence of a stable dispersion, as well as the gradual disappearance of diseases, in the system (2). The formula that follows is a representation of the stochastic adaptation relating to the model (2) is
27
Graph
in which indicate the variability in noise and are conventional one-dimensional autonomous Brownian movements. The additional parameters have the same relevance as they do in system (2).
In the sequel, let be a complete probability space and its filtration needs to fulfill the standard requirements (that is., it must be right continuous and comprise all -null sets), whilst are stated on the complete probability space. In addition, take For any matrix its transpose is indicated by
Utilizing as the solution of model (28) supplemented by ICs Furthermore, we utilize to represent and to show
Firstly, we assert an outcome about the existence–uniqueness of a global non-negative solution for system (28).
Theorem 7
Assume that there is a unique solution of structure (28) on for any starting value It stays in having probability 1 (a.s).
Proof
Here, we overlook the initial portion of the explanation just to display the essential Lyapunov function because it is comparable to Theorem 2.1 in[34].
Introducing a -functional on by
28
Graph
where the value of the non-negative constant will be obtained hereafter. When we implement Itô's algorithm[48] to , we obtain
29
Graph
where is determined by
30
Graph
Letting As a result, we have
31
Graph
where the constant is non-negative. According to Ref.[34], we similarly exclude the remaining portion of the explanation. The documentation is now complete.
Stationary distribution
Our primary concern with the stochastic outbreak framework is the virus's permanence. In this portion, we employ a novel method to demonstrate that structure (28) has a unique ESD, depending on the hypothesis of Khasminskii[49].
By developing appropriate Lyapunov functions, we will show adequate conditions for the development of a unique ESD. A key component of our major result's explanation is the lemma that follows.
Assume that is an ordinary time-homogeneous Markov phenomenon with . Its stochastic DE is as follows:
32
Graph
and the diffusion matrix is Consider the differential operator L connected to (36) as follows:
33
Graph
Lemma 1
([49]) Let us suppose the subsequent characteristics of a bounded open region with a regular boundary:
( ): In the region and some neighborhood therefore, the least significant eigenvalue of the diffusion matrix is bounded away from zero.
( ): a positive -function so that is negative for all
Then the Markov procedure has a stationary distribution Also, consider be a mapping which is positive in regard to the measure ones obtain
34
Graph
To begin with, we establish a few concepts for ease of use in later explanations. By resolving the subsequent (36) as
35
Graph
we find
36
Graph
Afterwards, by addressing the subsequent formula's:
37
Graph
which lead us
38
Graph
where with
Graph
Introduce
39
Graph
Theorem 8
If we suppose that , then structure (28) permits a unique ESD,
Proof
It is necessary to validate assumptions and in Lemma 1 for the purpose of establishing Theorem 8.
To begin with, we create an appropriate Lyapunov function and identify a closed set that ensures is negative in order to ensure the efficacy of in Lemma 1.
For this, let us suppose
Graph
Implementing Itô's technique to we find
40
Graph
Applying the variant we have
Again, considering (36), we have
41
Graph
Adopting the similar technique to and respectively, we have
42
Graph
Now, introducing a -mapping as
Graph
so that
43
Graph
where
Graph
where
Implementing the Itô's technique to and considering (42)–(44), we have
44
Graph
Then, we describe a -function as
45
Graph
which leads to
46
Graph
where
Graph
Graph
Implementing the Itô's technique to and considering (45)–(47), we have
47
Graph
Furthermore, we indicate
Graph
Implementing Itô's technique to one obtains
48
Graph
where
Graph
Analogously, implementing Itô's technique to and we find
49
Graph
Introducing
50
Graph
which leads to
51
Graph
where
Graph
Graph
Graph
As
Now, considering (49)–(52), we have
52
Graph
Furthermore, we describe
53
Graph
Thus, we conclude that and
In view of (45), (48) and (53), we have
54
Graph
Introducing
55
Graph
Again, implementing the Itô's technique to we have
56
Graph
Introducing
57
Graph
Implementing the Itô's technique to we have
58
Graph
Again, we describe
59
Graph
where fulfilling
Graph
Employing the Itô's technique to we have
60
Graph
where
61
Graph
and
62
Graph
Here, introducing a -function on
63
Graph
where is a sufficiently significant non-negative constant that satisfies
64
Graph
and
65
Graph
Examine that the minimal point of appears to exist, therefore we conclude
66
Graph
Merging (57), (59) and (61), we have
67
Graph
Next, we construct the following for a bounded closed set:
68
Graph
where is a non-negative constant that is small enough to meet the subsequent variants
69
Graph
having
70
Graph
For the sake of simplicity, we can split into the subsequent sixteen regions:
71
Graph
Evidently, Consequently, it is easy to demonstrate that
72
Graph
This verifies assumption of Lemma 1.
The diffusion matrix for model (28) is presented as follows:
73
Graph
It is evident that, matrix is positive definite This verifies assumption of Lemma 1. Thus, the model (28) has a unique stationary distribution and ergodic. This puts the proof to its conclusion.
Remark 2
For system (28), if the illness always endures. According to Theorem 8, if , the sickness will continue in the stochastic framework (28). In particular, in the absence of noise, that is, Observe that
Next, the equation is On the same instance, obtaining
74
Graph
Stochastic COVID-19 model without TB infection
By utilizing the identical technique from probabilistic framework (28) to incorporate random perturbation, we obtain the subsequent stochastic model:
75
Graph
Then, we state
Graph
The values of the parameters have similar significance within the system (75). Indicate
76
Graph
The two theorems that proceed are derived from "Codynamics model and preliminaries" and "Stochastic configuration of codynamics of TB-COVID-19 model" section using a similar methodology.
Theorem 9
Suppose there are initial values have unique solution of the model (75) with and the solution will exist in having probability 1 (a.s).
Theorem 10
Suppose that then model (75) possesses the ergodic functionality and yields a unique stationary distribution
Probability density function (P.D.F)
In what follows, we present a mathematical principle pertaining to the P.D.F associated with the subsequent probabilistic framework as:
77
Graph
where indicates the parameter whilst are some functions and is the Wiener technique.
Lemma 2
([50]) Suppose there is a mapping states the P.D.F associated to the formula (77):
Graph
Following that, we provide the prerequisites required to find the positive definite (P-D) 4D real symmetric matrix.
Lemma 3
Assume that there is a 4D real algebraic equation having while indicates the real symmetric matrix.
• If
-
Graph
- containing with and then
- 78
Graph
• If
-
Graph
- containing and then
- 79
Graph
• If
-
Graph
- containing and then
- 80
Graph
Proof
Indicate the -th significant main component of is which is expressed as
- Observe that then
-
Graph
- Furthermore, assertions (ii) and (iii) can be obtained in the same way.
Here, the precise representation of the density function of system (75) at a quasi-equilibrium point will be derived. In relation to analytical importance, it is important to note that the P.D.F can represent the majority of the unpredictable features of a probabilistic process.
Initially, we apply an analogous change to illustrate (75). For this, consider and Thus, system (75)'s corresponding expression is provided by
81
Graph
When we illustrate a quasi steady state where
82
Graph
Assume that Thus, system (81) can be expressed as
83
Graph
where
Furthermore, and
Define and we have
Graph
where
84
Graph
Next, we confirm that the real components of each of 's eigenvalues are negative. The characteristic polynomial of that corresponds to it is where
85
Graph
Following that, if then
86
Graph
Thus, and Subsequently it appears that A possesses every negative real-part eigenvalues that correspond to the Routh-Hurwitz stability condition[51].
With reference to Lemma 2, the Fokker–Planck equation below is satisfied by the relevant P.D.F to the Quasi-stationary condition of the system (81) can be expressed as
Graph
Given that is an invariant matrix, one can determine that is potentially identified as having a Gaussian distribution by incorporating the pertinent findings of Roozen[52]:
Graph
where justifying and is a real symmetric matrix fulfilling
87
Graph
If holds, we indicate can be found to possess the equivalent degree of positive definiteness. Following this, (87) has the structure that follows.
88
Graph
The mathematical structure of (88) is obtained by employing a finitely autonomous coherence theory, which gives us and then we have
Graph
where
Graph
and are decided upon thereafter.
Taking into account and the transformation between the frameworks (75) and (81) yields the following:
Graph
where
Theorem 11
Surmising that for any then the solution model (75) possess a log-normal P.D.F about as follows having is a positive definite matrix and the components and are described as
89
Graph
where
90
Graph
and the matrices and are illustrated in the subsequent result.
Proof
Case A: Considering
91
Graph
In view of the elimination matrix as
Graph
Consequently, we get
Graph
where
The subsequent sub-stages are then taken out of the appropriate evaluation.
Subphase AI Choose and then there is where and
Graph
Consequently, it is possible to write the appropriate formula of (91) as
Graph
By making the use of Lemma 3, we determine where
Graph
is a P-D symmetric matrix. Therefore, is also a P-D matrix.
Subcase AII Taking and also suppose that where
Graph
containing
Subcase AIII Taking and Moreover, suppose that where
Graph
containing and Hence, we get
Graph
By making the use of Lemma 3, we determine where
Graph
is a symmetric, semi P-D matrix. Thus,
Subcase AIV Taking and employing the analogous technique as we did in Subcasee AI. Suppose that so that Hence, we have where Thus, we conclude that is a P-D matrix.
Case B Considering
Graph
Assume that where
Graph
where and
Subcase BI When and suppose that where
Graph
containing and
In this way, we have
Graph
Taking into account Lemma 3, we have is a semi P-D matrix and
Graph
Finally,
Subcase BII When and applying the analogous approach as we did in the Subcase AI, then we attain so that hence, we get
Graph
where Thus,
92
Graph
is a P-D matrix.
Subcase BIII When and and suppose that where
Graph
Furthermore, we have so that Hence, we get
Graph
where Thus,
Graph
is a P-D matrix.
Subcase BIV When and suppose that where
Graph
Furthermore, we have so that Hence, we get
Graph
where Thus,
Graph
is a P-D matrix.
Subcase BV When and applying the identical technique from Subcase AIII), we obtain where
Graph
having
Graph
Moreover, and will be determined later.
Then, we obtain
Graph
where and
Graph
Therefore,
Subcase BVI If and applying the analogous approach as we did in Subcase AIV. Assume that so that Hence, we have
Graph
where Therefore, is a P-D matrix.
Case C Surmise that
Assume that where
Graph
having
Subcase CI When employing the identical approach as we applied in Subcase AIII and consider that where
93
Graph
containing and will be computed later.
Thus, we find
Graph
Using the fact of Lemma 3, we have where
Graph
Consequently,
Subcase CII When then applying Subcase AI with a similar technique, resulting in so that which leads to
Graph
where
Graph
Hence, we conclude that is a P-D matrix.
Case D Considering and also, we have where
Graph
Indicate so that Thus, we get
Graph
where
Graph
This concludes that is a P-D matrix. Finally, the expression is a P-D matrix. So, the solution of model (75) possess a log-normal P.D.F about as
Graph
This yields the intended result.
Numerical solutions of co-dynamics model using random perturbations
The computation methods of stochastic perturbations influence whenever differentiating expressions involve fractional differential compositions involving singular or nonsingular kernels, and classical prescriptions include this component. The fractional notions have an order corresponding to 0 and 1.
Caputo fractional derivative operator
The main objective of this study is to investigate the co-infection of the TB and COVID-19 models involving integer-order (2), power-law (3) and stochastic strategy for (28). This scheme incorporates substantial pulmonary inflammation, which makes the circulatory mechanism a key battleground for numerous ailments. In the situation where is chosen as the final propagation period, the mathematical framework will be built using the classical-order expression in the beginning, the power-law memory considered in the next step, and the stochastic configuration in the stages that follow. After the fact that the subsequent number pattern is provided to explain the incidence.
Specifically, we analyze the sectionally divided frameworks (2), (3) and (28) quantitatively by using the procedure given in[29] in the context of CFD. In order to outline the procedure, we conducted what follows:
Graph
Thus, it implies that
Graph
where
94
Graph
95
Graph
and
96
Graph
Caputo–Fabrizio fractional derivative operator
The aim of this research is to examine the co-infection of the TB and COVID-19 models using integer-order (2), exponential decay kernel (4) and the ensuing stochastic scheme (28). This plan includes significant pulmonary inflammation, which means that the circulatory system is a major site of disease combat for a variety of diseases. The mathematical structure will be constructed using the classical-order formulation at first, the exponential decay memory at a later stage, and the stochastic setting in the phases that proceed in the case when is selected as the ultimate dissemination time. Following this, the following numerical pattern is given to clarify this occurrence.
At this point, we examine the sequential configurations (2), (4) and (28) analytically by using the method outlined in Ref.[29] in the context of the CFFD. In order to lay out the procedure, we did what follows:
97
Graph
Thus, it implies that
98
Graph
Atangana–Baleanu–Caputo fractional derivative operator
The current research aims to investigate the co-infection of the stochastic technique (28) and the integer-order model (2) and the GML kernel TB and COVID-19 models (5). Significant pulmonary inflammation is present in this design, indicating that the circulatory system is a key area of illness defense for a number of illnesses. Initially, the classical-order interpretation will be used to build the computational framework; afterwards, the GML function will be implemented; and in the phases that follow, the stochastic configuration will be used in the scenario where is chosen as the eventual propagation time. After that, the subsequent numerical structure is provided to explain these instances. In particular, we analyze the sequential configurations (2), (5) and (28) numerically employing the algorithm defined in Ref.[29] in the framework of the ABCFD. In order to lay out the procedure, we did what follows:
Graph
Thus, it implies that
Graph
where the previous values of , and are found in (94)–(96).
Experimental outcomes and discussion
In order to support research ideas, we will demonstrate mathematical simulation techniques in the next part that make leverage of the Atangana and Araz approaches formerly mentioned in Ref.[29]. The appropriateness and usefulness of the planned TB-COVID-19 are demonstrated through a number of concrete instances, such as manpower reductions, delays in test result transformation, and limitations of analytical equipment. The accessibility and promptness of TB examinations have been severely impacted by these interruptions in the deterministic-probabilistic situation. Utilizing MATLAB 21, all quantitative and symbolic computations were performed.
Researchers are at present demonstrating a great deal of enthusiasm in the estimation of modeling characteristics from provided statistical information, and it is thought to be an essential component of quantitative disease investigations. The aforementioned section was added to the current investigation employing the popular nonlinear least squares method. Applying the previously described method, the settings were determined, and the structure was calibrated to actual codynamic situations found in Later research from the Philippines and South Africa revealed that, for a specific duration, COVID-19 patients having TB had a 2.17[53] and 2.7[54] worse probability of death, respectively, than COVID-19 individuals lacking TB[53]. Especially, the entire number of documented infections and fatalities over the time span between March 2020 (the initial incidence had been identified on March 12, 2020) and June 2022 were used to determine the characteristics of the model. Considering the implementation of (99), the Ordinary Least Square solution was employed to reduce the inaccuracy concepts, and the associated relative deviation is incorporated in assessing the quality of fit as
99
Graph
The documented accumulative infection rates are denoted by in this particular instance, while the total number of contaminated occurrences determined by modeling execution is denoted by . The people who are moved daily from the contaminated compartment to the confined compartment are added together to determine the estimated levels of progressive transmission. With the exception of , which is envisioned, all the parameters are estimates. When and the data in Fig. 6 has been fitted to the model.
Graph: Figure 6 The codynamics of TB-COVID-19 fitting outcomes considering the data obtained from WHO[55] for weekly reported cases.
The parameters' projected estimates are displayed in Table 2.
Table 2 Details on the system's characteristic.
Notations | Values | References |
---|
| 500 | Supposed |
| 0.0477 | 53 |
| 0.6 | 53 |
| 0.659 | 54 |
| 0.02 | Supposed |
| 0.01 | Estimated |
| 0.023 | 54 |
| 0.01 | 54 |
| 0.05 | Calculated |
| 0.03 | Calculated |
| 0.03 | Calculated |
| 0.003 | Calculated |
| 0.021 | Supposed |
| 0.09 | 53 |
| 0.048 | Estimated |
| 0.25 | 54 |
| 0.01 | Estimated |
| 0.002 | Estimated |
| 0.2 | Assumed |
| 0.05 | 53 |
| 0.056 | Estimated |
| 0.95 | Estimated |
| 0.9 | Estimated |
| 0.25 | 53 |
Example 1
To illustrate our results, we quantitatively generate the paths for the probabilistic sickness structures (2), (3) and (28) and their corresponding deterministic components. The starting points are and the time range is [0, 100] units. Table 2 allows us to re-select the parameters to represent the piecewise methodology assessment of the naturally occurring factor process for (2), (3) and (28), respectively.
Here, we calculate the fundamental reproductive quantity for the deterministic framework (2), which suggests that co-infection of TB and COVID-19 will continue to exist in the average situation in both submodels. To observe how noise concentration affects the behavior of the probabilistic framework (28), we select random perturbations This yields The existence of an ESD for the probabilistic model (28) is demonstrated by Theorem 8.
When a power-law-type kernel with a FO is employed on (3), Fig. 7a–h illustrates how incorporating two propagation interprets enhances the occurrence of ailments in comparison to a single procedure. In view of the CFD operator and biological-nature strategy, we also find that certain combinations prove more deadly than others. A transmission surge is produced by all interactions comprising the effective collaboration rate connecting and . This is followed by any coupling via the effective interface rate between and , and finally other procedures. The respective two ICs and are depicted in Figs. 8a–h and 9a–h, with various population schemes.
Graph: Figure 7 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (3) and (28) with the impacts of latent and active TB outbreaks using CFD operator having FO ω=0.98, low intensities and ICs (1000, 100, 10, 1000, 10, 7, 5, 1).
Graph: Figure 8 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (3) and (28) with the impacts of latent and active TB outbreaks using CFD operator having FO ω=0.98, low intensities and ICs (4000, 100, 40, 1000, 10, 7, 5, 4).
Graph: Figure 9 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (3) and (28) with the impacts of latent and active TB outbreaks using CFD operator having FO ω=0.98, low intensities and ICs (5000, 100, 50, 1000, 10, 7, 5, 5).
We employ the identical factors as in (3), when implementing the identical methodology to the CFFD operator of DEs (2), (4) and (28), respectively. However, we modify the ICs as previously mentioned. It is simple to compute the threshold factors and As seen in Fig. 10a–h, the co-infections are expected to continue in a typical way, supporting the result of Theorem 8 (see Figs. 11a–h, 12a–h). According to this research, co-infection will spread throughout the body and develop ineffective causative agents, whereas mycobacterium TB will go obsolete.
In light of the prevalent concentrations and parameterization fluctuations discussed above, an intervention plan based on the computational findings for (2), (5) and (28) seems to be effective. There is an ESD of a probabilistic framework (28), as shown by Theorems 8. These results suggest that co-infection will become increasingly permanent while TB will go extinct for ABCFD case. These are corroborated by Figs. 13a–h, 14a–h and 15a–h, respectively.
Example 2
For probabilistic co-infection systems (28) involving community propagation, it is challenging to define appropriate criteria for virus extermination considering the limits of statistical approaches. Nonetheless, we provide a numerical model of the disappearance of illnesses where the noise is high for a thorough explanation. For instance, in the actual environment, individuals haphazardly raise vaccination or exterminating rates to stop co-infection from spreading. This successfully removes contamination.
To illustrate that high levels of environmental disturbance will eventually cause TB to disappear, we set with the identical setting off rate as well as additional factors as in the aforesaid discussion. Following this, as Fig. 16a–h illustrates, co-infection will become extinct.
Example 3
For the probabilistic COVID-19 model in the absence of TB (75), the white noise and the IC and the remaining arguments are the similar as in Example 1. Thus, we determine and the quasi-equilibrium In view of Theorem 11.
100
Graph
As a result, the following is the relevant P.D.F where
Graph
Consequently, the four marginal D. Fs of are as follows:
Graph
Finally, population concentrations oscillate according to the quasi-stable equilibrium , as seen in Fig. 17.
Conclusions
In this article, a deterministic-stochastic model is being suggested to investigate the potential transmission of the codynamics of COVID-19 and TB. Taking into account the deterministic fractional model and stochastic approach, we have provided the qualitative characteristics such as positivity and boundedness, reproduction number and their allied outcomes for co-infection model (2), global positive solution and unique erogdicity for the co-dynamics of (28). Besides that, applying the Khasminskii notion and a suitable Lyapunov function, the existence of a stationary distribution in model (28) was analytically verified. Additionally, an accurate representation of the P.D.F regarding a quasi-equilibrium point of the random-perturbed COVID-19 model constitutes one of this research's particularly noteworthy discoveries. In fact, it has been determined that the validity and strength of our numerical outcomes and modeled estimates have been provided in a piecewise fractional DEs context. Furthermore, the outcomes of this study shed a spotlight on the P.D.F and stationary distribution of the probabilistic multidimensional framework at its quasi-equilibrium point. Although the ABCFD, CFFD and CFD have been demonstrated to be efficient in documenting various interaction practices, we contend that their ability to accomplish this effectively may be hindered by the vastness of biological systems. It follows that oscillation might eliminate signals that are widely dispersed, despite leaving infectious diseases uncontrolled.
Graph: Figure 10 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (3) and (28) with the impacts of latent and active TB outbreaks using CFFD operator having FO ω=0.98, low intensities and ICs (1000, 100, 10, 1000, 10, 7, 5, 1).
Graph: Figure 11 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (4) and (28) with the impacts of latent and active TB outbreaks using CFFD operator having FO ω=0.98, low intensities and ICs (4000, 100, 40, 1000, 10, 7, 5, 4).
Graph: Figure 12 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (4) and (28) with the impacts of latent and active TB outbreaks using CFFD operator having FO ω=0.98, low intensities and ICs (5000, 100, 50, 1000, 10, 7, 5, 5).
Graph: Figure 13 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (5) and (28) with the impacts of latent and active TB outbreaks using ABCFD operator having FO ω=0.98, low intensities and ICs (1000, 100, 10, 1000, 10, 7, 5, 1).
Graph: Figure 14 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (5) and (28) with the impacts of latent and active TB outbreaks using ABCFD operator having FO ω=0.98, low intensities and ICs (4000, 100, 40, 1000, 10, 7, 5, 4).
Graph: Figure 15 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 models (2), (5) and (28) with the impacts of latent and active TB outbreaks using ABCFD operator having FO ω=0.98, low intensities and ICs (5000, 100, 50, 1000, 10, 7, 5, 5).
Graph: Figure 16 Time evaluation plots for deterministic-probabilistic co-infection TB-COVID-19 model (2) and (28) with the impacts of latent and active TB outbreaks with low intensities and ICs (5000, 100, 50, 1000, 10, 7, 5, 5) using CFD when ω=0.95.
Graph: Figure 17 Numerical modeling of the outcome (S(τ),EC(τ),IC(τ),R(τ)) in system (75) is displayed in the upper portion row. The P.D.Fs and marginal D.Fs of S,EC,IC and R are displayed in the lower portion row, respectively, with ℘ȷ=0.02(ȷ=1,...,8) and R0κ=2.83212>1.
Predicting how TB is propagated by population mobility and random disturbances was challenging until this research was conducted. The research advances our knowledge of why TB still exists around the globe. Regarding stochastic TB systems, including community propagation, it is challenging to define adequate requirements for infection eradication considering the restrictive nature of computational approaches. On the other hand, we also offer a simulation of the disease's disappearance. To determine the necessary requirements for TB's endurance and extermination, additional investigation needs to be performed.
Numerous fascinating and open-ended, high-dimensional models deserve further consideration. It is vital for inquiry into phenomena that are impacted by additional factors, such as neural networking with stochastic resonance or oscillatory spectrum disruption, while examining the unpredictable nature of this form of contention. Such studies may include certain specific but complex concepts, including evaluating the effects of Lévy and Poisson noise or Markov processes. These pertinent issues might be covered in the upcoming analysis.
Acknowledgements
The authors extend their appreciation to Taif University, Saudi Arabia, for supporting this work through project number (TU-DSPP-2024-127).
Author contributions
All authors read and approved the final manuscript.
Funding
This research was funded by Taif University, Saudi Arabia, Project No.(TU-DSPP-2024-127).
Data availability
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Competing interests
The authors declare no competing interests.
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
[
References
1
Wu Z, McGoogan JM. 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; 323: 1239-1242. 1:CAS:528:DC%2BB3cXnsFCqsrs%3D. 10.1001/jama.2020.2648. 32091533
2
Chen N. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study. The Lancet. 2020; 395: 507-513. 1:CAS:528:DC%2BB3cXhvFOmsb8%3D. 10.1016/S0140-6736(20)30211-7
3
Pai M, Kasaeva T, Swaminathan S. Covid-19's devastating effect on tuberculosis care—A path to recovery. N. Engl. J. Med. 2022; 386: 1490-1493. 1:CAS:528:DC%2BB38XhtFGrsLjF. 10.1056/NEJMp2118145. 34986295
4
Ong CW, Goletti D. Impact of the global COVID-19 outbreak on the management of other communicable diseases. Int. J. Tuberc. Lung Dis. 2020; 24: 547-548. 1:STN:280:DC%2BB38vlvFGksw%3D%3D. 10.5588/ijtld.20.0140. 32398214
5
Loveday M. Opportunities from a new disease for an old threat: Extending COVID-19 efforts to address tuberculosis in South Africa. S. Afr. Med. J. 2020; 110: 1160-1167. 1:CAS:528:DC%2BB3MXislWrtL4%3D. 10.7196/SAMJ.2020.v110i12.15126. 33403958
6
Pai, M. Tuberculosis and COVID-19: Fighting a Deadly Syndemic. Forbes. https://www.forbes.com/sites/madhukarpai/2020/09/26/tuberculosis-and-covid-19-fighting-a-deadly-syndemic/ (Accessed 12 June 2023) (2020).
7
Hopewell PC, Reichman LB, Castro KG. Parallels and mutual lessons in tuberculosis and COVID-19 transmission, prevention, and control. Emerg. Infect. Dis. 2021; 27: 681-686. 1:CAS:528:DC%2BB3MXhtlWgtbjL. 10.3201/eid2703.203456. 33213689. 7920655
8
Wu F. A new coronavirus associated with human respiratory disease in China. Nature. 2020; 579: 265-269. 2020Natur.579.265W. 1:CAS:528:DC%2BB3cXksFKlsLc%3D. 10.1038/s41586-020-2008-3. 32015508. 7094943
9
Hadler SC. Epidemic Intelligence Service investigations of respiratory illness, 1946–2005. Am. J. Epidemiol. 2011; 174: S36-S46. 10.1093/aje/kwr309. 22135392
Yates A. The transmission of Mycobacterium tuberculosis in high burden settings. Lancet Infect. Dis. 2016; 16: 227-238. 10.1016/S1473-3099(15)00499-5. 26867464
Jayaweera M, Perera H, Gunawardana B, Manatunge J. Transmission of COVID-19 virus by droplets and aerosols: A critical review on the unresolved dichotomy. Environ. Res. 2020; 188: 109819. 1:CAS:528:DC%2BB3cXht1KisLbO. 10.1016/j.envres.2020.109819. 32569870. 7293495
World Health Organization. Tracking SARS-CoV-2 Variants. https://www.who.int/health-topics/typhoid/tracking-SARS-CoV-2-variants (Accessed 12 June 2023) (WHO, 2022).
Kline SE, Hedemark LL, Davies SF. Outbreak of tuberculosis among regular patrons of a neighborhood bar. N. Engl. J. Med. 1995; 333: 222-227. 1:STN:280:DyaK2Mzht1Sgtw%3D%3D. 10.1056/NEJM199507273330404. 7791838
Jassat W. Risk factors for COVID-19-related in-hospital mortality in a high HIV and tuberculosis prevalence setting in South Africa: A cohort study. Lancet HIV. 2021; 8: e554-e567. 1:CAS:528:DC%2BB3sXis1Onu7k%3D. 10.1016/S2352-3018(21)00151-X. 34363789. 8336996
Western Cape Department of Health in collaboration with the National Institute for Communicable DiseasesAfrica South. Risk factors for coronavirus disease 2019 (COVID-19) death in a population cohort study from the Western Cape Province, South Africa. Clin. Infect. Dis. 2021; 73: e2005-e2015. 10.1093/cid/ciaa1198
Muñiz-Salazar R. Impact of COVID-19 on tuberculosis detection and treatment in Baja California, Mexico. Front. Public Health. 2022; 10: 921596. 10.3389/fpubh.2022.921596. 35942259. 9356343
Subbaraman R. Digital adherence technologies for the management of tuberculosis therapy: Mapping the landscape and research priorities. BMJ Glob. Health. 2018; 3: e001018. 10.1136/bmjgh-2018-001018. 30364330. 6195152
Atangana A. Extension of rate of change concept: From local to nonlocal operators with applications. Results Phys. 2021; 19: 103515. 10.1016/j.rinp.2020.1
Atangana A, Gomez-Aguilar JF. Fractional derivatives with no-index law property: Application to chaos and statistics. Chaos Solit. Fract. 2018; 114: 516-535. 2018CSF.114.516A. 3856675. 10.1016/j.chaos.2018.07.033
Jarad F, Abdeljawad T, Hammouch Z. On a class of ordinary differential equations in the frame of Atangana–Baleanu fractional derivative. Chaos Solit. Fract. 2018; 117: 16-20. 2018CSF.117.16J. 3872475. 10.1016/j.chaos.2018.10.006
Caputo M. Linear models of dissipation whose Q is almost frequency independent-II. Geophys. J. Int. 1967; 13: 529-539. 1967GeoJI.13.529C. 10.1111/j.1365-246X.1967.tb02303.x
Caputo M, Fabrizio M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 2015; 1: 73-85
Atangana A, Baleanu D. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016; 20: 763-769. 10.2298/TSCI160111018A
Abdeljawad T, Al-Mdallal QM. Discrete Mittag–Leffler kernel type fractional difference initial value problems and Gronwall's inequality. J. Comput. Appl. Math. 2018; 339: 218-230. 3787689. 10.1016/j.cam.2017.10.021
Kumar S, Chauhan RP, Aly AA, Momani S, Hadid S. A study on fractional HBV model through singular and non-singular derivatives. Eur. Phys. J. Spec. Top. 2022; 231: 1885-1904. 10.1140/epjs/s11734-022-00460-6. 35251498. 8889534
Mekkaoui T, Atangana A, Araz SÍ. Predictor-corrector for non-linear differential and integral equation with fractal–fractional operators. Eng. Comput. 2021; 37: 2359-2368. 10.1007/s00366-020-00948-6
Atangana A. & Araz S. I. A Successive Midpoint Method for Nonlinear Differential Equations with Classical and Caputo–Fabrizio Derivatives (2023).
Atangana A, Araz Sİ. Deterministic–stochastic modeling: A new direction in modeling real world problems with crossover effect. Math. Biosci. Eng. 2022; 19: 3526-3563. 4488341. 10.3934/mbe.2022163. 35341263
Atangana A, Araz SI. New concept in calculus: Piecewise differential and integral operators. Chaos Solit. Fract. 2021; 145: 110638. 4233320. 10.1016/j.chaos.2020.110638
Qurashi MA, Rashid S, Jarad F. A computational study of a stochastic fractal-fractional hepatitis B virus infection incorporating delayed immune reactions via the exponential decay. Math. Biosci. Eng. 2022; 19: 12950-12980. 4477005. 10.3934/mbe.2022605. 36654030
Rashid S, Iqbal MK, Alshehri AM, Ashraf R, Jarad F. A comprehensive analysis of the stochastic fractal–fractional tuberculosis model via Mittag–Leffler kernel and white noise. Results Phys. 2022; 39: 105764. 10.1016/j.rinp.2022.105764
Chu YM, Sultana S, Rashid S, Alharthi MS. Dynamical analysis of the stochastic COVID19 model using piecewise differential equation technique. Comput. Model. Eng. Sci. 2023; 137: 2427-2464. 10.32604/cmes.2023.028771
Rashid S, Jarad F, El-Marouf SAA, Elagan SK. Global dynamics of deterministic–stochastic dengue infection model including multi specific receptors via crossover effects. AIMS Math. 2022; 8: 6466-6503. 4550604. 10.3934/math.2023327
Gafel HS, Rashid S, Elagan SK. Novel codynamics of the HIV-1/HTLV-I model involving humoral immune response and cellular outbreak: A new approach to probability density functions and fractional operators. AIMS Math. 2023; 8: 28246-28279. 4663027. 10.3934/math.20231446
Silva CJ, Torres DFM. A TB-HIV/AIDS co-infection model and optimal control treatment. Discret. Contin. Dyn. Syst. 2015; 35: 4639-4663. 10.3934/dcds.2015.35.4639
Duan B, Zheng Z, Cao W. Spectral approximation methods and error estimates for Caputo fractional derivative with applications to initial-value problems. J. Comput. Phys. 2016; 319: 108-128. 2016JCoPh.319.108D. 3506937. 10.1016/j.jcp.2016.05.017
Khader MM, Gómez-Aguilar JF, Adel M. Numerical study for the fractional RL, RC, and RLC electrical circuits using Legendre pseudo-spectral method. Int. J. Circ. Theor. Appl. 2021; 49: 3266. 10.1002/cta.3103
Bashiri T, Vaezpour SM, Nieto JJ. Approximating solution of Fabrizio–Caputo Volterra's model for population growth in a closed system by homotopy analysis method. J. Funct. Spaces. 2018; 2018: 3152502. 3748657
Dokuyucu MA, Celik E, Bulut H, Baskonus HM. Cancer treatment model with the Caputo–Fabrizio fractional derivative. Eur. Phys. J. Plus. 2018; 133: 1-6
Alkahtani BST. Chua's circuit model with Atangana–Baleanu derivative with fractional-order. Chaos Solit. Fract. 2016; 89: 547-551. 2016CSF.89.547A. 10.1016/j.chaos.2016.03.020
Algahtani OJJ. Comparing the Atangana–Baleanu and Caputo–Fabrizio derivative with fractional-order: Allen Cahn model. Chaos. Solit. Fract. 2016; 89: 552-559. 2016CSF.89.552A. 3507065. 10.1016/j.chaos.2016.03.026
Atangana A, Koca I. On the new fractional derivative and application to nonlinear Baggs and Freedman model. J. Nonlinear Sci. Appl. 2016; 9: 2467-2480. 3470409. 10.22436/jnsa.009.05.46
Kilbas AA, Srivastava HM, Trujillo JJ. Theory and Application of Fractional Differential Equations. 2006; Elsevier
Van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002; 180: 29-48. 1950747. 10.1016/S0025-5564(02)00108-6. 12387915
Martcheva M. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002; 180: 29. 1950747. 10.1016/S0025-5564(02)00108-6
Castillo-Chavez C, Song B. Dynamical models of tuberculosis and their applications. Math. Biosci. Eng. 2004; 1: 361. 2130673. 10.3934/mbe.2004.1.361. 20369977
Lipsitch M, Murray MB. Multiple equilibria: Tuberculosis transmission require unrealistic assumptions. Theor. Popul. Biol. 2003; 63: 169-170. 10.1016/S0040-5809(02)00037-0. 12615499
Mao X. Stochastic Differential Equations and Applications. 1997; Horwood Publishing
Khasminskii R. Stochastic Stability of Differential Equations. 1980; Sijthoff and Noordhoff. 10.1007/978-94-009-9121-7
Gardiner CW. Handbook of Stochastic Methods. 1983; Springer. 10.1007/978-3-662-02377-8
Ma Z, Zhou Y. Qualitative Analysis and Stability Method of Ordinary Differential Equations. 2001; Science Press
Roozen H. An asymptotic solution to a two-dimensional exit problem arising in population dynamics. SIAM J. Appl. Math. 1989; 49: 1793-1810. 1025962. 10.1137/0149110
Sy KTL, Haw NJL, Uy J. Previous and active tuberculosis increases risk of death and prolongs recovery in patients with COVID-19. Infect. Dis. 2020; 52: 902-907. 10.1080/23744235.2020.1806353
Boulle A. Risk factors for coronavirus disease 2019 (COVID-19) death in a population cohort study from the Western Cape Province, South Africa. Clin. Infect. Dis. 2020; 73: e2005-e2015
World Health Organization. WHO Information Note: COVID19: Considerations for Tuberculosis. 2021; World Health Organization
]
By Saima Rashid; Sher Zaman Hamidi; Saima Akram; Muhammad Aon Raza; S. K. Elagan and Beida Mohsen Tami Alsubei
Reported by Author; Author; Author; Author; Author; Author