我试着生成随机点并保留曲线内的点，这有点困难。下面是我的代码：

import numpy as np
import matplotlib.pyplot as plt
N=eval(input("Enter the number of points used for integration:"))
a,b=eval(input("Enter the range of integration:"))
def f(x):
    return x/(x**2 + 1)

x=np.random.uniform(a,b,N)
y=np.random.uniform(a,b,N)
a=np.linspace(a,b,N)
A=[f(j) for j in x]
xx=x.tolist()
yy=y.tolist()
xx.sort()
yy.sort()
X=[]
Y=[]
for i in range(N):
    if yy[i]<A[i]:
        Y.append(yy[i])
        X.append(xx[i])
plt.plot(a,f(a),color='black',label="given function")
plt.grid()
plt.scatter(X,Y)
plt.legend()
plt.show()

但是我得到的只是一个分布在整个网格中的分布，我在for循环中特别使用了if来表示我的条件，这是怎么回事？
我试着生成随机数，它们保留曲线下的那些，我所期望的正是这样，我试着用sort（）对随机数列表排序，但不知何故，这让事情变得更糟。