In the classical Heisenberg model each spin can take any orientation in space (recall that the classical partition function of a single spin in an external field B is zclassical=∫e−βBScosθsinθdθdϕ). Show that, in mean field theory,
⟨Sz⟩=Scoth[β(B^+B)S]−[β(B^+B)S]S