It has been an open problem whether the reduction of a matrix to Hessenberg (almost triangular) form by Gaussian similarity transformations is numerically stable [2, p. 364]. We settle this question ...