key: cord-0689944-gmve512q
authors: Panchal, Jitendra; Acharya, Falguni; Joshi, Kanan
title: A noninteger order SEITR dynamical model for TB
date: 2022-03-26
journal: Adv Cont Discr Mod
DOI: 10.1186/s13662-022-03700-0
sha: 87682bade1621a4a3d349eb1287e5a94b0083903
doc_id: 689944
cord_uid: gmve512q

This research paper designs the noninteger order SEITR dynamical model in the Caputo sense for tuberculosis. The authors of the article have classified the infection compartment into four different compartments such as newly infected unrecognized individuals, diagnosed patients, highly infected patients, and patients with delays in treatment which provide better detail of the TB infection dynamic. We estimate the model parameters using the least square curve fitting and demonstrate that the proposed model provides a good fit to tuberculosis confirmed cases of India from the year 2000 to 2020. Further, we compute the basic reproduction number as [Formula: see text] of the model using the next-generation matrix method and the model equilibria. The existence and uniqueness of the approximate solution for the SEITR model is validated using the generalized Adams–Bashforth–Moulton method. The graphical representation of the fractional order model is given to validate the result using the numerical simulation. We conclude that the fractional order model is more realistic than the classical integer order model and provide more detailed information about the real data of the TB disease dynamics.

Tuberculosis is a disease instigated by Mycobacterium tuberculosis bacteria (MTB), which commonly affects the lungs; however, it can also harm other parts of the body. It is an infectious disease; around ten percent of the latent infection develops into a highly infected disease and results in death. The transmission of TB from one person to another occurs through air droplets when an infectious person sneezes or coughs [2] . Generally, indications include fever, feeling of coldness, sweats at night, lost appetite, loss of weight, and weakness [2] . Tuberculosis has existed in the world for ages, which is apparent from the Egyptian mummy stored in the British Museum that reveals the infection in its spine due to tuberculosis; however, in the nineteenth and early twentieth centuries, the threat of prevalence in public increased, especially among urban [2] . Visibly the number of cases eventually degraded since the year 2000 [6] . By the year 2018, one-fourth population of the world was assumed to have a latent infection of tuberculosis, and there were around ten million active cases with one and half million deaths in the same year. This history ranks TB as the number one infectious disease causing death [6] . Eighty percent of the positive tested patients are from Asia and Africa, whereas only five to ten percent are tested positive in the USA by the tuberculin test [28] . For five decades, rampant and efficacious treatment for actively infected and inertly infected people has been made accessible. Medicines like Streptomycin combined with Pyrazinamide have been used to cure TB, whereas the two most impactful medications believed to combat Mycobacterium tuberculosis are Isoniazid and Rifampin [39] .

In the wake of the COVID-19 pandemic, India has observed a drop in detecting new TB cases by 1.3 million from 2019 to 2020 [1] . Additionally, according to the global 2021 tuberculosis report released by WHO, the number of deaths due to TB has increased thereafter, ranking India among the top countries contributing towards the deduction of diagnosis of TB during the COVID-19 crisis [1]. Approximately 4.1 million undiagnosed individuals are suffering from TB presently. Due to lockdown, the number of people visiting for treatments has also fallen notably [1] . India shared 34% of the total predicted 1.48 million deaths occurred by TB globally, whereas it surged to 13% within the country as compared to 2019 [1] .

Mathematical modelling for infectious diseases has been developed and studied thoroughly by many researchers for decades now. In comparison to the traditional integerorder models, the fractional-order models offer more precise and in-depth knowledge related to the complex patterns of numerous diseases owing to their inherited properties and explanation of memory [10, 27, 38] . For instance, Khan et al. [25] developed a novel fractional model for TB, [37] analysed the tuberculosis-HIV system, and [20] discussed a fractional model for fever named dengue. Furthermore, [17] developed a fractional order discrete differential model with a time delay to study the effects of endogenetic and exogenic recurrence in the qualitative behavior of tuberculosis. Moreover, Yang explained the impacts of being dormant for a long time and numerous contagions in the changing patterns of tuberculosis [51] . Bowong [14] examined the slow-moving and speedy evolution of the TB model that includes slow and fast development. The conclusion was if the number of elementary reproductions is smaller than one, the steady-state of disease-free is constant in contradiction to when it is larger than one. The optimal control problem for fractional TB disease model incorporating the effect of diabetes and stress was studied in [42] ; besides, its analysis using Atangana-Baleanu derivative was carried out in [10, 45] .

