Transcendental equation zj1=cj0 involving Bessel functions

Version 1.0.1 (2.27 KB) by Jie Zhang
Find N non-negative roots of transcendental equation z*J1(z) - c*J0(z) = 0, where J1(z) and J0(z) are Bessel functions of the first kind.
4 Downloads
Updated 2 Apr 2025

View License

ZJ1CJ0 - Find the first N non-negative roots of the transcendental equation z*J1(z) - c*J0(z) = 0, where J1(z) and J0(z) are Bessel functions of the first kind of order 1 and 0, respectively.
Syntax:
  • roots = zj1cj0(c, N)
  • roots = zj1cj0(c) % Default N: 10
Description:
This function finds the first N non-negative roots of the transcendental equation z*J1(z) = c*J0(z), where J0 and J1 are Bessel functions of the first kind.
Inputs:
  • c - A non-negative constant, typical the Biot number Bi in heat transder or the Sherwood number Sh in mass transfer.
  • N - (Optional) Number of positive roots to find. Default is 10.
Outputs:
  • roots - A 1xN array of the first N non-negative roots
Example:
  • roots = zj1cj0(2, 5) % Finds 5 roots when c = 2
  • roots = zj1cj0(2) % Finds 10 roots when c = 2 (default N)
See also: besselj, fzero
Copyright © 2025 Jie Zhang (https://orcid.org/0000-0002-8023-8144)

Cite As

Jie Zhang (2026). Transcendental equation zj1=cj0 involving Bessel functions (https://www.mathworks.com/matlabcentral/fileexchange/180578-transcendental-equation-zj1-cj0-involving-bessel-functions), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2021a
Compatible with any release
Platform Compatibility
Windows macOS Linux
Version Published Release Notes
1.0.1

For small Bi or Sh numbers, e.g., c < 0.04, the first root might be ignored in the previous version. Simply changing the start_z to a smaller value will solve the problem. -- 2025.04.02 V1.01

1.0.0