def METHOD_NAME(self, epsilon=1e-4, color_scheme=None):
    """Return pyplot of band structure calculation results."""
    import matplotlib.pyplot as plt
    fig, axarr = plt.subplots(2, 1)
    for ax in axarr:
        ax.xaxis.set_ticks_position("both")
        ax.yaxis.set_ticks_position("both")
        ax.xaxis.set_tick_params(which="both", direction="in")
        ax.yaxis.set_tick_params(which="both", direction="in")
        self._band_structure.plot(axarr, epsilon=epsilon, color_scheme=color_scheme)
    return plt