根据对称性,可以把积分化成在区域

0≤x≤1,0≤y≤1上积分的4倍

4∫∫(x^2+y^2)dxdy=4∫(0,1)dx∫(0,1)[x^2+y^2]dy=4∫(0,1)[x^2+1/3]dx=8/3