% F4_6
% Carol Lucas
% Script to call Glucose-Dependent Insulin Secretion Models
% Four state model, two state model and normalized model
%     for a range of external glucose values
%     GlucInsModel calls glucfun4 (4 states)
%       glucfun2 (2 states) , glucfunnorm (two states normalized)
mess=' Program pauses after each figure, hit return to continue.'
mess = ' Note that simulations can take a long time.'
GRAY=[8:2:20];
for k=1:length(GRAY)
    GlucInsModel_Init
    GNew=GRAY(k)
    GlucInsModel(GNew);
end