bessely#

Purpose#

Computes a Bessel function of the second kind (Weber’s function), \(Y_n(x)\).

Format#

y = bessely(n, x)#

Parameters:

n (NxK matrix or P-dimensional array where the last two dimensions are NxK) – the order of the Bessel function. Nonintegers will be truncated to an integer.
x (LxM matrix or P-dimensional array) – where the last two dimensions are LxM, ExE conformable with n.

Returns:

y (max(N,L) by max(K,M) matrix or P-dimensional array) – where the last two dimensions are max(N,L) by max(K,M).

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