Dit script bevat twee functies voor het berekenen van binomiale kansen: binomialpdf(k,N,p) geeft p(X=k) en binomialcdf(k,N,P) geeft p(X≤k) voor N-herhalingen en een succeskans p.
from math import * def factorial(n): n=int(n) if n<0: return None res=1 for i in range(1,n+1): res=res*i return res def n_choose_k(n,k): return int(factorial(n)/(factorial(k)*factorial(n-k))) def binomialpdf(k,N,p): return n_choose_k(N,k)*(p**k)*((1-p)**(N-k)) def binomialcdf(k,N,p): res=0 for i in range(k+1): res=res+binomialpdf(i,N,p) return res