opt108

qq = y0;
pp = x0;
sim('opt108s');
%print -deps -sopt108s opt108.eps

t = data.time;
u = data.signals.values(:,1);
yy = data.signals.values(:,2);
[tt,uu]=stairs(t,u);
plot(tt,uu,'k--', t,yy, 'k-')
v = axis;
legend('u','y')
title('b')
xlabel('t')
ylabel('u,y')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'opt109b.eps'

pause
qq = y;
pp = x;
sim('opt108s');
t = data.time;
u = data.signals.values(:,1);
yy = data.signals.values(:,2);
[tt,uu]=stairs(t,u);
plot(tt,uu,'k--', t,yy, 'k-')
axis(v);
legend('u','y')
title('a')
xlabel('t')
ylabel('u,y')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'opt109a.eps'
