function r=gtor(gamma,epsilon)
%function r=gtor(gamma,epsilon)

p = length(gamma);
aa=gamma(1);
r=[1 -gamma(1)];
for j = 2:p
  aa=[aa;0]+gamma(j)*[conj(flipud(aa));1];
  r=[r -fliplr(r)*aa];
end
if nargin == 2
  r = r*epsilon/prod(1-abs(gamma).^2);
end

