Dataset of Bessel function Jn maxima and minima to 600 orders and 10000 extrema

Bessel functions of the first kind are ubiquitous in the sciences and engineering in solutions to cylindrical problems including electrostatics, heat flow, and the Schrödinger equation. The roots of the Bessel functions are often quoted and calculated, but the maxima and minima for each Bessel function, used to match Neumann boundary conditions, have not had the same treatment. Here we compute 10000 extrema for the first 600 orders of the Bessel function J. To do this, we employ an adaptive root solver bounded by the roots of the Bessel function and solve to an accuracy of 10−19. We compare with the existing literature (to 30 orders and 5 maxima and minima) and the results match exactly. It is hoped that these data provide values needed for orthogonal function expansions and numerical expressions including the calculation of geometric correction factors in the measurement of resistivity of materials, as is done in the original paper using these data.