Todo puede calcularse a partir del $m=n=1$ diferenciando con respecto a $a$ y $b$ : $$ \frac{1}{(x+a)^m} = \frac{(-1)^{m-1}}{m!} \frac{\partial^{m-1}}{\partial a^{m-1}} \frac{1}{x+a}, $$ por lo que si escribimos $I_{m,n}(a,b)$ para su integral, $$ I_{m,n}(a,b) = \frac{(-1)^{m-1}}{m!n!} \frac{\partial^{m+n-2}}{\partial a^{m-1}\partial b^{n-1}} I_{1,1}(a,b). $$ $I_{1,1}$ es fácil de calcular: tenemos $$ I_{1,1}(a,b) = \frac{\log{a}-\log{b}}{a-b}, $$ escribiéndolo simétricamente. Ahora podemos aplicar la fórmula de Leibniz $$ D^n(fg) = \sum_{k=0}^n \binom{n}{k} (D^kf)(D^{n-k}g) $$ para calcular las derivadas. En particular, $$ \frac{(-1)^{m-1}}{m!} \frac{\partial^{m-1}}{\partial a^{m-1}} \frac{\log{a}-\log{b}}{a-b} \\ = \frac{\log{a}-\log{b}}{(a-b)^m} + \frac{(-1)^{m-1}}{m!} \sum_{k=1}^{m-1} \binom{m-1}{k} \frac{(-1)^{m-1-k} (m-k)!}{(a-b)^{m-k}} \frac{ (-1)^{k-1}(k-1)! }{a^k} \\ = \frac{\log{a}-\log{b}}{(a-b)^m} - \sum_{k=1}^{m-1} \frac{(m-1)!}{k!(m-k-1)!}\frac{ (m-k)!(k-1)!}{m!} \frac{1}{a^k(a-b)^{m-k}} \\ = \frac{\log{a}-\log{b}}{(a-b)^m} - \sum_{k=1}^{m-1} \left(\frac{1}{k}-\frac{1}{m} \right)\frac{1}{a^k(a-b)^{m-k}}. $$ Ahora se puede diferenciar con respecto a $b$ para obtener un par de series finitas y dobles para $I_{m,n}(a,b)$ .
Alternativamente, se puede combinar esta expresión, que es efectivamente $I_{m,1}(a,b)$ con una relación de recurrencia a partir de la integración por partes, a saber $$ \frac{1}{a^mb^n} + m I_{m+1,n}(a,b) + n I_{m,n+1}(a,b) = 0, $$ que da la respuesta en términos de una sola serie. Lamentablemente, parece que hemos perdido la simetría entre $a$ y $b$ en este proceso.
0 votos
functions.wolfram.com/FuncionesElementarias/Power/26/02/01/0004 escribe cada binomio como una función de meijer después de factorizar un a^m y b^n del denominador. Luego busca la función g de Meijer. Existe un teorema para la integral del producto de estas funciones sobre el eje real positivo. La solución puede reducirse a funciones elementales.