p <- seq(0,1,0.001)
v <- p*(1-p)
plot(p,v,type="l",xlab="Probabilita'",ylab="Varianza di una v.c. B(p)")
lines(c(0.5,0.5),c(0,0.25),lty=3,lwd=2)