R本身好像没有，但是有些包里有，比如rmutil包里的rlaplace函数

或者直接用inverse transform sampling自己写一个：

rdoubleexp <- function(n, mu, b) {
  y <- runif(n)
  x <- rep(NA, n)
  x[y <= 1/2] <- mu + b * log(2*y[y <= 1/2])
  x[y > 1/2] <- mu - b * log(2*(1-y[y > 1/2]))
  return(x)
}