Las respuestas anteriores no mencionar que también podemos encontrar la forma cerrada expresiones mediante el uso armónico de la suma de las fórmulas.
Introducir la suma
$$S_m(x) = \sum_{n\ge 1} \frac{(nx)^{2m-1}}{e^{2\pi nx}-1}$$
así que estamos buscando para encontrar $S_m(1).$
La suma plazo es armónico y puede ser evaluado por la inversión de sus Mellin transformar.
Recordar que el armónico suma de identidad
$$\mathfrak{M}\left(\sum_{k\ge 1} \lambda_k g(\mu_k x);s\right) =
\left(\sum_{k\ge 1} \frac{\lambda_k}{\mu_k^s} \right) g^*(s)$$
donde $g^*(s)$ es la transformada de Mellin $g(x).$
En el presente caso tenemos
$$\lambda_k = 1, \quad \mu_k = k \quad \text{y} \quad
g(x) = \frac{x^{2m-1}}{e^{2\pi x}-1}.$$
Necesitamos la transformada de Mellin $1/(e^{2\pi x}-1)$ que es
$$\int_0^\infty \frac{1}{e^{2\pi x}-1} x^{m-1} dx =
\int_0^\infty \frac{e^{-2\pi x}}{1-e^{-2\pi x}} x^{m-1} dx.$$
Ampliando el interior de la suma obtenemos
$$\int_0^\infty \sum_{q\ge 1} e^{-2\pi qx} x^{m-1} dx
= \sum_{q\ge 1} \int_0^\infty e^{-2\pi qx} x^{m-1} dx
= \Gamma(s) \sum_{q\ge 1} \frac{1}{(2\pi q)^s}
= \frac{1}{(2\pi)^s} \Gamma(s) \zeta(s).$$
De ello se desprende que la Mellin transformar $g^*(s)$ $g(x)$ está dado por
$$\frac{1}{(2\pi)^{s+2m-1}} \Gamma(s+2m-1) \zeta(s+2m-1).$$
Por lo tanto, la Mellin transformar $Q_m(s)$ $S_m(x)$ está dado por
$$\frac{1}{(2\pi)^{s+2m-1}} \Gamma(s+2m-1) \zeta(s)\zeta(s+2m-1).$$
El Mellin de inversión integral aquí es
$$\frac{1}{2\pi i} \int_{3/2-i\infty}^{3/2+i\infty} Q_m(s)/x^s ds$$
de la cual evaluamos desplazando a la izquierda por una expansión de alrededor de cero.
Necesitamos examinar que los polos de la función gamma plazo son cancelados por los dos zeta función de los términos. Los polos de la función gamma plazo son todos los números enteros menores o iguales a $-(2m-1).$ La llanura zeta plazo cancela todos los incluso queridos y también cancela el uno al $-(2m-1)+1$ desde el compuesto zeta término al $m>1.$ El compuesto zeta plazo cancela todos los impares, comenzando desde los $-(2m-1)-2.$ Esto deja sólo dos polos en $s=1$ y a las $s=-(2m-1).$
Por lo tanto tenemos los siguientes residuos que contribuyen:
$$\mathrm{Res}(Q_m(s)/x^s; s=1) =
\frac{1}{(2\pi)^{2m}} \Gamma(2m) \zeta(2m) \frac{1}{x}
\\= \frac{1}{(2\pi)^{2m}} (2m-1)! \frac{(-1)^{m+1} B_{2m} (2\pi)^{2m}}{2\times (2m)!}
\frac{1}{x} =
\frac{(-1)^{m+1} B_{2m}}{2\times 2m} = \frac{(-1)^{m+1} B_{2m}}{4m} \frac{1}{x}$$
y
$$\mathrm{Res}(Q_m(s)/x^s; s=-(2m-1)) = \zeta(-(2m-1)) \zeta(0) x^{2m-1}
\\= -\frac{1}{2} \times -\frac{B_{2m}}{2m} x^{2m-1}= \frac{B_{2m}}{4m} x^{2m-1}.$$
Al $m=1$ hay un poste en$s=0$, lo que contribuye
$$ \mathrm{Res}(Q_m(s)/x^s; s=0) = -\frac{1}{4\pi}.$$
Ahora nos especializamos en el caso de $m$ es impar y por lo menos tres. A continuación, podemos observar un hecho que es de vital importancia, a saber, que $Q(s)$ es impar en la línea $\Re(s) = -(m-1),$ que ahora íbamos a probar.
Recordar la fórmula para $Q(s),$, lo que fue
$$\frac{1}{(2\pi)^{s+2m-1}} \Gamma(s+2m-1) \zeta(s)\zeta(s+2m-1).$$
Ahora pon $s = -(m-1)+it$ para obtener
$$\frac{1}{(2\pi)^{m+it}} \Gamma(m+it) \zeta(-m+it+1)\zeta(m+it).$$
Utilizando la ecuación funcional de la Riemann zeta función en el primer zeta función plazo, esto se convierte en
$$\frac{2}{(2\pi)^{2m}} \Gamma(m+) \Gamma(m -) \zeta(m+) \zeta(m -)
\cos\left(\frac{1}{2}(m -)\pi\right).$$
Ahora el uso de la serie de la representación de la función zeta en la mitad de avión $\Re(s)>1$ es fácil ver que los dos zeta función de los términos son los conjugados de cada uno de los otros y, por tanto, su producto es un número real. Lo que es más importante el valor se mantiene igual al $t$ cambia de signo. Del mismo modo el límite de la definición de la función gamma muestra que el producto de las dos de la función gamma términos es un número real, que también se mantiene igual al $t$ cambia de signo. El factor de en frente, no depende del $t.$ Que deja
$$\cos\left(\frac{1}{2}(m -)\pi\right)
= \cos\left(\frac{\pi}{2}m - \frac{\pi}{2}\right)
= (-1)^{\frac{m-1}{2}} \sin\left(\frac{\pi}{2}\right).$$
Este término es, en efecto extraño en $t$ con el que concluye el argumento.
Volver a la original de cálculo ahora sabemos a cambio de la integral de $\Re(s) = 3/2$ $\Re(s) = -(m-1),$recogiendo sólo el residuo de a $s=1$ y se benefician del hecho de que la contribución de la izquierda segmento vertical es cero. (Set $x=1$ antes del turno.) Hemos establecido que para $m$ impar y $m>1,$
$$\sum_{n\ge 1} \frac{n^{2m-1}}{e^{2\pi n}-1} = (-1)^{m+1} \frac{B_{2m}}{4m}
= \frac{B_{2m}}{4m}.$$
Hay un término de corrección de $1/\pi/8$ al $m=1$ originarios como la mitad de la contribución de la pole en cero, que es en el lado izquierdo de la curva de nivel.