Abstract
The diabetes mellitus model (DMM) is explored in this study. Many health issues are caused by this disease. For this reason, the integer order DMM is converted into the time delayed fractional order model by fitting the fractional order Caputo differential operator and delay factor in the model. It is proved that the generalized model has the advantage of a unique solution for every time t. Moreover, every solution of the system is positive and bounded. Two equilibrium states of the fractional model are worked out i.e. disease free equilibrium state and the endemic equilibrium state. The risk factor indicator, R(0) is computed for the system. The stability analysis is carried out for the underlying system at both the equilibrium states. The key role of R(0) is investigated for the disease dynamics and stability of the system. The hybridized finite difference numerical method is formulated for obtaining the numerical solutions of the delayed fractional DMM. The physical features of the numerical method are examined. Simulated graphs are presented to assess the biological behavior of the numerical method. Lastly, the outcomes of the study are furnished in the conclusion section.