All published worksheets from http://sagenb.org
def main(): x=1.0 for n in range(50): x=math.sqrt((1.0-math.sqrt(1.0-(x/2.0)**2.0))**2.0+(x/2.0)**2.0) x=x*6.0*2.0**49.0 y=math.pi print '%1.50f' % x print '%1.50f' % y main()