where D
0
is the matrix formed by deleting the last column of
the coefficient matrix D of system (1) and b is the column
matrix containing the constants of the system. We employ
Dodgson’s condensation to condense S1 until we finally
arrive at a row matrix Sn consisting of n + 1 elements. Let
us denote these elements by D, D1, D2, D3, . . . , Dn, so that