I am currently working on finding a closed-form equation for the variance of the ordered data of the three-parameter Weibull distribution. One method has been proposed, but the published solution is using numerical integration. A much faster solution would need to simplify the expression, which lead me to this expression.
Is there a solution in closed-form for this equation in which $n \in \mathbb{N} >1$?
$$ \int_0^1 (1-u)^n (\log(-\log(1-u)))^2 \mathrm{d}u $$
If not, is there a good approximation?