where D is the first-quadrant part of the disk with center (0,0) and radius 1.

Solution:

Maple Check:

int( int(x*y, x = 0 ..sqrt(1-y^2)), y = 0..1);