sim('obr6awas');
vys1 = res;
sim('obr6awbs');
vys2 = res;

vys=vys1;
t=vys(:,1);
w=vys(:,2);
u=vys(:,3);
y=vys(:,4);
integ=vys(:,5);
subplot(2,1,1)
plot(t,w,'k:',t,y,'k-',t,u,'k--',t,integ,'k-.');
v=axis;
axis([v(1) v(2) v(3) v(4)]);
legend(' w', 'y', 'u', 'integ',4);
%xlabel('t');
%ylabel('');
title('(a)')

vys=vys2;
t=vys(:,1);
w=vys(:,2);
u=vys(:,3);
y=vys(:,4);
integ=vys(:,5);
subplot(2,1,2)
plot(t,w,'k:',t,y,'k-',t,u,'k--',t,integ,'k-.');
v=axis;
axis([v(1) v(2) v(3) v(4)]);
legend(' w', 'y', 'u', 'integ',4);
xlabel('t');
%ylabel('');
title('(b)')

print -deps 'obr6awab.eps'
