> root2on2 = sqrt(2)/2 root2on2 = 0.707 > M = [ root2on2, root2on2,0; root2on2, root2on2, 0; 0,0,1] M = 0.707 0.707 0 0.707 0.707 0 0 0 1 > R = schur(M) R = t z > new= schord(R,1,3) new = t z > new.t -1.11e-16 0 0 0 1 0 0 0 1.41 > new.z -0.707 0 0.707 0.707 0 0.707 0 1 0 > new.z*new.t*new.z' 0.707 0.707 0 0.707 0.707 0 0 0 1See also: Page
schur()
.