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