function I = cplxcomp(p1,p2)%  I = cplxcomp(p1,p2)% Compares two complex pairs which contain the same scalar elements%  but (possibly) at differrent indices.  This routine should be%  used after CPLXPAIR routine for rearranging pole vector and its%  corresponding residue vector.%      p2 = cplxpair(p1)%I=[];for j=1:1:length(p2)    for i=1:1:length(p1)	if (abs(p1(i)-p2(j)) < 0.0001)	   I=[I,i];        end    endendI=I';