Taylor = function (F,x0,n)
{x=x0; coeff=c(eval(F)); fatt=1
for (k in 1:n)
{F=D(F,"x"); fatt=fatt*k
coeff=c(coeff,eval(F)/fatt)}
coeff}