messyscript/binomial.py

#二项分布
from fractions import Fraction
from math import comb
n=3
p=Fraction(1,5)
print("EX=",n*p)
print("DX=",n*p*(1-p))
for i in range(n+1):
    print("P(X=",i,")=",comb(n,i)*pow(p,i)*pow(1-p,n-i))