from math import * global a, b a = [3] b = [2*sqrt(3)] i=0 while b[i]-a[i] > 0.0000000000000001: b_n = (2*a[i]*b[i])/(a[i] + b[i]) b.append(b_n) a_n = sqrt(a[i]*b[i+1]) a.append(a_n) i=i+1 print(i) print(a[i]) print(b[i]) print(3.141592653589793)