@@ 368-377 (lines=10) @@ | ||
365 | $xden = $xden * $xm1 + $lg_q1[$i]; |
|
366 | } |
|
367 | $res = $corr + $xm1 * ($lg_d1 + $xm1 * ($xnum / $xden)); |
|
368 | } else { |
|
369 | $xm2 = $y - 1.0; |
|
370 | $xden = 1.0; |
|
371 | $xnum = 0.0; |
|
372 | for ($i = 0; $i < 8; ++$i) { |
|
373 | $xnum = $xnum * $xm2 + $lg_p2[$i]; |
|
374 | $xden = $xden * $xm2 + $lg_q2[$i]; |
|
375 | } |
|
376 | $res = $corr + $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
377 | } |
|
378 | } elseif ($y <= 4.0) { |
|
379 | // --------------------- |
|
380 | // 1.5 .LT. X .LE. 4.0 |
|
@@ 378-390 (lines=13) @@ | ||
375 | } |
|
376 | $res = $corr + $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
377 | } |
|
378 | } elseif ($y <= 4.0) { |
|
379 | // --------------------- |
|
380 | // 1.5 .LT. X .LE. 4.0 |
|
381 | // --------------------- |
|
382 | $xm2 = $y - 2.0; |
|
383 | $xden = 1.0; |
|
384 | $xnum = 0.0; |
|
385 | for ($i = 0; $i < 8; ++$i) { |
|
386 | $xnum = $xnum * $xm2 + $lg_p2[$i]; |
|
387 | $xden = $xden * $xm2 + $lg_q2[$i]; |
|
388 | } |
|
389 | $res = $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
390 | } elseif ($y <= 12.0) { |
|
391 | // ---------------------- |
|
392 | // 4.0 .LT. X .LE. 12.0 |
|
393 | // ---------------------- |