Abstract
In this article, we consider a fractal-fractional tobacco mathematical model with generalized kernels of Mittag-Leffler functions for qualitative and numerical studies. From qualitative point of view, our study includes; existence criteria, uniqueness of solution and Hyers-Ulam stability. For the numerical aspect, we utilize Lagrange's interpolation polynomial and obtain a numerical scheme which is further illustrated simulations. Lastly, a comparative analysis is presented for different fractal and fractional orders. The numerical results are divided into four figures based on different fractal and fractional orders. We have found that the fractional and fractal orders have a significant impact on the dynamical behaviour of the model.