<s>[INST] Partial Fraction Decomposition of Exponential Generating Functions I want to see if it is possible to write
$$ \left(\frac{x}{e^x-1}\right) \left(\frac{x^2/2! }{e^x-1-x}\right) \left(\frac{x^3/3!}{e^x-1-x-x^2/2}\right)$$
as a linear combination of the factors
$$p(x)\frac{x}{e^x-1}+q(x)\frac{x^2/2}{e^x-1-x}+r(x)\frac{x^3/3!}{e^x-1-x-x^2/2}$$
I'm specifically looking for the functions $p(x), q(x)$ and $r(x)$ so I thought I would try a partial fraction decomposition method and this would imply that
$$\frac{A}{e^x-1}+\frac{B}{e^x-1-x}+\frac{C}{e^x-1-x-x^2/2}$$
where $A=xp(x), B=\frac{1}{2}x^2q(x), C=\frac{1}{6}x^3r(x)$.  So this means that
$$A(e^x-1-x)(e^x-1-x-x^2/2)+B(e^x-1)(e^x-1-x-x^2/2)+C(e^x-1)(e^x-1-x)=\frac{x^6}{12}$$
If we multiply through
$$A\left(e^{2x}-2e^x-2xe^x-\frac{x^2}{2}e^x+1+2x+\frac{3}{2}x^2+\frac{x^3}{2}\right)$$
$$+B\left(e^{2x}-2e^x-xe^x-\frac{x^2}{2}e^x+1+x+\frac{x^2}{2}\right)$$
$$+C\left(e^{2x}-2e^x-xe^x+1+x\right)$$
Thus we get
$$(A+B+C)(e^{2x}-2e^x)-(2A+B+C)xe^x-(A+B)\frac{x^2}{2}e^x+(A+B+C)+(2A+B+C)x+\left(\frac{3}{2}A+\frac{1}{2}B\right)x^2+\frac{1}{2}Ax^3=\frac{x^6}{12}$$
To me, now it seems that there is no solution since both $A+B+C=2A+B+C=0$ which means $A=0$ which then implies $B=C=0$.  Is this method incorrect for such a problem and if so is there another way to approach this problem?  Or did I make an error?
 [/INST] I think that your idea to decompose in simple fractions is good. You can put $K=\mathbb{Q}(x)$, replace $\exp(x)$ by $y$, and decompose the following fraction in $K(y)$:
$$\left(\frac{x}{y-1}\right) \left(\frac{x^2/2! }{y-1-x}\right) \left(\frac{x^3/3!}{y-1-x-x^2/2}\right)=\frac{A}{y-1}+\frac{B}{y-1-x}+\frac{C}{y-1-x-x^2/2}$$
To get $A$, you multiply by $y-1$ both side, simplify and put $y=1$, to get $B$, you multiply by $y-1-x$ both side, simplify and put $y=1+x$, etc .
When you are done, simply replace $y$ by $\exp(x)$.  
</s>