#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