python 割圆法求π的问题

为什么会报错

import math
def f(x):
    h = 1 - math.sqrt(1-(x/2)**2)
    return math.sqrt(h**2+(x/2)**2)

a = 1
k = 6
for i in range(100):#这里把n改成你想要设定的次数
    a = f(a)
    k*=2
print(a*k/2)

这样试试