在R中绘制密度函数的积分的正确方法是什么？

Question

我正在寻找一种在R中集成正态分布密度函数的方法，而我的方法是给我一个错误，即“ from”必须是有限的。

我尝试使用function()绘制分布，但这是我发现集成实际函数的最佳方法。

t <- seq(from = 87, to = 95, 0.1)

dicht <- function(x){ (1/(sqrt(2*pi)) * exp((-1/2)*(x-91)^2)) }

dichtI <- function(x){integrate(dicht, lower = 0, upper = 95)}

test_vect = Vectorize(dichtI, vectorize.args = 'x')
plot(t, dicht(t), type = "l")
curve(t, test_vect(t))

理想情况下，我希望能够绘制分布密度函数的积分函数。感谢您的提前帮助！