Block tridiagonal solver

Solves block tridiagonal systems of equations.

This is an implementantion of the Gauss algorithm for block tridiagonal systems. The blocks must all be the same size. This approach requires for loops so is of value when the matrix is large. Ideal for solving 1D BVPs using finite element or spectral methods.