$${n=1+\frac{Ne^2}{2\varepsilon_0 m_e \omega_0^2} \frac{1}{1 - \omega^2 / \omega_0^2}= A_\omega+B_\omega \omega^2}$$