Block Tridiagonal Matrix Solver
The Solver:
TRDBLK.f