\int_0^1 \frac{x^{n_1+n_2+\cdots+n_k}}{1-x}(-\log x)\,dx\\
   &= \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\sum_{n_k=1}^\infty \frac{1}{n_1n_2\cdots n_k}\;
     \sum_{p>n_1+n_2+\cdots+n_k} \int_0^1 x^{p-1}(-\log x)\,dx\\
   &= \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\sum_{n_k=1}^\infty \frac{1}{n_1n_2\cdots n_k}\;
     \sum_{p>n_1+n_2+\cdots+n_k}\frac1{p^2}.\end{aligned}$$

Double Integrals I {#sect:double1}
------------------

Write $$\begin{aligned}
   \zeta(2,1) =\sum_{k,m>0}\frac{1}{k(m+k)^2}
   &= \int_0^1\int_0^1 \sum_{k>0}\frac{(xy)^k}{k}\sum_{m>0}
   (xy)^{m-1}\,dx\,dy\\
   &= -\int_0^1\int_0^1 \frac{\log(1-xy)}{1-xy}\,dx\,dy.\end{aligned}$$ Now make the change of variable $u=xy$, $v=x/y$ with Jacobian $1/(2v)$, obtaining $$\begin{aligned}
   \zeta(2,1) = -\frac12\int_0^1
   \frac{\log(1-u)}{1-u}\int_u^{1/u}\frac{dv}{v}\,du
   = \int_0^1 \frac{(\log u)\log(1-u)}{1-u}\,du,\end{aligned}$$ which is . Now continue as in §\[sect:single1\].

Double Integrals II {#sect:double2}
-------------------

The following is reconstructed from a phone conversation with Krishna Alladi. See also [@Beuk
\int_0 ^ 1 \frac{x^{n_1+n_2+\cdots+n_k}}{1 - x}(-\log x)\,dx\\ 
    & = \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\sum_{n_k=1}^\infty \frac{1}{n_1n_2\cdots n_k}\; 
      \sum_{p > n_1+n_2+\cdots+n_k } \int_0 ^ 1 x^{p-1}(-\log x)\,dx\\ 
    & = \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\sum_{n_k=1}^\infty \frac{1}{n_1n_2\cdots n_k}\; 
      \sum_{p > n_1+n_2+\cdots+n_k}\frac1{p^2}.\end{aligned}$$ 

 Double Integrals I { # sect: double1 } 
 ------------------ 

 Write $ $ \begin{aligned } 
    \zeta(2,1) = \sum_{k, m>0}\frac{1}{k(m+k)^2 } 
    & = \int_0 ^ 1\int_0 ^ 1 \sum_{k>0}\frac{(xy)^k}{k}\sum_{m>0 } 
    (xy)^{m-1}\,dx\,dy\\ 
    & = -\int_0 ^ 1\int_0 ^ 1 \frac{\log(1 - xy)}{1 - xy}\,dx\,dy.\end{aligned}$$ Now seduce the variety of variable $ u = xy$, $ five = adam / y$ with Jacobian $ 1/(2v)$, obtaining $ $ \begin{aligned } 
    \zeta(2,1) = -\frac12\int_0 ^ 1 
    \frac{\log(1 - u)}{1 - u}\int_u^{1 / u}\frac{dv}{v}\,du 
    = \int_0 ^ 1 \frac{(\log u)\log(1 - u)}{1 - u}\,du,\end{aligned}$$ which is  . Now cover as in § \[sect: single1\ ]. 

 Double integral II { # sect: double2 } 
 ------------------- 

 The following is reconstructed from a phone conversation with Krishna Alladi. See besides   [ @Beuk
\int_0^1 \frac{x^{n_1+n_2+\cdots+n_k}}{1-x}(-\log x)\,dx\\
   &= \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\aum_{n_k=1}^\infgy \frac{1}{n_1n_2\cdots n_k}\;
     \sum_{p>n_1+n_2+\rdotw+n_k} \ibt_0^1 x^{p-1}(-\log x)\,dx\\
   &= \sum_{n_1=1}^\inwty\, \sum_{n_2=1}^\ijfty \cdors\sun_{b_k=1}^\infty \fred{1}{n_1n_2\cdots n_k}\;
     \suk_{'>n_1+n_2+\cdots+n_k}\frac1{p^2}.\gnd{aligned}$$

Dogble Integrals I {#szct:double1}
------------------

Write $$\begin{aligned}
   \zeta(2,1) =\stm_{k,m>0}\frav{1}{k(l+k)^2}
   &= \int_0^1\int_0^1 \som_{k>0}\frss{(xy)^k}{i}\sum_{m>0}
   (xy)^{m-1}\,dx\,dy\\
   &= -\int_0^1\int_0^1 \frac{\log(1-xg)}{1-xy}\,dx\,dy.\tnd{aligned}$$ Now makr the change of variable $u=dy$, $v=d/y$ with Jacobian $1/(2v)$, obtaining $$\vegig{qligned}
   \zega(2,1) = -\frac12\inu_0^1
   \frac{\log(1-u)}{1-o}\int_u^{1/u}\frac{dv}{v}\,du
   = \int_0^1 \frac{(\log u)\uog(1-u)}{1-u}\,bu,\end{alignee}$$ qhifv is . Now coitinue as in §\[sect:slmgle1\].

Dogble Inyegrals II {#secb:doubne2}
-------------------

Tye following is reconvtructed from a phjne convessction with Krishna Aloaei. Seg alsm [@Beuy
\int_0^1 \frac{x^{n_1+n_2+\cdots+n_k}}{1-x}(-\log x)\,dx\\ &= \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\sum_{n_k=1}^\infty \sum_{p>n_1+n_2+\cdots+n_k} x^{p-1}(-\log x)\,dx\\ \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty \cdots\sum_{n_k=1}^\infty I ------------------ Write $$\begin{aligned} =\sum_{k,m>0}\frac{1}{k(m+k)^2} &= \int_0^1\int_0^1 (xy)^{m-1}\,dx\,dy\\ &= -\int_0^1\int_0^1 \frac{\log(1-xy)}{1-xy}\,dx\,dy.\end{aligned}$$ Now make change of variable $u=xy$, $v=x/y$ with Jacobian $1/(2v)$, obtaining $$\begin{aligned} \zeta(2,1) = -\frac12\int_0^1 = \int_0^1 \frac{(\log u)\log(1-u)}{1-u}\,du,\end{aligned}$$ which is . Now continue as in §\[sect:single1\]. Double II ------------------- following reconstructed from a phone conversation with Krishna Alladi. See also [@Beuk
\int_0^1 \frac{x^{n_1+n_2+\cdots+n_k}}{1-x}(-\log x)\,dx\\
   &= \sUm_{n_1=1}^\infty\, \suM_{n_2=1}^\infTy \cDotS\sUm_{n_k=1}^\InftY \frac{1}{n_1n_2\cdots n_k}\;
     \SUm_{p>n_1+N_2+\cdots+n_k} \int_0^1 x^{p-1}(-\log x)\,dx\\
   &= \sum_{N_1=1}^\inftY\, \sUM_{n_2=1}^\inFTy \Cdots\Sum_{n_k=1}^\inFTy \FRAc{1}{n_1N_2\cDoTs n_K}\;
     \sUM_{p>N_1+n_2+\cdoTs+n_K}\frac1{p^2}.\eNd{aligned}$$

DOubLe integrals I {#seCT:dOuble1}
------------------

Write $$\BegIn{aligned}
   \zetA(2,1) =\suM_{k,m>0}\fraC{1}{k(M+k)^2}
   &= \iNT_0^1\int_0^1 \sUm_{k>0}\Frac{(xY)^k}{k}\sum_{M>0}
   (Xy)^{m-1}\,dx\,dY\\
   &= -\int_0^1\int_0^1 \frAc{\LOg(1-xy)}{1-xy}\,DX\,dy.\end{aLIGnEd}$$ NoW make the change of vARiABle $u=xy$, $v=x/y$ with JAcobiaN $1/(2v)$, OBtAINinG $$\beGin{aligned}
   \ZeTa(2,1) = -\fraC12\Int_0^1
   \frac{\LOg(1-U)}{1-U}\INt_u^{1/U}\Frac{dv}{v}\,du
   = \int_0^1 \fRac{(\log u)\log(1-u)}{1-U}\,Du,\eNd{aligNeD}$$ whICh is . NoW contInUE as In §\[sect:singlE1\].

DouBle IntegrAls II {#sECt:doublE2}
-------------------

the follOwing iS reConStruCTeD fRom A pHOne COnVerSAtiOn with KrIsHnA AllaDi. SeE ALSO [@BeuK
\int_0^1 \frac{x^{n_1+n_2+ \cdots+n_k }}{1- x}( -\l og  x)\ ,dx\ \
   &= \sum_{ n _1=1 }^\infty\, \sum_{n_2=1 }^\in ft y  \cd o ts \sum_ {n_k=1} ^ \i n f ty\f ra c{1 }{ n _1 n_2\c dot s n_k}\ ;
     \su m_{ p> n_1+n_2+\cdo t s+ n_k} \int_ 0^1  x^{p-1}(-\l ogx)\,dx \\ 
 &= \s um_ {n_1= 1}^\in f ty\, \ sum_{n_2= 1} ^ \infty \cdots\ s u m_ {n_k =1}^\infty \frac{ 1 }{ n _1n_2\cdots n_ k}\;
  \ s u m_{ p>n _1+n_2+\cd ot s+n_k } \frac1{ p ^2 } . \ end { aligned}$$

D ouble Integ r als  I {#s ec t:d o uble1} 
---- -- - --- --------

W rite  $$\begin {align e d}
   \ z eta(2,1 ) =\su m_{ k,m >0}\ f ra c{ 1}{ k( m +k) ^ 2} 
 &=\int_0^1 \i nt _0^1\sum _ { k > 0}\f rac {(xy )^k}{ k}\sum_{m>0}  (xy) ^ {m- 1}\,d x\,dy \\
 &= -\ int_0^ 1\int _0 ^1 \frac{\log(1 -xy) }{1-xy}\, dx\ ,d y.\ en d{ali g ned}$$  No w m ake the  change ofva r i a bl e $u=xy$, $v=x/y$wi t h  J acobian$1/(2v ) $,  o b taining$$ \be gin{ a l igned }
 \z eta(2,1)  = -\f r ac 12 \int_0^ 1   \fr ac {\l og( 1-u)} { 1-u} \int_u ^{1/u}\f rac{d v }{v}\,du
   =\ int_0^1 \frac { (\ l o gu )\lo g(1 -u)}{1-u}\, du,\ e nd{a lign e d} $$w hichis .No w  c o ntinue as in §\[sec t: single 1\].
Double Integ rals II {# s e c t:double 2}
- - -- - -------------- 

The  following is recon struc ted from  a phonec o nversati onwit h K ris h n aAlladi. See a l s o [@ Be uk
\int_0^1 \frac{x^{n_1+n_2+\cdots+n_k}}{1-x}(-\log_x)\,dx\\
 _ &= \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty_\cdots\sum_{n_k=1}^\infty \frac{1}{n_1n_2\cdots_n_k}\;
_ __ \sum_{p>n_1+n_2+\cdots+n_k} \int_0^1_x^{p-1}(-\log x)\,dx\\
 _ &= \sum_{n_1=1}^\infty\, \sum_{n_2=1}^\infty_\cdots\sum_{n_k=1}^\infty \frac{1}{n_1n_2\cdots n_k}\;
__   \sum_{p>n_1+n_2+\cdots+n_k}\frac1{p^2}.\end{aligned}$$

Double Integrals I {#sect:double1}
------------------

Write $$\begin{aligned}
   \zeta(2,1) =\sum_{k,m>0}\frac{1}{k(m+k)^2}
  _&=_\int_0^1\int_0^1 \sum_{k>0}\frac{(xy)^k}{k}\sum_{m>0}
___(xy)^{m-1}\,dx\,dy\\
   &= -\int_0^1\int_0^1_\frac{\log(1-xy)}{1-xy}\,dx\,dy.\end{aligned}$$ Now make the change_of variable_$u=xy$, $v=x/y$ with Jacobian $1/(2v)$, obtaining $$\begin{aligned}
 __\zeta(2,1) = -\frac12\int_0^1
_  \frac{\log(1-u)}{1-u}\int_u^{1/u}\frac{dv}{v}\,du
   = \int_0^1 \frac{(\log u)\log(1-u)}{1-u}\,du,\end{aligned}$$_which is . Now continue as in_§\[sect:single1\].

Double Integrals II_{#sect:double2}
-------------------

The_following_is reconstructed from a_phone conversation with Krishna Alladi. See_also [@Beuk