Poisson kernels on the half-plane are bell-shaped

Consider a second-order elliptic operator $L$ in the half-plane $\mathbb{R}\times (0,\mathrm{\infty })$ with coefficients depending only on the second coordinate. The Poisson kernel for $L$ is used in the representation of positive $L$-harmonic functions, that is, solutions of $Lu=0$. In probabilistic terms, the Poisson kernel is the density function of the distribution of the diffusion in $\mathbb{R}\times (0,\mathrm{\infty })$ with generator $L$ at the hitting time of the boundary. We prove that the Poisson kernel for $L$ is bell-shaped: its $n$th derivative changes sign $n$ times. In particular, it is unimodal and it has two inflection points (it is concave, then convex, then concave again).