plt.plot(x, mu, 'g', x, plo, 'r:', x, pup, 'r:')