Recently I have been working on the following question:
Let \begin{equation*} I_n:= \int\int_{[0,1]^2}\frac{(3x)^n(1-x^3)^n(3y)^n(1-y^3)^n}{(1+xy+x^2 y^2)^{2n+1}}dxdy. \end{equation*} Define \begin{equation*} L(\chi_{-3}, 2) := \sum_{n=1}^{\infty}\frac{\chi_{-3}(n)}{n^2} = \frac{1}{1^2}-\frac{1}{2^2}+\frac{1}{4^2}-\frac{1}{5^2}+... \end{equation*} I want to show that \begin{equation*} I_0=L(\chi_{-3}, 2), \qquad I_1=3L(\chi_{-3}, 2)-2. \end{equation*}
The reason I am working on this is because these $I_n$ satisfy an Apéry-like recursion: \begin{equation*} (n+1)^2 I_{n+1}-(10n^2+10n+3)I_n+9n^2 I_{n-1}=0, \end{equation*} from which it would follow that $I_n \in \mathbb{Q}+\mathbb{Q}L(\chi_{-3}, 2)$ for all $n$. However, I cannot find an original source on this, which is why I tried showing these initial values on my own. I can only show the case where $n=0$, which is as follows:
Let us recall the identity \begin{equation*} 1+z+z^2 = \frac{1-z^3}{1-z} \implies \frac{1}{1+z+z^2}=\frac{1-z}{1-z^3}. \end{equation*} Let $z=xy$, which results in \begin{equation*} I_0 = \int\int_{[0,1]^2}\frac{1}{1+xy+x^2 y^2}dxdy = \int\int_{[0,1]^2}\frac{1-xy}{1-(xy)^3}dxdy = \int\int_{[0,1]^2}\frac{1}{1-(xy)^3}-\frac{xy}{1-(xy)^3}dxdy. \end{equation*} Let us note that on the open unit square we have $|xy| < 1$, so we can write \begin{equation*} \frac{1}{1-(xy)^3} = \sum_{n=0}^{\infty}(xy)^{3n}, \qquad \frac{xy}{1-(xy)^3}= \sum_{n=0}^{\infty}(xy)^{3n+1}. \end{equation*} Plugging this back in our integral gives \begin{equation*} I_0 = \int\int_{[0,1]^2}\sum_{n=0}^{\infty}(xy)^{3n}-\sum_{n=0}^{\infty}(xy)^{3n+1}dxdy = \int_0^1\int_0^1\sum_{n=0}^{\infty}\left((xy)^{3n}-(xy)^{3n+1}\right)dxdy. \end{equation*} Now, using Fubini-Tonelli, we get that \begin{equation*} I_0 = \sum_{n=0}^{\infty} \int_0^1\int_0^1\left((xy)^{3n}-(xy)^{3n+1}\right)dxdy = \sum_{n=0}^{\infty} \frac{1}{(3n+1)^2}-\frac{1}{(3n+2)^2} = \sum_{n=1}^{\infty}\frac{\chi_{-3}(n)}{n^2}=L(\chi_{-3}, 2). \end{equation*}
Does have anyone have an idea for how to handle $n=1$? The same methods as in $n=0$ do not seem to work out if we expand it all (or maybe I do not see how all these terms combined to yield the right result?). Does anyone have a source which talks about this integral and what it is called? I could not find it online either.