In this article, we investigate a mathematical fractional model of tuberculosis that takes into account vaccination as a possible way to treat the disease. We use an in-host tuberculosis fractional model that shows how Macrophages and Mycobacterium tuberculosis interact to knowledge of how vaccination treatments affect macrophages that have not been infected. The existence of optimal control is proven. The Hamiltonian function and the maximum principle of the Pontryagin are used to describe the optimal control. In addition, we use the theory of optimal control to develop an algorithm that leads to choosing the best vaccination plan. The best numerical solutions have been discovered using the forward and backward fractional Euler method.
