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

exp_esatto = [cos(30) sin(30); -sin(30) cos(30)];

X = eye(2); % somma parziale
T = eye(2); % 1/d! * A^d, aggiornato ad ogni passo
norme = [];
errori = [];
for k = 1:d
    T = 1/k * A * T;
    X = X + T;
    errori(k) = norm(X - exp_esatto);
    norme(k) = norm(T);
end
figure(1);
semilogy(errori);
figure(2);
semilogy(norme);