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