A three-step method for obtaining asymptotically maximum likelihood ARMA spectral estimates using only linear transformations of the data is presented. The computational efficiency of the algorithm is comparable to that of Yule-Walker algorithms, but the three-step method gives asymptotically statistically efficient estimates. The implementation of the algorithm is discussed in detail, and numerical examples are presented to illustrate its performance.