目標：
從一個分布服從指數分布的隨機變量中去抽取1000個變量，問這1000個變量中至少有20個變量數值大於18的概率
其中指數分布的 lamb = 0.2
import numpy as np
'''
從一個分布服從指數分布的隨機變量中去抽取c個變量，問這c個變量中至少有t個變量數值大於h的概率
其中指數分布的 lamb = 0.2
'''
def ele_function(n, lamda, t, h, c): # 迭代次數
sum = 0
for i in range(n):
x = np.random.exponential(1 / lamda, c)
x_n = float(len(x[x >= h]))
if(x_n > t):
sum += 1
return sum / n
print(ele_function(100000, 0.2, 20, 18, 1000))
