clear
[t,x]=sim('simssl');

y = x(:,2);
plot(t,y,'k-')
xlabel('t')
ylabel('y')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 9])
print -deps 'obr422.eps'
