syms x
y=(exp(x)+exp(-x))/2;
taylor(f,x,0,'order',5)
ans =
x^4/24 + x^2/2 + 1