Abstract
This paper presents a systematic study of a mathematical model of glucose and insulin interaction with two time delays, with a focus on analytical studies, bifurcation analysis, and very well numerical simulations. This model based on the Intra-Venous Glucose Tolerance Test (IVGTT) and is presented with two time delays. One delay is the insulin response time to an increase in glucose concentration, and the hepatic glucose production time delay is the other. Then, we establish results on positivity, boundedness, and persistence. We also provide sufficient stability analysis conditions for both local and global asymptotic stability of the proposed models. For the latter, two different strategies are used: stability bifurcation analysis and Lyapunov-Krasovskii functionals. We investigate different regions of parameter space using two approaches, that yield different sets of sufficient conditions for global stability. The bifurcation graphs generated from our extensive and carefully designed simulations complement and confirm these analytical results. The insulin concentration level peaks after the glucose concentration level, according to the numerical simulations.