第一次实验

代码：

import numpy as np

import matplotlib.pylab as plt

n = int(input('输入一个整数'))

m = 0

import random

for i in range(n):

    x = random.random()

    y = random.random()

    if x ** 3 > y:

        m = m + 1

R = m / n

print(R)

x=np.linspace(0,1,num=150)

y=x**3

plt.plot(x,y)

plt.show()

结果：