In this paper we propose a way to solve the problem of singularities in model reference adaptive control of linear multi-input-multi-output (MIMO) systems using a parameter modification procedure based on the least squares covariance matrix inverse. The scheme does not require any explicit prior knowledge about the leading coefficient matrix associated with the control input and secures a uniform lower bound for the determinant of the estimate of this matrix. A global convergence analysis is presented.