К сожалению, Ваш рисунок не открывается, но, как решать задачу, понятно.
Сложность этой задачи не в вычислении интегралов, а в правильном составлении интегралов (подынтегральных выражений).
Да, применяем формула Гаусса-Остроградского. Вычисляем интеграл по объему:
[$8747$] div(A) dV,
где div(a) = [$8706$]Ax/[$8706$]x + [$8706$]Ay/[$8706$]y + [$8706$]Az/[$8706$]z
В данном случае div(A) = 3*x2 + 3*y2 + 3*z2 = 3*r2, где r - расстояние до начала координат.
То есть точки с одинаковой дивергенцией вектора А лежат на сферах с центром в начале координат. Поэтому интегрировать придется тоже "сферами".
Теперь смотрим какая у нас поверхность.
x2 + y2 + z2 = x или (x -1/2)2 + y2 + z2 = (1/2)2
То есть это сфера радиуса 1/2 с центром в точке (1/2; 0; 0), то есть ее "южный полюс" находится в начале координат, а "северный" в точке (1; 0; 0).
Значит, находясь в начале координат, будем через эту сферу проводить сферы радиусом от 0 до 1, центр которых в начале координат. В дальнейшем, чтобы не было путаницы, о какой сфере идет речь, сферу, уравнение которой дано в условиях задачи, будем называть телом или сферическим телом.
Опять же, для удобства будем обозначать радиус тела как R (R=1/2).
Тогда искомый интеграл [$8747$] div(A) dV равен:
[$8747$] (3*r2)*S(r)*dr с пределами от 0 до 2*R,
где r- радиус какой-либо сферы с центром в начале координат (на каждой такой сфере радиуса r одинаковое значение дивергенции = 3*r2);
S(r) - площадь части сферы радиуса r, находящейся внутри тела.
Значит, теперь ищем S(r).
Если смотреть из начала координат, то часть сферы радиуса r, находящаяся внутри тела, представляет собой поверхность сферы, вырезанную конусом, ось которого проходит через "северный" и "южный" полюса сферического тела (через первоначальную ось Х), а образующая составляет с осью некий угол [$966$], зависящий от r.
Возможно, Вам известно, как вычислить площадь этой части сферы.
У меня получилась формула S([$966$])=2*[$960$]*r2*(1- cos([$966$])).
Кстати, при [$966$]=2*[$960$] получаем хорошо знакомую формулу полной поверхности сферы.
При этом cos( [$966$] )= r/(2*R), где R - радиус сферического тела.
Проверив эти две формулы, можно их затем подставить в вышеуказанный интеграл и получить решение задачи. Ниже кратко излагается, как эти формулы получены.
1) Для S([$966$]). Желательно сделать рисунок.
По сути, сначала мы ищем площадь очень тонкого кольца, вырезаемого на сфере для угла между [$629$] и [$629$]+d[$629$]. "Ширина" такого кольца r*d[$629$], где [$629$] - это угол, который, по сути, тот же угол [$966$], но, поскольку в пределах интегрирования будет стоять [$966$], то переменную интегрирования назвали [$629$]. При этом [$629$]MAX=[$966$]. Длина окружности кольца равна 2*[$960$]*r*sin([$629$]).
А потом площади всех этих бесконечно тонких колец складываем (интегрируем).
То есть [$8747$] (2*[$960$]*r*sin([$629$]) )*r*d[$629$] с пределами от 0 до [$966$].
2) Как находим [$966$]=[$966$](r). Делаем рисунок окружности радиуса R (сферическое тело). Из центра окружности проводим в "южный полюс" (в нижнюю точку окружности) радиус, а затем из "южного полюса" окружности проводим хорду длиной r. Необходимо найти угол между хордой и радиусом, вернее косинус этого угла.
У меня получилось cos( [$966$] )= r/(2*R)
Проверьте мои расчеты. Вроде бы, все.