import random import numpy as np from matplotlib import pyplot as plt plt.rcParams['font.sans-serif'] = ['SimHei'] plt.rcParams['axes.unicode_minus'] = False m = 1000000 n = 0 for i in range(m): x = random.uniform(0,1) y = random.uniform(0,1) if x**3>y: n += 1 R = n/m print('y=x**3在[0,1]上的定积分是:',R) x = np.arange(0, 1, 0.1) y = x * x *x plt.xlabel('x') plt.ylabel('y') plt.plot(x, y) plt.show()