% Script P3
clc
T0 = [7 4 3 ;...
0 2 8 ;...
0 -9 5 ]
[Q,T] = Swap3(T0);
fprintf('errQ = norm(Q''*Q - I) = %10.3e\n',norm(Q'*Q-eye(3,3)))
fprintf('errT = norm(Q''*T0*Q - T) = %10.3e\n',norm(Q'*T0*Q-T))
format long
Q = Q
T = T