 crout with pivot algorithm 16 