solve

solve, a MATLAB code which solves a linear system of equations A*x=b using Gauss elimination.

In MATLAB, of course, one can simply type "x = A\b" and get an answer, even if the matrix is ill-conditioned, singular, complex, rectangular (over-determined or under-determined).

Nonetheless, this library is useful because:

• it shows the steps involved in a simple Gauss-elimination procedure;
• it can be used to count the number of operations in Gauss-elimination;
• it can be used as a guide for writing a corresponding program in a non-MATLAB environment;
• its performance can be compared to corresponding versions in C or Fortran;
• it can be used as a starting point for exploring band storage, sparse storage, iterative solutions and other topics in linear algebra.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

solve is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.

Related Data and Programs:

solve_test

r8lib, a MATLAB code which contains many utility routines using double precision real (R8) arithmetic.

Source Code:

• r8mat_fs.m, factors and solves a linear system.
• r8vec_print.m, prints an R8VEC.