td = [0 0.5 0.75 1 1.5 2 2.5 3 4 6];
gn = [70 150 165 145 90 75 65 75 80 75];
gd = [100 185 210 220 195 175 105 100 85 90];

xlab = '$t$ (hr)'; 		% X-label
ylab = 'Glucose (mg/dl)'; 			% Y-label
mytitle = 'GTT Model'; 	% Title

xx = linspace(0,6,200);
pn = [79.1814 171.5474 0.99272 1.81274 0.90056];
pd = [95.2124 263.1528 0.63349 1.03037 1.51604];
yn = pn(1)+pn(2)*exp(-pn(3)*xx)...
    .*cos(pn(4)*(xx-pn(5)));
yd = pd(1)+pd(2)*exp(-pd(3)*xx)...
    .*cos(pd(4)*(xx-pd(5)));

plot(xx,yn,'b-');
hold on;
plot(xx,yd,'r-');
plot(td,gn,'ob',td,gd,'or');
grid
xlim([0 6]);
ylim([0 250]);
fontlabs = 'Times New Roman'; 
xlabel(xlab,'FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');  
ylabel(ylab,'FontSize',14,'FontName',fontlabs,...
    'interpreter','latex');  
title(mytitle,'FontSize',16,'FontName',...
    'Times New Roman','interpreter','latex'); 
set(gca,'FontSize',12);

print -depsc GTT_modelA.eps		
