1 Matrix4d A = Matrix4d::Random(4,4); 2 cout << "Here is a random 4x4 matrix:" << endl << A << endl; 3 HessenbergDecomposition<Matrix4d> hessOfA(A); 4 Matrix4d pm = hessOfA.packedMatrix(); 5 cout << "The packed matrix M is:" << endl << pm << endl; 6 cout << "The upper Hessenberg part corresponds to the matrix H, which is:" 7 << endl << hessOfA.matrixH() << endl; 8 Vector3d hc = hessOfA.householderCoefficients(); 9 cout << "The vector of Householder coefficients is:" << endl << hc << endl; 10