I'm doing a project on random matrices and its applications. I have the joint probability density and want to calculate the probability of $s=\sum_{j=1}^N\lambda_j^2$. So we have
$$P(s)=C_{N,K}\int...\int\delta\left(s-\sum_{j=1}^N\lambda_j^2\right)\delta\left(1-\sum_{j=1}^N\lambda_j\right)\prod_i\lambda_i^{K-N}\prod_{i<j}(\lambda_i-\lambda_j)^2\prod_id\lambda_j$$
The limits being from $-\infty$ to $\infty$ for all the integrals and $N,K$ are constants. Now I tried expanding the Vandermonde determinant after getting rid of the delta functions and tried to use Selberg integral but unfortunately this is turning out to be very tedious and cumbersome. Even if its possible we'll end up sums of functions of gamma functions. Is there any other way to do this?