-
Notifications
You must be signed in to change notification settings - Fork 0
/
bessel.h
59 lines (56 loc) · 1.61 KB
/
bessel.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
#pragma once
/* bessel.c
Copyright (c) 1998
Kapteyn Institute Groningen
All Rights Reserved.
*/
/*
#> bessel.dc2
Function: BESSEL
Purpose: Evaluate Bessel function J, Y, I, K of integer order.
Category: MATH
File: bessel.c
Author: M.G.R. Vogelaar
Use: See bessj.dc2, bessy.dc2, bessi.dc2 or bessk.dc2
Description: The differential equation
2
2 d w dw 2 2
x . --- + x . --- + (x - v ).w = 0
2 dx
dx
has two solutions called Bessel functions of the first kind
Jv(x) and Bessel functions of the second kind Yv(x).
The routines bessj and bessy return the J and Y for
integer v and therefore are called Bessel functions
of integer order.
The differential equation
2
2 d w dw 2 2
x . --- + x . --- - (x + v ).w = 0
2 dx
dx
has two solutions called modified Bessel functions
Iv(x) and Kv(x).
The routines bessi and bessk return the I and K for
integer v and therefore are called Modified Bessel
functions of integer order.
(Abramowitz & Stegun, Handbook of mathematical
functions, ch. 9, pages 358,- and 374,- )
The implementation is based on the ideas from
Numerical Recipes, Press et. al.
This routine is NOT callable in FORTRAN.
Updates: Jun 29, 1998: VOG, Document created.
#<
*/
//#> bessel.h
#if !defined(_bessel_h_)
#define _bessel_h_
extern float bessj(int, double);
extern float bessy(int, double);
extern float bessi(int, double);
extern float bessk(int, double);
//double bessj( int, double );
//double bessy( int, double );
//double bessi( int, double );
//double bessk( int, double );
#endif