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