Caro Luiz O Cesar Rabak já lhe deu uma excelente dica de leitura Compreenda qual o limite de representação numérica da sua máquina Experimente os limites da sua máquina Estude as bibliotecas gmp e Rmpfr https://en.wikipedia.org/wiki/Floating_point options(digits = 22) .Machine ## as.numeric("14472334024676221") # [1] 14472334024676220 ## se sua máquina representasse esse número, não seria necessário usar a biblioteca Rmpfr as.numeric("-0.999999999999999999999999998975999999999999914663505297831951000646311908486498") # [1] -1 ## epsilon da maquina eps = 1 while(eps + 1 > 1){eps = eps/2} (eps = eps*2) ## 1 - 0.00000000000000001 1 - 0.0000000000000001 p <- mpfr(c(500,600,700)^(-10), precBits=255) - mpfr(1,precBits=255) class(p) str(p) bons estudos!