The evolution of TB was studied mathematically using epidemiological modeling [13, 21] for a long time; however, Waaler bagged the initial effort to model the disease epidemiology of TB [48] . Later, Kermack and McKendrick bifurcated various stages of disease into diverse compartments, which were originated by [24] and expounded by Baile [33] , Anderson and May [9] . In [46] , the C-F fractional-order derivative is applied to analyze the TB model. Agusto and Cook analyzed a deterministic model to study transmission after separation, cases that forgot to turn up for checkups, and cases with drug resistance [7] . Melsew, Adekunle, and Cheng worked to prove the presence of widespread heterogeneity among infections of people with active TB [31] .

Furthermore, researchers developed mathematical models to study the treatments and ways of cure like vaccination. The mathematical model in [32] demonstrates the efficiency of the vaccination to cure an infection of TB. Egonmwan et al. [36] established a newfangled TB model to inspect the effect of diagnosis and treatment for the infectious mass of TB. The study in [12] shows that diagnosis, proper vaccination drives, and treatments can restrain uncertain aspects of TB; however [40] mentions that even deficient vaccination of tuberculosis plays a role in inhibiting the spread of the disease. Robert studied a TB system with the consequence of reversion of infection [49] .

Some authors have also analyzed the relation between TB and coronavirus looking forward to the recent pandemic. The interpretations given by Marimuthu, Nagappa, Sharma, Basu, and Kishore Chopra. suggest a considerable hike of COVID infection among the patients with TB [30] . Also, Iyengar and Jain predicted that due to distributed attention during COVID-19 pandemic the tracking of TB patients and the vital activities for TB may be affected [23] . Besides, the work of Cilloni and Fu suggested to act aptly towards combating TB to avoid the adverse situation as a result of COVID-19 pandemic [18] .

The integer-order derivative is local in nature, which means that the integer-order derivative is useful in understanding nature of the function at one point only and in the neighborhood of that point. The fractional-order derivatives can overcome the challenges imposed by traditional integer-order derivatives, and it also covers the aspects of integer order derivatives. There are various fractional order derivatives, the most notable of which are the Riemann-Liouville fractional derivative and the Caputo fractional derivative in applications. The Caputo fractional derivative is more suitable than the Riemann-Liouville fractional derivative as it requires initial conditions containing limit values at t = 0. The former gives more accurate and precise information regarding the complexities of numerous diseases. Fractional order models are better than conventional integer-order models because of hereditary properties and memory description [42, 43] . However, in reality, any phenomena of nature can be expressed better by a biological system of nonintegerorder rather than integer-order models [17] . The fractional-order derivatives can apprehend nonlocal relationships in time and space, and they also offer more degree of freedom and a precise representation of nonlinear phenomena [17] . Saeed Ahmad et al. stated the use of the fusion method to develop a series solution for the analytical solution differential equations of fractional order which includes nonsingular derivative [8] . For TB disease, a three-strained noninteger order model is studied in [41] . A fractional order model is also used by [16] to analyze the dynamics of TB. Weronika used Lypanove theory to check whether a noninteger order system for tuberculosis is globally stable or not [50] . Moreover, [44] modeled a fraction order TB model using Caputo derivative from the factual data of Khyber Pakhtunkhwa.

