Les dejo un script de Matlab en el que nos vamos a basar para la explicación de hoy: ising2D0.