正太分布图 求面积 保留4位有效数字

/**
 * 正太分布图 求面积  保留4位有效数字
 * User : FHD
 * Time : 2020/9/1 0001 17:56
 * @param $a
 * @return float
 */
public function normdist($a)
{
    $p = 0.2316419;
    $b1 = 0.31938153;
    $b2 = -0.356563782;
    $b3 = 1.781477937;
    $b4 = -1.821255978;
    $b5 = 1.330274429;

    $x = abs($a);
    $t = 1 / (1 + $p * $x);

    $val = 1 - (1 / (sqrt(2 * M_PI)) * pow(M_E, (-1 * pow($a, 2) / 2))) *
        ($b1 * $t + $b2 * pow($t, 2) + $b3 * pow($t, 3) + $b4 * pow($t, 4) + $b5 * pow($t, 5));
    if ($a < 0) {
        $val = 1 - $val;
    }
    return round($val, 4);

}