
# comandi per analizzare le prime 1000 cifre decimali di pi greco

# file generato con 
#  echo "scale=10000; 4*a(1)" | BC_LINE_LENGTH=0 bc -l > pi1000.txt 

pi <- read.table('pi1000.txt')   # Non funziona!

file.info('pi1000.txt')
fi <- file.info('pi1000.txt')
str(fi)

fileName = 'pi1000.txt'
readChar(fileName, file.info(fileName)$size)

pi1000 <- readChar(fileName, file.info(fileName)$size - 1)
ch.pi1000 <- strsplit(pi1000, split='')
str(ch.pi1000)

dec.pi = as.numeric( ch.pi1000[[1]][-(1:2)] ) 
barplot(table(dec.pi), col='cyan')


