目标：
从一个分布服从指数分布的随机变量中去抽取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))
