271 	double	x, z, s, c, ss, cc, r, t, p0, q0;  in __k_j0f()  local
293 			z = t * t;  in __k_j0f()
294 			p0 = (pr[0] + z * (pr[1] + z * (pr[2] + z * (pr[3] +  in __k_j0f()
295 			    z * (pr[4] + z * (pr[5] + z * pr[6])))))) /  in __k_j0f()
296 			    (ps[0] + z * (ps[1] + z * (ps[2] + z * (ps[3] +  in __k_j0f()
297 			    z * (ps[4] + z * (ps[5] + z))))));  in __k_j0f()
298 			q0 = ((qr[0] + z * (qr[1] + z * (qr[2] + z * (qr[3] +  in __k_j0f()
299 			    z * (qr[4] + z * (qr[5] + z * qr[6])))))) /  in __k_j0f()
300 			    (qs[0] + z * (qs[1] + z * (qs[2] + z * (qs[3] +  in __k_j0f()
301 			    z * (qs[4] + z * (qs[5] + z))))))) * t;  in __k_j0f()
311 	z = x * x;  in __k_j0f()
314 		r = r0[0] + z * (r0[1] + z * (r0[2] + z * r0[3]));  in __k_j0f()
315 		s = s0[0] + z * (s0[1] + z * (s0[2] + z * s0[3]));  in __k_j0f()
316 		return (one + z * (r / s));  in __k_j0f()
321 		r = r * z + r1[i];  in __k_j0f()
322 		s = s * z + s1[i];  in __k_j0f()
358 	double	x, z, s, c, ss, cc, t, p0, q0, u, v;  in __k_y0f()  local
380 			z = t * t;  in __k_y0f()
381 			p0 = (pr[0] + z * (pr[1] + z * (pr[2] + z * (pr[3] +  in __k_y0f()
382 			    z * (pr[4] + z * (pr[5] + z * pr[6])))))) /  in __k_y0f()
383 			    (ps[0] + z * (ps[1] + z * (ps[2] + z * (ps[3] +  in __k_y0f()
384 			    z * (ps[4] + z * (ps[5] + z))))));  in __k_y0f()
385 			q0 = ((qr[0] + z * (qr[1] + z * (qr[2] + z * (qr[3] +  in __k_y0f()
386 			    z * (qr[4] + z * (qr[5] + z * qr[6])))))) /  in __k_y0f()
387 			    (qs[0] + z * (qs[1] + z * (qs[2] + z * (qs[3] +  in __k_y0f()
388 			    z * (qs[4] + z * (qs[5] + z))))))) * t;  in __k_y0f()
394 	z = x * x;  in __k_y0f()
397 		u = u * z + u0[i];  in __k_y0f()
398 	v = v0[0] + z * (v0[1] + z * (v0[2] + z * (v0[3] + z * v0[4])));  in __k_y0f()
437 	double	x, z, s, c, ss, cc, r, t, p1, q1;  in __k_j1f()  local
461 			z = t * t;  in __k_j1f()
462 			p1 = (pr0[0] + z * (pr0[1] + z * (pr0[2] + z *  in __k_j1f()
463 			    (pr0[3] + z * (pr0[4] + z * pr0[5]))))) /  in __k_j1f()
464 			    (ps0[0] + z * (ps0[1] + z * (ps0[2] + z *  in __k_j1f()
465 			    (ps0[3] + z * (ps0[4] + z * (ps0[5] + z))))));  in __k_j1f()
466 			q1 = ((qr0[0] + z * (qr0[1] + z * (qr0[2] + z *  in __k_j1f()
467 			    (qr0[3] + z * (qr0[4] + z * qr0[5]))))) /  in __k_j1f()
468 			    (qs0[0] + z * (qs0[1] + z * (qs0[2] + z *  in __k_j1f()
469 			    (qs0[3] + z * (qs0[4] + z * (qs0[5] + z))))))) * t;  in __k_j1f()
482 	z = x * x;  in __k_j1f()
485 		r = a0[0] + z * (a0[1] + z * (a0[2] + z * a0[3]));  in __k_j1f()
486 		s = b0[0] + z * (b0[1] + z * (b0[2] + z * b0[3]));  in __k_j1f()
487 		t = x * half + x * (z * (r / s));  in __k_j1f()
491 			r = r * z + a1[i];  in __k_j1f()
492 		s = b1[0] + z * (b1[1] + z * (b1[2] + z * (b1[3] + z * b1[4])));  in __k_j1f()
526 	double	x, z, s, c, ss, cc, u, v, p1, q1, t;  in __k_y1f()  local
548 			z = t * t;  in __k_y1f()
549 			p1 = (pr0[0] + z * (pr0[1] + z * (pr0[2] + z *  in __k_y1f()
550 			    (pr0[3] + z * (pr0[4] + z * pr0[5]))))) /  in __k_y1f()
551 			    (ps0[0] + z * (ps0[1] + z * (ps0[2] + z *  in __k_y1f()
552 			    (ps0[3] + z * (ps0[4] + z * (ps0[5] + z))))));  in __k_y1f()
553 			q1 = ((qr0[0] + z * (qr0[1] + z * (qr0[2] + z *  in __k_y1f()
554 			    (qr0[3] + z * (qr0[4] + z * qr0[5]))))) /  in __k_y1f()
555 			    (qs0[0] + z * (qs0[1] + z * (qs0[2] + z *  in __k_y1f()
556 			    (qs0[3] + z * (qs0[4] + z * (qs0[5] + z))))))) * t;  in __k_y1f()
562 	z = x * x;  in __k_y1f()
565 		u = c0[0] + z * (c0[1] + z * (c0[2] + z * c0[3]));  in __k_y1f()
566 		v = d0[0] + z * (d0[1] + z * (d0[2] + z * (d0[3] + z * d0[4])));  in __k_y1f()
570 			u = u * z + c1[i];  in __k_y1f()
571 		v = d1[0] + z * (d1[1] + z * (d1[2] + z * (d1[3] + z * d1[4])));  in __k_y1f()
610 	double	a, b, temp, x, z, w, t, q0, q1, h;  in jnf()  local
701 		z = w + h;  in jnf()
702 		q1 = w * z - one;  in jnf()
706 			z += h;  in jnf()
707 			temp = z * q1 - q0;  in jnf()