All published worksheets from http://sagenb.org
import math def main(): x=math.sqrt(20000.0) i=-1.0 for n in range(50): x=math.sqrt(((100.0-math.sqrt(10000.0-(x/2.0)**2.0))**2.0 + (x/2.0)**2)) y=(x*2**51.0)/100.0 k=math.pi print '%1.50f' % y print '%1.50f' % k main()