function X = exptaylor(A, d)
% esponenziale tramite la serie di MacLaurin, fermandoci a grado d

X = eye(2); % somma parziale
T = eye(2); % 1/d! * A^d, aggiornato ad ogni passo
for k = 1:d
    T = 1/k * A * T;
    X = X + T;
end