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(2,2); 11*bf2c3715SXin Li A << cos(pi/3), -sin(pi/3), 12*bf2c3715SXin Li sin(pi/3), cos(pi/3); 13*bf2c3715SXin Li std::cout << "The matrix A is:\n" << A << "\n\n"; 14*bf2c3715SXin Li std::cout << "The matrix square root of A is:\n" << A.sqrt() << "\n\n"; 15*bf2c3715SXin Li std::cout << "The square of the last matrix is:\n" << A.sqrt() * A.sqrt() << "\n"; 16*bf2c3715SXin Li } 17