def factorielle(n): if n==0: return 1 else: return n*factorielle(n-1) def C(n,k): return factorielle(n)//(factorielle(k)*factorielle(n-k)) def binomFdp(n,p,k): return C(n,k)*p**k*(1-p)**(n-k)