polygon()

Un truco para rellenar regiones a destacar en un gráfico es utilizar la función polygon().
El gráfico anterior se consigue con el siguiente código:

curve(exp,0.5,3.5,lwd=2)

t<-seq(1,3,by=0.01)
x <- c(1,t,3)
y <- c(0,exp(t),0)
polygon(x,y,density=20)

title("Exponencial entre 1 y 3")

Observemos la utilización de los vértices adecuados para reseguir la curva exponencial.
El parámetro density=20 controla la densidad de las linias (por pulgada). Tambíen es posible rellenar el polígono con un color sólido como en el siguiente ejemplo con la densidad normal:

El código para dibujarlo es:

curve(dnorm,-3,3,axes=F) # dibujamos la curva pero no los ejes
axis(1) # ahora los ejes
axis(2)

t <- seq(0,1,by=0.01)
x <- c(0,t,1)
y <- c(0,dnorm(t),0)
polygon(x,y,col="grey") # la región en gris

curve(dnorm,-3,3,lwd=2,add=T) # la curva encima
abline(h=0)
title("Probabilidad entre 0 y 1")