Stable matrix inverse method for array DF using trapezoidal factorisation

An approach to finding the inverse of the covariance matrix required for DF using maximum-likelihood methods (MLM) is to factorise the matrix into LD (lower triangular-diagonal) form. An alternative factorisation, trapezoidal-diagonal (TD) is outlined which is both more stable and, for a small number of signals, more economical.