Descubrí la siguiente fórmula mientras trabajaba en la secuencia A141697 de la OEIS. No tengo ni idea de si es algo trivial o no. Me encantaría saber más al respecto.
$$ \textrm{Li}_\nu\left(z\right)= \frac { 6 (1+z)^{-\nu-1} + \displaystyle\sum_{k=0}^{-\nu-1} \displaystyle\left( -6 \displaystyle{{-\nu-1}\displaystyle\choose k}+7\sum_{j=0}^{k+1}(-1)^j (k-j+1)^{-\nu} {{-\nu+1}\ \choose j} \right) z^k } { 7 (1 - z)^{-\nu+1} }z $$ cuando $\nu$ es un número entero negativo.
En caso de que haya escrito mal la fórmula, aquí están mis códigos en dos idiomas diferentes:
Para Pari-GP:
mypolylog(n, x) = { ( 6*(x+1)^(-n-1)
+ sum(k=0,-n-1, (-6*binomial(-n-1,k)
+ 7*sum(j=0,k+1, (-1)^j * (k-j+1)^(-n) * binomial(-n+1,j)))*x^k) ) * x
/ (7*(1-x)^(-n+1) ) }
Para Mathematica:
mypolylog[n_, x_] := (6*(x+1)^(-n-1)
+ Sum[(x^k*(-6*Binomial[-n-1, k]
+ 7*Sum[(-1)^j*(k-j+1)^(-n)*Binomial[-n+1, j], {j, 0, k+1}])),
{k, 0, -n-1}]) / ( 7*(1 - x)^(-n+1) ) * x
0 votos
¿Cómo se obtiene la fórmula? ¿Has probado a reformular las funciones racionales utilizando los números de Stirling o Eulerian dados en es.wikipedia.org/wiki/Polilogaritmo#Valores_particulares ?
2 votos
@gammatester Descubrí empíricamente esta fórmula utilizando una herramienta que escribí para descubrir identidades: github.com/baruchel/oeis que detectó A141697(n)=3*A168524(n)-2*A154337(n). La fórmula anterior proviene de eso.
0 votos
Buen trabajo. Puedes usar $6=t$ y $7=t+1$ donde $t\neq -1$ en su fórmula.
1 votos
Esto también podría aplicarse a A141696.