k[ll]=kEpshiftinv_gelesen[ll+1];//allocate and calculate gaudin matrix:gsl_matrix * m = gsl_matrix_alloc (N, N);for(intii=0; ii<N; ii++)for(intjj=0; jj<N; jj++) gsl_matrix_set (m, ii, jj, matrixelementenorm(k, c, ii, jj));//calculate determinant via LU decompositioni...