1#include <unsupported/Eigen/MatrixFunctions>
2#include <iostream>
3
4using namespace Eigen;
5
6int main()
7{
8  Matrix4cd A = Matrix4cd::Random();
9  MatrixPower<Matrix4cd> Apow(A);
10
11  std::cout << "The matrix A is:\n" << A << "\n\n"
12	       "A^3.1 is:\n" << Apow(3.1) << "\n\n"
13	       "A^3.3 is:\n" << Apow(3.3) << "\n\n"
14	       "A^3.7 is:\n" << Apow(3.7) << "\n\n"
15	       "A^3.9 is:\n" << Apow(3.9) << std::endl;
16  return 0;
17}
18