function x=solvhank(c,r,b)
% solvhank(c,r,b) == hankel(c,r)\b
% solvhank(c,b) == hankel(c)\b

% Dirk Laurie 2000-08-11
n=size(c); k=n:-1:1;
if nargin<3, b=r; r=c(k); end
x=solvtoep(c(k),r,b(k));
 
