double dbesj1(double);
void dbmp1(double,double*,double*);
double dcsevl(double,double[],long);
void dinits(double[],long,double,long*);