j1f
float j1f(float);
PROTO_NORMAL(j1f);
return(x*(u/v) + tpi*(j1f(x)*logf(x)-one/x));
DEF_NONSTD(j1f);
if(n==1) return(j1f(x));
b = j1f(x);