Abstract
A general discussion on applying the cluster variation method to non-uniform Ising and Heisenberg models is given. In particular, a general formula is provided for obtaining the linear response to increases of the external fields and the parameters of interactions for the uniform Ising and Heisenberg models in the cluster variation method. As an illustration, it is shown that the expressions of the spin-pair correlation function in the pair and the square approximation are readily derived by the formula, for the Ising models on the square and simple cubic lattices.