naginterfaces.library.specfun.bessel_​y_​complex

naginterfaces.library.specfun.bessel_y_complex(fnu, z, n, scal)

bessel_y_complex returns a sequence of values for the Bessel functions $Y_{\nu +n}\left(z\right)$ for complex $z$, non-negative $\nu$ and $n=0,1,\dots ,N-1$, with an option for exponential scaling.

For full information please refer to the NAG Library document for s17dc

https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/s/s17dcf.html

Parameters

fnufloat

$\nu$, the order of the first member of the sequence of functions.

zcomplex

$z$, the argument of the functions.

nint

$N$, the number of members required in the sequence $Y_{\nu }\left(z\right),Y_{\nu +1}\left(z\right),\dots ,Y_{\nu +N-1}\left(z\right)$.

scalstr, length 1

The scaling option.

$scal=\text{'U'}$

The results are returned unscaled.

$scal=\text{'S'}$

The results are returned scaled by the factor $e^{-\left|Im\left(z\right)\right|}$.

Returns

cycomplex, ndarray, shape $\left(n\right)$

The $N$ required function values: $cy[i-1]$ contains $Y_{\nu +i-1}\left(z\right)$, for $\text{i}=1,2,\dots ,N$.

nzint

The number of components of $cy$ that are set to zero due to underflow. The positions of such components in the array $cy$ are arbitrary.

Raises

NagValueError

(errno $1$)

On entry, $z=(0.0,0.0)$.

(errno $1$)

On entry, $fnu=⟨value⟩$.

Constraint: $fnu\ge 0.0$.

(errno $1$)

On entry, $scal$ has an illegal value: $scal=⟨value⟩$.

(errno $1$)

On entry, $n=⟨value⟩$.

Constraint: $n\ge 1$.

(errno $2$)

No computation because $\left|z\right|=⟨value⟩<⟨value⟩$.

(errno $2$)

No computation because $Re\left(z\right)=⟨value⟩>⟨value⟩$, $scal=\text{'U'}$.

(errno $3$)

No computation because $fnu+n-1=⟨value⟩$ is too large.

(errno $5$)

No computation because $\left|z\right|=⟨value⟩>⟨value⟩$.

(errno $5$)

No computation because $fnu+n-1=⟨value⟩>⟨value⟩$.

(errno $6$)

No computation – algorithm termination condition not met.

Warns

NagAlgorithmicWarning

(errno $4$)

Results lack precision because $\left|z\right|=⟨value⟩>⟨value⟩$.

(errno $4$)

Results lack precision because $fnu+n-1=⟨value⟩>⟨value⟩$.

Notes

bessel_y_complex evaluates a sequence of values for the Bessel function $Y_{\nu }\left(z\right)$, where $z$ is complex, $-\pi <argz\le \pi$, and $\nu$ is the real, non-negative order. The $N$-member sequence is generated for orders $\nu$, $\nu +1,\dots ,\nu +N-1$. Optionally, the sequence is scaled by the factor $e^{-\left|Im\left(z\right)\right|}$.

Note: although the function may not be called with $\nu$ less than zero, for negative orders the formula $Y_{-\nu }\left(z\right)=Y_{\nu }\left(z\right)\cos \left(\pi \nu \right)+J_{\nu }\left(z\right)\sin \left(\pi \nu \right)$ may be used (for the Bessel function $J_{\nu }\left(z\right)$, see bessel_j_complex()).

The function is derived from the function CBESY in Amos (1986). It is based on the relation $Y_{\nu }\left(z\right)=\frac{H_{\nu }^{\left(1\right)}\left(z\right)-H_{\nu }^{\left(2\right)}\left(z\right)}{2i}$, where $H_{\nu }^{\left(1\right)}\left(z\right)$ and $H_{\nu }^{\left(2\right)}\left(z\right)$ are the Hankel functions of the first and second kinds respectively (see hankel_complex()).

When $N$ is greater than $1$, extra values of $Y_{\nu }\left(z\right)$ are computed using recurrence relations.

For very large $\left|z\right|$ or $(\nu +N-1)$, argument reduction will cause total loss of accuracy, and so no computation is performed. For slightly smaller $\left|z\right|$ or $(\nu +N-1)$, the computation is performed but results are accurate to less than half of machine precision. If $\left|z\right|$ is very small, near the machine underflow threshold, or $(\nu +N-1)$ is too large, there is a risk of overflow and so no computation is performed. In all the above cases, a warning is given by the function.

References

NIST Digital Library of Mathematical Functions

Amos, D E, 1986, Algorithm 644: A portable package for Bessel functions of a complex argument and non-negative order, ACM Trans. Math. Software (12), 265–273