Bivariate Gumbel copula
\begin{equation} C_{\theta}^{Gu}(u_{1},u_{2})=\exp\{-((-\ln{u}_{1})^{\theta}+(-\ln{u}_{2})^{\theta})^{\frac{1}{\theta}}\},1\leq\theta<\infty\\ \end{equation}