\begin{Schunk}
\begin{Sinput}
 x <- seq(0,1,by=.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}
