  The critical points of the 3-states two-layer Potts model on square lattice
for different interlayer couplings (Kx, Ky,and Kz) are calculated with high
precision using probabilistic cellular automata with Glauber algorithm, where
Kx and Ky are the nearest-neighbor interactions within each layer in the x and
y directions, respectively and Kz is the interlayer coupling. The obtained
results are 0.726, 0.807, 0.928, 0.987 and 1.00 for the cases Kz per Kx = 1,
0.5, 0.1, 0.01 and 0.001, respectively. Then, the critical curve has been
constructed for this model.
