The main objective of this research is to investigate a new fractional mathematical model involving a nonsingular derivative operator to discuss the clinical implications of diabetes and tuberculosis coexistence. The new model involves two distinct populations, diabetics and nondiabetics, while each of them consists of seven tuberculosis states: susceptible, fast and slow latent, actively tuberculosis infection, recovered, fast latent after reinfection, and drug-resistant. The fractional operator is also considered a recently introduced one with Mittag–Le er nonsingular kernel. The basic properties of the new model including non-negative and bounded solution, invariant region, and equilibrium points are discussed thoroughly. To solve and simulate the proposed model, a new and e cient numerical method is established based on the productintegration rule. Numerical simulations are presented, and some discussions are given from the mathematical and biological viewpoints. Next, an optimal control problem is dened for the new model by introducing four control variables reducing the number of infected individuals. For the control problem, the necessary and su cient conditions are derived and numerical simulations are given to verify the theoretical analysis
