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

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

pause(1)

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

pause(1)

plot(t,x(:,3),'k-', t,x(:,6),'k--')
xlabel('t')
ylabel('x_3')
legend('nonlinear model', 'linearised model')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 9])
print -deps 'obr428.eps'
