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