Added lu(), fixed roundoff errors