Matrix H in Subspace 26 with Dimension 1

H is the Hamiltonian. Here we set µ=0. For grand-canonical calculations add -6µ to the main diagonal.
-2*(h - U)