用Riemann和求积分

Question

如何用Riemann和求Z~N(0,1)的J=P(0<=Z<=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