function y = system_1(x) N = length(x); z(1) = x(1); z(2:N+1) = x(1:N) z(N+2) = x(N); for n = 1:N y(n) = median(z(n:n+2)); end