function pout = refroots(pin)
 ord = length(pin) - 1;
 zers = roots(pin);
 for i=1:ord
   if real(zers(i))>0
     zers(i) = -zers(i);
   end
 end
 pout = real(pin(1)*poly(diag(zers)));
%
% Copyright MUSYN INC 1991,  All Rights Reserved
