python-模拟math.pi

编写循环控制代码用下面公式逼近圆周率(精确到小数点后15位)，并且和math.pi的值做比较。

重点：

• 阶乘写成一个函数更加方便使用

import mathdef factorial(k):    result = 1    while k >= 1:        result *= k        k -= 1    return resultk = 0sum = 0while k <= 15:    a = factorial(4*k)    b = 1103 + 26390*k    c = (factorial(k))**4    d = 396**(4*k)    sum += a*b/c/d    k += 1ans = (2*(2**0.5))/9801*sumans = 1/ansprint "My answer: %.15f"%ansprint "math.pi: ",math.pi运算结果My answer: 3.141592730013306math.pi:  3.14159265359