library(VGAM)
theta <-2;
n <- 100
x <-rlaplace(n,scale=theta)
fn <- function(X,theta=2){(1/2*theta)^n*exp((-1/theta)*sum(abs(X)))}
x = sort(x)
output = sapply(x,fn)
plot(x,output,type="l")