picomath

sub expm1 {
    my $x = $_[0];

    if (abs($x) < 1e-5) {
        return $x + 0.5*$x*$x;
    } else {
        return exp($x) - 1.0;
    }
}

1;