function [H, Q] = sym_hessenberg4(A)
% converts 4x4 A=A' to upper Hessenberg form using
% Givens transformations
%    / a11 a12 a13 a14 \         /h11 h12  0   0  \
%A = | a12 a22 a23 a24 | --> H = |h12 h22 h23  0  |
%    | a13 a23 a33 a34 |         | 0  h23 h33 h34 |
%    \ a14 a24 a34 a44 /         \ 0   0  h34 h44 /
% (1) zero first column