
fun=inline('exp(x)');
interval=[-1,1];
order=3;

order_vett=1:10;

% Punti tests e valutazione funzione.
x=linspace(interval(1),interval(2),10000);
fx=feval(fun,x);

for index=1:length(order_vett)
    order=order_vett(index);
    
    % Calcola i coefficienti di Taylor, in ordine crescente per grado.
    v=(0:order)';
    coeffs=1./gamma(v+1);
    
    % Valuta l'errore in norma infinito.
    px=polyval(flipud(coeffs),x);
    inferr=norm(fx-px,inf);
    
    fprintf('\n \t [DEGREE]: %3.0f [ERR. INFTY NORM]: %2.2e',order,inferr);
end


