Este script contém 2 funções para calcular probabilidades binomiais: binomialpdf(k,N,p) devolve p(X=k) e binomialcdf(k,N,P) devolve p(X≤k) para N repetições e uma probabilidade de sucesso 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