This research paper establishes a fractional-order mathematical model in the Caputo sense to analyze the dynamical behavior of the spread of tuberculosis infection. The study represents various compartments of TB transmission infection for the dormant individuals (not diagnosed as infected) and individuals identified as infected with TB. The following system considers more intermediate compartments as diagnosis, high infection, and cure, to study different treatment and recovery aspects. Basic definitions are given in Sect. 2. Section 3 contains model formulation. The basic reproduction number and endemic equilibria are given in Sect. 4 and Sect. 5. In Sect. 6, parameter estimation and model fitting with the real data of the reported cases in India during the year 2000 to 2020 are given. The existence and uniqueness of the solution for the model is proved in Sect. 7. The numerical results of the approximate solution for the fractional-order model in the Caputo sense are given in Sect. 8. Numerical simulations with different scenarios for the model are given in Sect. 9, which are useful for the Ministry of Health and Family Welfare, Government of India to implement an action plan to control TB infection in India. Finally, conclusion is provided in Sect. 10.

This section discusses the prerequisites and elementary notions used throughout this paper.

Definition 2.1 ( [35, 36] ) For an integrable function f, the Caputo derivative of fractional order α ∈ (0, 1) is given by

Also, the fractional integral of order α with Re(α) > 0 is given by

Definition 2.2 ( [11, 15] ) For f ∈ H 1 (c, d) and d > c, the C-F derivative of fractional order α ∈ (0, 1) for f is given by

where t ≥ 0, M(α) is a normalization function that depends on α and M(0) = M(1) = 1. If f / ∈ H 1 (c, d) and 0 < α < 1, this derivative for f ∈ L 1 (-∞, d) is given by

Also, the C-F fractional integral is presented by

Definition 2.3 ( [35, 36] ) The L. T. of Caputo fractional differential operator of order α is given by

which can also be written as

The compartmental model for the transmission dynamics of TB is represented graphically in Fig. 1 . The list of state variables and their description with values are given in Table 1 . Rate at which I(t) join I 1 (t) or I 2 (t) compartment 0.6102 Fitted k 2 Rate at which I(t) is classified as I 1 (t) or I 2 (t) compartment 0.4213 Fitted a 7 Rate at which E 1 (t) join I(t) compartment 0.5654 Fitted a 8 Rate at which

The population is separated into nine compartments which are classified based on the infection status. Here, λ is the contact rate. At time t, S(t) represents the number of susceptible individuals at risk of getting infected. a 1 λ is the rate of the susceptible individual classified as infected and joining I(t) class. E(t) is the number of the new latently infected individual being in sufficient contact with an infected person but not contagious. a 2 σ is the rate at which the new latently infected individual joins the diagnosed latently infected E 1 (t) and others associated with the undiagnosed latently infected E 2 (t) class based on the TB test result. The undiagnosed actively infected people I(t) are those who become contagious and can spread infection to others. a 7 λ and a 4 λ + k 1 are the rates at which individuals from E 1 (t) and E 2 (t) become actively infected. The individual from I(t) class diagnosed as actively infected with rate a 6 k 2 joins the class I 1 (t), which represents the individuals diagnosed as actively infected with speedy treatment. The people diagnosed as actively infected but delayed in treatment are joining I 2 (t) at rate (1a 6 )k 2 . r 0 , r 1 , and r 2 are the rates at which infected individuals from E 1 (t), I 1 (t), and I 2 (t) started treatment and moved to treatment compartment T(t). R(t) represents the number of recovered or removed individuals as well as those who stopped treatments after partial recovery. A person who has completed treatment and recovered moves to R(t) at rate r 3 . In MDR and XDR compartments of TB, an individual from T(t) class is moved to I(t) compartment at rate a 1 a 5 λ. The recovered people can become susceptible at the rate a 8 . The integer-order model is unable to represent the dynamics of the real-world problem between two points. The fractional-order mathematical model is a more reliable and powerful tool for minimizing error created by neglected parameters in modeling [27, 34] . The fractional order dynamical system of the TB infection model using Caputo derivative of order α ∈ (0, 1) is given by

with the nonnegative initial condition

Therefore, we can observe that the solution of system (1) is nonnegative and bounded [29] .

The overall dynamics of the population attained by adding all six equations of model (1) is as follows:

