# Задаем параметры
y = w2

# Создаем массив для аргумента x
x = np.linspace(-40, 40, 1000)

# Вычисляем функцию z
z = np.sqrt(1/2*(x**2 + y**2 + np.sqrt((y**2 - x**2)**2 + 4*a**2)))

# Строим график
plt.figure(figsize=(12, 6))
plt.plot(x, z)
plt.xlabel('w1')
plt.ylabel('Ω')
plt.title('График функции Ω')
plt.show()