Matlab - Serie de mclaurin

function s=mlacurint1(x,n)

s=0;

for i=1:n

    s=s+((-1)^(i+1))*x^i/factorial(i);

end

sol_matlab=log(1+x)

>> s=mlacurint1(0.5,50)

sol_matlab =

    0.4055

s =

    0.3935

>> s=mlacurint1(-0.5,50)

sol_matlab =

   -0.6931

s =

   -0.6487