El siguiente método se debe a Euler, también explicado aquí . Escribamos:
$$S_r = x_1^r + x_2^r +x_3^r$$
para $r\in \mathbb{N}$
Considera la función:
$$f(u) = -\sum_{j=1}^3\log\left(1-\frac{x_j}{u}\right)$$
La expansión en serie alrededor del infinito de esta función viene dada por:
$$f(u) = \sum_{k=1}^{\infty}\frac{S_k}{k u^k}\tag{1}$$
También tenemos que $\exp\left[-f(u)\right]$ es un polinomio de tercer grado en $u^{-1}$ . Esto significa que podemos calcular $S_4$ exponenciando la serie (1), calculando el término de cuarto orden en $u^{-1}$ y poniendo el resultado a cero.
En este caso particular podemos simplificar esto, escribiendo:
$$f(u) = 12 \sum_{k=1}^{\infty}\frac{1}{k u^k} + \sum_{k=4}^{\infty}\frac{S_k'}{k u^k} = -12 \log\left(1-\frac{1}{u}\right) + \sum_{k=4}^{\infty}\frac{S_k'}{k u^k}$$
donde $S'_k = S_k - 12$ .
Por lo tanto, tenemos:
$$\exp\left[-f(u)\right] = \left(1-u^{-1}\right)^{12}\left[1- \frac{S_4 - 12}{4 u^4}\right] +\mathcal{O}(u^{-5})$$
Siendo éste un polinomio de 3er grado en $u^{-1}$ significa que el coeficiente de $u^{-4}$ es igual a cero. Por lo tanto, tenemos:
$$S_4 = 4 \left[\binom{12}{4} + 3\right] = 1992$$