atan_f.m

A function to find necessary angle

function theta=atan_f(a,b)
    theta=zeros(1,length(a));
    for n=1:length(a)
        if (a(n)>=0) && (b(n)>=0)
            theta(n)=atan(b(n)/a(n));
        elseif (a(n)<=0) && (b(n)>=0)
            theta(n)=atan(b(n)/a(n))+pi;
        elseif (a(n)>=0) && (b(n)<=0)
            theta(n)=atan(b(n)/a(n))+2*pi;
        else
            theta(n)=atan(b(n)/a(n))+pi;
        end
    end
end