Сколько чисел после десятичной точки вы можете показать, используя R?

Question

Мне было интересно, как далеко вы можете напечатать это:

http://www.iheartchaos.com/post/16393143676/fun-with-math-dividing-one-by-998001-yields-a

В R.

например: 1/998001

Answer 1

Вы можете использовать пакет mpc ,

# 3000 is the precision in bits
> mpc(1, 3000) / mpc(998001, 3000)
[1] "(1.0020030040050060070080 ...
...

Answer 2

Вот решение, которое не требует каких-либо модных пакетов, только база R и воспоминание об арифметике начальной школы:

longdiv <- function(num,den, limit=3000, spaces=TRUE) {
    i <- 0
    if( num < den ) {
        cat('0.')
        num <- num*10
    } else {
        cat( num %/% den, '.', sep='')
        num <- (num %% den) * 10
    }
    while( i < limit ) {
        i <- i + 1
        if( num < den ){
            cat('0')
            num <- num*10
        } else {
            cat( num %/% den )
            num <- (num %% den) * 10
        }
        if(spaces){
            if( (i %% 30) == 0 ) {
                cat('\n')
            } else if( (i %% 3) == 0 ) {
                cat(' ')
            }
        }
        flush.console()
    }
    cat('\n')
}   


longdiv(1,998001)