bday <-
function(n){
    bday = 1
    nm1 = n - 1
    for(j in 1:nm1){
       bday = bday*((365-j)/365)
    }
    bday <- 1 - bday
    return(bday)
}