Background: I'm an undergraduate at an institution with no researchers in analytic number theory, and no ties to the analytic number theory community. I believe I have found what is, as far as I can tell after some googling, a new family of integral representations of $\zeta(2n+1)$. I was told I should post this here by an algebraic number theorist at my university, to see if it was a known result or not.


Feel free to look at Equations (4), (5), and (7), to see if you recognize them before reading the entire text. Most of this text is for the special case of $\zeta(3)$, but all results generalize quite readily to the case of $\zeta(2n+1)$, which I discuss at the end.


**Beginning of setup**


As a preliminary step, we will split the sum $\zeta(3)$ into even and odd parts. Observe that $$
\sum\_{n=1}^\infty \frac{1}{(2n)^3} = \frac{1}{8} \sum\_{n=1}^\infty \frac{1}{n^3} = \frac{1}{8}\zeta(3)
$$
Therefore $$
\sum\_{n=1}^\infty \frac{1}{(2n+1)^3} = \frac{7}{8}\zeta(3) \tag{1}
$$
We now write $$
\frac{1}{(2n+1)^3} = \left(\int\_0^1 x^{2n} dx\right)^3 = \int\_0^1 \int\_0^1 \int\_0^1 x^{2n} y^{2n} z^{2n} \tag{2} dx dy dz
$$
and plug this triple integral into Equation (1), obtaining $$
\zeta(3) = \frac{8}{7} \sum\_{n=0}^\infty \left[ \int\_0^1 \int\_0^1 \int\_0^1 x^{2n} y^{2n} z^{2n} dxdydz \right]
$$
Using the absolute convergence of the geometric series on $(0,1)$, we can interchange the sum and integrals, obtaining $$
\zeta(3) = \frac{8}{7} \int\_0^1 \int\_0^1 \int\_0^1 \frac{1}{1-x^2y^2z^2} dxdydz \tag{3}
$$


**End of setup**


With this part out of the way, I will now describe the integral representations that form the meat of my question.


Consider the $u$-substitution for the integral given by Equation (3) $$
x = \frac{\sinh{u}}{\cosh{v}}, \qquad y = \frac{\sinh{v}}{\cosh{w}}, \qquad z = \frac{\sinh{w}}{\cosh{u}}
$$
Some computation shows that $$
dxdydz = [1-(\tanh{u}\tanh{v}\tanh{w})^2]dudvdw
$$
But dividing both sides by $1-(\tanh{u}\tanh{v}\tanh{w})^2$ and rewriting the left-hand side in terms of $(x,y,z)$ gives $$
\frac{1}{1-x^2y^2z^2}dxdydz = dudvdw
$$
After changing limits, the transformed integral then reads $$
\zeta(3) = \frac{8}{7} \int\_0^\infty \int\_0^{\sinh^{-1}(\cosh(w))} \int\_0^{\sinh^{-1}(\cosh(v))} dudvdw \\
= \frac{8}{7} \int\_0^\infty \int\_0^{\sinh^{-1}(\cosh(w))} \sinh^{-1}(\cosh(v)) dv dw \tag{4}
$$


But the most interesting part of the above is the following generalization. Let $f : (0,\infty) \to (0,\infty)$ be a function satisfying the following three properties.


1. $f$ is surjective, with $\lim\limits\_{x\to 0} f(x) = 0$ and $\lim\limits\_{x \to \infty} f(x) = \infty$
2. $f$ is invertible
3. $f$ is differentiable


Now consider the $u$-substitution $$
x = \frac{\sinh{f(u)}}{\cosh{f(v)}}, \qquad y = \frac{\sinh{f(v)}}{\cosh{f(w)}}, \qquad z = \frac{\sinh{f(w)}}{\cosh{f(u)}}
$$
We then find that $$
dxdydz = f'(u)f'(v)f'(w)[1-(\tanh{f(u)}\tanh{f(v)}\tanh{f(w)})^2]dudvdw
$$
and hence $$
\zeta(3) = \frac{8}{7} \int\_0^\infty \int\_0^{g(w)} \int\_0^{g(v)} f'(u)f'(v)f'(w) dudvdw \tag{5}
$$
where $g(x) = f^{-1}(\sinh^{-1}(\cosh(f(x))))$.


All of the above can be generalized to $\zeta(2n+1)$. Equation (3) becomes $$
\zeta(2n+1) = \frac{2^{2n+1}}{2^{2n+1}-1} \int\_0^1 \int\_0^1 \cdots \int\_0^1 \frac{1}{1-x\_1^2 x\_2^2 \cdots x\_{2n+1}^2} dx\_1 dx\_2 \cdots dx\_{2n+1} \tag{6}
$$
Our $u$-substitution becomes (where $u\_{(2n+1)+1} = u\_1$) $$
x\_i = \frac{\sinh(f(u\_i))}{\cosh(f(u\_{i+1}))}
$$
which transforms Equation (6) to $$
\zeta(2n+1) = \frac{2^{2n+1}}{2^{2n+1}-1} \int\_0^\infty \int\_0^{g(u\_{2n})} \int\_0^{g(u\_{2n-1})} \cdots \int\_0^{g(u\_1)} f'(u\_1)f'(u\_2)\cdots f'(u\_{2n+1}) du\_1 du\_2 \cdots du\_{2n+1} \tag{7}
$$


My questions are then: 


1. Are Equations (4), (5), and (7) known results?
2. Equation (4) is a volume integral. Exploring the region of integration in Mathematica numerically, it looks like an octant of a hyperbolic cube with vertices at infinity. Is this in fact the case?



Is this integral representation of $\zeta(2n+1)$ known?
As indicated in my comment, some of these integrals are essentially known, and involve the hyperbolic "Beukers-Kolk-Calabi" change of variables.


In particular, in [this paper](http://arxiv.org/pdf/1003.3602v1.pdf), Z. Silagadze shows in (27) that
\begin{equation\*}
 \zeta(n) = \frac{2^n}{2^n-1}\int\_0^1\cdots\int\_0^1 \frac{dx\_1\cdots dx\_n}{1-x\_1^2\cdots x\_n^2},
\end{equation\*}
which clearly yields (7) above. 


The hyperbolic change of variables is used later in that paper to obtain a related formula over the $2n$-simplex $\Delta\_{2n}$, namely
\begin{equation\*}
 \zeta(2n+1) = -\frac{1}{2n}\frac{2^{2n+1}}{2^{2n+1}-1}\int\_{\Delta\_{2n}} \log(\tan x\_1 \tan x\_2\ldots \tan x\_{2n})dx\_1\cdots dx\_{2n}.
\end{equation\*}
This integral is amenable to further interesting modifications.


*EDIT:*


1. The arXiv paper cited above has been published [in the journal "Resonance"](http://www.ias.ac.in/resonance/php/cissue.php)
2. [Another paper](http://arxiv.org/pdf/1207.2055.pdf) by the same author (Z. Silagadze) comments on these integrals.
3. The "Beukers-Kolk-Calabi" change of variables was also considered in [this paper](http://arxiv.org/abs/1003.2170).
4. [This paper](http://arxiv.org/pdf/1007.3436.pdf) (cited below by Benjamin Dickman) is also relevant here.
5. The relation of these representations to [Amoebas](https://en.wikipedia.org/wiki/Amoeba_%28mathematics%29) is quite interesting, and worthy of further exploration.


