\begin{Schunk}
\begin{Sinput}
 x <- seq(0, 1, by = 0.005)
 b1 <- c(3, 0.7, 1.2)
 b2 <- c(5.6, 0.58, 0.2)
 pbeta1 <- dbeta(x, b1[1], b2[1])
 pbeta2 <- dbeta(x, b1[2], b2[2])
 pbeta3 <- dbeta(x, b1[3], b2[3])
 plot(x, pbeta1, type = "n", xlab = "x", ylab = "Probability Density", ylim = c(0, 
      4))
 lines(x, pbeta1, lty = 1)
 lines(x, pbeta2, lty = 2)
 lines(x, pbeta3, lty = 3)
 legend(0.55, 3.5, legend = c("Beta(3,5.6)", "Beta(0.7, 0.58)", "Beta(1.2,0.2)"), 
      lty = 1:3)
\end{Sinput}
\end{Schunk}
