dbenford <- function(x){
log10(1 + 1/x)
}
pbenford <- function(q){
cumprobs <- cumsum(dbenford(1:9)) return(cumprobs[q])
[Validate]
Generated with the Darkfish Rdoc Generator 2.