System (1) is moderated by substituting the Caputo fractional time-derivative. In this moderate system, the dimension of the system will not remain the same for the right and left sides. We use -an auxiliary parameter to resolve this problem. As per the discussion, the fractional-order model for t > 0 and α ∈ (0, 1) with the same initial conditions is given by

The feasible region for model (1) is given by

We prove that the closed set is the feasible region of system (3).

Proof We add all the terms in system (3) to obtain the overall population in the fractional order, i.e., α-

To obtain the population size, we use the Laplace transform as follows:

where N(0) is the initial population size. After simplifying we get

The following table represents the parametric values used for model (1) fitted or estimated using data available and reported cases of TB infection in India.

Consequently, the closed set is positive invariant concerning fractional-order model (3). The solutions of system (1) are known as equilibrium points, and they have two equilibrium points as follows:

(1) Infection-free equilibrium points;

(2) Endemic equilibrium points, which are discussed in the following sections.

The basic reproduction number 0 (also known as threshold parameter) is used to analyze the average number of secondary infected users rising from an average primary infected users in entirely susceptible population for the stability of system (1). 0 is obtained by using the next generation matrix [47] . The new infectious rates represented using matrix F, and other transferred rates within compartments by matrix V are given as follows: 

The basic reproduction number of system (3) is given by the spectral radius

λa 5 (k 2 + μ)(r 2 + μ)(a 3 r 0 σ (1a 2 ) + r 0 a 2 σ (a 4 λ + k 1 + a 3 + μ)) + λa 5 k 2 (r 1 a 6 (r 2 + μ) + r 2 (1a 6 ))(σ (1a 2 ))((a 4 λ + μ)(a 7 λ + r 0 + μ) + a 3 a 7 λ) + (a 2 a 7 λσ )(a 4 λ + k 1 + a 3 + μ) ⎞ ⎟ ⎟ ⎟ ⎠ (k 2 + μ)(r 2 + μ)(σ + μ)(a 7 λ + r 0 + μ)(a 4 λ + k 1 + a 3 + μ)(a 1 a 5 λ + r 3 + μ) .

The basic reproduction number is obtained as 0 = 1.7307 using the parameter values given in Table 1 for India.

If 0 > 1 then the disease will persist in the community since one diseased individual will infect more than one susceptible individual on average. This is possible if (k 2 + μ)(r 2 + μ)(σ + μ)(a 7 λ + r 0 + μ)(a 4 λ + k 1 + a 3 + μ)(a 1 a 5 λ + r 3 + μ) < 1, i.e., the values of parameters k 1 , k 2 , a 1 , a 4 , a 5 ,and a 7 decrease, which means that the latently or actively infected individual remains in the population and may spread the infection to others. Also, the decrement of the parameters r 2 , r 3 , σ , r 0 , and a 3 means that the infected individuals and relapse cases are not taking proper treatment or left treatments before recovery and may propagate infection to others.

If 0 < 1 then on average one diseased individual can only infect one other person, and the disease will eventually die out. Proof To obtain the TB-infection free equilibrium at a point E 0 , the Jacobian matrix [19, 22] is as follows:

Define s(J) = max{ λ : λ is an eigenvalue of J}, where s(J) is the simple eigenvalue of matrix J with positive eigenvector, then we have 0 < 1 ⇔ s(J) < 0, see [47] for more details. Therefore,

-(a 4 λ + k 1 + a 3 + μ), -(k 2 + μ), -(r 1 + μ), -(r 2 + μ), -(a 1 a 5 λ + r 3 + μ), -(a 8 + μ)

< 0.

Hence, if 0 < 1, then s(J) < 0, E 0 is locally asymptotically instable.

We use the following equations to determine the equilibrium points for the fractionalorder model (3):

The algebraic solution of the equation provides equilibrium points of the system, and if 0 > 1, then system (3) has a positive endemic equilibrium E * 1 = (S * , E * , E * 1 , E * 2 , I * , I * 1 , I * 2 , T * , R * ) and the Jacobian is given as

