The attached script performs step by step the Cordic algorithm for circular case in vectoring mode.
Input: x and y coordinates, z=0 (initial phase).
Output: in x(n)=module multiply by An - gain constant. And in z(n)=phase (atan(x/y)).
Diego Barragán (2020). Using CORDIC for phase and magnitude computation (https://www.mathworks.com/matlabcentral/fileexchange/35843-using-cordic-for-phase-and-magnitude-computation), MATLAB Central File Exchange. Retrieved .