Your scenario describes the extension of the radius from 0 to 1, with theta wrapping one half of the x-y plane. The reason why you're then getting 0, is that you're going -pi/2 to pi/2. I.E. a negative quadrant to a positive quadrant, and because you're not allowing for symmetry, the negative and positive are cancelling to produce the nullity. Also, for phi, the range is 0 to pi dropping DOWN from the z pole.
I was just discussing the correct way to do this with the esteemed Ahmad (read: he essentially told me how it is

)
Phi makes the angle with the vertical, yeah? So to drop down from the vertical z-pole to the x-y plane and not cross it you need an angle of pi/2 only (0, pi/2). To sweep the volume of the ball you then need to rotate around the whole x,y plane 0 to 2pi.
So it should be
0 < r < 1
0 < theta < 2pi
0 < phi < pi/2
I think they (solutions) have their phi and their theta swapped around.
P.S. Suitcase packing is coming along nicely

I have a lotttt of material to take in, I've realised.