x <- seq(-3,3,0.01)
y1 <- dnorm(x,0,1)
y2 <- 0.5*dnorm(x,0,1)
plot(x,y1,type=\"l\",bty=\"L\",xlab=\"X\",ylab=\"dnorm(X)\")
points(x,y2,type=\"l\",col=\"red\")
p
x <- seq(-3,3,0.01)
y1 <- dnorm(x,0,1)
y2 <- 0.5*dnorm(x,0,1)
plot(x,y1,type="l",bty="L",xlab="X",ylab="dnorm(X)")
points(x,y2,type="l",col="red")
polygon(c(x,rev(x)),c(y2,rev(y1)),col="skyblue")
idx <- which(x %in% -3:3)
text(x[idx], y1[idx], labels = round(y1[idx], 2))
text(x[idx], y2[idx], labels = round(y2[idx], 2), col="red")