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