dbenford <- function(x){
log10(1 + 1/x)
}
pbenford <- function(q){
cumprobs <- cumsum(dbenford(1:9)) return(cumprobs[q])
}
dbenford <- function(x){
log10(1 + 1/x)
}
pbenford <- function(q){
cumprobs <- cumsum(dbenford(1:9)) return(cumprobs[q])
}