rkwbs(n, alpha, beta, a, b)$$F(x; \alpha, \beta, a, b) = 1 - {1 - \Phi(\nu)^a}^b$$ where $\beta > 0$ is a scale parameter and the other positive parameters $\alpha$, $a$, and $b$ are shape parameters in that $\nu = \alpha^{-1}\rho(x/\beta)$, $\rho(z) = z^{1/2} - z^{-1/2}$ and $\Phi(\cdotp)$ denotes the standard normal distribution function.
Cordeiro, G. M., and M. Castro. 2011. A new family of generalized distributions. J. Stat. Comput. Simulation, 81(7), 883-898.
Cordeiro, G.M. Lemonte, A.J. and Ortega, E.M.M. (2013). An extended fatigue life distribution. Statistics 47, 626-653.
a = 2.0
b = 2.5
alpha = 1.0
beta = 1.0
rkwbs(5, alpha, beta, a, b)Run the code above in your browser using DataLab