In: Advanced Math
Solve the linear system Az = b using the following methods:
I. The GE+PP algorithm for sparse (banded) linear systems, which is the default algorithm used by Matlab’s “\” operator when the matrix (call it Asparse) is of sparse type. You may find it easiest to set up the matrix using the spdiags command.