b[n_, p_, i_]:= Binomial[n, i]p^i(1-p)^(n-i)(*Probability of i heads if a biased coin (P(heads)=p) is flipped n times.*)

bd[n_, p_]:=Table[{i, b[n,p, i]}, {i, 0, n}];(*Binomial probability mass function, displayed as a table.*)

bdn[n_, p_]:=Table[{i, b[n,p, i], N@b[n,p, i]}, {i, 0, n}];(*Same as previous, but with a column showing probabilities in decimal notation.*)S

TableForm[bdn[3, .6]]

TableForm[bdn[100, .6]]

Manipulate[ListPlot[bd[n, p], Filling->Axis, AxesOrigin->{-1,0}, PlotRange->{0,1}], {n, 1, 100,1}, {p, .01, .99}]