function x=glev(r,b);
%GLEV  	General Levinson recursion.
%----
%USAGE: x = glev(r,b)
%
%	The Levinson recursion solves the Toeplitz equations
%		R x = b
%	where R=toeplitz(r) and b is an arbitrary vector.
%	
%
%---------------------------------------------------------------
% copyright 1996, by M.H. Hayes.  For use with the book 
% "Statistical Digital Signal Processing and Modeling"
% (John Wiley & Sons, 1996).
% page 268
%---------------------------------------------------------------

r=r(:);
p=length(b);
a=1;
x=b(1)/r(1);
epsilon=r(1);
for j=2:p;
	g=r(2:j)'*flipud(a);
	gamma(j-1)=-g/epsilon;
	a=[a;0] + gamma(j-1)*[0;conj(flipud(a))];
	epsilon=epsilon*(1 - abs(gamma(j-1))^2);
	delta=r(2:j)'*flipud(x);
	q=(b(j)-delta)/epsilon;
	x=[x;0] + q*[conj(flipud(a))];
end;