Ce script contient deux fonctions pour le calcul des probabilités binomiales : binomialpdf(k,N,p) renvoie p(X=k) et binomialcdf(k,N,P) renvoie p(X≤k) pour N répétitions et une probabilité de succès 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