Abstract
An efficient finite difference approach is adopted to analyze the solution of a novel fractional-order mathematical model to control the co-circulation of double strains of dengue and COVID-19. The model is primarily built on a non-integer Caputo fractional derivative. The famous fixed-point theorem developed by Banach is employed to ensure that the solution of the formulated model exists and is ultimately unique. The model is examined for stability around the infection-free equilibrium point analysis, and it was observed that it is stable (asymptotically) when the maximum reproduction number is strictly below unity. Furthermore, global stability analysis of the disease-present equilibrium is conducted via the direct Lyapunov method. The non-standard finite difference (NSFD) approach is adopted to solve the formulated model. Furthermore, numerical experiments on the model reveal that the trajectories of the infected compartments converge to the disease-present equilibrium when the basic reproduction number ([Formula: see text]) is greater than one and disease-free equilibrium when the basic reproduction number is less than one respectively. This convergence is independent of the fractional orders and assumed initial conditions. The paper equally emphasized the outcome of altering the fractional orders, infection and recovery rates on the disease patterns. Similarly, we also remarked the importance of some key control measures to curtail the co-spread of double strains of dengue and COVID-19.