b=1; a=[1 -1.2 .9]; n=400; x=randn(n,1); y=filter(b,a,x); plot(y); title('Time Series of y_k = (-1.2)y_{k-1} + .9y_{k-2} + e_k'); print -deps ts7