restart; f:=sin(x^2); S1:=series(f,x=0,12); P:=convert(S1,polynom); A1:=subs(x=1/2,f); A2:=subs(x=1/2,P); evalf(A1); evalf(A2);