aaa(A, B, C)

function [a1, b1, c1, a2, b2, c2] = aaa(A, B, C)
%AAA  gives both solutions to the angle-angle-angle problem, in radians.
%
%   AAA(A, B, C) will result in NaNs if the existence condition 
%   |pi - |A|-|B|| <= |C| <= pi - ||A| - |B||
%   is not met. 
%
%   See also AAAD.

% Rody P.S. Oldenhuis
% Delft University of Technology
% Last edited: 23/Feb/2009
    
    % first solution
    a1 = acos2((cos(A) + cos(B).*cos(C)) ./ (sin(B).*sin(C)), A);
    b1 = acos2((cos(B) + cos(A).*cos(C)) ./ (sin(A).*sin(C)), B);
    c1 = acos2((cos(C) + cos(A).*cos(B)) ./ (sin(A).*sin(B)), C);
    
    % second solution
    a2 = 2*pi - a; 
    b2 = 2*pi - b;
    c2 = 2*pi - c;

    % check constraints
    indices = ( abs(pi - abs(A)-abs(B)) <= abs(C) <= pi - abs(abs(A)-abs(B)) );
    a1(indices) = NaN;    a2(indices) = NaN;
    b1(indices) = NaN;    b2(indices) = NaN;
    c1(indices) = NaN;    c2(indices) = NaN;
 
end