shape
  and scale.dEV(x, shape = 1, scale = 1, log = FALSE)
pEV(q, shape = 1, scale = 1, lower.tail = TRUE, log.p = FALSE)
qEV(p, shape = 1, scale = 1, lower.tail = TRUE, log.p = FALSE)
hEV(x, shape = 1, scale = 1, log = FALSE)
HEV(x, shape = 1, scale = 1, log.p = FALSE)
rEV(n, shape = 1, scale = 1)length(n) > 1, the length
    is taken to be the number required.dEV gives the density,
  pEV gives the distribution function,
  qEV gives the quantile function,
  hEV gives the hazard function,
  HEV gives the cumulative hazard function, and
  rEV generates random deviates.  Invalid arguments will result in return value NaN, with a warning.
scale parameter $a$ and
  shape parameter $\sigma$ has hazard function given by
  $$h(x) = (b/\sigma)(x/\sigma)^(b-1)\exp((x / \sigma)^b)$$
  for $x \ge 0$.