function z = sinint(x)
%SININT Sine integral function.
%   SININT(x) = int(sin(t)/t,t,0,x).
%
%   See also COSINT.

%   Copyright (c) 1993-98 by The MathWorks, Inc.
%   $Revision: 1.5 $  $Date: 1998/06/05 21:36:04 $

if findstr(char(x),'Inf')
   z = sym(pi/2);
   stat = 0;
else
   [z,stat] = maple('map','sinint',x);
end

if stat
   % Singularities
   z = x;
   for k = 1:prod(size(x))
      [z(k),stat] = maple('sinint',x(k));
      if stat
         z(k) = Inf;
      end
   end
end
