cálculo do aniversário

aniversario <- function(n)   {  
  return(p <- 1 - choose(365, 365 - n) * factorial(n) / 365 ^ n)
}

resp23 = aniversario(23)
resp23

[1] 0.5072972

resp50 = aniversario(50)
resp50

[1] 0.9703736

resp10 = aniversario(10)
resp10

[1] 0.1169482