Точки пересечения графиков

import numpy as np
import matplotlib.pyplot as plt
import math
t_obs = 5
lambda_val = 0.2
n_values = np.arange(1, 21)
alpha = lambda_val / t_obs
factorials = [math.factorial(n) for n in n_values]
factorials = [1 / f for f in factorials]
 
P0 = np.sum([(alpha ** k) * factorials[k - 1] for k in n_values])
Pn = (alpha ** n_values) * factorials
P_obs = 1 - Pn
Nk = alpha * (1 - Pn)
Kz = Nk / n_values
графики P_obs и Kz.
intersections = np.argwhere(np.diff(np.sign(P_obs - Kz)))
if intersections.size > 0:
    n_optimal_values = n_values[intersections]
    print("Оптимальные количества процессоров:", n_optimal_values)
else:
    print("Не найдено точек пересечения графиков.")
plt.figure(figsize=(10, 5))
plt.plot(n_values, P_obs, label='P_obs')
plt.plot(n_values, Kz, label='Kz')
if intersections.size > 0:
    plt.scatter(n_optimal_values, P_obs[intersections], color='red', marker='x',
                label='Оптимальные количества процессоров')
plt.xlabel('Количество процессоров (n)')
plt.ylabel('Значение')
plt.title('Графики P_obs и Kz от количества процессоров')
plt.legend()
plt.grid()
plt.show()