Current implementation uses ``kwant.physics.Bands``. This example fails
when one substitute ``k_x**2`` with for example ``k_x**4`` due to farther
then nearest neighbours hoppings. Using wraparound provides also easier way
to extend it to multidimensional band structures (except of plotting).