The method of least squares is a well known technique of curve fitting, giving estimates for the values of the coefficients in equations such as y = ax, y = ax + b, etc, in terms of observed values of the variables x, y. It is less widely realised that the technique can be extended to estimating matrix coefficients.