#ifndef bessel2_h
#define bessel2_h

#include <default.h>

static int imaxarg1, imaxarg2;
#define IMAX(a,b) (imaxarg1 = (a), imaxarg2 = (b),(imaxarg1) > (imaxarg2) ? (imaxarg1) : (imaxarg2) )
#define SIGN(a,b) ((b) >= 0.0 ? fabs(a) : -fabs(a))

#define EPS 1.0e-16
#define FPMIN 1.0e-30
#define MAXIT 10000
#define XMIN 2.0
#define PI 3.141592654
#define NUSE1 7
#define NUSE2 8
#define RTPIO2 1.2533141

void sphbes(int n, double x, double *sj, double *sy,
	    double *sjp, double *syp);
void bessjy(double x, double xnu, double *rj, double *ry,
	    double *rjp, double *ryp);
void beschb(double x, double *gam1, double *gam2, double *gampl,
	    double *gammi);
double chebev(double a, double b, double c[], int m, double x);

#endif