Fast Pentadiagonal System Solver

Solves symmetric and asymmetric pentadiagonal systems.

Solves the problem Ax=b when A is pentadiagonal (5-banded) and strongly nonsingular. This is much faster than x=A\b for large matrices. The algorithm will check to see if A is symmetric and use a more efficient algorithm if it is. Users are encouraged to improve and redistribute this script.