%Second order system
t=0:0.1:16;
subplot(221)
y=exp(-0.25*t).*sin(t);
yu=exp(-0.25*t);
yl=-yu;
y0=zeros(size(t));
y1=ones(size(t));
p1=plot(t,y1,'b-');
set(p1,'linewidth',2)
axis([0 16 -0 1.2])
xlabel('t')
ylabel('y')
title('integrator')
%
y=sin(t);
subplot(223)
p2=plot(t,y,'b-');
set(p2,'linewidth',2)
axis([0 16 -1.2 1.2])
xlabel('t')
ylabel('y')
title('oscillatory');
print -depsc int-osc.eps