% RI2FQ Convert complex pole/zero into frequency/Q pole/zero representation.
%
%  [f0, q]= ri2fq(c)
%
% M Hewitson 26-01-07
%
% $Id: ri2fq.m,v 1.3 2008/10/24 09:40:19 hewitson Exp $
%

function [f0, q]= ri2fq(c)
  
  if(nargin==0)
    disp('usage: [f0, q]= ri2fq(c)');
    return
  end
  
  a = real(c(1));
  b = imag(c(1));

  f0 = sqrt(a^2 + b^2)/(2*pi);
  q = 0.5*sqrt(1 + b^2/a^2);
  
end
% END