fixed multiplication, implemented LDU decomposition without partial pivoting