## Run script twice.  First time is okay, second time segfaults.
A=[-3.6 0 9.72;3.4 1.9 8.67;3.8 1.7 0;-1.6 4 0;-1.6 4 0];
[Q, R, p] = qr (A, 'vector')
