def jnjnp_zeros(nt):
"""Compute nt (<=1200) zeros of the bessel functions Jn and Jn'
and arange them in order of their magnitudes.
Returns
-------
zo[l-1] : ndarray
Value of the lth zero of of Jn(x) and Jn'(x). Of length `nt`.
n[l-1] : ndarray
Order of the Jn(x) or Jn'(x) associated with lth zero. Of length `nt`.
m[l-1] : ndarray
Serial number of the zeros of Jn(x) or Jn'(x) associated
with lth zero. Of length `nt`.
t[l-1] : ndarray
0 if lth zero in zo is zero of Jn(x), 1 if it is a zero of Jn'(x). Of
length `nt`.
See Also
--------
jn_zeros, jnp_zeros : to get separated arrays of zeros.
"""
if not isscalar(nt) or (floor(nt)!=nt) or (nt>1200):
raise ValueError("Number must be integer <= 1200.")
nt = int(nt)
n,m,t,zo = specfun.jdzo(nt)
return zo[1:nt+1],n[:nt],m[:nt],t[:nt]
def jnjnp_zeros(nt):
"""Compute nt (<=1200) zeros of the bessel functions Jn and Jn'
and arange them in order of their magnitudes.
Returns
-------
zo[l-1] : ndarray
Value of the lth zero of of Jn(x) and Jn'(x). Of length `nt`.
n[l-1] : ndarray
Order of the Jn(x) or Jn'(x) associated with lth zero. Of length `nt`.
m[l-1] : ndarray
Serial number of the zeros of Jn(x) or Jn'(x) associated
with lth zero. Of length `nt`.
t[l-1] : ndarray
0 if lth zero in zo is zero of Jn(x), 1 if it is a zero of Jn'(x). Of
length `nt`.
See Also
--------
jn_zeros, jnp_zeros : to get separated arrays of zeros.
"""
if not isscalar(nt) or (floor(nt) != nt) or (nt > 1200):
raise ValueError("Number must be integer <= 1200.")
nt = int(nt)
n, m, t, zo = specfun.jdzo(nt)
return zo[1:nt + 1], n[:nt], m[:nt], t[:nt]
def jnjnp_zeros(nt):
"""Compute nt (<=1200) zeros of the bessel functions Jn and Jn'
and arange them in order of their magnitudes.
Outputs (all are arrays of length nt):
zo[l-1] -- Value of the lth zero of of Jn(x) and Jn'(x)
n[l-1] -- Order of the Jn(x) or Jn'(x) associated with lth zero
m[l-1] -- Serial number of the zeros of Jn(x) or Jn'(x) associated
with lth zero.
t[l-1] -- 0 if lth zero in zo is zero of Jn(x), 1 if it is a zero
of Jn'(x)
See jn_zeros, jnp_zeros to get separated arrays of zeros.
"""
if not isscalar(nt) or (floor(nt)!=nt) or (nt>1200):
raise ValueError, "Number must be integer <= 1200."
nt = int(nt)
n,m,t,zo = specfun.jdzo(nt)
return zo[:nt],n[:nt],m[:nt],t[:nt]
def jnjnp_zeros(nt):
"""Compute nt (<=1200) zeros of the bessel functions Jn and Jn'
and arange them in order of their magnitudes.
Outputs (all are arrays of length nt):
zo[l-1] -- Value of the lth zero of of Jn(x) and Jn'(x)
n[l-1] -- Order of the Jn(x) or Jn'(x) associated with lth zero
m[l-1] -- Serial number of the zeros of Jn(x) or Jn'(x) associated
with lth zero.
t[l-1] -- 0 if lth zero in zo is zero of Jn(x), 1 if it is a zero
of Jn'(x)
See jn_zeros, jnp_zeros to get separated arrays of zeros.
"""
if not isscalar(nt) or (floor(nt) != nt) or (nt > 1200):
raise ValueError, "Number must be integer <= 1200."
nt = int(nt)
n, m, t, zo = specfun.jdzo(nt)
return zo[:nt], n[:nt], m[:nt], t[:nt]
