Evolutionary computational method for tuberculosis model with fuzziness

This work investigates the computational study of a six-compartmental mathematical model of tuberculosis disease dynamics with the impact of vaccination. Traditional mathematical models presume that all variables are precise and can be measured or calculated precisely. However, in many real-world scenarios, variables may need to be more accurate or easier to quantify, resulting in model uncertainty. Considering this, fuzziness is introduced into the model by taking the contact, recovery, and death rates due to disease as fuzzy membership functions. Two numerical computational schemes, forward Euler and nonstandard finite difference (NSFD), are designed to solve the model. The positivity and convergence for the developed method are investigated, which are significant characteristics of these dynamical models, and it is revealed that these features are preserved in the extended scheme. Numerical computations are performed to support the analytical results. The numerical and computational results indicate that the proposed NSFD method adequately represents the dynamics of the disease despite the uncertainty and heterogeneity. Moreover, the obtained method generates plausible predictions that regulators can use to design and develop control strategies to support decision-making.