double dbesj0(double);
void dbmp0(double,double*,double*);
double dcsevl(double,double[],long);
void dinits(double[],long,double,long*);