1Matrix3d m = Matrix3d::Ones();
2m(0,2) += 1e-4;
3cout << "Here's the matrix m:" << endl << m << endl;
4cout << "m.isOnes() returns: " << m.isOnes() << endl;
5cout << "m.isOnes(1e-3) returns: " << m.isOnes(1e-3) << endl;
6