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