plt.gca().set_ylim(ymin - margin * ydiff, ymax + margin * ydiff)