Lines Matching defs:j1
33 static char sccsid[] = "@(#)j1.c 8.2 (Berkeley) 11/30/93";
67 /* double j1(double x), y1(double x)
69 * Method -- j1(x):
70 * 1. For tiny x, we use j1(x) = x/2 - x^3/16 + x^5/384 - ...
71 * 2. Reduce x to |x| since j1(x)=-j1(-x), and
73 * j1(x) = x/2 + x*z*R0/S0, where z = x*x;
74 * (precision: |j1/x - 1/2 - R0/S0 |<2**-61.51 )
76 * j1(x) = sqrt(2/(pi*x))*(p1(x)*cos(x1)-q1(x)*sin(x1))
89 * j1(nan)= nan
90 * j1(0) = 0
91 * j1(inf) = 0
97 * y1(x) = 2/pi*(j1(x)*(ln(x/2)+Euler)-1/x-x/2+5/64*x^3-...)
98 * therefore y1(x)-2/pi*j1(x)*ln(x)-1/x is an odd function.
100 * y1(x) = x*U(z)/V(z) + (2/pi)*(j1(x)*ln(x)-1/x), z= x^2
147 j1(double x)
171 * j1(x) = 1/sqrt(pi) * (P(1,x)*cc - Q(1,x)*ss) / sqrt(x)
277 return (x*(u/v) + tpi*(j1(x)*log(x)-one/x));
Indexes created Mon Sep 29 13:09:55 GMT 2025