where S * = ( +a 8 R * ) (λ+μ) , E * = λ(1-a 1 )(S * +a 5 T * ) (σ +μ)

, E 1 * = 1 (a 7 λ+r 0 +μ) (a 2 σ E * +a 3 E 2 * ), E 2 * = σ (1a 2 ) (a 4 λ + k 1 + a 3 + μ) E * , I * = 1 (k 2 + μ) a 1 λ(S * +a 5 T * ) + a 7 λE 1 * +(a 4 λ + k 1 )E 2 * , I 1 * = a 6 k 2 r 1 + μ I * , I 2 * = (1a 6 )k 2 r 2 + μ I * , T * = (a 8 + μ) r 3 × R * or (r 0 E 1 * +r 1 I 1 * +r 2 I 2 * ) (a 1 a 5 λ + r 3 + μ) R * = ((r 0 ) a 2 (a 4 λ + k 1 + μ) + a 3 r 0 λσ (r 2 + μ) k 2 r 1 a 6 + (1a 6 )r 2 (a 7 λ + r 0 + μ)

× (1a 2 )(a 4 λ + k 1 + a 3 + μ) × ( 0 -1) / r 1 a 6 (r 2 + μ)(σ + μ)(a 7 λ + r 0 + μ)λσ (a 4 λ + k 1 + a 3 + μ)(1a 1 )(1a 2 ) × 0 .

From the above equations, we observe that the endemic equilibrium E * 1 = (S * , E * , E * 1 , E * 2 , I * , I * 1 , I * 2 , T * , R * ) only exists if 0 > 0. Figure 3 represent the model fitting to the real data for various fractional order.

From Fig. 4 , we can see the long-term trend in the number of reported cases of tuberculosis infection in India, which demonstrates a significant increase in the number of recorded cases and gives a terrible sign to the Ministry of Health and Family Welfare, Government of India. The model parameters given in Table 1 are calculated from the reported data and will be employed in numerical simulations. 

In this section, we use the generalized Adams-Bashforth-Moulton technique given in [35] to prove the uniqueness of the solution for system (3) as follows:

By using Lemma 3.1, the system is given as

In the ensuing theorem the kernels Q i , i = 1, 2, 3, 4, 5, 6, 7, 8, 9, satisfy the Lipschitz condition and contraction.

Proof For S and S * we have

Thus, for Q 1 the Lipchitz condition is obtained, and if 0 ≤ λ + μ < 1 then Q 1 is a contraction.

Similarly, the Lipschitz condition for Q i , i = 2, 3, 4, 5, 6, 7, 8, 9, is given as follows:

Q 6 (t, I 1 ) -Q 6 (t, I 1 * ) ≤ d 6 I 1 (t) -I 1 * (t) ,

where d 2 = σ + μ, d 3 = a 7 λ + r 0 + μ, d 4 = a 4 λ + k 1 + a 3 + μ, d 5 = k 2 + μ, d 6 = r 1 + μ, d 7 = r 2 + μ, d 8 = a 1 a 5 λ +r 3 +μ,and d 9 = a 8 +μ are bounded functions, if 0 ≤ d i < 1, i = 2, 3, 4, 5, 6, 7, 8, 9, then Q i , i = 2, 3, 4, 5, 6, 7, 8, 9, are contractions. According to system (5), consider the following recursive forms:

P 5n (t) = I n (t) -I n-1 (0) = 1-α α t 0 Q 5 (τ , I n-1 ) -Q 5 (τ , I n-2 ) (tτ ) α-1 dτ , P 6n (t) = (I 1 ) n (t) -(I 1 ) n-1 (0) = 1-α α t 0 Q 6 τ , (I 1 ) n-1 -Q 6 τ , (I 1 ) n-2 (tτ ) α-1 dτ , P 7n (t) = (I 2 ) n (t) -(I 2 ) n-1 (0) = 1-α α t 0 Q 7 τ , (I 2 ) n-1 -Q 7 τ , (I 2 ) n-2 (tτ ) α-1 dτ ,

(E 2 ) 0 (t) = (E 2 )(0), I 0 (t) = I(0), (I 1 ) 0 (t) = (I 1 )(0), (I 2 ) 0 (t) = (I 2 )(0), T 0 (t) = T(0) and R 0 (t) = R(0).

We take the norm of the first equation in the above system, then

With Lipchitz condition (6), we have

As a similar way, we obtained

Thus, we can write that

T n (t) = n j=1 P 8j (t), R n (t) = n j=1 P 9j (t).

The existence of a solution is given in the next theorem.

Proof From the recursive technique and Eq. (7) and Eq. (8) we conclude that

Thus, the system has a continuous solution. To prove that the above functions construct a solution for model (2), we assume that

(I 1 )(t) -(I 1 )(0) = (I 1 ) n (t) -W 6n (t),

By repeating the method, we obtain

At t 1 , we get

As n approaches to ∞, this implies W 1n (t) → 0. Similarly, we can obtain W in (t) → 0, i = 2, 3, 4, 5, 6, 7, 8, 9. Hence the theorem is proved.

To prove the uniqueness of the solution, consider that the system has another solution such as S (t), E (t), E 1 (t), E 2 (t), I (t), I 1 (t), I 2 (t), T (t), and R (t), then we have

We take the norm of this equation

It follows from Lipschitz condition (3) that

(I 1 ) n+1 = (I 1 ) 0 + h α 1-α (α + 2) a 6 k 2 I p n+1 -(r 1 + μ)(I 1 ) p n+1

a i,n+1 a 6 k 2 I i -(r 1 + μ)(I 1 ) i ,

a i,n+1 r 0 (E 1 ) i + r 1 (I 1 ) i + r 2 (I 2 ) i -(a 1 a 5 λ + r 3 + μ)T i ,

(1a 6 )k 2 I i -(r 2 + μ)(I 2 ) i ,

φ i,n+1 r 0 (E 1 ) i + r 1 (I 1 ) i + r 2 (I 2 ) i -(a 1 a 5 λ + r 3 + μ)T i ,

i ≤ n, and j = 1, 2, 3, 4, 5.

The following are the graphical simulations for the fractional orders α = 0.5, 0.6, 0.7, 0.8, 0.9, 0.99, 1. In many countries, the BCG vaccine is given to a newly born baby. It can be observed from the figure that the number of susceptible individuals decreases over time as people are immunized through BCG vaccination. In Fig. 5 , we can observe that the number of susceptible individuals decreases over time.

In Fig. 6 , due to vaccination, the number of exposed individuals to TB also decreases over time. However, many people are still reported as TB-infected with different compartments. As NLI individuals do not have any symptoms and feel sick, they are not spreaders of TB bacteria to others; but can be identified using skin or blood tests. We can also observe from Fig. 3 that initially the number of NLI individuals increases, and then due to the decrease in the number of susceptible individuals concerning time, the number of NLI individuals decreases slowly over time.

Because of awareness regarding TB infection, people are taking precautions in the initial compartments, and if diagnosed as TB-infected, the person becomes diagnosed latently infected (DLI) and then starts taking treatment. In Fig. 7 , we can observe that the number of DLI individuals is increasing slowly over time as NLI individuals are diagnosed positive. In Fig. 8 , as discussed earlier, from the symptoms of NLI-TB we can observe that the number of ULI individuals is increasing over time. As the duration of the treatment is long, some people stop taking treatment as they feel better and hence do not recover completely, which may convert them into actively infected and spread the infection to other people.

In Fig. 9 , due to the awareness about TB disease, infected people start taking medication immediately. We can observe that the number of DAI with prompt treatment increases over the time after diagnosed as positively infected.

The individual who is not diagnosed in the initial compartment may become actively infected and also spread the infection to others. People having symptoms of TB infection but not diagnosed as TB-infected and people diagnosed as positive but not starting treatment may move to treatment class after reaching more severe conditions. In Fig. 10 , we can observe that the infection of some people who are diagnosed as infected becomes quite severe, and then they start taking medication, due to which the number of DAI individuals with delay in treatment increases after some time. Figure 11 represents the number of individuals under the treatment after diagnosed as infected. This includes the number of individuals in DLI, DLI with prompt treatment and DLI with delay treatment. Perhaps, during the treatment some people may become drugresistant and can be classified as MDR and XDR. As in all compartments, people are getting TB infection, the number of infected individuals increases initially and then decreases over time. Figure 12 represents the number of recovered/removed individuals. The duration of medication or treatment takes longer time for TB. So, after getting partial recovery, some people stop medication or treatment. The individual recovered from this stage may become susceptible for the infection and may have chance of getting infected again. From Fig. 13 we observe that if the contact rate λ of an infected individual is considered as lower than the recovery rate r 3 , the number of infected people rises at first, then drops significantly. Figure 14 represents the behavior of the infected individual when the contact rate λ, treatment rates r 1 and r 2 are considered as 1 and the recovery rate as r 3 = 0.1, then we can observe a rapid increase in the number of infected individuals and increased burden on treatment class. Moreover, the recovery rate is lower than the contact rate and the treatment rate, which decreases the number of infected individuals due to the mortality through infection.

For Fig. 15 , we have considered the values of the contact rate λ = 1, the rate at which an individual from the treatment class joins the infected class a 1 = a 5 = 1, the treatment rates r 1 and r 2 are considered as 0.5 and the recovery rate as r 3 = 0.1. As the recovery rate is 

In this paper, a dynamical fractional-order mathematical model in the Caputo sense is proposed for the behavioral investigation of TB infection in India. Using the nonlinear least square algorithm, we estimated the parameters used in this model using TB-infected cases reported in India from the year 2000 to 2020, which shows that the model curve provides a good fit to the real data. Moreover, the future trend of the model curve provides indication to the decision maker or policy maker for better devising disease prevention and control measures. Further, we present the endemic equilibria, locally asymptotically stable, and the TB infection-free equilibrium in terms of the basic reproduction number. The value of the basic reproduction is obtained as 0 = 1.7307 demonstrating the importance of ap- propriate control and prevention strategies. The existence and uniqueness of the approximate solution for the model is derived using the generalized Adams-Bashforth-Moulton method in the numerical results. The graphical representation is provided to demonstrate the flow of all compartments for α = 0.5, 0.6, 0.7, 0.8, 0.9, 0.99, 1. The integer order model represents the rate of infection at order 1, but fails to provide any information between order 0 and 1, which is important to identify the initial behavior of the infection. However, our model represents the fractional rate of infection at different order between 0 and 1 to capture the significant information and analyze the intricate dynamics of the infection of tuberculosis for better apprehension, which can be observed in all the above graphical representation. Tuberculosis infection is curable, but still the infection is transmitting rapidly between the community. Also, from the real data of TB infection in India [3] we can observe that the confirmed cases of TB infection decreased during the COVID-19 lockdown in India from March 2020 to December 2020. During the lockdown, people were using masks, sanitizer, maintained social distancing, and followed COVID-19 guidelines; and this action plan played a vital role inpreventing TB infection among the population of India. Further, we can extend this model by considering the regimen changed compartment or the exempted compartment for non-evaluate individual or both separately.

Mathematical Model of MDR-TB and XDR-TB with Isolation and Lost to Follow-up

Mathematical analysis of tuberculosis control model using nonsingular kernel type Caputo derivative

Infectious Diseases of Human: Dynamics and Control

New fractional derivatives with non-local and nonsingular kernel: theory and application to heat transfer model

New numerical approach for fractional differential equations

A mathematical model of the tuberculosis epidemic

The intrinsic transmission dynamics of tuberculosis epidemics

Global analysis of a dynamical model for transmission of tuberculosis with a general contact rate

A new definition of fractional derivative without singular kernel

Non-integer order analysis of the impact of diabetes and resistant strains in a model for TB infection

A fractional-order model with time delay for tuberculosis with endogenous reactivation and exogenous reinfections

The potential impact of the COVID-19 pandemic on the tuberculosis epidemic a modelling analysis

On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations

A fractional calculus-based model for the simulation of an outbreak of Dengue fever

Using a mathematical model to evaluate the efficacy of TB control measures

Backward bifurcations in Dengue transmission dynamics

Tuberculosis and COVID-19 in India-double trouble

A contribution to the mathematical theory of epidemics

A new fractional model for tuberculosis with relapse via Atangana-Baleanu derivative

Modeling the transmission dynamics of tuberculosis in Khyber Pakhtunkhwa Pakistan

A new analysis of the Fornberg-Whitham equation pertaining to a fractional derivative with Mittag-Leffler-type kernel

Robbins Basic Pathology, 8th edn

Modeling and forecasting of COVID-19 spreading by delayed stochastic differential equations

COVID-19 and tuberculosis: a mathematical model based forecasting in Delhi

Heterogeneous infectiousness in mathematical models of tuberculosis: a systematic review

Mathematical model on pulmonary and multidrug-resistant tuberculosis patients with vaccination

Bailey: the mathematical theory of epidemics

Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications

SEIR epidemic model for COVID-19 transmission by Caputo derivative of fractional order

Fractional Integrals and Derivatives: Theory and Applications

A TB-HIV/AIDS coinfection model and optimal control treatment

On the analysis of chemical kinetics system pertaining to a fractional derivative with Mittag-Leffler type kernel

Tuberculosis models with fast and slow dynamics: the role of close and casual contacts

An SVEIRE model of tuberculosis to assess the effect of an imperfect vaccine and other exogenous factors

Comparative study for multi-strain tuberculosis (TB) model of fractional order

Optimal control for a fractional tuberculosis infection model including the impact of diabetes and resistant strains

A mathematical model for COVID-19 transmission by using the Caputo fractional derivative

A fractional model for the dynamics of TB virus

A fractional model for the dynamics of tuberculosis (TB) using Atangana-Baleanu derivative

A fractional model for the dynamics of tuberculosis infection using Caputo-Fabrizio derivative

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

The use of mathematical models in the study of the epidemiology of tuberculosis

Mathematical models of tuberculosis reactivation and relapse

Uniform asymptotic stability of a fractional tuberculosis model

Assessing the effects of multiple infections and long latency in the dynamics of tuberculosis

Not applicable.

Not applicable.

The data used in this manuscript are available on the following websites: (1) https://tbfacts.org/tb-statistics-india/;(2) https://tbcindia.gov.in; (3) https://pib.gov.in/PressReleseDetail.aspx?PRID=1606209; (4) https://www.worldometers. info/world-population/india-population/#:~:text=India%202020%20population%20is%20estimated,(and%20 dependencies)%20by%20population; (5) https://www.thehindu.com/data/data-tb-deaths-on-a-seven-year-high-as-casenotifications-and-outpatient-visits-dipped-in-2020/article37098076.ece.

The authors declare that they have no competing interests.

All authors contributed equally in writing of this manuscript. All authors read and approved the final manuscript.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 

Theorem 7.3 The solution of model (3) is unique if the following condition holds:Proof Suppose that condition (9) holdsThen S(t) -S (t) = 0. Therefore, we get S(t) = S (t). Likewise, the same equality can be shown for E, E 1 , E 2 , I, I 1 , I 2 , T, and R.

In this section, we present the numerical results for the TB-model (3). We used the Adams-Bashforth-Moulton scheme [11] . Set h = T N , t n = nh, n = 0, 1, 2, . . . , N ∈ Z + , we can write system (3) as follows: + h α 1-α (α + 2) n i=0 a i,n+1 a 1 λ(S i + a 5 T i ) + a 7 λ(E 1 ) i + (a 4 λ + k 1 )(E 2 ) i -(k 2 + μ)I i ,