Abstract
In this study, we formulate a noninteger-order mathematical model via the Caputo operator for the transmission dynamics of the bacterial disease tuberculosis (TB) in Khyber Pakhtunkhwa (KP), Pakistan. The number of confirmed cases from 2002 to 2017 is considered as incidence data for the estimation of parameters or to parameterize the model and analysis. The positivity and boundedness of the model solution are derived. For the dynamics of the tuberculosis model, we find the equilibrium points and the basic reproduction number. The proposed model is locally and globally stable at disease-free equilibrium, if the reproduction number Script capital R0<1. Furthermore, to examine the behavior of the various parameters and different values of fractional-order derivative graphically, the most common iterative scheme based on fundamental theorem and Lagrange interpolation polynomial is implemented. From the numerical result, it is observed that the contact rate and treatment rate have a great impact on curtailing the tuberculosis disease. Furthermore, proper treatment is a key factor in reducing the TB transmission and prevalence. Also, the results are more precise for lower fractional order. The results from various numerical plots show that the fractional model gives more insights into the disease dynamics and on how to curtail the disease spread.