使用黎曼和求积分

Question

如何使用黎曼和求出J = P（0 <= Z <= 1），其中Z〜N（0,1）？

所以J是函数（1 / sqrt（2 * pi））* exp ^（（-x ^ 2）/ 2）从0到1的整数

这是我在R中实现此方法的方法

 m<-5000
a<-0
b<-1
w<-(b-a)/m
x<-seq(a+(w/2),b-(w/2),w)
h<-(1/sqrt(2*pi))*exp^((-x^2)/2)

# Error in exp^((-x^2)/2) : non-numeric argument to binary operator

sum(h*w)
#Error: object 'h' not found

我不知道为什么标记这种错误，我键入is.numeric(x)并返回TRUE，所以问题出在哪里，如果我只组合数字呢？

Answer 1

删除^功能后面的exp。试试：

m<-5000
a<-0
b<-1
w<-(b-a)/m
x<-seq(a+(w/2),b-(w/2),w)
h<-(1/sqrt(2*pi))*exp((-x^2)/2)

sum(h*w)
[1] 0.3413447