public static function upperIncompleteGamma($s, $x) { if ($s <= 0) { throw new Exception\OutOfBoundsException("S must be > 0. S = {$s}"); } return self::gamma($s) - self::lowerIncompleteGamma($s, $x); }