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