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 10*bf2c3715SXin Li MatrixXd A(3,3); 11*bf2c3715SXin Li A << 0, -pi/4, 0, 12*bf2c3715SXin Li pi/4, 0, 0, 13*bf2c3715SXin Li 0, 0, 0; 14*bf2c3715SXin Li std::cout << "The matrix A is:\n" << A << "\n\n"; 15*bf2c3715SXin Li std::cout << "The matrix exponential of A is:\n" << A.exp() << "\n\n"; 16*bf2c3715SXin Li } 17