Не по сфере, а по шару. Если бы тоже самое было с поверхностным интегралом по сфере, тоже ноль бы получился.
А ноль, потому что область интегрирования симметрична относительно плоскости z=0 и функция под интегралом нечётна по z.
Это и у Вас видно - в интеграле, который в середине, у Вас должно получиться

в пределах от

до

, что и даст этот заранее ясный ноль. Замена

на

в верхнем пределе интегрирования даст верный результат.
Маленькое замечание. А стоило ли вообще здесь переходить к сферическим координатам? Если сразу через повторные, не гораздо ли проще получится?
