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.
2 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 (2025). 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

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
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