From Pirani's lectures on General Relativity I got the following identity which he asks the reader to prove by induction which is easy
$$\partial_{\alpha_1} \partial_{\alpha_2} ... \partial_{\alpha_n} \frac{1}{r}= (-1)^n \frac{(2n)!}{2^n n!} \frac{x^{\alpha_1} x^{\alpha_2} ... x^{\alpha_n}}{r^{2n+1}}+ \text{kronecker delta terms}$$
However, I am not entirely satisfied by the proof by induction since any induction proof assumes knowledge of the given expression. The expression in question is often found out by trial and error method, which seems like a formidable task for this particular identity given its complexity (or I am just dumb).
Is there a systematic proof of this identity which does not use induction or if not how can we guess the expression after computing it for a few small values of $n$ or via some ansatz. I am guessing there must be an elegant method to arrive at this formula "without guessing" as the coefficient $\frac{(2n)!}{2^n n!}$ seems familiar from somewhere: probably Rodrigues formula for Legendre polynomials.