这篇文章就不多说了，直接上代码，大家自己理解

我仍然使用简单易懂的Python语言

# -*- coding: utf-8 -*- #声明编码
from __future__ import division
import time #导入时间模块，用来显示所用时间

程序根据马青公式计算圆周率，通过计算机计算更多位数以减少误差

在马青公式的运用过程中也感谢各位学霸的帮助

while True: #开始循环，此为死循环，为程序永不停止所设计
  number = int(input('请输入想要计算到小数点后的位数:'))
# 多计算10位，防止尾数取舍的影响
  time1=time.time()###开始计时（取一时间），这句代码的位置我认为放在这里是最合适的
  number1 = number+10
  # 算到小数点后number1位
  b = 10**number1
# 求含4/5的首项
  x1 = b*4//5
# 求含1/239的首项
  x2 = b//-239
# 求第一大项
  he = x1+x2
#设置下面循环的终点，即共计算n项
  number *= 2
#循环初值=3，末值2n,步长=2
  for i in range(3,number,2): #再次开始循环
      # 求每个含1/5的项及符号
      x1 //= -25
      # 求每个含1/239的项及符号
      x2 //= -57121
      # 求两项之和
      x = (x1+x2) // i
      # 求总和
      he += x

  # 求出π
  pai = he*4
#舍掉后十位
  pai //= 10**10
# 输出圆周率π的值
  paistring=str(pai)
  result=paistring[0]+str('.')+paistring[1:len(paistring)]
  print (result)

  time2=time.time()###取另一个时间，用差值得出耗时
  print (u'总共耗时：' + str(time2 - time1) + 's')

  print()#输出空行，美观实用

这个程序在CMD下运行的效率还是不错的