The commutation of matrices in (1) allow them to be diagonalised, and thus eigenvalues can be found. The partition function is calculated from the maximal eigenvalue, resulting in a free energy per site of