% P1 
% Checks out SkewReduce
clc
randn('seed',0)
n = 6; A = randn(n,n); A = A - A';
[Q,T] = SkewReduce(A);
% Check Q...
errQ = norm(Q'*Q - eye(n,n))
% Check Residual...
errResid = norm(Q'*A*Q - T)
% Display T
format short
T = T

