Spontaneous magnetization of the Ising model on a layered square lattice

We have studied the Ising model on a layered square lattice with four different coupling constants and two different magnetic moments. The partition function at zero magnetic field is derived exactly. We propose a formula for the spontaneous magnetization which agrees with the exact low-temperature series expansion up to the 16th order and reduces to the exact result of Au-Yang and McCoy in a special case.