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