for N = 5 : 5 : 20 
  Delta = exp(1) - somma_eN(N) ; 
  fprintf('N = %2d, Delta = %.15f\n', N, Delta) ; 
end