compute power spectrum of an ARMA process defined by x(n)=a(1)x(n-1)+...+a(l)x(n-l)+e(n)+b(1)e(n-1)+...+b(k)e(n-k)