with plt.style.context(pubStyle):
    fig, ax = plt.subplots(1,1,figsize=(10,7))
    for h in hevolve:
        ax.plot(h[0], h[1], 'ko')
        
    line = lambda x, m, b:  m*x + b
    fit, covar = curve_fit(line, hevolve[:, 0], hevolve[:, 1])
    X = np.linspace(hevolve[:, 0].min(), hevolve[:, 0].max())
    ax.plot(X, line(X, *fit), 'k--')
    ax.set_xlabel("$\Delta$ Helium Mass Fraction", fontsize=20)
    ax.set_ylabel("Turn off offset from mean [mag]", fontsize=20)
    ax.annotate(f'$O(Y)={fit[0]:0.2f}Y+{fit[1]:0.2f}$', xy=(-0.04, 0.30), color='black', fontsize=20)

    fig.savefig("Figures/HeliumMeanOffset.pdf")