#! /usr/bin/env python # Print digits of pi forever. # # The algorithm, using Python's 'long' integers ("bignums"), works # with continued fractions, and was conceived by Lambert Meertens. # # See also the ABC Programmer's Handbook, by Geurts, Meertens & Pemberton, # published by Prentice-Hall (UK) Ltd., 1990. import sys def main(): k, a, b, a1, b1 = 2, 4, 1, 12, 4 while True: # Next approximation p, q, k = k*k, 2*k+1, k+1 a, b, a1, b1 = a1, b1, p*a+q*a1, p*b+q*b1 # Print common digits d, d1 = a//b, a1//b1 while d == d1: output(d) a, a1 = 10*(a%b), 10*(a1%b1) d, d1 = a//b, a1//b1 def output(d): # Use write() to avoid spaces between the digits sys.stdout.write(str(d)) # Flush so the output is seen immediately sys.stdout.flush() if __name__ == "__main__": main()
Name | Type | Size | Permission | Actions |
---|---|---|---|---|
.beer.pyo.40009 | File | 703 B | 0644 |
|
.eqfix.pyo.40009 | File | 4.52 KB | 0644 |
|
.fact.pyo.40009 | File | 1.14 KB | 0644 |
|
.find-uname.pyo.40009 | File | 1.47 KB | 0644 |
|
.from.pyo.40009 | File | 749 B | 0644 |
|
.lpwatch.pyo.40009 | File | 2.54 KB | 0644 |
|
.makedir.pyo.40009 | File | 732 B | 0644 |
|
.markov.pyo.40009 | File | 3.93 KB | 0644 |
|
.mboxconvert.pyo.40009 | File | 3.17 KB | 0644 |
|
.morse.pyo.40009 | File | 4.33 KB | 0644 |
|
.newslist.pyo.40009 | File | 7.55 KB | 0644 |
|
.pi.pyo.40009 | File | 921 B | 0644 |
|
.pp.pyo.40009 | File | 2.28 KB | 0644 |
|
.primes.pyo.40009 | File | 921 B | 0644 |
|
.queens.pyo.40009 | File | 2.95 KB | 0644 |
|
.script.pyo.40009 | File | 1.21 KB | 0644 |
|
.unbirthday.pyo.40009 | File | 2.93 KB | 0644 |
|
.update.pyo.40009 | File | 2.69 KB | 0644 |
|
README | File | 1009 B | 0644 |
|
beer.py | File | 459 B | 0755 |
|
beer.pyc | File | 703 B | 0644 |
|
beer.pyo | File | 703 B | 0644 |
|
eqfix.py | File | 6.16 KB | 0755 |
|
eqfix.pyc | File | 4.52 KB | 0644 |
|
eqfix.pyo | File | 4.52 KB | 0644 |
|
fact.py | File | 1.11 KB | 0755 |
|
fact.pyc | File | 1.14 KB | 0644 |
|
fact.pyo | File | 1.14 KB | 0644 |
|
find-uname.py | File | 1.18 KB | 0755 |
|
find-uname.pyc | File | 1.47 KB | 0644 |
|
find-uname.pyo | File | 1.47 KB | 0644 |
|
from.py | File | 874 B | 0755 |
|
from.pyc | File | 749 B | 0644 |
|
from.pyo | File | 749 B | 0644 |
|
lpwatch.py | File | 2.77 KB | 0755 |
|
lpwatch.pyc | File | 2.54 KB | 0644 |
|
lpwatch.pyo | File | 2.54 KB | 0644 |
|
makedir.py | File | 510 B | 0755 |
|
makedir.pyc | File | 732 B | 0644 |
|
makedir.pyo | File | 732 B | 0644 |
|
markov.py | File | 3.51 KB | 0755 |
|
markov.pyc | File | 3.93 KB | 0644 |
|
markov.pyo | File | 3.93 KB | 0644 |
|
mboxconvert.py | File | 3.11 KB | 0755 |
|
mboxconvert.pyc | File | 3.17 KB | 0644 |
|
mboxconvert.pyo | File | 3.17 KB | 0644 |
|
morse.py | File | 4.21 KB | 0755 |
|
morse.pyc | File | 4.33 KB | 0644 |
|
morse.pyo | File | 4.33 KB | 0644 |
|
newslist.doc | File | 2.36 KB | 0755 |
|
newslist.py | File | 11.17 KB | 0755 |
|
newslist.pyc | File | 7.55 KB | 0644 |
|
newslist.pyo | File | 7.55 KB | 0644 |
|
pi.py | File | 888 B | 0755 |
|
pi.pyc | File | 921 B | 0644 |
|
pi.pyo | File | 921 B | 0644 |
|
pp.py | File | 3.73 KB | 0755 |
|
pp.pyc | File | 2.28 KB | 0644 |
|
pp.pyo | File | 2.28 KB | 0644 |
|
primes.py | File | 603 B | 0755 |
|
primes.pyc | File | 921 B | 0644 |
|
primes.pyo | File | 921 B | 0644 |
|
queens.py | File | 2.19 KB | 0755 |
|
queens.pyc | File | 2.95 KB | 0644 |
|
queens.pyo | File | 2.95 KB | 0644 |
|
script.py | File | 962 B | 0755 |
|
script.pyc | File | 1.21 KB | 0644 |
|
script.pyo | File | 1.21 KB | 0644 |
|
unbirthday.py | File | 3.07 KB | 0755 |
|
unbirthday.pyc | File | 2.93 KB | 0644 |
|
unbirthday.pyo | File | 2.93 KB | 0644 |
|
update.py | File | 2.69 KB | 0755 |
|
update.pyc | File | 2.69 KB | 0644 |
|
update.pyo | File | 2.69 KB | 0644 |
|