\begin{Schunk}
\begin{Sinput}
  xseq <- seq(0.001,0.999, by=0.001)
  xd <- dbeta(xseq, 0.9, 0.9)
  plot(xd ~ xseq, type="l", xlab="A Beta Variate",ylab="Probability")
  lines(c(0,0),c(0,xd[1]),lty=4)
  lines(c(1,1),c(0,xd[1]),lty=4)
\end{Sinput}
\end{Schunk}
