Context
I am studying electrodynamics using Zangwill [1]. I have a narrow question regarding a proposition found there in. In the section on the Drude model of conducting matter that applies to particles with mean travel time $\tau$, a characteristic plasma frequency $\omega_p$, the Fourier transform of the permittivity, $\hat{\epsilon}(\omega)$, the Fourier transform of the permiability, $\hat{\mu}(\omega)$, and frequency parameter $\omega$. Adjusting from what Zangwill writes in [1], the proposition is that
In the high frequency limit when $$\frac{\hat{\epsilon}(\omega)}{\epsilon_o}\approx 1 - \frac{\omega_p^2}{\omega^2} \quad (\omega\,\tau \gg 1)$$ is valid, the dispersion relation for transverse waves $$k(\omega) = \omega\sqrt{\hat{\mu}(\omega)\,\hat{\epsilon}(\omega)} = \frac{\omega}{c}\,\hat{n}(\omega)$$ reads $$ \boxed{\omega^2 = \omega^2_p + c^2\,k^2 .}\tag{1} $$ I have made several attempts to maniupulate these and other equations to obtain (1). Yet, I come up short.
Question
The narrow question here is, how does one derive the high-frequency (or collisionless) dispersion relation, i.e., Eq.(1), for transverse waves using the Drude model?
Bibliography
Zangwill, Modern Electrodynaics, 2013 pp. 630-632 .