__ieee754_y1
if(n==1) return(sign*__ieee754_y1(x));
b = __ieee754_y1(x);
extern double __ieee754_y1 __P((double));
return __ieee754_y1(x);
z = __ieee754_y1(x);