syms x; >> fx=x^2*cos(x) fx = x^2*cos(x) >> [an,bn,f]=fseries(fx,x,16,-pi,pi) syms x; >> fx=x^2*cos(x)
f(x)=sum(4(1+n^2)*pi*(-1)^(n+1)*cos(n*x)/(n^2-1)^2,{n,0,infinity})