Let me give an example
Suppose I want to expand cos(x)4 as in a series of functions Ck=(1/2−cos(kx)). That is cos(4x)=k0C0+k1C1+k2C2⋯ and I am looking for coefficients k
Now, I can use trig_reduce() to express cos(x)4 as a sum cos(kx), but how do I get get cos(x)4 in terms of Ck.