function hh = hhat(A,n)
n = max(n,size(A,1));
[h,tha] = stricthh(A);
if length(tha<2)
  tha = repmat(tha,n-1,1);
end
hh = zeros(n);
for i=1:n
  for j=i+1:n
    x = 1;
    for k=i:j-1
      x = x * tha(k);
    end
    hh(i,j) = x;
  end
end
hh = hh + eye(n);
return