plt.plot(x, y0, 'k', label='True model')