plt.gca().set_xlim(xmin - margin * xdiff, xmax + margin * xdiff)