How do i compute a determinant using leibniz's (big and inefficient) formula?

For one, but that might be irrelevant to you, this formula is not the way you want to actually compute determinants; there are much more efficient methods for that. But I will try to help you understand the formula. We have an n×n n × n -matrix A A . Its elements are ai,j a i , j as usual, or a[i,j] a [ i , j ] in more C-like notation. The outer sum sums over all permutations (this set is Sn S n ), so that is an outer loop. A permutation is just a reordering of the indices, like (1,3,2) ( 1 , 3 , 2 ) for n=3 n = 3 ; for n=3 n = 3 there are 3! 3 ! permutations of the indices: (1,2,3),(1,3,2),(2,1,3),(2,3,1),(3,1,2),(3,2,1) ( 1 , 2 , 3 ) , ( 1 , 3 , 2 ) , ( 2 , 1 , 3 ) , ( 2 , 3 , 1 ) , ( 3 , 1 , 2 ) , ( 3 , 2 , 1 ) , and for n n indices there are n! n ! in general. Here we can interpret this as a function: (3,2,1) ( 3 , 2 , 1 ) is the function that sends 1 1 to 3 3 , 2 2 to 2 2 and 3 3 to 1 1 , e.g. Now, each permutation has a sign sgn(σ) sgn ⁡ ( σ ) , which is either +1 or -1. The identity has sign +1, while interchanging two elements of a permutation changes the sign to the opposite one. So (1,3,2),(3,2,1),(2,1,3) ( 1 , 3 , 2 ) , ( 3 , 2 , 1 ) , ( 2 , 1 , 3 ) are odd permutations on n=3 n = 3 , because we interchanged one pair. One can prove that half the permutations are even (have sign +1), and the other half has sign -1. Now for each permutation σ σ we take the multiplication of the elements of A A as specified by the permutation, interpreted as a map. So (1,3,2) ( 1 , 3 , 2 ) is the map 1→1,2→3,3→2 1 → 1 , 2 → 3 , 3 → 2 , so we get the product a1,1a2,3a3,2 a 1 , 1 a 2 , 3 a 3 , 2 , where the first index is the original, the second its image under σ σ . The formula says that det(A) det ( A ) is the sum of all such products where we consider all permutations of the index set, and the products from an odd permutation get a minus sign. So for n=3 n = 3 we get det(A)=a1,1a2,2a3,3−a1,1a2,3a3,2−a1,3a2,2a3,1−a1,2a2,1a3,3+a1,2a2,3a3,1+a1,3a2,1a3,2 det ( A ) = a 1 , 1 a 2 , 2 a 3 , 3 − a 1 , 1 a 2 , 3 a 3 , 2 − a 1 , 3 a 2 , 2 a 3 , 1 − a 1 , 2 a 2 , 1 a 3 , 3 + a 1 , 2 a 2 , 3 a 3 , 1 + a 1 , 3 a 2 , 1 a 3 , 2 where the minus signs correspond to the odd permutations from above. So you need algorithms to enumerate all permutations, and compute their signs, and then loop over all these products. It's quite complicated to get right, and not very efficient. But this formula allows one to prove some facts for determinants, so it has its uses. There's an example.