Consistent MFA correlation functions of Ising models for all temperatures

A correct calculation of the Ising model correlation function C(q) = 〈(S(q) ? 〈S(q)〉) (S(-q) ? 〈S(-q)〉)〉 in the MFA results in

$\text{C(q)=}\text{〈S}{}^2\text{〉?〈〉}{}^2\text{1?(〈S}{}^2\text{〉?〈S〉}{}^2\text{βJ(q}\text{1}\text{N}\text{∑}\text{q}\text{1}\text{1?〈S〉}{}^2\text{βJ(q)}{}^{\mathrm{?1}}\text{C(q) fulfills the exact sum rule N}{}^{-1}\text{Σ}{}_{\mathrm{q}}\text{C(q) = 〈S}{}^2\text{〉 ? 〈S〉}{}^2$

. Previous literature supposed a violation of this sum rule to be a characteristic disadvantage of this approximation.