(Scipy) 連続分布(確率密度)の総和の求め方

Question

Scipyで確率密度の総和を求めたいです．
このためには，積分が必要だと思いますが，どのように実行されるのでしょうか．
当然ではありますが，以下のコードを実行すると，1を超えます．

from scipy.stats import uniform

n = 11
xs = np.arange(0, n, 0.01)
a = 3
b = 8 - a

print(uniform.pdf(xs, loc=a, scale=b))

Answer 1

simpsあるいはrombを使って積分できます。
rombの方が精度はいいですが、データ数に制限があります。
https://org-technology.com/posts/integrate-function-fixed-sample.html

from scipy.stats import uniform
from scipy.integrate import simps

x = np.linspace(0, 1, 1000)
y = uniform.pdf(x)
simps(y, x)