Abstract
In this paper, we prove that for $s\in(0,\infty)$, the harmonic mean of $E_k(s)$ and $E_k(1/s)$ is always less than or equal to $\Gamma(1-k,1)$. Where $E_k(s)$ is the generalized exponential integral function, $\Gamma(u,s)$ is the upper incomplete gamma function and $k\in \mathbb{N}$.