(l - aa2) * factorial(l + bb2) * factorial(l - bb2) * factorial(l + cc2) * factorial(l - cc2) * (2*l + 1)); sfac1 = factorial(3*l + 1); sfac2 = factorial(l); dcg = sqrt(sfac2*sfac2*sfac2 / sfac1); cglist[idxcg_count] = sum * dcg * sfaccg; ...