function [ y] = myexp13(x,n) %calcolo dell'esponenziale %strumento approssimazione di Taylor %(sviluppo di centro 0 e raggio x arrestato all'ordine n) %f(x)=e^x app 1 +x +x^2/2 +....+x^n/factorial(n) %y= (....((1+x)+x^2/2)+x^3/factorial(3))...)+x^n/factorial(n) %ripetere l'operazione di somma accumulando il valore su y % % y=1; % % for k=1:n % % y=y+x^k/factorial(k); % % end %inefficente x^k=x^(k-1)*x; factorial(k)=.... %riduco il costo computazionale=numero di operazioni %dimostrate che nell'implementazione 1 il costo e' circa n^2 y=1; %p=1; m=1; z=1; for k=1:n % p=p*x; % m=m*k; z=z*x/k; % y=y+p/m; y=y+z; end %implementazione piu' efficiente %il costo diviene lineare in n (4n operazioni aritmetiche) end