Optimal estimates for the gradient of harmonic functions in the multidimensional half-space

A representation of the sharp constant in a pointwise estimate of the gradient of a harmonic function in a multidimensional half-space is obtained under the assumption that function’s boundary values belong to $L^p$. This representation is concretized for the cases $p=1,2,$ and $\mathrm{\infty }$.