1typedef Matrix<double,3,Dynamic> Matrix3Xd; 2Matrix3Xd M = Matrix3Xd::Random(3,5); 3Projective3d P(Matrix4d::Random()); 4cout << "The matrix M is:" << endl << M << endl << endl; 5cout << "M.colwise().homogeneous():" << endl << M.colwise().homogeneous() << endl << endl; 6cout << "P * M.colwise().homogeneous():" << endl << P * M.colwise().homogeneous() << endl << endl; 7cout << "P * M.colwise().homogeneous().hnormalized(): " << endl << (P * M.colwise().homogeneous()).colwise().hnormalized() << endl << endl;