Problem with Schur decomposition

[david allen]

What is supposed to be a upper quasi-triangular matrix has non-zero 
elements below the first sub-diagonal. This is demonstrated in the 
attached file. Computer and OS information are also in the file.

schur.cpp
Description: Text Data