numpy计算服从指数分布的概率

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