integrate ((cos^2(u) * cos^2(v) *sin^2(u) + cos^2(u) * sin^2(v) *cos^2(u) * cos^2(v) + sin^2(u) * cos^2(u) *sin^2(v))*cos(u) ) du dv, -pi/2<=u<=pi/2, 0<=v<= 2*pi