__ieee754_j1
return(x*(u/v) + tpi*(__ieee754_j1(x)*__ieee754_log(x)-one/x));
b = __ieee754_j1(x);
w = __ieee754_j1(x);
if(n==1) return(__ieee754_j1(x));
extern double __ieee754_j1 __P((double));
return __ieee754_j1(x);
z = __ieee754_j1(x);