Function File: a = lauchli (n)

Function File: a = lauchli (n,mu)

Creates the matrix [ ones(1,n); mu*eye(n) ] The value mu defaults to sqrt(eps). This is an ill-conditioned system for testing the accuracy of the QR routine.

                A = lauchli(15);
                [Q, R] = qr(A);
                norm(Q*R - A)
                norm(Q'*Q - eye(rows(Q)))