1538.sagews

var('r,theta,phi')
integrate(integrate(r^2*sin(theta),theta,(90-30)/180*pi,(90+30)/180*pi),phi,0,2*pi)/(4*pi*r^2)

1/2