https://github.com/sheetjs/bessel
Bessel Functions in JS
https://github.com/sheetjs/bessel
Last synced: 10 months ago
JSON representation
Bessel Functions in JS
- Host: GitHub
- URL: https://github.com/sheetjs/bessel
- Owner: SheetJS
- License: apache-2.0
- Created: 2013-03-16T02:07:44.000Z (almost 13 years ago)
- Default Branch: master
- Last Pushed: 2019-10-08T06:14:22.000Z (over 6 years ago)
- Last Synced: 2025-04-19T08:42:40.881Z (10 months ago)
- Language: JavaScript
- Homepage: http://oss.sheetjs.com/bessel/
- Size: 61.5 KB
- Stars: 35
- Watchers: 8
- Forks: 6
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# bessel
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
## Installation
With [npm](https://www.npmjs.org/package/bessel):
```bash
$ npm install bessel
```
In the browser:
```html
```
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`
## Usage
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:
```js
// 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
```
## Testing
`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:
```mathematica
(* 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.
## License
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.
## Badges
[](https://saucelabs.com/u/bessel)
[](https://travis-ci.org/SheetJS/bessel)
[](https://coveralls.io/r/SheetJS/bessel?branch=master)
[](https://npmjs.org/package/bessel)
[](https://david-dm.org/sheetjs/bessel)
[](https://ghit.me/repo/sheetjs/bessel)
[](https://github.com/SheetJS/bessel)