Task

Solve   Ax=b   using Gaussian elimination then backwards substitution.

A   being an   n by n   matrix.

Also,   x and b   are   n by 1   vectors.

To improve accuracy, please use partial pivoting and scaling.



See also

   the Wikipedia entry:   Gaussian elimination


[[   1.00   0.00   0.00   0.00   0.00   0.00  -0.01]
 [   1.00   0.63   0.39   0.25   0.16   0.10   0.61]
 [   1.00   1.26   1.58   1.98   2.49   3.13   0.91]
 [   1.00   1.88   3.55   6.70  12.62  23.80   0.99]
 [   1.00   2.51   6.32  15.88  39.90 100.28   0.60]
 [   1.00   3.14   9.87  31.01  97.41 306.02   0.02]]→[A]
Matr>List(rref([A]),7,L1)
L1