The task is to solve a system of simultaneous linear equations:

                    Ax = b,
where A is the coefficient matrix, b is the right-hand side, and x is the solution. If there are several right-hand sides, then
                    AX = B,  
where the columns of B are the individual right-hand sides, and the columns of X are the corresponding solutions.

Two types of routines, Driver and Computational, are provided for solving systems of linear equation, including in the forms of A^TX=B and A^HX=B.