Matlab - Series de Taylor en Matlab

function a = ealax(x, nterm)

if length(x) ~= 1 them

    error('x debe ser un escalar')

end % if

if length(nterm) < 1 them

    error('nterm debe ser uno o mayor')

end %if

suma = 0 %suma inicial es cero

for i = 0 : nterm-1

    suma = suma + (x^i/factorial(i))

end %for

a = suma;

function y=taylor1(Y,n)

syms x

Y=eval(Y);

y=0;

x=0;

for l=0:n-1

y=y+eval(diff(Y,l))*(sym('x')^l)/factorial(l);

end

end

>> taylor1('sin(x)',3)

 ans =

- x^3/6 + x

>> taylor1('exp(x)',3)

 ans =

x^3/6 + x^2/2 + x + 1

taylor1('log(x+1)',3)

ans =

x^3/3 - x^2/2 + x