function S=calcola_In(n)

% Calcolo di \sum_{k=0}^n \frac{1}{k!} \int_0^{1/4} (x^2)^k

% NB: La funzione e' scalare e non vettoriale, ovvero "n" e' un numero!

k=0:n;
argomento_somma=((1/4).^(2*k+1))./((2*k+1).*gamma(k+1));
S=sum(argomento_somma);

% v=(1/4).^(2*k+1);
% u=1./gamma(k+1);
% S=v*u';

% NB: si poteva fare pure
% S=0;
% for k=0:n
%     argomento=(1/4)^(2*k+1)/((2*k+1)*gamma(k+1));
%     S=S+argomento;
% end