uemean=x1(1);
unmean=x2(1);
nc=length(x1-1)/2;
A=x2(2:2:2*nc);
B=x2(3:2:2*nc+1);
C=x1(2:2:2*nc);
D=x1(3:2:2*nc+1);

uemag=sqrt(C.^2+D.^2);
unmag=sqrt(A.^2+B.^2);
uephase=atan2(D,C);
unphase=atan2(B,A);

A1=.5*(A+D);
B1=.5*(C-B);
A2=.5*(A-D);
B2=.5*(C+B);           
R1=sqrt(A1.^2.+B1.^2.);
R2=sqrt(A2.^2.+B2.^2.);
RMA=R1+R2;
RMI=R1-R2;
G1=atan2(B1,A1);
G2=atan2(B2,A2);
OREN= 0.5*(G1+G2); 
PHASE=0.5*(G2-G1);