 trdiag algorithm 17 