sx = linspace(1, 10, 100);
c = ones(1, 10);
z = zeros(1,10);
y = [1 zeros(1, 39) 1];
stem(y)
axis([-0 42 -0.1 1.5])
xlabel('Numer próbki');
ylabel('Amplituda');
title('Impuls jednostkowy');