De Robert Z de la respuesta, tenemos
m2+n2mn(m2−n2)2=12mn(m+n)2+12mn(m−n)2.
Por lo tanto, la suma
S=∑m,n∈N1m≠n1mn(m2−n2)2
se puede dividir en dos partes:
S1=∑m,n∈N1m≠n12mn(m+n)2=∑m,n∈N112mn(m+n)2−∞∑n=118n4
y
S2=∑m,n∈N1m≠n12mn(m−n)2=∑m,n∈N1m>n12mn(m−n)2+∑m,n∈N1m<n12mn(n−m)2.
Veamos S2 más de cerca. El primer sumando satisface
∑m,n∈N1m>n12mn(m−n)2=∑m,k∈N112m(m+k)k2,
donde k=m−n. El segundo sumando satisface
∑m,n∈N1m>n12mn(m−n)2=∑n,l∈N112(n+l)nl2,
donde l=n−m. Por lo tanto, se puede cambiar el nombre de las variables ficticias y obtener
S2=∑m,k∈N112m(m+k)k2+∑n,l∈N112(n+l)nl2=∑m,n∈N112mn2(m+n)+∑m,n∈N112m2n(m+n).
Así,
S2=∑m,n∈N1(12mn2(m+n)+12m2n(m+n))=∑m,n∈N112m2n2.
Esto demuestra que
S2=12∑m∈N11m2∑n∈N11n2=12(ζ(2))(ζ(2))=12(π26)2=π472.
Ahora, para S1, vemos que el último término es
∞∑n=118n4=18ζ(4)=18(π490)=π4720.
Nos quedamos con el primer término, que es
∑m,n∈N112mn(m+n)2=ζ(3,1)=14ζ(4)=14(π490)=π4360,
por Robert Z de la respuesta. Por lo tanto,
S1=π4360−π4720=π4720.
La combinación de los dos resultados anteriores, tenemos
S=S1+S2=π4720+π472=11π4720=118ζ(4).
Aquí es otro gran recurso para aprender acerca de las múltiples zeta función. Corolario 1.64 en la página 26 y discusiones que llevan a esta conclusión son importantes para este problema.
Aquí está una elemental forma de mostrar que la ζ(3,1)=14ζ(4). Recordemos que
ζ(s1,s2,…,st)=∑n1,n2,…,nt∈N1n1>n2>…>nt1ns11ns22⋯nstt.
En primer lugar, observar que
1m2n2=1(m+n)2n2+1(m+n)2m2+2(m+n)3m+2(m+n)3n.
Tomando la suma de la expresión de arriba más de m,n∈N1, obtenemos
(ζ(2))2=ζ(2,2)+ζ(2,2)+2ζ(3,1)+2ζ(3,1)=2(ζ(2,2)+2ζ(3,1)),
así
ζ(2,2)+2ζ(3,1)=12(ζ(2))2=54ζ(4).
Ahora, tenga en cuenta que
52ζ(4)=(ζ(2))2=∑m,n∈N11m2n2=∑n∈N11n4+∑m,n∈N1m>n1m2n2+∑m,n∈N1m<n1m2n2.
Por lo tanto,
52ζ(4)=ζ(4)+ζ(2,2)+ζ(2,2),
hacer
ζ(2,2)=34ζ(4).
La combinación de (1) y (2), obtenemos
ζ(3,1)=12(54ζ(4)−ζ(2,2))=14ζ(4).