aminit;

%Illustration of stability of solution
t=0:0.01:6;
x0=sin(t)+0.8*cos(2.2*t)+2.5;

clf; subplot(311)
pl=plot(t,x0, 'k');
set(pl,'Linewidth',AM_ref_linewidth);

hold on
pl=plot(t,x0+0.5,'r', t,x0-0.5, 'r');
set(pl,'Linewidth',AM_ref_linewidth);

x=x0-0.4*sin(t-0.3);
pl=plot(t,x,'b--');
set(pl,'Linewidth',AM_data_linewidth);

amaxis([0 6 0 4.5]);
set(gca,'Ytick',[0 2 4]);
xlabel('time');
ylabel('state');

amprint('stable-lyap.eps');
