Pure-JS Bessel Functions

  • bessel

Downloads in past


3101.0.26 years ago11 years agoMinified + gzip package size for bessel in KB


Pure-JS implementation of Bessel functions J,Y,I,K (for the browser and nodejs). Emphasis on correctness and performance for integer order.
The standard notation is used here:
- J is the Bessel function of the first kind - Y is the Bessel function of the second kind - I is the modified Bessel function of the first kind - K is the modified Bessel function of the second kind


With npm:
$ npm install bessel

In the browser:
<script src="bessel.js"></script>

The browser exposes a variable BESSEL
The script will manipulate module.exports if available (e.g. in a CommonJS require context). This is not always desirable. To prevent the behavior, define DO_NOT_EXPORT_BESSEL


In all cases, the relevant function takes two arguments (value, order).
The return value is a JS number. NaN signals an error in calculation.
  • BESSEL.besselj(x, n) computes J_{n}(x)

  • BESSEL.bessely(x, n) computes Y_{n}(x)

  • BESSEL.besseli(x, n) computes I_{n}(x)

  • BESSEL.besselk(x, n) computes K_{n}(x)

For example:
// var BESSEL = require('bessel'); // uncomment this line if in node
BESSEL.besselj(1.5,0)              // 0.5118276712499389
BESSEL.bessely(1.5,0)              // 0.38244892476502895
BESSEL.besseli(1.5,0)              // 1.6467232021476754
BESSEL.besselk(1.5,0)              // 0.2138055693236539

var Y = BESSEL.bessely
Y(Math.PI, 5) + Y(Math.PI, 3) - (2 * 4 / Math.PI) * Y(Math.PI, 4) // 0


make test will run the nodejs-based test.
To generate the excel.tsv test cases, make a 6-column Excel Sheet:
| x | n | BESSELI | BESSELJ | BESSELK | BESSELY | |---|---|:--------------:|:--------------:|:--------------:|:--------------:| | x | n |BESSELI(A1,B1)|BESSELJ(A1,B1)|BESSELK(A1,B1)|BESSELY(A1,B1)|
To generate the mma.tsv test cases, use the Mathematica Bessel Functions:
(* Bessel_[value, order] *)
F[x_,n_]:={x/2,n,BesselI[n,x/2], BesselJ[n,x/2], BesselK[n,x/2], BesselY[n,x/2]}
Do[Print[ExportString[N[F[x,n],10],"csv"]], {n,1,3}, {x,1,20} ]

Note: Each function follows Excel semantics (value, order). Other platforms like Mathematica reverse the argument order.


Please consult the attached LICENSE file for details. All rights not explicitly granted by the Apache 2.0 License are reserved by the Original Author.


Build Status
Build Status
Coverage Status
NPM Downloads
Dependencies Status