YsA = Ys[:, 0]
sortedYsA = np.argsort(YsA)
plt.plot(YsA[sortedYsA], chi2[sortedYsA])
plt.xlim(