BESSELJ is a MATLAB library which evaluates Bessel J functions of noninteger order.
The computer code and data files made available on this web page are distributed under the GNU LGPL license.
BESSELJ is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version and a Python version.
FN, a MATLAB library which approximates elementary and special functions using Chebyshev polynomials; functions include Airy, Bessel I, J, K and Y, beta, confluent hypergeometric, error, gamma, log gamma, Pochhammer, Spence; integrals include hyperbolic cosine, cosine, Dawson, exponential, logarithmic, hyperbolic sine, sine; by Wayne Fullerton.
You can go up one level to the MATLAB codes.