function [ a ] = myarg( z )
% Returns the polar angle of a complex number
% Input
%   z = a complex number
% Output
%   a = polar angle
x = real(z); % real part of z
y = imag(z); % imaginary part of z
a = atan( y/x ) + pi/2 * sign(y) * (1 - sign(x));
end

