Maximum likelihood estimation of the GLS model with unknown parameters in the disturbance covariance matrix