Y = tf([1 0],[1 -0.9],-1);
covar(Y,1)