picomath

def expm1(x)
    if x.abs < 1e-5
        return x + 0.5*x*x
    else
        return Math.exp(x) - 1.0
    end
end