Three-loop β-function for the bosonic non-linear sigma model