sys = tf(1, [1 1]);
[y,t]=impulse(sys);
plot(t,y)
xlabel('t / T_1');
ylabel('T_1 g');

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

