 Solution

Use source spherical coordinates with radius $s$, polar angle $\alpha$ measured from the positive $z$-axis, and azimuth $\theta$. The lower half of the spherical shell has $a\le s\le 2a$, $\pi /2\le \alpha \le \pi$, and $0\le \theta \le 2\pi$. Its volume element is $s^2\sin \alpha dsd\alpha d\theta$. At the observation point $(0,0,h)$, $h>0$, the separation is $\sqrt{s^2+h^2-2sh\cos \alpha }$, soThe angular integral is elementary: put $q=\cos \alpha$ to obtainWe used $\sqrt{(s+h{)}^2}=s+h$ since both $s,h$ are positive. HenceThus the axial potential of a uniform hemispherical shell, at every point of the positive $z$-axis, isThere is no need to split into $h<a$, $a<h<2a$, and $h>2a$: the observation point is outside the lower hemisphere in all three cases. As checks, $\varphi (h)\to -3\pi G\rho a^2$ as $h\downarrow 0$, and at large $h$,The last factor is exactly the total mass of the lower half of the spherical shell.