Updated heisenberg.cpp

N=1 is bulk now.
Correctly implemented default longrange 1/x^3 coupling