Hi,

I have this taylor series expansion:

> taylor(exp(-x)*sin(x), x, 11);

Now I need to do a legendre economization. So I want to write x, x², x³ etc (from the taylor series) as linear combinations of the Legendre polynomials. My idea was to simply write something like:

>solve({aP(0,x)+bP(1,x)+cP(2,x) = x²}, {a,b,c});

This is supposed to give coefficients for the legendre expansion of x². This doesn't work because it gives solutions in terms of x, which is useless to me. Is there any way?