n:20;

p1:0.1;

p2:0.3;

B(p,k):=binomial(n,k)*p**k*(1-p)**(n-k);

wxplot2d([B(p1,x),B(p2,x)], [x,0,n], [y,0,1])$