Stability analysis and numerical simulations of IVGTT glucose-insulin interaction models with two time delays

    Sayed Saber   Affiliation
    ; Ahmad Alalyani Affiliation


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.

Keyword : glucose-insulin regulatory system, insulin secretion, ultradian oscillation, delay differential equation model

How to Cite
Saber, S., & Alalyani, A. (2022). Stability analysis and numerical simulations of IVGTT glucose-insulin interaction models with two time delays. Mathematical Modelling and Analysis, 27(3), 383–407.
Published in Issue
Aug 12, 2022
Abstract Views
PDF Downloads
Creative Commons License

This work is licensed under a Creative Commons Attribution 4.0 International License.


E. Ackerman, L.C. Gatewood, J.W. Rosevear and G.D. Molnar. Model studies of blood-glucose regulation. Bulletin of Mathematical Biophysics, 27:21–37, 1965.

E. Ackerman, J.W. Rosevear and W.F. McGuckin. A mathematical model of the glucose-tolerance test. Physics in Medicine and Biology, 9(2):203–213, 1964.

R.N. Bergman and C. Cobelli. Minimal modeling, partition analysis, and the estimation of insulin sensitivity. Fed Proc., 39(1):110–115, 1980.

R.N. Bergman, Y.Z. Ider, C.R. Bowden and C. Cobelli. Quantitative estimation of insulin sensitivity. American Journal of Physiology-Endocrinology and Metabolism, 236(6):E667, 1979.

V.W. Bolie. Coefficients of normal blood glucose regulation. Journal of Applied Physiology, 16(5):783–788, 1961.

A. Caumo, R.N. Bergman and C. Cobelli. Insulin Sensitivity from Meal Tolerance Tests in Normal Subjects: A Minimal Model Index. The Journal of Clinical Endocrinology & Metabolism, 85(11):4396–4402, 2000.

A. De Gaetano and O. Arino. Mathematical modelling of the intravenous glucose tolerance test. Journal of Mathematical Biology, 40:136–168, 2000.

L.C. Gatewood, E. Ackerman, J.W. Rosevear, G.D. Molnar and T.W. Burns. Tests of a mathematical model of the blood-glucose regulatory system. Computers and Biomedical Research, 2(1):1–14, 1968.

T.A. Gresl, R.J. Colman, T.C. Havighurst, L.O. Byerley, D. B. Allison, D.A. Schoeller and J. W. Kemnitz. Insulin sensitivity and glucose effectiveness from three minimal models: effects of energy restriction and body fat in adult male rhesus monkeys. American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, 285(6):R1340–R1354, 2003.

Y. He, M. Wu and J.-H. She. Delay-dependent stability criteria for linear systems with multiple time delays. IEE Proceedings - Control Theory and Applications, 153(4):447–452, 2006.

W.M. Hirsch, H. Hanisch and J.-P. Gabriel. Differential equation models of some parasitic infections: Methods for the study of asymptotic behavior. Communications on Pure and Applied Mathematics, 38(6):733–753, 1985.

Y. Kuang. Delay differential equations: With applications in population dynamics. Boston, Academic Press, 1993.

C.P. Li and F.R. Zhang. A survey on the stability of fractional differential equations. The European Physical Journal Special Topics, 193:27–47, 2011.

J. Li and Y. Kuang. Analysis of a model of the glucoseinsulin regulatory system with two delays. SIAM Journal on Applied Mathematics, 67(3):757–776, 2007.

J. Li, Y. Kuang and B. Li. Analysis of IVGTT glucose-insulin interaction models with time delay. Discrete and Continuous Dynamical Systems - B, 1(1):103–124, 2001.

J. Li, Y. Kuang and C.C. Mason. Modeling the glucoseinsulin regulatory system and ultradian insulin secretory oscillations with two explicit time delays. Journal of Theoretical Biology, 242(3):722–735, 2006.

J. Li, M. Wang, A. De Gaetano, P. Palumbo and S. Panunzi. The range of time delay and the global stability of the equilibrium for an IVGTT model. Mathematical Biosciences, 235(2):128–137, 2012.

G. Pacini and R.N. Bergman. MINMOD: a computer program to calculate insulin sensitivity and pancreatic responsivity from the frequently sampled intravenous glucose tolerance test. Computer Methods and Programs in Biomedicine, 23(2):113–122, 1986.

S. Saber, E.B.M. Bashier, S.M. Alzahrani and I.A. Noaman. A mathematical model of glucose-insulin interaction with time delay. Journal of Applied & Computational Mathematics, 3(7):416–421, 2018.

G.M. Steil, A. Volund, S.E. Kahn and R.N. Bergman. Reduced Sample Number for Calculation of Insulin Sensitivity and Glucose Effectiveness From the Minimal Model: Suitability for Use in Population Studies. Diabetes, 42(2):250–256, 1993.