3f56f1f9
 dmayerich
initial commit
|
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
|
#ifndef bessH
#define bessH
#define _USE_MATH_DEFINES
#include <math.h>
#include <complex>
using namespace std;
#define eps 1e-15
#define el 0.5772156649015329
int msta1(double x,int mp);
int msta2(double x,int n,int mp);
int bessjy01a(double x,double &j0,double &j1,double &y0,double &y1,
double &j0p,double &j1p,double &y0p,double &y1p);
int bessjy01b(double x,double &j0,double &j1,double &y0,double &y1,
double &j0p,double &j1p,double &y0p,double &y1p);
int bessjyna(int n,double x,int &nm,double *jn,double *yn,
double *jnp,double *ynp);
int bessjynb(int n,double x,int &nm,double *jn,double *yn,
double *jnp,double *ynp);
int bessjyv(double v,double x,double &vm,double *jv,double *yv,
double *jvp,double *yvp);
int bessik01a(double x,double &i0,double &i1,double &k0,double &k1,
double &i0p,double &i1p,double &k0p,double &k1p);
int bessik01b(double x,double &i0,double &i1,double &k0,double &k1,
double &i0p,double &i1p,double &k0p,double &k1p);
int bessikna(int n,double x,int &nm,double *in,double *kn,
double *inp,double *knp);
int bessiknb(int n,double x,int &nm,double *in,double *kn,
double *inp,double *knp);
int bessikv(double v,double x,double &vm,double *iv,double *kv,
double *ivp,double *kvp);
int cbessjy01(complex<double> z,complex<double> &cj0,complex<double> &cj1,
complex<double> &cy0,complex<double> &cy1,complex<double> &cj0p,
complex<double> &cj1p,complex<double> &cy0p,complex<double> &cy1p);
int cbessjyna(int n,complex<double> z,int &nm,complex<double> *cj,
complex<double> *cy,complex<double> *cjp,complex<double> *cyp);
int cbessjynb(int n,complex<double> z,int &nm,complex<double> *cj,
complex<double> *cy,complex<double> *cjp,complex<double> *cyp);
int cbessik01(complex<double>z,complex<double>&ci0,complex<double>&ci1,
complex<double>&ck0,complex<double>&ck1,complex<double>&ci0p,
complex<double>&ci1p,complex<double>&ck0p,complex<double>&ck1p);
int cbessikna(int n,complex<double> z,int &nm,complex<double> *ci,
complex<double> *ck,complex<double> *cip,complex<double> *ckp);
int cbessiknb(int n,complex<double> z,int &nm,complex<double> *ci,
complex<double> *ck,complex<double> *cip,complex<double> *ckp);
int cbessjyva(double v,complex<double> z,double &vm,complex<double>*cjv,
complex<double>*cyv,complex<double>*cjvp,complex<double>*cyvp);
int cbessikv(double v,complex<double>z,double &vm,complex<double> *civ,
complex<double> *ckv,complex<double> *civp,complex<double> *ckvp);
#endif
|
