1#include <unsupported/Eigen/MatrixFunctions>
2#include <iostream>
3
4using namespace Eigen;
5
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