function A=invhank(c,r)
% invhank(c,r) == inv(hankel(c,r))
 
% Dirk Laurie 2000-08-11 
n=size(c); k=n:-1:1;
if nargin<2, r=c(k); end
A=invtoep(c(k),r); A(:,k)=A;
