docs: function gam1(x) {
        if (x <= -1) { return NaN; }
        if (x < -0.5 || x > 1.5) { return 1 / gamma(x + 1) - 1; }
        if (x <= 0) { return x * (1 + w.evalAt(x)); }
        if (x <= 0.5) { return x * w1.evalAt(x); }
        if (x <= 1) { return (x - 1) / x * w.evalAt(x - 1); }
        return (x - 1) / x * (w1.evalAt(x - 1) - 1);
     }