bessely
==============================================

Purpose
----------------

Computes a Bessel function of the second kind (Weber's function), :math:`Y_n(x)`.

Format
----------------
.. function:: y = bessely(n, x)

    :param n: the order of the Bessel function. Nonintegers will be truncated to an integer.
    :type n: NxK matrix or P-dimensional array where the last two dimensions are NxK

    :param x: where the last two dimensions are LxM, ExE conformable with n.
    :type x: LxM matrix or P-dimensional array

    :return y: where the last two dimensions are max(N,L) by max(K,M).

    :rtype y: max(N,L) by max(K,M) matrix or P-dimensional array

Examples
----------------

::

    // Create the sequence 0.1, 0.2, 0.3, 0.4, 0.5
    x = seqa(0.1, 0.1, 5);

    // Create the sequence 1, 1.1, 1.2, 1.3, 1.4
    x2 = seqa(1, 0.1, 5);

    /*
    ** Calculate a first order bessel function against 'x' and
    ** calculate a third order bessel function agains 'x2'
    ** NOTE: The '~' provides horizontal concatenation
    */
    ord = { 1 3 };
    y = bessely(n, x~x2);

After the code above:

::

        -6.459 -5.822         0.100  1.000
        -3.324 -4.507         0.200  1.100
    y = -2.293 -3.590  x~x2 = 0.300  1.200
        -1.781 -2.930         0.400  1.300
        -1.471 -2.442         0.500  1.400

.. seealso:: Functions :func:`besselj`, :func:`mbesseli`, :func:`besselk`