código matlab aproximação taylor

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function y = aprox_taylor(f,a,x,k)
soma = 0;
for n=0:k-1
 fl = feval(f,n);
 soma = soma + fl*(x-a).^n/factorial(n);
 %f^n(a)= {n! se n é par / -n! se n é impar
 %x - pt investigacao
 
 % k - numero termos utilizados/ a - 
end
y= soma;
end