1 2 3 4 5 6 7 8
m=augment([1;2], [3;4]) transpose(m) det(m) lu(m, l, u) l, u inverse(m) 3/2*m 3.0/2*m