请问如何用python求圆周率 （并且用numpy的方法）

#python#
要求求误差，希望运行结果如图 并且希望可以对应一亿个点
谢谢

img


pi = 0
N = 100
for k in range(N):
  pi += 1/pow(16, k) * (4 / (8 * k + 1) - 2 /(8 * k + 4) - 1/(8 * k + 5) - 1 /(8 * k + 6))
print("圆周率值是:{}".format(pi))

使用numpy计算圆周率，可参考下这里：
https://www.cnblogs.com/Flat-White/p/14807002.html