Abstract
The present study aims to investigate a new fractional model describing the dynamical behaviour of the tuberculosis infection. In this new formulation, we use a recently introduced fractional operator with Mittag-Leffler nonsingular kernel. To solve and simulate the proposed model, a new and efficient numerical method is developed based on the product-integration rule. Simulation results are provided and some discussions are given to verify the theoretical analysis. The results indicate that employing the nonsingular operator can extract the hidden aspects of the model under consideration while these features are invisible when we use the ordinary time-derivatives. Therefore, the non-integer calculus supplies more flexible models describing the asymptotic behaviours of the real-world phenomena and helps us to better understand their complex dynamics.