|
static const static double double | one = 1.00000000000000000000e+00 |
|
static const static double double | C1 = 4.16666666666666019037e-02 |
|
static const static double double | C2 = -1.38888888888741095749e-03 |
|
static const static double double | C3 = 2.48015872894767294178e-05 |
|
static const static double double | C4 = -2.75573143513906633035e-07 |
|
static const static double double | C5 = 2.08757232129817482790e-09 |
|
static const static double double | C6 = -1.13596475577881948265e-11 |
|
double y double | x |
|
double y double | y |
|
int | ix |
|
| z = x * x |
|
| r = z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * C6))))) |
|
| else |
|
| a = one - qx |
|