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