# (c) Rory Clune
# UNFINISHED

from numpy import *

nbits = 100
bonds = random.randn(nbits)
min_energy = -sum(fabs(bonds))
nsteps = 300
population = 100
beta = array([10])#, 1, 0, 0.1, 0.01])

spins = zeros((population, nbits))
# energy = zeros((
for temp in range (0,beta.size):
	for n in range (0, population):
		spins[n,:] = (random.rand(nbits) > 0.5)*2-1
	for t in range(0,nsteps):
		energy