MATLAB function that computes $A^n$, being $A$ an square matrix and $n$ an natural number?

Question

Do any of you know a matlab function that computes $A^n$ just like in the following example?

This question arises in response of discrete time state space systems.

Answer 1

Try this:

syms n;
A = sym([0 1; -1/6 5/6]);
[V,D] = eig(A);
An = V*(D^n)*inv(V);
pretty(An);