function Y = bessely(nu,Z)
%BESSELY Symbolic Bessel function, Y(nu,z).

%   Copyright (c) 1993-98 by The MathWorks, Inc.
%   $Revision: 1.6 $  $Date: 1997/11/29 01:05:19 $

nu = sym(nu);
Z = sym(Z);

if all(size(nu) == 1)
   [Y,stat] = maple('map2','bessely',nu,Z);
elseif all(size(Z) == 1)
   [Y,stat] = maple('map','bessely',nu,Z);
else
   error('One argument must be a scalar.')
end


if stat
   % Singularity encountered.
   if all(size(nu) == 1)
      nu = nu(ones(size(Z)));
   else
      Z = Z(ones(size(nu)));
   end
   Y = Z;
   for k = 1:prod(size(Z));
      [Y(k),stat] = maple('bessely',nu(k),Z(k));
      if stat
         Y(k) = inf;
      end
   end
end
