Bessel function of first kind
J = besselj(nu,Z)
J = besselj(nu,Z,1)
The differential equation
where ν is a real constant, is called Bessel's equation, and its solutions are known as Bessel functions.
Jν(z) and J–ν(z) form a fundamental set of solutions of Bessel's equation for noninteger ν. Jν(z) is defined by
where Γ(a) is the gamma function.
Yν(z) is a second solution of Bessel's equation that is linearly independent of Jν(z). It can be computed using bessely.
J = besselj(nu,Z) computes the Bessel function of the first kind, Jν(z), for each element of the array Z. The order nu need not be an integer, but must be real. The argument Z can be complex. The result is real where Z is positive.
If nu and Z are arrays of the same size, the result is also that size. If either input is a scalar, it is expanded to the other input's size.
format long z = (0:0.2:1)'; besselj(1,z) ans = 0 0.09950083263924 0.19602657795532 0.28670098806392 0.36884204609417 0.44005058574493
The Bessel functions are related to the Hankel functions, also called Bessel functions of the third kind,
where is besselh, Jν(z) is besselj, and Yν(z) is bessely. The Hankel functions also form a fundamental set of solutions to Bessel's equation (see besselh).
 Abramowitz, M., and I.A. Stegun, Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series #55, Dover Publications, 1965, sections 9.1.1, 9.1.89, and 9.12, formulas 9.1.10 and 9.2.5.
 Carrier, Krook, and Pearson, Functions of a Complex Variable: Theory and Technique, Hod Books, 1983, section 5.5.
 Amos, D.E., "A Subroutine Package for Bessel Functions of a Complex Argument and Nonnegative Order," Sandia National Laboratory Report, SAND85-1018, May, 1985.
 Amos, D.E., "A Portable Package for Bessel Functions of a Complex Argument and Nonnegative Order," Trans. Math. Software, 1986.