“…In the other case, the model must be solved by approximation method or numerical method, such as the molecular field approximation (MFA), the Bethe-Peierls approximation (B-PA) [3], the Bethe lattice approximation (BLA) [4], in which the lattice is simplified as an infinite Cayley tree with the coordination number z, the effective field theory (EFT) [5], the coherent anomaly method (CAM) [6], the double-chain approximation (DCA) [7], which is a natural extension of the pair approximation, the cluster variation method (CVM) [8], the series expansion methods (SE) [1,9], the Monte Carlo technique (MC) [10], and so on. The transition temperatures obtained by the former seven approximation methods are higher than the exact value, the B-PA and BLA give the same transition temperatures for S = 1/2, and the DCA with the large cluster gives the transition temperature closing mostly to the exact value.…”