function L25_gaussian

fs = 1;
tr = 2;
sigma = tr * fs / 2 / 1.28;

b = firgauss(10, 'minorder', sigma)'
b = b / sum(b);
a = zeros(1, length(b));
a(1) = 1;

H = tf(b, a, 1/fs);

figure(4);
step(H);
