# 2006 年 05 月 23 日作成
xx<-seq(-6,6,by=0.01)
n=10000000
x<-rnorm(n,0,1)
# col=2 は赤，col=3 は緑，col=4 は青，col=5 は空色，col=6 はピンク，col=7 は黄色
plot(xx,dnorm(xx,0,1),xlim=c(-5,5),ylim=c(0,0.5),col=2,type="l")
par(new=T)
plot(density(x),xlim=c(-5,5),ylim=c(0,0.5))
p=sum(x<=1/2)
p/n
[1] 0.6914858
pnorm(1/2,0,1)
[1] 0.6914625