Calculation of The Critical Temperature for Anisotropic Two-Layer Ising   Model Using The Transfer Matrix Method
  A new finite-size scaling approach based on the transfer matrix method is developed to calculate the critical temperature of anisotropic two-layer Ising ferromagnet, on strips of r wide sites of square lattices. The reduced internal energy per site has been accurately calculated for the ferromagnetic case, with the nearest neighbor couplings Kx, Ky (where Kx and Ky are the nearest neighbor interactions within each layer in the x and y directions, respectively) and with inter-layer coupling Kz, using different size-limited lattices. The calculated energies for different lattice sizes intersect at various points when plotted versus the reduced temperature. It is found that the location of the intersection point versus the lattice size can be fitted on a power series in terms of the lattice sizes. The power series is used to obtain the critical temperature of the unlimited two-layer lattice. The results obtained, are in good agreement with the accurate values reported by others. 
M. Ghaemi, M. Ghannadi, B. Mirza