function [connections] = initialize_hamiltonian(N)
connections = rand(N,N)-1/2;  % Random matrix, each entry between -1/2 and 1/2
connections = (connections + connections');% Make the matrix symmetric
                                           % Not necessary, but physically
                                           % a very sensible thing...
    