If A is a complex and quadratic matrix there exists a unitary Q such that Q*AQ = T = D + N with D being the diagonal matrix consisting of the eigenvalues and N being strictly upper tridiagonal.
A = [3 6 1
23 13 1
0 3 4];
T = schur(A);
We also display the runtime of schur
dependent on the square root of matrix elements: