% RP.m
% Rectangular Coordinates to Polar Coordinates transform
function [in_an,in_pn]=RP(ich3,qch3)
n=length(ich3);
for i=1:n
    in_an(i)=sqrt(ich3(i)^2+qch3(i)^2);
    if(ich3(i)>=0)
        if(ich3(i)==0)
            in_pn(i)=pi/2;
        else
        in_pn(i)=atan(qch3(i)/ich3(i));
        end 
    else
        in_pn(i)=pi+atan(qch3(i)/ich3(i));
    end
end