\begin{Schunk}
\begin{Sinput}
 library(VGAM)
 xvals <- seq(0.1,10,len=1000)
 plot(xvals, dpareto(xvals, location=1, shape=1), type="l", xlab="possible values of x", ylab="probability of x", main="Pareto Probability Density")
 text(2, .6, "shape=1, location=1", pos=4, col=1)
 
\end{Sinput}
\end{Schunk}
