Erf.php

Go to the documentation of this file.

<?php
 
namespace PhpOffice\PhpSpreadsheet\Calculation\Engineering;
 
use PhpOffice\PhpSpreadsheet\Calculation\Functions;
 
class Erf
{
    private static $twoSqrtPi = 1.128379167095512574;
 
    public static function ERF($lower, $upper = null)
    {
        $lower = Functions::flattenSingleValue($lower);
        $upper = Functions::flattenSingleValue($upper);
 
        if (is_numeric($lower)) {
            if ($upper === null) {
                return self::erfValue($lower);
            }
            if (is_numeric($upper)) {
                return self::erfValue($upper) - self::erfValue($lower);
            }
        }
 
        return Functions::VALUE();
    }
 
    public static function ERFPRECISE($limit)
    {
        $limit = Functions::flattenSingleValue($limit);
 
        return self::ERF($limit);
    }
 
    //
    //    Private method to calculate the erf value
    //
    public static function erfValue($value)
    {
        if (abs($value) > 2.2) {
            return 1 - ErfC::ERFC($value);
        }
        $sum = $term = $value;
        $xsqr = ($value * $value);
        $j = 1;
        do {
            $term *= $xsqr / $j;
            $sum -= $term / (2 * $j + 1);
            ++$j;
            $term *= $xsqr / $j;
            $sum += $term / (2 * $j + 1);
            ++$j;
            if ($sum == 0.0) {
                break;
            }
        } while (abs($term / $sum) > Functions::PRECISION);
 
        return self::$twoSqrtPi * $sum;
    }
}