from scipy.stats import hypergeom

n1 = 21; n2 = 6
k1 = 15; k2 = 1
alpha1 = 0.05

kl = int(hypergeom.ppf(alpha1, n1 + n2, k1 + k2, n2))
kr = int(hypergeom.ppf(1.0 - alpha1, n1 + n2, k1 + k2, n2))

print("kl =", kl, " k2 =", k2, " kr =", kr)
print("k2<=kl =", k2<=kl)
print("k2>=kr =", k2>=kr)

#2020.06.23 10:13:49.184	Python	kl = 2;  k2 = 2;  kr = 6
#2020.06.23 10:13:49.184	Python	k2<=kl = True
#2020.06.23 10:13:49.184	Python	k2>=kr = False