function Y = atan(X)
% ATAN    Inverse tangent.
% (Quaternion overloading of standard Matlab function.)

% Copyright © 2006 Stephen J. Sangwine and Nicolas Le Bihan.
% See the file : Copyright.m for further details.

error(nargchk(1, 1, nargin)), error(nargoutchk(0, 1, nargout))

if isreal(X)
    
    % X is a real quaternion, and we compute the inverse tangent of an
    % isomorphic complex number using the standard Matlab atan function,
    % then construct a quaternion with the same axis as the original
    % quaternion.
    
    Y = isoquaternion(atan(isocomplex(X)), X);
else
    
    % X is a complex quaternion, and therefore we cannot use the method
    % above for real quaternions, because it is not possible to construct
    % an isomorphic complex number.
    
    error('quaternion/atan is not yet implemented for complex quaternions');
end;