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