Physics and Mathematics

Create a 1-D array of 20 elements. At each element generate a +1 or a -1 spin. Keep track of how many +1 spins you have. Determine the total energy of the system using Eij = SiSj. Assume periodic boundary conditions, so the first and last element are assumed to be adjacent.

Run this simulation many times (each time randomly generating spin), and keep track of energy vs number of +1 spins.

Create a short report of what you did, what you found and what it means. The report should include your (well commented) code.