Ecco come studiare la legge di distribuzione del "n° di biscotti difettosi in una confezione" considerata alla voce "Leggi di distribuzione (variabili discrete)": la poligonale che la rappresenta graficamente, il calcolo della media (#m) e della varianza (#v), con i comandi eseguiti automaticamente (copiare le seguenti righe e clicca [CLP]). Cambiando #n e #p si può studiare qualunque binomiale.

n
punt=2
f(x)=binom
#n=6 _ #p=1/8
plot x:0..#n n=#n y:f|14
oo
-
h(x)=x*f(x) _ #m=[0,#n] h sum
h(x)=(x-#m)^2*f(x) _ #v=[0,#n] h sum

Ecco cosa battere in R:
n <- 6; p <- 1/8
binom <- function(m) choose(n,m)*p^m*(1-p)^(n-m)
plot(binom,0,6,n=7,type="p")
abline(h=seq(0.1,0.4,0.1),lty=3, col="grey")
abline(v=0,h=0,col="blue")
plot(binom,0,6,n=7,add=TRUE,type="h",col="red")
plot(binom,0,6,n=7,add=TRUE,lty=3,col="red")
x <- seq(0,6,1); media <- sum(x*binom(x)); media
  0.75