17faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez#include <unsupported/Eigen/MatrixFunctions> 27faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez#include <iostream> 37faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez 47faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandezusing namespace Eigen; 57faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez 67faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandezint main() 77faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez{ 87faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez Matrix4cd A = Matrix4cd::Random(); 97faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez MatrixPower<Matrix4cd> Apow(A); 107faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez 117faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez std::cout << "The matrix A is:\n" << A << "\n\n" 127faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez "A^3.1 is:\n" << Apow(3.1) << "\n\n" 137faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez "A^3.3 is:\n" << Apow(3.3) << "\n\n" 147faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez "A^3.7 is:\n" << Apow(3.7) << "\n\n" 157faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez "A^3.9 is:\n" << Apow(3.9) << std::endl; 167faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez return 0; 177faaa9f3f0df9d23790277834d426c3d992ac3baCarlos Hernandez} 18