v <- function(p) p*(1-p)
s <- function(p) sqrt(p*(1-p))
h <- function(p) -(p*log(p) + (1-p)*log(1-p))
p <- (0:1000)/1000
plot(s(p) ~ p, type="l", ylab="standard deviation")