HW5

Matrix Decompositions

1)
As discussued a long time ago it is frequenctly usefull
to have a program that decomposes a matrix A into
	A = L U
with L lower triangular with 1's on the diagonal
and U upper triangular.  

Write a Mathematica or Matlab function to do this.
The input should be a square matrix A.
The output should be two square matrices L and U
of the appropriate form.
It does not need to test for zero pivots or
perform any pivoting.
Test your program on a small matrix.

2) Repeat the process for (1) for a QR factorization
using Householder rotations.