clear; clc;

A = [0 1; 1 0];

trial1 = exp(A)
trial2 = eye(size(A));
for n = 1 : 10
    trial2 = trial2 + A ^ n / factorial(n);
end
trial2
trial3 = exp(1) ^ A