plt.plot(inds, mu1o, '*', label='mu1o')