|
| 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 |
| |