diff --git a/TBit.cabal b/TBit.cabal
--- a/TBit.cabal
+++ b/TBit.cabal
@@ -10,7 +10,7 @@
 -- PVP summary:      +-+------- breaking API changes
 --                   | | +----- non-breaking API additions
 --                   | | | +--- code changes with no API change
-version:             0.4.2.1
+version:             0.4.2.2
 
 -- A short (one-line) description of the package.
 synopsis:            Utilities for condensed matter physics tight binding calculations.
@@ -77,7 +77,7 @@
   -- other-modules:       
   
   -- Other library packages from which modules are imported.
-  build-depends:       base ==4.8.*, containers ==0.5.*,
+  build-depends:       base >=4.7 && < 4.9, containers ==0.5.*,
                        hmatrix ==0.16.*,
                        mtl ==2.2.*, parallel ==3.2.*, deepseq ==1.4.*,
                        integration ==0.2.*, numeric-tools ==0.2.*, fgl ==5.5.*,
diff --git a/TBit/Systems/HoneycombLattice.hs b/TBit/Systems/HoneycombLattice.hs
--- a/TBit/Systems/HoneycombLattice.hs
+++ b/TBit/Systems/HoneycombLattice.hs
@@ -57,8 +57,30 @@
                        ky = (k!1) :+ 0.0
                        y  = realPart ky
                     in (4 >< 4) 
-            [ hz + v - 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) + 2*soc*sin(kx) , 0, t + (2*t*cos(kx/2))/cis(0.5*sqrt(3)*y), r*(negate iC + (iC*(cos(kx/2) + sqrt(3)*sin(kx/2)))/cis(0.5*sqrt(3)*y))
-            , 0, -hz + v + 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) - 2*soc*sin(kx), r*(negate iC + (iC*(cos(kx/2) - sqrt(3)*sin(kx/2)))/cis(0.5*sqrt(3)*y)), t + (2*t*cos(kx/2))/cis(0.5*sqrt(3)*y) 
-            , t + 2*cis(0.5*sqrt(3)*y)*t*cos(kx/2), r*(iC - (negate iC)*cis(0.5*sqrt(3)*y)*(cos(kx/2) - sqrt(3)*sin(kx/2))), -hz - v + 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) - 2*soc*sin(kx), 0
-            , r*(iC - iC*cis(0.5*sqrt(3)*y)*(cos(kx/2) + sqrt(3)*sin(kx/2))), t + 2*cis(0.5*sqrt(3)*y)*t*cos(kx/2), 0, hz - v - 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) + 2*soc*sin(kx) ]
+
+            -- NBs: cis(t) === exp(i*t)
+            --    : iC === i === sqrt(-1)
+            -- Row 1:
+            [ hz + v - 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) + 2*soc*sin(kx) 
+            , 0
+            , t + (2*t*cos(kx/2))/cis(0.5*sqrt(3)*y)
+            , r*(negate iC + (iC*(cos(kx/2) + sqrt(3)*sin(kx/2)))/cis(0.5*sqrt(3)*y))
+
+            -- Row 2:
+            , 0
+            , -hz + v + 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) - 2*soc*sin(kx)
+            , r*(negate iC + (iC*(cos(kx/2) - sqrt(3)*sin(kx/2)))/cis(0.5*sqrt(3)*y))
+            , t + (2*t*cos(kx/2))/cis(0.5*sqrt(3)*y) 
+
+            -- Row 3:
+            , t + 2*cis(0.5*sqrt(3)*y)*t*cos(kx/2)
+            , r*(iC - iC*cis(0.5*sqrt(3)*y)*(cos(kx/2) - sqrt(3)*sin(kx/2)))
+            , -hz - v + 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) - 2*soc*sin(kx)
+            , 0
+
+            -- Row 4:
+            , r*(iC - iC*cis(0.5*sqrt(3)*y)*(cos(kx/2) + sqrt(3)*sin(kx/2)))
+            , t + 2*cis(0.5*sqrt(3)*y)*t*cos(kx/2)
+            , 0
+            , hz - v - 4*soc*cos((sqrt(3)*ky)/2)*sin(kx/2) + 2*soc*sin(kx) ]