In this paper, the mathematical model which describes the glucose-insulin homeostasis in healthy rats is
considered. The model is discretized by constructing a nonstandard finite difference (NSFD) scheme to
obtain the numerical solutions. The equilibrium point of the discretized model is determined and stability
analysis of the discretized model is discussed. The effect of time step sizes on 4th order Runge-Kutta
method and NSFD method is presented. Also, comparison of NSFD scheme solution, Runge-KuttaFehlberg Method (RKF45) solution and analytical solution are presented in graphical form. The
effectiveness of the proposed method in the solution and interpretation of the model is observed.