x <- seq(0,5,length=1000)
plot(x = x,
     y = dnorm(x, 1.5, 0.4),
     type = "l",
     col = "white",
     axes = FALSE,
     mgp = c(2, 2, 2),
     ylim=c(0,2),                         # Set limit of y-axis
     frame.plot=TRUE,
     xlab = "theta",
     ylab = "plausibility",
     font.main = 1,
     main=paste("Distributions"),
     lwd=2
)
polygon(x,dnorm(x, 1.5, 0.4),col=1,border = NULL)
polygon(x,dnorm(x, 1, 0.5),col=2,border = NULL)