% Script P1
  clc
% This generates a random A that has a complex eigenvalue...
  n=6;
  randn('seed',0)
  A = randn(n,n);
  [X,D] = eig(A);
  lambda = D(1,1)
  z = X(:,1)
 
  [Q,H] = Hess32(A,z,lambda);
  
  format short
  H = H
  
  disp(sprintf('norm(Q''*Q - eye(n,n) = %10.3e',norm(Q'*Q - eye(n,n))))
  disp(sprintf('norm(Q''*A*Q - H)     = %10.3e',norm(Q'*A*Q - H)))
  