Function definition and center of expansion

f:log(x); x0:1;

All Taylor polynomials up to degree n

n:2;
create_list(sum(at(diff(f,x,j),x=x0)/j!*(x-x0)^j,j,0,k),k,0,n)$
append([f],%)$
plot2d(%,[x,x0-0.7,3])$