I came on the following multiple integral while renormalizing elliptic multiple zeta values: $$\int_0^1\cdots \int_0^1\int_1^\infty {{1}\over{t_n(t_{n-1}+t_n)\cdots (t_1+\cdots+t_n)}} dt_n\cdots dt_1.$$ Only the variable $t_n$ goes from $1$ to $\infty$, the others all go from $0$ to $1$. Numerically, I seem to be getting a rational multiple of $\pi^n$. I would like to prove this. Has anyone ever seen an integral like this before?
Edit: After reworking out why I thought it would be a power of $\pi$, I now have a different integral, which numerically really does seem to give a rational multiple of $\zeta(n)$ for each $n$ (though I have only been able to go up to n=4 numerically). I want to integrate $1/(z_1\cdots z_n)$ over the part of the simplex $0\le z_n\le \cdots \le z_1\le 1$ in which $|z_i -z_{i+1}|>\varepsilon$, then let $\varepsilon\rightarrow 0$. The integral is $$\int_{n\epsilon}^{1-\varepsilon}\int_{(n-1)\varepsilon}^{z_n-\epsilon}\cdots \int_{\varepsilon}^{z_2-\varepsilon} {{1}\over{z_1\cdots z_n}} dz_n\cdots dz_1.$$$$\int_{n\epsilon}^{1-\varepsilon}\int_{(n-1)\varepsilon}^{z_1-\epsilon}\cdots \int_{\varepsilon}^{z_{n-1}-\varepsilon} {{1}\over{z_1\cdots z_n}} dz_n\cdots dz_1.$$ It actually diverges when $\varepsilon\rightarrow 0$, but it can be regularised like ordinary multizeta values by computing the integral as a power series in $ln(\epsilon)$$ln(\varepsilon)$ and $\epsilon$$\varepsilon$ and then taking its constant term. It's related to the previous integral by the variable change $z_1=\epsilon(t_1+\cdots+t_n),\ldots,z_n=\epsilon t_n$$z_1=\varepsilon(t_1+\cdots+t_n),\ldots,z_n=\varepsilon t_n$, but the new bounds of the integral form an $(n+1)$-angled polyhedron, not a cube.