Neste trabalho simulamos computacionalmente o modelo ± J numa rede quadrada para o cálculo dos mínimos e máximos locais de energia. Para tal, fizemos o uso de variáveis Booleanas para armazenar as configurações variáveis de spins e as configurações fixas de ligações ferro e antiferromagnéticas entre sítios vizinhos. Definindo uma dinâmica conveniente dos spins na rede, estudamos o gráfico da Energia × S, onde S representa um elemento do espaço de estados de dimensão N = L². Utilizamos redes de dimensões L = 8, 16, 32 e 64

In this work we have simulated numerically the ± J Ising model in a square lattice to obtain the minimum and maximum local energies. We used Boolean variables to store the variable configuration of the spins and the fixed configuration of ferro and antiferromagnetic bonds linking nearest neighbors. We have defined a convenient dynamic of the spins on the lattice and we analyzed the behavior of this dynamic by plotting Energy × S, where S represents an element of the associated space of states with dimension N = L². We take lattice with linear sizes L = 8, 16, 32 and 64