Arithmetic-geometric mean/Calculate Pi: Difference between revisions

Content added Content deleted
Line 1,211: Line 1,211:
# Compute pi to one million decimal digits:
# Compute pi to one million decimal digits:
p <- 3322000
p <- 3322000
x <- agm(20, p)
roundMpfr(agm(20, p) - Const("pi", p), 64)

# Check
roundMpfr(x - Const("pi", p), 64)
# 1 'mpfr' number of precision 64 bits
# 1 'mpfr' number of precision 64 bits
# [1] 4.90382361286485830568e-1000016</lang>
# [1] 4.90382361286485830568e-1000016

# Save to disk
f <- file("pi.txt", "w")
writeLines(formatMpfr(x, 1e6), f)
close(f)</lang>


=={{header|Racket}}==
=={{header|Racket}}==