利用python的numpy和matplotlib库
import matplotlib.pyplot as plt
from matplotlib.patches import Circle
import numpy as np
n = 5000 # 投掷10000个点
blue_x = []
blue_y = []
red_x = []
red_y = []
for i in range(1, n):
x = np.random.uniform(0, 1, 1)
y = np.random.uniform(0, 1, 1)
distance = pow(x ** 2 + y ** 2, 0.5)
if distance < 1:
blue_x.append(x), blue_y.append(y)
else:
red_x.append(x), red_y.append(y)
fig = plt.subplots(1, 1)
plt.scatter(blue_x, blue_y, color='b')
plt.scatter(red_x, red_y, color='r')
plt.xlabel('X')
plt.ylabel('Y')
plt.xlim((0, 1)) # 修改坐标轴范围
plt.ylim((0, 1))
fig = plt.gca()
fig.set_aspect(1)
# 画圆
circle = Circle((0, 0), radius=1, facecolor='none', edgecolor='black', linewidth=3, alpha=1)
fig.add_patch(circle)
# 显示图像
plt.show()