clear
sim('sin1');
plot(t,u,'k:',t,y,'k-')
axis([0 180 -5 5])
legend('u','y')
xlabel('t [min]')
ylabel('u [°C]  y[°C]')
grid on

pause (1)
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 9])
print -deps 'obr431.eps'
