Integrate[(x^2 + y^2), {x, 0, r}, {y, 0, Sqrt[r^2 - x^2]}]