Написать код на питон

Написать код на пито
Исследовать средний удельный доход в единицу времени ? в системе массовог
обслуживания М|М|4|2, где интенсивность входящего потока ?=9 час-
, интенсивност
обслуживания на одном приборе µ=2 час-
, доход от обслуживания одной заявки составляе
с1=6 у.е., плата в единицу времени за простой одного прибора с2=2 у.е./час, плата в единиц
времени ожидания заявки в очереди равна с3=1 у.е./час, штраф за потерю заявки с4=3 у.е
случайную величину моделировать как
a = [9, 11, 13, 15, 17, 17, 8
где i = 0, если все каналы свободны
1, если 1 занят
2, если 2 занято
3, если 3 занято
4, если 4 занято
5, если 4 занято, 1 в очеред
6, если 4 занято, 2 в очеред
ответ после реализации кода должен быть ~32,3
def ksi(i)
gamma = np.random.uniform(0, 1
return -1/a[i] * np.log(gamma
Оценить средний удельный доход в единицу времени с помощью компьютерног
моделирования. Для этого построить имитационную модель функционирования указанно
системы массового обслуживания как реализацию траектории введенного случайног
процесса ?(t). Смоделировать n=1000 реализаций на интервале t?(0;1000), для каждо
реализации найти удельный дохода в единицу времени за 1000 часов и по полученны
данным построить оценку для ?.