// Updated(18-7-07)
// 11.8

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