Para justificar la representación integral, se observa que la
$$
\frac{-\ln(1-x)}{x}=\sum_{n\geq 1} \frac{x^{n-1}}{n}
$$
converge normalmente (por lo tanto de manera uniforme absolutamente) en $[0,1/\phi^2]$.
Así que término por término integración es posible y los rendimientos de la fórmula deseada.
Ahora para calcular la integral, tenemos que saber que una antiderivada de $-\ln(1-x)/x$ es el llamado dilogarithm. De acuerdo a Wikipedia, Don Zagier dijo una vez que el dilogarithm es la única función matemática con sentido del humor.
Así
$$
\sum_{n\geq 1} \frac{1}{\phi^{2n}n^2}=\mbox{Li}_2(1/\phi^2)=\mbox{Li}_2\left( \frac{3-\sqrt{5}}{2}\right).
$$
He encontrado aquí: http://en.wikipedia.org/wiki/Polylogarithm que
$$
\mbox{Li}_2\left( \frac{3-\sqrt{5}}{2}\right)=\frac{\pi^2}{15}-\ln^2\left(\frac{\sqrt{5}-1}{2}\right).
$$
Esa es su fórmula, ya $\ln^2(1/\phi)=\ln^2(\phi)$.
El dilogarithm fórmula que se utilizó anteriormente es más probable demostrado aquí: http://books.google.fr/books?id=u_UVn_iquj0C&pg=PA141&lpg=PA141&dq=trilogarithm+identities&source=bl&ots=5lp7VWoN1p&sig=xkvQm5jZFj_bcb0CpPKqsmIa1Fw&hl=fr&sa=X&ei=d5EOUd_iPLK-0QH8k4CAAQ&ved=0CDoQ6AEwAQ#v=onepage&q=trilogarithm%20identities&f=false