theta denotes the parameter vector (alpha, beta, delta, lambda,
mu). Either provide the parameters individually OR provide theta.
The function searches for a root between qmin and qmax with
uniroot.
Boundaries can either be supplied by the user or a built-in approach using
the stable distribution is used.