我正在为 Ising 模型使用 2D 到 5D 的 Metropolis 算法，并且我想计算自旋-自旋相关函数。

$$c(r)=<s_is_r>−〈s_i〉〈s_r 〉$$ 但我不确定如何估算

$$<s_is_r> = \frac{1}{N}\sum s_is_r$$

我能做什么$s_is_r$数字上？