LU Factorizations

[Nick Allen]

I have a question about using Octave to determine an LU matrix
factorization.  Straight out of the manual for the "lu" function I run:

[l, u, p] = lu (a)

l =
  1.00000  0.00000
  0.33333  1.00000

u =
  3.00000  4.00000
  0.00000  0.66667

p =
  0  1
  1  0

but if i then do l*u I get:
  3  4
  1  2

when it should result in the original matrix a, which was:
  1  2
  3  4

Could someone please explain to me what is going on?  Is this a problem
with Octave, or is it a fault in my basic understanding of LU
factorizations?  By the way I am running Octave 2.1.36 on an ia32 Linux
platform.

Thanks
Nick

-- 
_____________________________
Nick Allen <address@hidden>





-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------