for mu = 1:Mu

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

end