@@ 355-364 (lines=10) @@ | ||
352 | $xden = $xden * $xm1 + $lg_q1[$i]; |
|
353 | } |
|
354 | $res = $corr + $xm1 * ($lg_d1 + $xm1 * ($xnum / $xden)); |
|
355 | } else { |
|
356 | $xm2 = $y - 1.0; |
|
357 | $xden = 1.0; |
|
358 | $xnum = 0.0; |
|
359 | for ($i = 0; $i < 8; ++$i) { |
|
360 | $xnum = $xnum * $xm2 + $lg_p2[$i]; |
|
361 | $xden = $xden * $xm2 + $lg_q2[$i]; |
|
362 | } |
|
363 | $res = $corr + $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
364 | } |
|
365 | } elseif ($y <= 4.0) { |
|
366 | // --------------------- |
|
367 | // 1.5 .LT. X .LE. 4.0 |
|
@@ 365-377 (lines=13) @@ | ||
362 | } |
|
363 | $res = $corr + $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
364 | } |
|
365 | } elseif ($y <= 4.0) { |
|
366 | // --------------------- |
|
367 | // 1.5 .LT. X .LE. 4.0 |
|
368 | // --------------------- |
|
369 | $xm2 = $y - 2.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 = $xm2 * ($lg_d2 + $xm2 * ($xnum / $xden)); |
|
377 | } elseif ($y <= 12.0) { |
|
378 | // ---------------------- |
|
379 | // 4.0 .LT. X .LE. 12.0 |
|
380 | // ---------------------- |