From Cramer’s rule, it is known that |D| is the determinant
of the coefficient matrix D and |Dk| is the determinant
of the matrix formed from the original coefficient matrix
by replacing the kth column of D by the column with the
elements b1,. . . ,bn [10],[11], [12], [14]. Thus, the solution of
the system (1) is given by the formulas