1 #include <unsupported/Eigen/MatrixFunctions> 2 #include <iostream> 3 4 using namespace Eigen; 5 main()6int main() 7 { 8 const double pi = std::acos(-1.0); 9 Matrix3d A; 10 A << cos(1), -sin(1), 0, 11 sin(1), cos(1), 0, 12 0 , 0 , 1; 13 std::cout << "The matrix A is:\n" << A << "\n\n" 14 "The matrix power A^(pi/4) is:\n" << A.pow(pi/4) << std::endl; 15 return 0; 16 } 17