Abstract
Modeling glucose-insulin regulatory system plays a key role for treating diabetes, a serious health problem for numerous patients. The effect of the incommensurate fractional-order derivatives on a glucose-insulin regulatory model is studied in this work. It has been shown that the model exhibits some interesting dynamics, such as chaos and coexisting attractors, in response of a specific change in such derivatives’ values, even if it was slight. When comparing such model with some previous models, we have deduced a clear presence of wider chaotic regions once the values of these incommensurate-orders are changed.