Quanti numeri dopo la virgola decimale puoi mostrare usando R?

Question

Mi chiedevo quanto lontano si può stampare questa:

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

In R.

esempio: 1/998001

Answer 1

è possibile utilizzare il pacchetto mpc,

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

Answer 2

Ecco una soluzione che non richiede alcun pacchetto di fantasia, basta basare R e un flashback per eleme aritmetica della scuola media:

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)