First Submit

Made comments, tested code.  Only works for 2x2 matrices?  Get a weird
error otherwise: “Lapack routine dgesv: system is exactly singular:
U[10,10] = 0 “ or some variation of U[10,10]