Made Trdiagonal matrix calculation little more efficient
Added ConvergenceException, to be thrown when tridiagonal function, doesn't converge