%
% Function to return a Jordan block of a given size sz
% with eigenvalue z.
%
function J=jorblk(z,sz)

   J = diag(ones(sz,1)*z)+...
       [zeros(sz-1,1) diag(ones(sz-1,1)); zeros(1,sz)];
