syms x;

x_0 = 0;
n=4;

f(x) = log(1+x);

T(x) = taylor(f,x,x_0);

disp(T);