syms xf = 3*x^2 + x*exp(x) -cos(x^2);f0 = subs(f,x,0);f1 = subs(f,x,1);dfx = diff(f);df0 = subs(dfx,x,0);
