The question is that given in the title;
Prove that all roots of the equation $z\tan z = k$ lie in $\Bbb R$, where $k$ is a positive non-zero real number.
I attempted a somewhat "brute force" approach by simply letting $z = a + bi$, with $a, b \in \Bbb R$. Under this assumption, if $a = 0$ (that is, $z$ is purely imaginary), we come to the conclusion that $k$ must be of the form
$$-b\tanh b = k.$$
If $b > 0$ then $\tanh b \in (0,1]$ and so $-b\tanh b < 0$. If $b < 0$ then $\tanh b \in [-1, 0)$ since $\tanh$ is odd, and so again, $-b\tanh b < 0$, and so for $k \in \Bbb R_{>0}$, we cannot have purely imaginary $z$.
Consider then the case that $z = a + bi$ where $a, b \neq 0$. Then after some manipulation, we end up with
$$z\tan z = \frac{\tan a + i\tanh b}{1 - i\tan a \tanh b}.$$
Since both $\tan$ and $\tanh$, with domain restricted to $\Bbb R$, also have codomain $\Bbb R$, this ratio must be complex, since $\tanh b = 0$ iff $b = 0$, which we excluded under assumption, and so $k$ is complex if $z$ is complex.
The only remaining possibility is that $z \in \Bbb R$, the result has been proven.
I think this is somewhat hand-wavy at best; can anyone give me a more satisfactory/elegant proof of this result?