Tuberculosis (TB) is a deadly contagious disease that affects vital organs of the body, especially the lungs. Although the disease is preventable, there are still concerns about its continued spread. Without effective prevention or appropriate treatment, TB infection can be fatal to humans. This paper presents a fractional-order TB disease (FTBD) model to analyze TB dynamics and a new optimization method to solve it. The method is based on the basis functions of generalized Laguerre polynomials (GLPs) and some new operational matrices of derivatives in the Caputo sense. Finding the optimal solution to the FTBD model is reduced to solving a system of nonlinear algebraic equations with the aid of GLPs using the Lagrange multipliers method. A numerical simulation is also carried out to determine the impact of the presented method on the susceptible, exposed, infected without treatment, infected with treatment, and recovered cases in the population.