integrate(cos theta*sin theta,{theta,0,pi/2},{phi,0,2*pi})