| @@ 357-366 (lines=10) @@ | ||
| 354 | $xden = $xden * $xm1 + $lg_q1[$i]; |
|
| 355 | } |
|
| 356 | $res = $corr + $xm1 * ($lg_d1 + $xm1 * ($xnum / $xden)); |
|
| 357 | } else { |
|
| 358 | $xm2 = $y - 1.0; |
|
| 359 | $xden = 1.0; |
|
| 360 | $xnum = 0.0; |
|
| 361 | for ($i = 0; $i < 8; ++$i) { |
|
| 362 | $xnum = $xnum * $xm2 + $lg_p2[$i]; |
|
| 363 | $xden = $xden * $xm2 + $lg_q2[$i]; |
|
| 364 | } |
|
| 365 | $res = $corr + $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
| 366 | } |
|
| 367 | } elseif ($y <= 4.0) { |
|
| 368 | // --------------------- |
|
| 369 | // 1.5 .LT. X .LE. 4.0 |
|
| @@ 367-379 (lines=13) @@ | ||
| 364 | } |
|
| 365 | $res = $corr + $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
| 366 | } |
|
| 367 | } elseif ($y <= 4.0) { |
|
| 368 | // --------------------- |
|
| 369 | // 1.5 .LT. X .LE. 4.0 |
|
| 370 | // --------------------- |
|
| 371 | $xm2 = $y - 2.0; |
|
| 372 | $xden = 1.0; |
|
| 373 | $xnum = 0.0; |
|
| 374 | for ($i = 0; $i < 8; ++$i) { |
|
| 375 | $xnum = $xnum * $xm2 + $lg_p2[$i]; |
|
| 376 | $xden = $xden * $xm2 + $lg_q2[$i]; |
|
| 377 | } |
|
| 378 | $res = $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
| 379 | } elseif ($y <= 12.0) { |
|
| 380 | // ---------------------- |
|
| 381 | // 4.0 .LT. X .LE. 12.0 |
|
| 382 | // ---------------------- |
|