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