| 488 | */ |
| 489 | |
| 490 | double p1evl(double x, double coef[],int N) |
| 491 | { |
| 492 | double ans; |
| 493 | double *p; |
| 494 | int i; |
| 495 | |
| 496 | p = coef; |
| 497 | ans = x + *p++; |
| 498 | i = N-1; |
| 499 | |
| 500 | do |
| 501 | ans = ans * x + *p++; |
| 502 | while( --i ); |
| 503 | |
| 504 | return( ans ); |
| 505 | } |
| 506 | |
| 507 | double TWOOPI = 6.36619772367581343075535E-1; /* 2/pi */ |
| 508 | double THPIO4 = 2.35619449019234492885; /* 3*pi/4 */ |