cdesante

#The birthday Problem.

#Simple computation:
N <- 1:50
pshare <- 1- exp(-(N*(N-1)/2)/365)
cbind(N, pshare)

#just want one value:
bday <- function(X){
N <- X