for
mu = 1:Mu
uybar(mu,:) = filter(1,[1, -Az(mu,mu)],T*uexcite(mu,:));
end