Nastavení okrajových podmínek rovnice – Python – Fórum – Programujte.com

Dobrý den, pomohl by mi někdo prosím s nastavením hodnot pro tuto rovnici? Výsledek by měl vypadat podobně jako na obrázku. Každý graf je pro jiné u_0. Jeden by měl mít 2 kořeny a druhý 4 (průsečík v 1). Moc děkuji

Připojen obrázek.

import matplotlib.pyplot as plt
import math
import numpy as np
fig, (ax1, ax2) = plt.subplots(nrows=2, ncols=1, figsize=(5,5))

# Equation - dispersion relation 
omega_p = 500
k       = 1
u_0     = 5
u_02    = 10
m       = 9
M       = 2
ksi     = np.arange(0.00001, 0.001, 0.00003)
omega   = 1/(ksi*omega_p)
ksi_0   = k*u_0/omega_p
ksi_02  = k*u_02/omega_p
F       = (m/M)/ksi**2 + 1/(ksi-ksi_0)**2
F2      = (m/M)/ksi**2 + 1/(ksi-ksi_02)**2

# Fig 1
ax1.plot(ksi, F)
# Fig 2
ax2.plot(ksi, F2)
plt.show()