Working tri-diagonal solver
Working tri-diagonal solver
Fills  diagonal with random
Fills E successfully
Moved common functions to matrix.cpp