% filter of system output and input

function yf=filt1(X,fil,R,m)

R=inv(R');

for i=1:m
     R(i,:)=R(i,:)/(abs(R(i,:))*ones(m,1));
     end
Z=R*X';
yf=Z(fil);