import matplotlib.pyplot as plt
import numpy as np


print("Введите чему равна лямбда:")
lymbda = int(input())
print("Введите степень эпсилона(10^(-x):")
epsilon = 10 ** (-int(input()))
print(epsilon)
print("Введите чему равен мю:")
mu = int(input())
t = 0
e =  2.718281828459045
deltat = 0.001
#Вычисление чему равно время, когда система входит в стационарный режим
p = 1 / (lymbda + mu) * (mu + lymbda * e ** (-(lymbda + mu) * t))
while abs(1 / (lymbda + mu) * (mu + lymbda * e ** (-(lymbda + mu) * (t + deltat))) - 1 / (lymbda + mu) * (mu + lymbda * e ** (-(lymbda + mu) * t))) / (1 / (lymbda + mu) * (mu + lymbda * e ** (-(lymbda + mu) * t))) >= epsilon:
    t = t + deltat
print("Ответ:",t)



t = np.linspace(0, 4, 100)                     # от 0 до 4 сделать 100 точек
p0 = (1 / (lymbda + mu)) * (mu + lymbda * e ** (-(lymbda + mu) * t))                         # y1 - тоже много точек
p1 = (lymbda / (lymbda + mu) * (1 - e ** (-t * (lymbda + mu))))
fig, ax = plt.subplots()                        # будет 1 график, на нем:
ax.plot(t, p0, color="blue", label="p(t)")      # функция p0(t), синий, надпись p(t)
ax.plot(t, p1, color="blue", label="p(t)")      # функция p0(t), синий, надпись p(t)
ax.set_xlabel("t")                              # подпись p горизонтальной оси t
ax.set_ylabel("p0")                              # подпись p вертикальной оси p
ax.legend()                                     # показывать условные обозначения

plt.show()                                      # показать рисунок