On répète n fois une expérience de Bernoulli.
La probabilité de succès est p.
X est le nombre de succès. X prend les valeurs k qui vont de 0 à n.
n = 5
p = 0.6
k = 0
from scipy.stats import binom
binom.pmf(k,n,p)
binom.cdf(k,n,p)
from scipy.stats import binom
from pylab import plt
import numpy as np
fig,ax = plt.subplots(1,1)
x = np.arange(n+1)
ax.plot(x,binom.pmf(x,n,p),'ro',ms=0)
ax.vlines(x,0,binom.pmf(x,n,p),colors='r',lw=10,alpha=0.5)
plt.show()
print("k","P(X=k)")
for k in range(n+1):
print(k,binom.pmf(k,n,p))