用matplotlib实现

    1. import matplotlib.pyplot as plt
    2. import numpy as np
    4. # 支持中文
    5. plt.rcParams['font.sans-serif'] = ['SimHei']  # 用来正常显示中文标签
    6. plt.rcParams['axes.unicode_minus'] = False  # 用来正常显示负号
    7. count = 10000
    8. plt.figure(figsize=(11,11))  # 设置画布长宽
    9. a,hits = 5,0            # 设置随机数范围
    10. lsx = []                     # 落在圆外的点的列表x
    11. lsy = []                     # 落在圆外的点的列表y
    12. lsPx = []                    # 落在圆内的点的列表x
    13. lsPy = []                    # 落在圆内的点的列表y
    14. for i in range(count):
    15.     x,y = np.random.uniform(-a,a),np.random.uniform(-a,a)
    16.     pos = (x ** 2 + y ** 2) ** 0.5  # 计算坐标（x,y）到原点的距离
    17.     if pos <= 5.0:                  # 如果距离小于等于5，点在圆内
    18.         lsPx.append(x)
    19.         lsPy.append(y)
    20.         cl = 'red'
    21.     else:                           # 如果距离大于5.0，点在圆外
    22.         lsx.append(x)
    23.         lsy.append(y)
    24.         cl = 'blue'
    25. plt.scatter(lsx,lsy,facecolors='blue')  # 绘制散点，落在圆外在点颜色用蓝色
    26. plt.scatter(lsPx,lsPy,facecolors='red') # 绘制散点，落在圆内在点颜色用红色
    27. PI = 4 * len(lsPx) / count              # 落在圆内的点数与总数的比值为面积比值，可计算PI值
    28. plt.title('模拟{}次时PI值为{}'.format(count,PI))  # 用算得的PI值做图的标题
    29. plt.show()                              # 显示绘制结果

    image.png
    turtle实现

    1. import random                   # 导入随机数库
    2. import turtle                   # 导入turtle库
    4. n = int(input())
    5. # random.seed(10)
    6. # turtle.tracer(10)              刷新间隔,提高绘制速度
    8. turtle.pensize(4)
    9. turtle.penup()
    10. turtle.goto(-300, -300)
    11. turtle.pendown()
    12. for i in range(4):
    13.     turtle.forward(600)
    14.     turtle.left(90)
    15. turtle.forward(300)
    16. turtle.pencolor('green')
    17. turtle.circle(300)
    18. hits = 0                                    # 落在圆内的计数器初值设为 0
    19. for i in range(1, n+1):
    20.     x, y = random.random(), random.random() # 生成两个随机数模拟一个点的坐标
    21.     signx,signy = random.choice('+-'),random.choice('+-')
    22. # 随机生成x和y的符号“+”或“-”
    23.     pos = (x ** 2 + y ** 2)**0.5  # 计算坐标（x,y）到原点的距离
    24.     if pos <= 1.0:                # 如果距离小于等于1，点在圆内
    25.         hits = hits + 1
    26.         turtle.pencolor('red')    # 落在圆内在点颜色用红色
    27.     else:                         # 如果距离大于1，点在圆外
    28.         turtle.pencolor('blue')   # 落在圆内在点颜色用蓝色
    29.     x = float(signx+str(x))       # 将 x 值与符号拼接并转为浮点数
    30.     y = float(signy+str(y))       # 将 y 值与符号拼接并转为浮点数
    31.     turtle.penup()
    32.     turtle.goto(x * 300, y * 300) # 画笔抬起并移动到数值放大400倍的x,y处
    33.     turtle.pendown()
    34.     turtle.dot(5)                 # 画一个半径为 3 的圆点
    35.     if i % 10000 == 0:             # 实验为10000的倍数次时
    36.        turtle.pencolor('black')
    37.        pi = 4 * (hits / i)        # 根据落在圆内外的点数量估算PI值
    38.        turtle.penup()             # 画笔抬起
    39.        turtle.goto(320, 150- i // 1000 * 30)  # 移动到区域外记录当前PI值
    40.        turtle.pendown()           # 画笔抬起
    41.        turtle.write("{}次时PI的值是{:.4f}".format(i,pi),font=("宋体", 18, "normal"))
    42. turtle.update()                   # 刷新
    43. turtle.hideturtle()               # 隐藏光标
    44. turtle.update()                   # 刷新
    45. turtle.done()                     # 结束绘制

    TIM截图20200419214406.png