gy.plot(hom.solution["epsilon"]["y"].real, geometry=lattice, ax=ax[1])