Three-loop β-function for the two-dimensional non-linear σ-model with a Wess-Zumino-Witten term