BOX(x,y):=[[x-0.5,0],[x+0.5,0],[x+0.5,y],[x-0.5,y],[x-0.5,0]]

BWF(n,p,k):=COMB(n,k)*p^k*(1-p)^(n-k)

BVF(n,p,k):=SUM(BWF(n,p,i),i,0,k)

BVI(n,p,a,b):=SUM(BWF(n,p,i),i,a,b)

GRAPH(n,p):=VECTOR(BOX(k,BWF(n,p,k)),k,0,n)

mu(n,p):=n*p

sigma(n,p):=SQRT(n*p*(1-p))

INTERV(k,n,p):=[[[mu(n,p)-k*sigma(n,p),0],[mu(n,p)-k*sigma(n,p),1]],[[mu(n,p)~
,0],[mu(n,p),1]],[[mu(n,p)+k*sigma(n,p),0],[mu(n,p)+k*sigma(n,p),1]]]

N(x,m,s):=1/(s*SQRT(2*pi))*#e^(-0.5*((x-m)/s)^2)