1
2/* @(#)e_asin.c 1.3 95/01/18 */
3/*
4 * ====================================================
5 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 *
7 * Developed at SunSoft, a Sun Microsystems, Inc. business.
8 * Permission to use, copy, modify, and distribute this
9 * software is freely granted, provided that this notice
10 * is preserved.
11 * ====================================================
12 */
13
| 1
2/* @(#)e_asin.c 1.3 95/01/18 */
3/*
4 * ====================================================
5 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 *
7 * Developed at SunSoft, a Sun Microsystems, Inc. business.
8 * Permission to use, copy, modify, and distribute this
9 * software is freely granted, provided that this notice
10 * is preserved.
11 * ====================================================
12 */
13
|
17
18/* __ieee754_asin(x)
19 * Method :
20 * Since asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ...
21 * we approximate asin(x) on [0,0.5] by
22 * asin(x) = x + x*x^2*R(x^2)
23 * where
24 * R(x^2) is a rational approximation of (asin(x)-x)/x^3
--- 89 unchanged lines hidden --- | 16
17/* __ieee754_asin(x)
18 * Method :
19 * Since asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ...
20 * we approximate asin(x) on [0,0.5] by
21 * asin(x) = x + x*x^2*R(x^2)
22 * where
23 * R(x^2) is a rational approximation of (asin(x)-x)/x^3
--- 89 unchanged lines hidden --- |