function y=ang_az_cal(XA,YA,XB,YB,XC,YC)
%Calculo do angulo azimutal entre tres pontos, medido no sentido
%retrogado de B para C
%A=ponto estacao
%B=ponto visado 1
%C=ponto visado 2
y=rumo_cal(XA,YA,XC,YC)-rumo_cal(XA,YA,XB,YB);
if (y<0)
    y=2*pi+y;
else
    if (y>pi)
        y=2*pi-y;
    end
end
end

