Раз область
- единичный круг с центром в начале координат, то данный интеграл будем вычислять в полярных координатах.
Заменяя:
и дифференцируя 

Из определения
,т.е. получим что
и угол изменяется в пределах 
Подставим в подынтегральное выражение в полярной системе координат

Осталось подставить все данные в двойной интеграл и перейти к повторному интегралу, то есть
