Ach, guckst du. Jetzt ist alles klar. Das Integral kannst du in Polarkoordinaten sofort lösen. In Polarkoordinaten ist:
$$\binom{x}{y}=\binom{r\cos\varphi}{r\sin\varphi}$$Die Bedingung \(0<x^2+y^2\le1\) bedeutet, dass \(r\in]0;1]\) liegt. Die Bedingungen \(0\le x\le1\) und \(0\le y\le1\) bedeuten, dass du dich im ersten Quadranten befindest, also \(\varphi\in[0;\frac{\pi}{2}]\). Das Flächenelement transformiert sich beim Übergang zu Polarkoordinaten gemäß: \(dx\,dy=r\,dr\,d\varphi\). Damit lautet das Integral:
$$I=\int\limits_0^{\pi/2}d\varphi\int\limits_0^1dr\frac{1}{1+r^2}r=\int\limits_0^{\pi/2}d\varphi\left[\frac{1}{2}\ln\left(1+r^2\right)\right]_0^1$$$$\phantom{I}=\frac{1}{2}\ln(2)\int\limits_0^{\pi/2}d\varphi=\frac{1}{2}\ln(2)\frac{\pi}{2}=\frac{\pi}{4}\ln(2)$$