Differentiating the mean-field free energy and imposing stationarity gives
The hyperbolic-function identity turns this into the self-consistency equation