besselzero {CircularDDM}

Find First k Positive Zeros for the Bessel Functions

Description

Find first k positive zeros of the Bessel function J(n,x) or Y(n,x) using Halley's method.

Usage

besselzero(nu, k, kind)

Arguments

nu

The order of the corresponding Bessel function.

k

an integer for first k positive zeros.

kind

0, 1, or 2. A switch selects besselI, besselJ or besselY

Value

a vector

References

besselzero.m

Examples

nu <- seq(0, 5, length.out=10)
output <- matrix(numeric(5*length(nu)), nrow=5)
  for(i in 1:length(nu)) {
    output[,i] <- besselzero(nu[i], 5, 1)
  }
output

output <- matrix(numeric(5*length(nu)), nrow=5)
for(i in 1:length(nu)) {
    output[,i] <- besselzero(nu[i], 5, 2)
}
output