import numpy as np import matplotlib.pylab as plt class MC: def __init__(self): self.a=2000 def MC(self): n=0; for i in range(self.a): x=np.random.uniform(0,1) y=np.random.uniform(0,1) if y<=x**3: n=n+1 plt.scatter(x,y,c='green') else: plt.scatter(x, y, c='red') s=n/self.a; print(s) def draw(self): x = np.linspace(0, 1, 200) y = x ** 3 plt.plot(x, y) plt.show() MC().MC() MC().draw()