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