Лучше в ЦСК
x=r*cos(fi);
y=r*sin(fi);
3=<z=<sqrt[16-r^2]
0=<fi=<2*pi;
0=<r=<sqrt[7];
dxdydz=rdrd(fi)dz
Посчитаете, получится должно V1=11*pi/3, так как на отрицательных зет, -sqrt[16-x^2-y^2]=<z=<-3
то V=2*V1