j1f
extern float j1f(float x);
return x*(u/v) + tpi*(j1f(x)*logf(x)-1.0f/x);
w = j1f(x);
return j1f(x);
b = j1f(x);