theta denotes the parameter vector (alpha, delta, lambda).
Either provide the parameters alpha, delta, lambda
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.