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