function y = system_3(x) % clip outside of range [-1,1] range = 1.0; N = length(x); n = 1:1:N; y = x; ul = range*ones(1,N); ll = -range*ones(1,N); sup_ul = x > ul; sup_ll = x < ll; y(sup_ul) = ul(sup_ul); y(sup_ll) = ll(sup_ll);