Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import random
- from math import log
- l_1 = 0.1
- l_2 = 0.2
- T = 10
- N = 1000
- def get_random(l_):
- return - (log(random.uniform(0, 1)) / l_)
- time = []
- m = 0
- for _ in range(0, N):
- tmp = min(get_random(l_1), get_random(l_2))
- if tmp < T:
- time.append(tmp)
- else:
- m += 1
- print("Среднее время выживания: {}".format(float(sum(time)) / max(len(time), 1)))
- print("Вероятность что не сбит: {}".format(m / N))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement