n:10;

p:0.3;

q:1-p;

m:n*p;

s:sqrt(n*p*q);

B(x):=binomial(n,x)*p**x*q**(n-x);

load(distrib)$

N(x):=pdf_normal(x,m,s);

wxplot2d([B(x),N(x)], [x,0,n])$

b:makelist(floor(B(x)*1000+0.5)/1000.0,x,0,n);

n:makelist(floor(N(x)*1000+0.5)/1000.0,x,0,n);

print("Vergleich =",transpose(b),transpose(n))$