This is the Blume–Capel model. In the mean-field approximation, write and replace
Because every site has coordination number , the resulting energy is
The one-site partition function is therefore
With , the mean-field partition function is . Taking gives