\begin{Schunk}
\begin{Sinput}
 leftedge <- c(1,3,5)
 x1 <- seq(leftedge[1],20,len=1000)
 y1 <- dpareto(x1, location = leftedge[1], shape = 1)
 x2 <- seq(leftedge[2], 20, len = 1000)
 y2 <- dpareto(x2, location = leftedge[2], shape = 1)
 x3 <- seq(leftedge[3], 20, len = 1000)
 y3 <- dpareto(x3, location = leftedge[3], shape = 1)
 plot(x1, y1, type="l", xlab = "possible values of x", ylab = "probability of x", ylim = c(0,1))
 lines(x2, y2, lty = 2, col=2)
 lines(x3, y3, lty = 3, col=4)
 legend("topright", c("location=1","location=3","location=5"), lty=1:3, col=c(1,2,4))
\end{Sinput}
\end{Schunk}
