1MatrixXd ones = MatrixXd::Ones(3,3);
2cout << "The operator norm of the 3x3 matrix of ones is "
3     << ones.operatorNorm() << endl;
4