w = np.linspace(1.0,5.0,100) alpha = 0.01
y = np.power((w-3),2)
plt.plot(w,y,'r') plt.title(r'$y=(\omega-3)^2 $') plt.xlabel(r'$\omega$') plt.ylabel('y') plt.grid();