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

plot(t,x(:,1),'k-', t,x(:,2),'k--')
xlabel('t [s]')
ylabel('x_1 [dm]')
legend('nonlinear model', 'linearised model')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 9])
print -deps 'obr424.eps'
