function y = bds1_5a(a, x, yn1)

N = length(x);
y = zeros(1,N);

y(1) = a*yn1+x(1);

for j = 2:N
    y(j) = a*y(j-1) + x(j);
end

