schur-1.m

来自「GNU Octave is a high-level language, pri」· M 代码 · 共 4 行

a = [1, 2, 3; 4, 5, 9; 7, 8, 6];
[u, s] = schur (a);
size (u) == [3, 3] && size (s) == [3, 3] && abs (s - u' * a * u) < sqrt (eps)