plt.plot(w, S1(w))