function eA = taylor_exp(A,N)
eA = zeros(size(A));
for n=0:N
   eA = eA + A^n/factorial(n);
end
